数字信号处理有限冲击响应数字滤波器设计.docVIP

数字信号处理教程 课程设计报告书 设计题目：有限冲击响应数字滤波器设计 姓 名： 金燕 学 号： 班 级： 电信101班 信息与电子工程学院 2012年12月24日 数字信号处理课程设计 ——有限冲激响应数字滤波器设计 实验目的：加深对数字滤波器的常用指标和设计过程的理解。 实验原理：低通滤波器的常用指标： 通带边缘频率，阻带边缘频率 ，通带，最小阻带衰减 通带峰值起伏， 最小阻带衰减。 数字滤波器有IIR和FIR两种类型，它们的特点和设计方法不同。 根据阻带衰减选择窗函数 根据过渡带宽度确定N值 求所设计的FIR滤波器的单位抽样响应 计算频率响应，验算指标是否满足要求 窗系数需要实现用窗函数blackman（N）, hamming(N),hanning(N) kaiser(N)产生。 wd=boxcar(N);%数组wd返回N点矩形窗函数 wd=triang(N);%数组wd返回N点三角窗函数 wd=hanning(N);%数组wd返回N点汉宁窗函数 wd=hamming(N);%数组wd返回N点汉明窗函数 wd=blackman(N);%数组wd返回N点布莱克曼窗函数 wd=kaiser(N,beta);%数组wd返回给定β值的 N点凯泽窗函数 以书上P343例7-1为例 wn=hamming(33);%产生窗函数，N=33 nn=[0:1:32]; alfa=(33-1)/2; hd=sin(0.4*pi*(nn-alfa+eps))./(pi*(nn-alfa+eps));%eps为一个非常小的数，防止出%现零为除数 h=hd.*wn; %hd为一个行向量，wn为一个列向量，需要转置 [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1)));%通过求对数得到以db为单位的幅度特性 axis([0,1,-100,10]); grid; xlabel(归一化频率/pi); ylabel(幅度/dB); 三、设计内容： 设计FIR低通滤波器，通带边界频率 ，阻带边界频率 ，阻带衰减不小于50dB。dB，则应该选择海明窗，其阻带最小衰减-53dB满足要求。 △W=Wst-Wp=0.1*pi N=6.6*pi/△W=6.6*pi/(0.1*pi)=66 Wc=0.5*(Wst+WP)=0.25*pi t=(N-1)/2 hd(n)= sin(Wc*(n-t))./(pi*(n-t)),n≠t Wc/pi, n=t 程序： wn=hamming(66); nn=[0:1:65]; alfa=(66-1)/2; hd=sin(0.25*pi*(nn-alfa+eps))./(pi*(nn-alfa+eps)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1))); axis([0,1,-100,10]); grid; xlabel(归一化频率/pi); ylabel(幅度/dB); 汉宁窗： wn=hamming(69); nn=[0:1:68]; alfa=(69-1)/2; hd=sin(0.25*pi*(nn-alfa+eps))./(pi*(nn-alfa+eps)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1))); axis([0,1,-100,10]); grid; xlabel(归一化频率/pi); ylabel(幅度/dB); 设计FIR高通滤波器，通带边界频率 ，阻带边界频率 ，阻带衰减不小于0dB。t=(N-1)/2 hd(n)= （sin(pi*(n-t))-sin(Wc*(n-t))）./(pi*(n-t)),n≠t 程序： wn=hamming(111); nn=[0:1:110]; alfa=(111-1)/2; hd=(sin(pi*(nn-alfa+eps))-sin(0.25*pi*(nn-alfa+eps)))./(pi*(nn-alfa+eps)); h=hd.*wn; [h1,w1]=freqz(h,1); plot(w1/pi,20*log10(abs(h1))); axis([0,1,-100,10]); grid; xlabel(归一化频率/pi); ylabel(幅度/dB); 设计FIR带通滤波器，边界频率，] ，，]，阻带衰减不小于0dB。t=(N-1)/2 hd(n)=