有时候我们需要的滤波器响应并不是简单的带通带阻,我们需要的是一些特定频率响应的滤波器,这时候我们就需要设计任意响应的滤波器设计,用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,逼近的滤波器再高频段响应就比较好,但是低频段拟合的不够好,谁有更好的方法可以给我留言。