matlab 任意频率响应滤波器的设计

有时候我们需要的滤波器响应并不是简单的带通带阻，我们需要的是一些特定频率响应的滤波器，这时候我们就需要设计任意响应的滤波器设计，用matlab以很好地帮我们解决这个问题，当然这样的滤波器是用数字滤波器来实现的。

采样率60K，对于20Hz~20KHz的信号，信号通过一个带阻网络，在某些频率会衰减的比较厉害，这个带阻网络模拟的就是我们通信的环境，这样我们通常需要一个补偿滤波器，这个滤波器的响应和带阻网络的响应正好成倒数关系，这样来构成一个均衡网络。

代码实现：

a=load('data_60K.txt');
a=a/1961;
fs=60000;
semilogx(20*(1:3000),20*log10(a));grid; %原始网络曲线
legend('带阻网络');
a=1./a;
f=1:3000;f=f*20;f=f/fs;   %归一化
N=1000;     %FIR滤波器阶
b=firpm(N,f,a);    %计算N阶FIR逼近的系数
fidrin=fopen('coe_60K_firpm.txt','w');

fprintf(fidrin,'%d\n',b);
[h,w] = freqz(b,1, length(a));
figure;
subplot(2,1,1),
semilogx(f*fs,20*log10(a),'r',f*fs,20*log10(abs(h)));grid; %绘制均衡网络曲线
xlabel('a');
legend('均衡网络');
subplot(2,1,2);
semilogx(f*fs,20*log10(abs(h)));grid;  %FIR逼近函数幅频特性曲线;
semilogx(f*fs,20*log10(abs(h))-20*log10(a'));grid;
xlabel('b');
legend('FIR逼近');

将滤波器系数导入quartus的fir IP核就可以实现我们所需要的滤波器设计了。

此方法是基于matlab firpm函数设计的，可以参考matlab中的firpm说明。

不过也可以尝试另外的设计方法，firpm设计的滤波器从图中可以看出，高频段响应不好，会发生振荡，但是用fir2，逼近的滤波器再高频段响应就比较好，但是低频段拟合的不够好，谁有更好的方法可以给我留言。