数字信号处理报告2010.docVIP

数字信号处理实验报告 专业：通信工程 指导老师：吴国平 姓名:张仁琴 学号：20081001863 班级序号：075081-15 FIR数字滤波器设计 一．实验目的 1、设计和应用FIR低通滤波器； 2、掌握FIR数字滤波器的窗函数设计法； 3、了解设计参数（窗型、窗长）的影响。 二、实验内容 （一）、利用 “窗函数设计法”设计FIR数字滤波器。 选择各种窗型、各种窗长、布莱克曼窗的参数的不同数值。观察（时域）窗序列形状、（频域）窗谱形状、主瓣和副瓣的宽度与高度。 1、了解各种窗函数及其主要参数： 1）bartlett：巴特利特窗 2）blackman:布莱克曼窗 3）rectangle:矩形窗 4）hamming:汉明窗 5）hanning:汉宁窗 6）triang:三角窗 7）chebwin:切比雪夫窗 8）kaiser:凯瑟窗 窗函数 窗谱性能指标 加窗后滤波器性能指标 旁瓣峰值（dB） 主瓣宽度（×2π/N） 过渡带宽（×2π/N） 阻带最小衰减 矩形窗 -13 2 0.9 -21 巴特利特窗 -25 4 3.05 -25 汉宁窗 -31 4 3.1 -44 汉明窗 -41 4 3.3 -53 布莱克曼窗 -57 6 5.5 -74 凯瑟窗 -57 5 5 -80 2、学习MATLAB中设计FIR滤波器的函数： 1）fir1:采用窗函数法设计FIR滤波器，调用格式为： b=fir1(N,Wn);b=fir1(N,Wn,’hifh’);b=fir1(N,Wn,’stop’); 其中N为滤波器的阶次，因此滤波器长度为N+1；Wn是通带截止频率，其值在0~1之间，1对应抽样频率的一半；b式设计好的滤波器系数h(n)。窗函数默认值为汉明窗。 2）fir2:采用窗函数法设计具有任意幅频特性的FIR滤波器，调用格式为：b=fir2(N,F,M); 其中F是频率向量，值在0~1之间，M是与F相对应的所希望的幅频响应，默认汉明窗。 3）remez：采用切比雪夫最佳一致逼近FIR数字滤波器。同时，还可以用来设计希尔伯特变换器和差分器。调用格式： b=remez(N,F,A); b=remez(N,F,A,W)等。 其中N为滤波器的阶次，b为设计的滤波器的系数，因此滤波器长度为N+1；F是频率向量，其值在0~1之间；A是对应的F的各频段上的理想幅频响应；W是各频段上的加权向量。 4）remexord：采用切比雪夫一致逼近设计FIR数字滤波器时所需要的滤波器阶次。调用格式: [N,Fo,Ao,W]=remexord(F,A,DEV,Fs); 其中F、A的含义同上。为通带和阻带上的偏差；该函数输出的是符合要求的滤波器阶次N,频率向量Fo，幅度向量Ao和加权向量W. 5）firls：用最小平方法设计线性相位FIR数字滤波器。 6）fircls：用约束的最小平方法设计线性相位FIR数字滤波器。 7）firclsl：用带约束的最小平方设计法设计线性相位FIR低通和高通滤波器。 8）firrcos：用来设计低通线性相位FIR数字滤波器，其过渡带为余弦函数形状。 （二）、编写MATLAB程序，完成以下工作。 1、用窗函数设计一个线性相位的FIR低通滤波器，并满足性能指标；通带边界频率wp=0.5π，阻带边界频率为ws=0.66π，阻带衰减不小于40dB,通带波纹不大于3dB。 1）MATLAB程序： wp=0.5*pi;ws=0.66*pi; %性能指标 wdelta=ws-wp; %通带宽度 N=ceil(8*pi/wdelta) %滤波器长度 if rem(N,2)==0 %类似mod N=N+1; end Nw=N; %窗口长度 wc=(ws+wp)/2; %截止频率 n=0:N-1; %理想滤波器的冲激响应 alpha=(N-1)/2; m=n-alpha+0.00001; hd=sin(wc*m)./(pi*m); win=hamming(Nw); %汉宁窗的时域响应 h=hd.*win; %滤波器实际冲激响应 b=h; figure; freqz(b,1,512); 上述FIR的MATLAB运行结果图如下： 2、用矩形窗和汉明窗设计一个低通滤波器，所希望的频率响应截止频率之间为1，在之间为0，分别取N=11,21,41， 观察并分析N不同时其频谱响应特点。 1）MATLAB程序： N=10; b1=fir1(N,0.25,boxcar(N+1)); %用矩形窗作为冲激响应的窗函数 b2=fir1(N,0.25,hamming(N+1)