高通有限冲激响应滤波器FIR于DSP上的实现.docVIP

高通有限冲激响应滤波器FIR于DSP上的实现 IIR数字滤波器能保留一些典型模拟滤波器优良的幅度特性，但设计中只考虑了幅度特性，没考虑相位特性，所设计的滤波器相位特性一般是非线性的。为了得到线性相位特性，对IIR滤波器必须另外加相位校正网络，是滤波器设计变的复杂，成本也高。FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性。稳定和线性相位特性是FIR滤波器突出的优点。它的系统函数为 本设计利用窗函数法来设计FIR滤波器。由于理想的滤波器的单位取样响应是无限长的，且是非因果的，为了构造一个长度为N的线性相位滤波器，只有将单位取样响应截取一段，这里就用到了窗函数。窗函数设计法的基本原理是用有限长的单位脉冲响应序列h(n)逼近，为期望理想滤波器的单位脉冲响应序列。方法是用窗函数w(n)将截断，并进行加权处理，即 h(n)= w(n) 用h(n)作为实际设计的FIR数字滤波器的单位脉冲响应序列，相应的频率响应函数为用窗函数设计的实际滤波器的频响函数。调整窗口长度N可以有效的控制过渡带的宽度，减少带内波动以及加大阻带的衰减只能从窗函数的形状上找解决方法。这里利用了哈明窗，它的旁瓣峰值幅度为-41/dB，过渡带宽为8pi/N，阻带最小衰减为-53/dB。 利用matlab产生一个80阶的高通滤波器，采用哈明窗，采样频率为6kHz，截止频率为3kHz，matlab程序如下： % n： 阶数 % fp：高通截止频率 % fsample：采样率 % w：窗函数值 % h：返回设计好的FIR滤波器系数 n=120; fsample=7000; fp=3000; w=hamming(n+1); wp=fp/(fsample/2); %ws=fs/(fsample/2); h=fir1(n,wp,high,w); m=0:n; figure(1) subplot(2,1,1);plot(m,w);grid; ylabel(w(n));xlabel(n); subplot(2,1,2);plot(m,h);grid; ylabel(h(n));xlabel(n); figure(2) freqz(h) xlabel(\omega^pi) ylabel(|H(e^j\omega|dB)) subplot(2,1,2) xlabel(\omega^pi) 程序产生的图形如下： 图一 图一中上图是窗函数参数，下图是滤波器参数 图二 图二中上图是相对幅值响应，下图是相位响应 2. 利用matlab求的h的参数，编写dsp的程序，相应程序如下： /* fir_test.c */ #include math.h #define L 80 /* Number of FIR filter coefficients */ #define Fs 6000 /* 9000 Hz sampling frequency */ #define T 1/Fs #define f1 400 /* 800 Hz frequency */ #define f2 1800 /* 1800 Hz frequency */ #define f3 3300 /* 3300 Hz frequency */ #define PI 3.1415926 #define w1 (2*PI*f1*T) /* 2*pi*f1/Fs */ #define w2 (2*PI*f2*T) /* 2*pi*f2/Fs */ #define w3 (2*PI*f3*T) /* 2*pi*f3/Fs */ #define a1 0.333 /* Magnitude for wave 1 */ #define a2 0.30 /* Magnitude for wave 2 */ #define a3 0.30 /* Magnitude for wave 3 */ extern float firC(float *,float *,unsigned int); /* High-pass FIR filter coefficients */ float hh[]={-1.14361737749628e-020 ,2.81071199533121e-005, -0.000104762148808631, 0.000190622649920012 ,-0.000188360514955522, 5.31949187585617e-018, 0.00038632360096985, -0.000834992