次のように、片面パワースペクトルを正常に生成したとします。
X_mags = abs(fft(signal));
bin_vals = [0 : N-1];
fax_Hz = bin_vals*fs/N;
N_2 = ceil(N/2);
plot(fax_Hz(1:N_2), 20*log10(X_mags(1:N_2)));`
ここで、最初のグラフの上に 2 番目のグラフをプロットします。
hold on;
最後に、信号に対して LPC 分析を実行し、周波数応答を計算します。周波数応答は、パワー スペクトルの上にプロットする必要があるため、次のようになります。
[a, g] = lpc(signal,N);
[h,f] = freqz(b,a,N,fs);
plot(?);
簡単にするために、パラメータがすべて正しく与えられていると仮定しましょう。周波数応答を正しく表示するには、プロット関数をどのように記述すればよいでしょうか? 単純な plot(f) は機能しません。
誰かが理由を説明できますか？ありがとう