《基于频率采样法FIR带通滤波器设计》由会员分享，可在线阅读，更多相关《基于频率采样法FIR带通滤波器设计（15页珍藏版）》请在金锄头文库上搜索。
1、*大学电子信息工程学院DSP课程设计报告题 目：基于频率采样法的 FIR带通滤波器的设计专业班级：通信工程组 员：*指导老师：*2011年 6月 13日2目 录一、设计目的3二、设计要求3三、设计原理4四、程序代码6五、调试分析12六、设计总结13七、设计心得13八、参考文献3一、 设计目的：1、掌握用频率取样法设计 FIR DF的方法，并掌握该方法的MATLAB编程。2、熟悉频率取样理论，熟悉内插函数及其应用。3、了解 FIR DF的频率特性和相位特性，观察过渡带取样点对滤波器幅频特性的影响。 二、 设计内容:基于 MATLAB结合 FFT和 IFFT，利用频率采样法设计 FIR数字带通滤波
2、器，然后用自己设计的滤波器对采集的加噪后的语音信号进行滤波，并将滤波前后的信号进行比较，回放语音信号。三、设计原理：FIR滤波器的单位脉冲响应 是有限长的（ ） ，其nh10Nnz变换为 的（N-1）阶多项式：1znNnhzXYH0可得 FIR滤波器的系统差分方程：nxbmnxbNnxbnxbny m10110因此，FIR 滤波器又称为卷积滤波器。FIR 滤波器的频率响4应表达式为： 信号通过 FIR滤波器不失ejwnNnjwhH10真条件是在通带内具有恒定的幅频特性和相位线性特性。理论上可以证明：当 FIR滤波器的系数满足下列中心对称条件：时，滤波器设计在逼近平直nNhnNhn11或 者幅频
3、特性的同时，还能获得严格的线性相位特性。线性相位 FIR滤波器的相位滞后和群延迟在整个频带上是相等且不变的。对于一个 N阶的线性相位 FIR滤波器，群延迟为常数，即滤波后的信号简单地延迟常数个时间步长。这一特性使通带频域内信号通过滤波器后仍保持原有波形形状而无相位失真。FIR滤波器设计的方法很多，如窗函数法、频率采样法以及其他的各种优化设计方法，本次课程设计使用频率采样法设计 FIR带通滤波器。频率采样法是从频域出发。因为有限长序列 h(n)又可用其离散傅立叶变换 H(k)来唯一表示，H(k)与所要求的 FIR 滤波器系统函数 Hd(z)之间存在着频率取样关系。即 Hd(z)在 Z 平面单位圆
4、上按角度等分的取样值等于 Hd(k)的各相应值，就以此 Hd(k)值作为实际 FIR 数字滤波器频率特性的取样值 H(k),或者说 H(k)正是所要求的频率响应 H(ejw)的 N 各等间隔的取样值。频率采样法就是根据频域采样理论，由滤波特性指标构造希望逼近的滤波器频响函数 Hd(ej)，对其在0,2上采样得到:5然后，就可求出单位脉冲响应 h(n)，或是系统函数 H(z)。这样，h(n)或是 H(z)就是滤波器的设计结果。频率取样法设计的基本思想：把给出的理想频率响应进行取样，通过 IDFT 从频谱样点直接求得有限脉冲响应。其设计过程如下频率取样法的关键是正确确定数字频域系统函数 H(k)在
5、0，2 内的 N 个样点，其约束条件为H(k)=H(N-k)（m）=- (N-m) 0kN-1频率采样法的优点是可以在频域直接设计，并且适合最优化设计；缺点是采样频率只能等于 的整数倍，因而不能确保截N/2止频率 的自由取值，要想实现自由地选择截止频率，必须增wc1102 NkeHkkNjdd ,L内 插 公 式 1,1,011 ,101 NkzWkHNzzHnIDFTnhkNL)()()()(H2 jDTFIDTFNkjd eHnhkHe6加采样点数 N，但是这又使计算量加大。四、程序代码:语音程序：filename=111;s,fs,nbits=wavread(filename);soun
6、d(s,fs,nbits); %回放语音信号n = length (s) ; %求出语音信号的长度t=0:1/fs:(n-1)/fs;Y=fft(s); %傅里叶变换figure(1)subplot(2,1,1);plot(s);title(原始信号波形);grid onsubplot(2,1,2);plot(abs(Y);title(原始信号频谱)sound(s,fs);grid on仿真结果：图一带通滤波器程序：N=40;alfa=(40-1)/2;k=0:N-1;w1=(2*pi/N)*k;T1=0.109021; T2=0.59417456;hrs=zeros(1,5),T1,T2,o
7、nes(1,7),T2,T1,zeros(1,9),T1,T2,ones(1,7),T2,T1,zeros(1,4);hdr=0,0,1,1,0,0; wd1=0,0.2,0.35,0.65,0.8,1;k1=0:floor(N-1)/2); k2=floor(N-1)/2)+1:N-1;angH=-alfa*(2*pi)/N*k1,alfa*(2*pi/N*(N-k2);7H=hrs.*exp(j*angH);h=real(ifft(H);db,mag,pha,grd,w = freqz_m(h,1);Hr,ww,a,L =Hr_Type2(h);subplot(2,2,1)plot(w1(
8、1:21)/pi,hrs(1:21),o,wd1,hdr)axis(0,1,-0.1,1.1);title(带通：N=40，T1=0.109021， T2=0.59417456)ylabel(Hr(k);set(gca,XTickMode,manual,XTick,0,0.2,0.35,0.65,0.8,1)set(gca,YTickMode,manual,YTick,0,0.059,0.109,1);grid %绘制带网格的图像subplot(2,2,2); stem(k,h);axis(-1,N,-0.4,0.4)title(脉冲响应); ylabel(h(n); text(N+1,-0.
9、4,n)subplot(2,2,3); plot(ww/pi,Hr,w1(1:21)/pi,hrs(1:21),o);axis(0,1,-0.1,1.1);title(振幅响应)xlabel(频率 （单位：pi）); ylabel(Hr(w)set(gca,XTickMode,manual,XTick,0,0.2,0.35,0.65,0.8,1);set(gca,YTickMode,manual,YTick,0,0.059,0.109,1);gridsubplot(2,2,4); plot(w/pi,db); axis(0,1,-100,10);gridtitle(幅度响应);xlabel(频
10、率 （单位：pi）); ylabel(分贝)set(gca,XTickMode,manual,XTick,0,0.2,0.35,0.65,0.8,1)set(gca,YTickMode,manual,YTick,-60;0);set(gca,YTickLabelMode,manual,YTickLabels,60;0);仿真结果：8图二总程序：%设计一个通带为 0.35 到 0.65，阻带下截止为 0.2，上截止为 0.8Clear;N=40;alfa=(40-1)/2;k=0:N-1;w1=(2*pi/N)*k;T1=0.109021; T2=0.59417456;hrs=zeros(1,5
11、),T1,T2,ones(1,7),T2,T1,zeros(1,9),T1,T2,ones(1,7),T2,T1,zeros(1,4);hdr=0,0,1,1,0,0; wd1=0,0.2,0.35,0.65,0.8,1;k1=0:floor(N-1)/2); k2=floor(N-1)/2)+1:N-1;angH=-alfa*(2*pi)/N*k1,alfa*(2*pi/N*(N-k2);H=hrs.*exp(j*angH);h=real(ifft(H);db,mag,pha,grd,w = freqz_m(h,1);Hr,ww,a,L =Hr_Type2(h);figure(1)subpl
12、ot(2,2,1)plot(w1(1:21)/pi,hrs(1:21),o,wd1,hdr)axis(0,1,-0.1,1.1);title(带通：N=40，T1=0.109021， T2=0.59417456)ylabel(Hr(k);set(gca,XTickMode,manual,XTick,0,0.2,0.35,0.65,0.8,1)set(gca,YTickMode,manual,YTick,0,0.059,0.109,1);grid %绘制带网格的图像subplot(2,2,2); stem(k,h);axis(-1,N,-0.4,0.4)title(脉冲响应); ylabel(h
13、(n); text(N+1,-0.4,n)subplot(2,2,3); plot(ww/pi,Hr,w1(1:21)/pi,hrs(1:21),o);axis(0,1,-0.1,1.1);title(振幅响应)xlabel(频率 （单位：pi）); ylabel(Hr(w)set(gca,XTickMode,manual,XTick,0,0.2,0.35,0.65,0.8,1);set(gca,YTickMode,manual,YTick,0,0.059,0.109,1);gridsubplot(2,2,4); plot(w/pi,db); axis(0,1,-100,10);gridtit
14、le(幅度响应);xlabel(频率 （单位：pi）); ylabel(分贝)set(gca,XTickMode,manual,XTick,0,0.2,0.35,0.65,0.8,1)set(gca,YTickMode,manual,YTick,-60;0);set(gca,YTickLabelMode,manual,YTickLabels,60;0);9s,fs,nbits=wavread(sj.wav);%信号 de 取样频率为 44100HZx=s(:,1);sound(x,fs);L=length(x);f=fs*(0:L-1)/L;t=0:1/fs:(L-1)/fs; %将所加噪声信号的点数调整到与原始信号相同d=0.03*abs(max(x)*cos(2*pi*22000*t); %噪声为 500 和 3300Hz 的余弦信号dz=cos(2*pi*11025*t);xd=x.*dz;xz=xd+d;sound(xz,fs); %播放加噪声后的语音信号X=fft(x); %求信号的频谱XD=fft(xd); %信号调制后的频谱XZ=fft(xz);figure(2)subplot(3,1,1);plot(
