スペクトル形状のグラフ化

ある音データを最大エントロピー法によってベクトル化されたtxtデータがあります。 -1以上1未満の数字が縦に4000個程並べられているのですが、このデータを縦軸dB/Hzで幅が0から-50,横軸Hzで0から2000のグラフを表示させるプログラムのアドバイスをいただきたいです。

Answer 1

[Tohru Kikawada]

最大エントロピー法の出力がマイナスの範囲で出てくる部分が理解できていません。複素数ということでしょうか。パワースペクトルなので0以上の値が出てくるという認識を持っています。

最大エントロピー法でスペクトル推定する方法としては こちら のページのユール・ウォーカー AR 法が分かりやすいかと思います。

ご参考になれば幸いです。

% 音声データのロード
load mtlb

% ユール・ウォーカー AR 法では、最大エントロピー推定器として同じ結果が生成される
order = 14;
[Pxx,F] = pyulear(mtlb,order,1024,Fs);

% 可視化
plot(F,10*log10(Pxx));
xlabel('周波数(kHz)');
ylabel('パワースペクトル密度(dB/Hz)');
grid on;

[Yusaku Ando]

回答ありがとうございます。 おかげさまで理想のグラフが表示できました。