r - Rでピンクノイズをシミュレートする方法

Question

の出力を時系列として扱うことで、ホワイトノイズを実現できることを私は知っていますrnorm()。ピンクノイズをシミュレートする方法について何か提案はありますか?

Answer 1

パッケージtuneRにnoiseは、ホワイト ノイズまたはピンク ノイズのいずれかである波オブジェクトを生成できる関数があります。

require(tuneR)
w <- noise(kind = c("white"))
p <- noise(kind = c("pink"))
par(mfrow=c(2,1))
plot(w,main="white noise")
plot(p,main="pink noise")

編集:上記の方法ではベクトルが生成されないことに気付きました(doh)。それをベクトルに変換する残忍な方法は、以下のコードを追加することです。

writeWave(p,"p.wav")#writes pink noise on your hard drive
require(audio)#loads `audio` package to use `load.wave` function
p.vec <- load.wave("path/to/p.wav")#this will load pink noise as a vector