1/fゆらぎの傾きについて
私は今音の1/f ゆらぎについて研究していています。
まず傾きが既知のホワイトノイズ(傾き0)やピンクノイズ(傾き-1)、ブラウンノイズ(傾き-2)の音をMATLABで高速フーリエ変換を行ってスペクトル解析を行って周波数ごとのパワーを算出してエクセルで両対数をとって直線近似をしています。
結果、ホワイトノイズは0と正しいのですが、ピンクノイズは傾き-0.48、ブラウンノイズは傾き-0.98と約1/2の値となってしまいます。
強度はMATLABで算出された強度Pそのままを使っているのですが、文献などを見ると、dBの値で強度が表されいるので、そこに原因があるかと思うのですが・・・
何かわかる方がいらっしゃればよろしくお願いいたします。
以下プログラム↓
[y,Fs]=wavread('pinknoize.wav');
y1 = y(:,1);
Y = fft(y1,Fs/10);
Y2 = abs(fftshift(Y));
Y3 = Y2(length(Y2)/2+1:length(Y2));
l = linspace(0,Fs/2,length(Y3));
