MATLAB实验四。【DOC精选】.doc

题目四： 一：实验题目： Using the bilinear transformation and a lowpass analog Butterworth prototype, design a highpass digital filter operating at a rate of 20kHz and having passband extending to 5kHz with maximum passband attenuation of 0.5dB,and stopband ending at 4kHz with a minimum stopband attenuation of 10dB. 二：实验目的： （1）熟悉MATLAB语言在数字信号处理中的应用； （2）掌握数字滤波器的计算机仿真方法； （3）掌握IIR数字滤波器设计的一般步骤； （4）熟悉用双线性变换法设计IIR数字滤波器的原理和方法； （5）运用MATLAB平台编辑和设计程序，熟悉MATLAB语言及其调试方法； （6）复习数字信号处理课程中的滤波器原理和IIR数字滤波器设计的原理和方法。 三：实验原理： 数字滤波器的幅度与相位性能指标与模拟滤波器的指标是基本一致的，IIR数字滤波器的设计和模拟滤波器的设计有着很紧密的关系，通常要先设计出适当的模拟滤波器在经过一定的频带变换把它转化为所需的数字IIR滤波器。模拟滤波器具有无限长的脉冲相应，IIR数字滤波器也具有无限长脉冲响应，因此二者相似。IIR滤波器设计的基本思想就是利用复变函数映射把模拟滤波器转换成数字滤波器。双线形变换法是利用s=2*（1-z-1）/T*（1+z-1）将s域转换到z域，从而得到系统函数H（Z）。由于巴特沃斯低通滤波器和高通滤波器只有在归一化频率时和去归一化时关系式互为倒数，因此只要在这两处变换系数，就可以利用沃斯低通滤波器来设计高通滤波器。 四：实验步骤简述： （1）首先将数字频率转换成角频率； （2）其次将数字角频率转换成模拟角频率； （3）计算模拟低通滤波器传输函数； （4）将模拟低通滤波器传输函数转换成模拟高通滤波器传输函数； （5）利用双线性变法将模拟高通滤波器转换成数字高通滤波器。 五：程序框图： 六：源程序： %高通滤波器指标 Fp=5000; Fs=4000; Ap=0.5; As=10; F=20000; T=1/F; wp=2*pi*Fp/F; ws=2*pi*Fs/F; Wp=2*tan(wp/2)/T; Ws=2*tan(ws/2)/T; W1=1/Wp; W2=1/Ws; %求巴特沃思滤波器阶数 [N,W1]=buttord(W1,W2,Ap,As,s) %设计巴特沃思滤波器阶数 [z,p,k]=buttap(N); %计算分子，分母多项式系数 num=k*real(poly(z)); den=real(poly(p)); %频带变换为模拟高通滤波器 [numT denT]=lp2hp(num,den,Wp); %利用双线性变换转换为数字高通滤波器 [numd,dend]=bilinear(numT,denT,F) w=0:pi/1024:pi; %计算系统的频率响应 [H,phw]=freqz(numd,dend,w); %画出系统的对数幅频特性 subplot(2,1,1);plot(w/pi,20*log10(abs(H))); axis([0,1,-30,2]);grid; xlabel(归一化频率); ylabel(幅度 /dB); subplot(2,1,2);plot(w/pi,(abs(H)));grid; xlabel(归一化频率); ylabel(|H(w)|); 七：程序结果及图表： 运行结果：N = 7 numd = 0.0284 -0.1989 0.5968 -0.9946 0.9946 -0.5968 0.1989 -0.0284 dend = 1.0000 -0.7512 1.1198 -0.4631 0.2699 -0.0566 0.0119 -0.0007 二阶基本节形式： H(Z)=(1+Z-1)7/(1.85-0.15Z-1)(2.1007-0.555Z-1+1.3443Z-2)(2.7825-0.555Z-1+0.6625Z-2)(3.2541-0.555Z-1+0.1909Z-2) 八：实验总结： 通过对数字高通滤波器的设计，使我对模拟滤波器和数字滤波器的原理有了更进一步的了解，并且熟悉了设计IIR数字滤波器的一般步骤。已知模拟低通滤波器的设计方法，通过模拟滤波器的频率变换转化为模拟高通，带通和带阻滤波器