使用具有保持功能的Matlab绘制频谱图

Question

假设我已经成功生成了单边功率谱，如下所示：

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)));`

现在我想在第一个图形上绘制第二个图形：

hold on;

最后，我对信号进行LPC分析并计算频率响应。 频率响应应绘制在功率谱的顶部，因此：

[a, g] = lpc(signal,N);
[h,f] = freqz(b,a,N,fs);
plot(?);

为了简单起见，假设所有参数均正确给出，如何编写绘图函数以正确显示频率响应？ 简单的绘图（f）不起作用。

有人可以解释为什么吗？ 谢谢

Answer 1

响应h很复杂，因此您需要通过乘以其复共轭来获得响应的大小。

plot(f, 10*log10(h.*conj(h)));

请注意使用10 * log10，因为上面的运算将保持在h的幅度响应平方。

或者，如果您对处理复数不太感兴趣，则可以取复数值的绝对值，并确保为20 * log10，因为abs不会使数值平方

plot(f, 20*log10(abs(h)));

Answer 2

一个简单的plot(f)试图绘制频率矢量，不是吗？

检查以下代码：

X_mags   = abs(fft(signal));
bin_vals = [0 : N-1];
fax_Hz   = bin_vals*fs/N;
N_2      = ceil(N/2);

[a, g]   = lpc(signal,N);
[h, f]   = freqz(b, a, N, fs);

figure,

hold on,
plot(f, 20*log10(abs(h)), 'r');
plot(fax_Hz(1:N_2), 20*log10(X_mags(1:N_2)));
title('Frequency Spectrum');
xlabel('Frequency (Hz)'); 
ylabel('Amplitude (dB)');
legend('Frequency Response', 'Single-sided Power spectrum')

顺便说一句，有一个用于分贝计算的MATLAB函数db() 。 那可能有用。