matlab fir2频率响应与幅度响应不对应

时间:2017-07-01 09:49:34

标签: matlab audio filter signal-processing

我试图在Matlab中设计一个带fir2函数的FIR,我正在努力得到滤波器的频率响应,这与我想要实现的幅度函数相对应。

plots of desired and computed response

以下是我用于滤镜设计的代码,其中x是幅度响应,如图所示

fs = 48000;    %sample rate
fny = fs/2;    %Nyquist frequency
tabs = 512;    %tabs

%frequency resolution 1/3 octave, length=28
f = [0 63 80 100 125 160 200 250 315 400 500 630 800 1000 1250 1600 
2000 2500 3150 4000 5000 6300 8000 10000 12500 16000 20000 fny];

fn = f/(fny);

coeffs = fir2(tabs,fn,x); 

figure(1);
[freq_response,fc] = freqz(coeffs,tabs);
H = abs(freq_response);

subplot(2,1,1);
semilogx(fc*fny,H, 'r');
hold on;
xlim([20 20000]);
title('frequency response of fir2 filter');
xlabel('frequency [Hz]');
ylabel('magnitude [dB]');
grid on;
set(gca,'XTick',[20 50 100 200 500 1000 2000 5000 10000 20000]);

subplot(2,1,2);
semilogx(f,20*log10(x),'b');
title('desired magnitude function');
xlabel('frequency [Hz]');
ylabel('fir coefficients');
grid on;
xlim([20 20000]);
set(gca,'XTick',[20 50 100 200 500 1000 2000 5000 10000 20000]);

我做错了什么?

1 个答案:

答案 0 :(得分:1)

首先我们应该观察到函数freqz需要滤波器传递函数的分子和分母多项式系数。然后,您还可以指定其他参数,例如频率响应中所需的点数,但这将强制分子和分母多项式系数之后。对于没有显式分母的FIR滤波器,隐式分母是常数1。因此,在您的案例中调用freqz的正确方法是:

[freq_response,fc] = freqz(coeffs,1,tabs); 

请注意,默认情况下,返回的频率是覆盖[0,pi]范围的角频率。要获得以Hz为单位的频率,您需要使用pi重新调整角频率(正如您尝试的那样,但缺少f*fny/pi因子)。或者,您可以简单地将采样率传递给freqz函数:

[freq_response,fc] = freqz(coeffs,1,tabs,fs);

然后,您对图的y轴刻度存在问题。假设您尝试以dB为单位绘制幅度响应,则应将计算出的响应H(线性标度)转换为dB:

[freq_response,fc] = freqz(coeffs,1,tabs,fs);
H = abs(freq_response);
...
semilogx(fc,20*log10(H), 'r');

在同一图表上绘制所需和计算的响应应该会给您类似以下的内容:

enter image description here