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

实验报告 一：实验题目： 题目三 Using the bilinear transformation and a lowpass analog Butterworth prototype, design a lowpass digital filter operating at a rate of 20kHz and having passband extending to 4kHz with maximum passband attenuation of 0.5dB,and stopband starting at 5kHz with a minimum stopband attenuation of 10dB. 二：实验目的 掌握用matlab设计数字低通滤波器的方法，学会利用matlab中的函数实现滤波器的设计。 三：实验原理 双线形变换法是利用s=2*（1-z-1）/T*（1+z-1）将s域转换到z域， 从而得到系统函数H（Z） 四：实验步骤简述 （1）复习有关巴特沃斯模拟滤波器设计和用双线形变换法设计数字滤波器的内容； （2）用双线形变换法设计出数字滤波器的系统函数H(z)； （3）先将数字频率转换成角频率，再将数字角频率转换成模拟角频率； （4）计算模拟低通滤波器传输函数； （5）利用双线性变法将模拟高通滤波器转换成数字滤波器。 五：程序框图 六：源程序 wp=4000*(2*pi)/20000; Wp=2*tan(wp/2)*20000; rp=0.5; ws=5000*(2*pi)/20000; Ws=2*tan(ws/2)*20000; rs=10; ksp=sqrt((10.^0.05-1)/(10.^1-1)); lsp=Ws/Wp; A=abs(-log10(ksp)/log10(lsp)); N=ceil(A) Fs=20000; FN=Fs/2; fc=Wp*(10^(0.1*rp)-1)^(-1/(2*N))/(2*pi); Fc=fc/FN; [b,a]=butter(N,Fc); [z,p,k]=butter(N,Fc); b a z p k subplot(3,1,1) [H,f]=freqz(b,a,512,Fs); plot(f,abs(H)) xlabel(frequcy(Hz)); ylabel(magnitude responsse) subplot(3,1,2) zplane(b,a) subplot(3,1,3) plot(f,20*log10(abs(H))) 七：程序结果及图表 N = 7 b =0.0247 0.1727 0.5180 0.8633 0.8633 0.5180 0.1727 0.0247 a =1.0000 0.5211 1.0260 0.3276 0.2331 0.0391 0.0099 0.0005 z =-1 -1 -1 -1 -1 -1 -1 p = -0.0963 + 0.7929i -0.0963 - 0.7929i -0.0621 + 0.2274i -0.0621 - 0.2274i -0.0726 + 0.4795i -0.0726 - 0.4795i -0.0590 k = 0.0247 八：实验总结 在本题中，我更一步的掌握了matlab，学会了在此软件中调用巴特沃兹函数进行低通滤波器的设计。其中传输函数 H(z)=(0.247+0.1727z-1+0.5180z-2+0.8633z-3+0.8633z-4+0.5180z-5+0.1727z-6+0.0247z-7)/ (1.000+0.5211z-1+1.0260z-2+0.3276z-3+0.2331z-4+0.0391z-5+0.0099z-6+0.0005z-7) 分布图如上图所示。 一：实验题目： 题目五： DTMF信号的生成是通过将两个有限长数字序列相加而实现；而对DTMF信号的检测是通过计算DTMF信号的DFT；然后测量在给定8个频率上的能量而实现。已知采样频率为8KHz，DTMF信号点数N＝205，对DTMF信号进行205点DFT。右表为DTMF数字。 1209Hz 1336Hz 1447Hz 1663Hz 1 697Hz 1 2 3 A 2 770Hz 4 5 6 B 3 852Hz 7 8 9 C 4 941Hz * 0 # D 要求：① 通过键盘任意输入16个键之