数字信号处理实验(吴镇扬版)matlab程序.docx

《数字信号处理实验(吴镇扬版)matlab程序.docx》由会员分享，可在线阅读，更多相关《数字信号处理实验(吴镇扬版)matlab程序.docx（31页珍藏版）》请在课桌文档上搜索。

1、(1)数组的加、减、乘、除和乘方运算。输入A=1234,B=3456,求C=A+B,D=A-B,E=A.*B,F=A./B,G=AJB并用Stem语句画出A、B、C、D、E、F、G。clearall;a=l234;b=3456;c=a+b;d=a-b;e=a.*b;f=a./b:g=a.b;n=l:4;subplot(4,2,1);stem(n,a);xlabelCn,);xlim(O5);ylabelCA,);subplot(4,2,2);Stem(I,b);xlabelCn);Xlim(O5);ylabel(B);subplot(4,2,3);stem(n,c);xlabelCn,);xl

2、im(05);ylabel(C);subplot(4,2,4);stem(n,d);xlabelCn,);xlim(O5);ylabel(D);subplot(4,2,5);Stem(I,e);xlabelCn);XIim(05);ylabelCE,);subplot(4,2,6);stem(n,f);xlabelCn,);xlim(05);ylabelCF,);subplot(4,2,7);stem(n,g);xlabelCn,);xlim(O5);ylabel(G);(2)用MATLAB实现下列序列:a) x(n)=0.8nOWnWI5b) x(n)=e(0.2+3j)n0n15c) x(

3、n)=3cos(0.125n+0.2)+2sin(0.25n+0.1)0n15d)将c)中的x(n)扩展为以16为周期的函数xl6(n)=x(n+16),绘出四个周期。e)将c)中的x(n)扩展为以10为周期的函数xl0(n)=x(n+10),绘出四个周期。clearall;N=0:15;xa=0.8.N;figure;subplot(2,1,1);stem(N,xa);xlabel(,n,)jxlim(O16);ylabel(,xa,);xb=exp(0.2+3*j)*N);subplot(2,1,2);stem(N,xb);xlabel(*n,);xlim(016);ylabel(,xb,

4、);figure;xc=3*cos(0.125*pi*N+0.2*pi)+2*sin(0.25*pi*N+0.l*pi);subplot(3,1,1);stem(N,xc);xlabel(,n*);xlim(016);ylabelCxc,);k=0:3;m=0;fori=l:4forj=l:16m=m+1;n(m)=N(j)+16*k;xl6(m)=3*cos(0.125*pi*n(m)+0.2*pi)+2*sin(0.25*pi*n(m)+0.l*pi);endendsubplot(3,1,2);stem(n,xl6);xlabel(,n,);ylabelCxl6,);forj=l:10xl

5、(j)=xl6(j);endfori=l:3form=l:10xl(i*l+m)=xl(m);endendn=l:40;subplot(3,1,3);stem(n,xl);xlabel(,n,);ylabelCxl,);TTTYY70QQC一2468 IO 121416246IO121416(3) x(n)=l,-l,3,5,产生并绘出下列序列的样本:a) Xi(n)=2x(n+2)-(n-l)-2x(n)5b) x2(n)=nx(n-k)k=lclearalln=l:4;T=4;x=l-135;x(5:8)=x(l:4);subplot(2,1,1);stem(l:8,x);grid;for

6、i=l:4ifi-l0xl(i)=2*x(i+2)-(i-l)-2*x(i);elsex1(i)=2*x(i+2)-x(i-1+T)-2*x(i);endendxl(5:8)=xl(1:4);subplot(2,1,2);stem(l:8,xl);grid;y轴以及图形上方均须加上适当的标注:(4)绘出下列时间函数的图形，对X轴、a)X(t)=sin(2t)0t10sb)X(t)=Cos(100t)sin(t)0t4sta=O:O.05:10;xa=sin(2*pi*ta);subplot(2,1,1);plot(ta,xa);xlabelCt,);ylabel(幅度);tb=O:O.01:4

7、;xb=cos(100*pi*tb).*sin(pi*tb);subplot(2,1,2);plot(tb,xb);xlabelCt,);ylabel(幅度);(5)编写函数StePShift(n,nl,n2)实现u(n-n),nln=O;stem(n,x);(6)给一定因果系统H(Z)=(1+2z-+1)/(1-O67z-+0.9Z=)求出并绘制H(Z)的幅频响应与相频响应。clearall;b=l,sqrt(2),1;a=l,-O.67,O.9;h,w=freqz(b,a);am=20*logl0(abs(h);subplot(2,1,1);plot(w,am);ph=angle(h);s

8、ubplot(2,1,2);plot(w,ph);（7）计算序列8-2-123和序列23-1-3）的离散卷积，并作图表示卷积结果。(1)、观察高斯序列的时域和幅频特性，固定信号xa(n)中参数p=8,改变q的值，使q分别等于2,4,8,观察它们的时域和幅频特性，了解当q取不同值时,对信号序列的时域幅频特性的影响；固定q=8,改变p，使P分别等于8,13,14,观察参数P变化对信号序列的时域及幅频特性的影响，观察P等于多少时，会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。解：程序见附录程序一：n=0:l:15;%p=8不变，q变化(2,4,8

9、);p=8;q=2;%p=8;q=2;xa1=exp(-(n-p).2)q);subplot(5,2,l);plot(n,xal,xlabel(,tT);ylabel(,xa(n);title(,p=8q=2,)xkl=abs(fft(xal);subplot(5,2,2);stem(n,xkl)xlabel(k);ylabel(,Xa(k),);title(,p=8q=2,)p=8;q=4;%p=8;q=4;xa1=exp(-(n-p).2)q);subplot(5,2,3);plot(n,xal,xlabel(,tT);ylabel(,xa(n)1);title(,p=8q=4,)xkl=

10、abs(fft(xal);subplot(5,2,4);stem(n,xkl)xlabel(,k,);ylabel(,Xa(k),);title(,p=8q=4)p=8;q=8;%p=8;q=8;xa1=exp(-(n-p).2)q);subplot(5,2,5);plot(n,xal,xlabel(tT);ylabel(,xa(n);xk1=abs(fft(xal);title(,p=8q=8,)subplot(5,2,6);stem(n,xkl)xlabel(,k);ylabel(Xa(k)1);title(,p=8q=8,)%q=8不变,p变化(8,13,14);p=8;q=8;%p=8

11、;q=8;xa1=exp(-(n-p).2)q);subplot(5,2,5);plot(n,xal,xlabel(,tT);ylabel(,xa(n)1);xkl=abs(fft(xal);title(,p=8q=8)subplot(5,2,6);stem(n,xkl)xlabel(,k,);ylabel(,Xa(k)1);title(,p=8q=8,)p=13；q=8；%p=13;q=8;xa1=exp(-(n-p).2)q);subplot(5,2,7);plot(n,xal,xlabel(,tT);ylabel(,xa(n),);xkl=abs(fft(xal);title(,p=13

12、q=8)subplot(5,2,8);stem(n,xkl)xlabel(,k);ylabel(Xa(k);title(,p=13q=8,)p=14;q=8;%p=14;q=8;xa1=exp(-(n-p).2)q);subplot(5,2,9);plot(n,xaI,xlabel(,tT);ylabel(,xa(n);title(,p=14q=8)xkl=abs(fft(xal);subplot(5,2J0);stem(n,xkl)xlabel(,k,);ylabel(,Xa(k);titlefp=14q=8,)幅频特性Q(?9qqqqq9?p=14 q=8t/Tk分析：由高斯序列表达式知n

13、=p为期对称轴；当P取固定值时，时域图都关于n=8对称截取长度为周期的整数倍，没有发生明显的泄漏现象；但存在混叠，当q由2增加至8过程中，时域图形变化越来越平缓，中间包络越来越大，可能函数周期开始增加，频率降低，渐渐小于fs2,混叠减弱；当q.固定不变，P变化时，时域对称中轴右移，截取的时域长度渐渐地不再是周期的整数倍，开始无法代表一个周期，泄漏现象也来越明显，因而图形越来越偏离真实值，p=14时的泄漏现象最为明显，混叠可能也随之出现；(2)、观察衰减正弦序列xb(n)的时域和幅频特性，a=0.1,f=0.0625,检查谱峰出现位置是否正确，注意频谱的形状，绘出幅频特性曲线，改变f,使f分别等

14、于0.4375和0.5625,观察这两种情况下，频谱的形状和谱峰出现位置，有无混叠和泄漏现象？说明产生现象的原因。et,nsin(2/?),0n5其他nl=0:l:15;xb1=exp(-0.1*n1).*sin(2*pi*0.0625*n1);subplot(3,2,l);plot(nl,xbl,-*);xlabel(,n);ylabel(,x(n)t);title(1f=0.0625);xkl=abs(fft(xb1);subplot(3,2,2);stem(nl,xkl)xlabel(,k);ylabel(X(k);title(1f=0.0625,);n2=0:1:15;xb2=exp(

15、-0.1*n2).*sin(2*pi*0.4375*n2);subplot(3,2,3);plot(n2,xb2,xlabel(,n,);ylabel(x(n);title(,f=0.4375,);xk2=abs(fft(xb2);subplot(3,2,4);stem(n2,xk2)xlabel(,k);ylabel(,X(k)1);title(,f=0.4375,);n3=0:1:15;xb3=exp(-0.1*n3).*sin(2*pi*0.5625*n3);subplot(3,2,5);plot(n3,xb3,-*);xlabel(,n);ylabel(,x(n)t);title(,f

16、=0.5625,);xk3=abs(fft(xb3);subplot(3,2,6);stem(n3,xk3)xlabel(,k);ylabel(,X(k)1);title(,f=0.5625,);时域特性幅频特性f=0.0625f=0.0625分析：当f=fl=0.0625时，谱峰位置出现正确，存在在混叠现象，时域采样为一周期，不满足采样定理。当40.4375和0.5625时，时域图像关于Y轴对称，频域完全相同。这是因为频域图是取绝对值的结果,所以完全相同。另外由于时域采样为6个半周期,满足采样定理，无混叠；但由于截取长度不是周期整数倍，出现泄漏。(3)、观察三角波和反三角波序列的时域和幅频特

17、性，用N=8点FFT分析信号序列xc(n)和xd(n)的幅频特性，观察两者的序列形状和频谱曲线有什么异同？绘出两序列及其幅频特性曲线。在xc(n)和xd(n)末尾补零，用N=16点FFT分析这两个信号的幅频特性，观察幅频特性发生了什么变化？两情况的FFT频谱还有相同之处吗？这些变化说明了什么？三角波序列：nn3xc(n)=-8-n,4nl0,其他,0 3 ,4 7 ,其他反三角波序列：4-nxc(n)=n-40clearall;n=0:3;k=l:8;%定义三角波序列Xc(n+1)=n;Xc(n+5)=4-n;%定义反三角波序列Xd(n+1)=4-n;Xd(n+5)=n;%三角波特性%subp

18、lot(2,2,1);plot(k-l,Xc);XlabelCnXylabel(时域特性);IeXKl,3；三角波);subplot(2,2,2);plot(k-ltabs(ffi(,Xc);Xlabdm;ylabel(幅频特性)exl(4,l:三角波);%反三角波特性%subplot(2,2,3);PlOKk-1,Xd);XlabeHn);ykibe(时域特性);iexl(3,3；反三角波上subplot(2,2,4);plot(k-l,ahs(fft(Xd);*1疝1(%);)，1而1(幅频特性)；1凶(4,10，反三角波*%末尾补O,计算32点FFTXc(9:32)=0;Xd(9:32)

19、=0;k=l:32;figure;%三角波特性%subplot(2,2,1);plot(k-l,Xc);XlabelCnXylabel(时域特性);IeXKl,3；三角波);subplot(2,2,2);plot(k-ltabs(ffi(,Xc);XIabem;ylabel(幅频特性exl(4,10,三角波);%fan%subplot(2,2,3);PlOKk-1,Xd);XIabelen)Nabel(时域特性);IeXi(3,3；反三角波*subplot(2,2,4);plot(k-l,ahs(fft(Xd);*1疝1(%);)，1而1(幅频特性)；1凶(4,10，反三角波*N=8时域和幅度

20、频谱图：时域特性频域特性分析：由图知，三角波序列和反三角波序列的时域图像成镜像关系，但频域图像完全一样，只是因为幅频图是对X(k)的值取绝对值。N=32时域和幅度频谱图：时域特性频域特性分析：由实验所得的图形知，N=32点时/()和/5)的幅频特性都更加密集，更多离散点的幅值显示，“栅栏效应”减小，分辨率提高，而对于/()来说变化更加明显。在原序列的末端填补零值，变动了DFT的点数，人为的改变了对真实频谱采样的点数和位置，相当于搬动了“尖桩栅栏”的位置，从而使得频谱的峰点和谷点暴露出来。N=32时，/()和/()的频谱差别较大，但总体趋势仍然都是中间最小，两侧呈对称。(4)、一个连续信号含两个

21、频率分量，经采样得x(n)=sin2*0.125n+cos2*(0.125+f)nn=0,l,N-I已知N=16,Af分别为1/16和1/64,观察其频谱;当N=128时,Af不变，其结果有何不同，为什么？clearall;%N=16%N=I6;delf=l/16;n=/0:N-1;xl(n+l)=sin(2*pi*0.125.*n)+cos(2*pi*(0.125+det.*n);detf=l/64;x2(n+l)=sin(2*pi*0.125.*n)+cos(2*pi*(0.125+detf).*n);%N=16,detf=1/16%subplot(2,2,1);stem(n,xl);ho

22、ld;plot(n,xl);XlabelCnXyIabeI(时域特性);lext(6,l：N=16,deif=l/16Xsbplot(2,2,2);stem(n,abs(fft(xl);xlabel(,n,):ylabel(幅值特性XteXI(6,4,N=16,detfr116,);%N=16,detf=1/64%sbplot(2,2,3);stem(n,x2);xlabel(,n);ylabel(，时域特性XIeXH6,1,N=16,detf=1/64,);SMbPIol(2,2,4);stem(n,abs(ffi(x2);XlabeICnXylabelC幅值特性*exl(6,4,N=I6,

23、detf=1/64%N=128%N=128;detf=l/16;n=0:N-1;x3(n+l)=sin(2*pi*O.I25.*n)+cos(2*pi*(0.125+detJ).*n);detf=1/64;x4(n+1)=sin(2*pi*0.125.*n)+cos(2*pi*(0.125+det.n);%N=128,detf=1/16%figure;subplot(2,2,1);stem(n,x3);XlabelcnxyIabeK时域特性XaXiS(Io128-22);text(6,L5,N=28,detf=l16);sbplot(2,2,2);stem(n,abs(fft(x3);Xlab

24、eKn)Nabel(幅值特性XaXiSao128-IO70);text(40,60,N=128,detf=1/16,);%N=128,detf=1/64%subplot(2,2,3);stefn(n,x3);XlabelnXylabel(时域特性XaXiSao128-22);text(6,L5,N=128tdetf=l16,);subplot(2,2,4);stem(n,abs(fft(x4);XlabeKn)NabeIe幅值特性XaXiSuO128-1070);text(40,60,N=I28,detf=116);510152 10N=16,det1B451015IR一rIqrsN=128l

25、det1116120 Oo 80 6040 2020406080100120CE)N=128,detf116120 g80604020100120(5)、用FFT分别实现xa(n)(p=8,q=2)和xb(n)(a=0.1,f=0.0625)的16点圆周卷积和线性卷积。-O51015clearall;N=16;n=O:N-l;P=8；q=2;Xa(n+l)=exp(-(np).2.q);a=0.1;f=0.0625;Xb(n+l)=exp(-a.*n).*sin(2*pi*f.*n);%16点循环卷积Fa=fft(Xa);Fb=fft(Xb);Fx=Fa.*Fb;X51=ifft(Fx);st

26、em(n,X51);%16点线性卷积Xa(N+k2*N-l)=0;Xb(N+l:2*N-l)=0;Fa=fft(Xa);Fb=fft(Xb);Fc=Fa.*Fb;X52=ifft(Fc);figure;stem(1:2*N-1,X52);(7)用FFT分别计算xa(n)(p=8,q=2)和Xb(n)(a=O.1,f=O.0625)的16点循环相关和线性相关，问一共有多少种结果，他们之间有何异同点。clearall;N=16;n=0:N-l;P=8；q=2;Xa(n+l)=exp(-(np).2.q);a=0.l;f=0.0625;Xb(n+l)=exp(-a.*n).*sin(2*pi*f.*

27、n);N=Iength(Xa);%16点循环相关Fa=fft(Xa,2*N);Fb=fft(Xb,2*N);Fx=conj(Fa).*Fb;X71=real(ifft(Fx);X71=X71(N+2:2*N)X71(1:N);n=(-N+l):(N-I);stem(n,X71);%16点线性相关Xa(N+l:2*N-l)=0;Xb(N+1:2*N-I)=0;Fa=fft(Xa);Fb=fft(Xb);Fc=conj(Fa).*Fb;X72=real(ifft(Fc);figure;stem(1:2*N-1,X72);(8)用FFT分别计算xa(n)(p=8,q=2)和xb(n)(a=0.1,f

28、=0.0625)的自相关函数。clearall;N=16;n=O:N-l;P=8;q=2;Xa(n+l)=exp(-(n-p).t.q);a=0.1;f=0.0625;Xb(n+l)=exp(-a.*n).*sin(2*pi*f.*n);%自然对数的底：e=:2.71828182845904523536N=Iength(Xa);%Xa(n)16点自相关Fa=fft(Xa,2*N);Fb=fft(Xb,2*N);Fl=conj(Fa).*Fa;X81=real(ifft(Fl);X81=X81(N+2:2*N)X81(1:N);n=(-N+l):(N-I);subplot(2,1,1);stem

29、(n,X81);xlabel(,n,);ylabcl(幅度)；%Xb(n)16点自相关Fb=fft(Xb,2*N);F2=conj(Fb).*Fb;X82=real(ifft(F2);X82=X82(N+2:2*N)X82(1:N);%n=(-N+l):(N-I);subplot(2,1,2);stem(n,X82);(1)fc=0.3kHz,b=0.8dB,fr=0.2kHz,A，=20dB,T=L3；设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。解：程序：clear;fc=300;fr=200;fs=1000;rp=0.8;rs=20;wc=2*fs*tan(2*pi*f

30、c(2*fs);wt=2*fs*tan(2*pi*fr(2*fs);N,wn=cheb1Ord(WC,wt,rp,rs,s);B,A=cheby1(N,rp,wn,high,s,);bz,az=bilinear(B,A,fs);h,w=freqz(bz,az);f=w*fs(2*pi);plot(f,20*log 10(abs(h);axis(0,fs2,-80,10);grid;Xlabel(频率/Hz);ylabel(幅度 dBt);100 10 20国0-40 50607080050100150200250300350400450500Hz分析：f=200Hz时阻带衰减大于3OdB,通过

31、修改axis(0,fs2,-80,10)axis(200,fs2,-l,l)发现通带波动rs满足0.8。bz=0.0262-0.10470.1570-0.10470.0262az=1.00001.52891.65370.94520.2796系统函数为：_0.0262OOWzT+0.1570z-2-0.1Q47Z-+0.0262z41+1.5289z,+1.6537z2+0.9452z30.2796z4(2) fc=0.2kHz, b=ldB, L=O3kHz, Ar = 25d&T = Ims;分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器，观察所设计数字滤波器的幅频特性曲 线

32、，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。解：程序：clear;fs= 1000;fc=200;fr=300;rp= 1 ;rs=25;%脉冲响应不变法wp=2*pi*fc;ws=2*pi*fr;N, wn = buttord(wp, ws, rp, rs, s);bl a 1 =butter(N,wn,s,);bz 1 ,az 1 =impinvar(b l,al,fs);h I ,w=freqz(bz 1 ,az 1);%双线性变换法wp=2*fs*tan(2*pi*fcfs2);ws=2*fs*tan(2*pi*frfc2);N, wn = buttord(wp, w

33、s, rp, rs, s);b2 a2=butter(N,wn,s,);bz2,az2=bilinear(b2,a2,fs);Ih2,w=freqz(bz2,az2);f=w(2*pi)*fs;figure; plot(f,abs(h 1 ),-.r,f,abs(h2),-b);grid; XlabeIC频率/Hz); ylabe(幅度)； Iegendc脉冲响应不变法；双线性变换法 title，巴特沃思低通滤波器，线性幅度谱上1.4巴特沃思低通滤波器，线性幅度谱bzl=0.00000.00020.01530.09950.14440.06110.00750.00020.00000azI=LOO

34、OO-1.91992.5324-2.20531.3868-0.63090.2045-0.04500.0060-0.0004因此脉冲响应不变法的系统函数为：H0.0002+0.0153z-2+0.0995z7+().1444-0.061Iz7-0.00752-0.00022-，=l-l.9199z1+2,5324z2-2.2053zs+1.3868z4-0.6309-5+0.2()45z*,-0.()4507+0.(X)60h-0.(XX)4zybz2=0.01790.10720.26810.35750.26810.10720.0179az2=1.0000-0.60190.9130-0.2989

35、0.1501-0.02080.0025因此双线性变换法的系统函数为：Ool79+09072zT+0.2681z-2+0.3575zT+0.2681z-+0.1072zT+0.0179z6hilZ-l-0.6019z-,+0.9130z-2-0.2989z3+0.1501z4-0.0208z5+0.0025z-6分析：脉冲响应不变法的N=9,双线性变换法的N=6,由图知它们都满足要求，但脉冲响应的衰减较快，双线性变换的过渡带窄一些，且阶数比脉冲小，容易实现。(3)利用双线性变换法分别设计满足下列指标的巴特沃思型、切比雪夫型和椭圆型数字低通滤波器，并作图验证设计结果：fc=2kHz,0.5dB,f

36、r=2kHz,f=w(2*pi)*fs;At40dB,fs = 8kHz .解：程序：clear;fs=8000;fc= 1200;fr=2000;rp=0.5;rs=40;%巴特沃思低通滤波器wp=2*fs*tan(2*pi*fcfs2);ws=2*fs*tan(2*pi *frfs2);IN, wn = buttord(wp, ws, rp, rs, s);bl al=butter(N,wn,s,);bzl,az 1 =bilinear(bl,al,fs);h I ,w=freqz(bz 1 ,az 1);Hl=20*logl0(abs(hl);figure;plot(f,Hl);%对数幅

37、度谱axis(0,fs2,-l00,10);grid;XIabele频率/Hz);ylabel(S);titleC巴特沃思低通滤波器，对数幅度谱)；%切比雪夫低通滤波器wc=2*fs*tan(23*cpi*fc(2*fs);wt=2*fs*tan(2*pi*fr(2*fs);IN,wn=cheb1ord(wc,wt,rp,rs,s,);b2,a2=chebyI(N,rp,wnJOW,s);bz2,az2=bilinear(b2,a2,fs);h2,w=freqz(bz2,az2);H2=20*logl0(abs(h2);f=w*fs(2*pi);figure;plot(f,H2);axis(0

38、,fs2,-100,10);grid;XIabelC频率/Hz);ylabelC幅度dB);IiUeC切比雪夫低通滤波器，对数幅度谱)；%椭圆型数字低通滤波器wp=2*fs*tan(2*pi*fcfs2);%双线性变换法ws=2*fc*tan(2*pi*frfs2);lN,wpJ=ellipord(wp,ws,rp,rs,s,);b3,a3=ellip(N,rp,rs,wplow,s,);bz3,az3=bilinear(b3,a3,fs);h3,w=freqz(bz3,az3);H3=20*logl0(abs(h3);f=w(2*pi)*f;figure;PIOt(f,H3);axis(0,

39、fs2!00,10);grid;XIabe1(频率/Hz);ylabel(幅度dB);title。椭圆型数字低通滤波器，对数幅度谱)；巴特沃思低通滤波器，对数幅度谱频率/Hzbzl=0.00040.00320.01290.03020.04530.04530.03020.01290.00320.0004azl=1.0000-2.79964.4582-4.54123.2404-1.63300.5780-0.13700.0197-0.0013系统函数为：O.OOO4+O.OO32zJ+O.OI29z-0.0302z-3+0.0453z7+().0453?-+0.0302z*+0.0129z=+0.0

40、032Z-X+().0(X)4二Z11-2.7996,+4.4582z-2-4.5412z3+3.24O4z4-1.6330z-s+0.5680z-6-0.1370-7+0.0197z-8-0.0013z-分析：N=9,为九阶巴特沃思低通滤波器，从图中可以看出通带波动和阻带衰减都满足设计要求。切比雪夫低通滤波器:p三三Ooo-1000X切比雪夫低通滤波器，对数幅度谱OOO-560OO-8-9bz2=0.0026az2=1.0000系统函数为:5001000150020002500频率/Hz30000.01320.02640.02640.01320.0026-2.97754.2932-3.512

41、41.6145-0.333435004000H2()=0.0026+0.0132ZT+0.0264z”+0.02642点+0.0132ZT+0.0026z51-2.9775Z-I+4,2932z2-3.5124z31.6154ZT-0.3334z5分析：N=5,为五阶切比雪夫低通滤波器，从图中可以看出通带波动和阻带衰减都满足设计要求。椭圆型数字低通滤波器:10-100 0/5001000150020002500300035004000频率/Hz椭圆型数字低通滤波器，对数幅度谱0.06650.03630.0389J2.3658 -1.32500.3332bz=10.03890.0363az=1.0000-2.1444系统函数为：H3(z)=0.0389-0.0363z - 2. 1444z-l + 2. 3658z-2 - 1. 3250z-3 + O. 3332z-4分析：N=4,为四阶椭圆型数字低通滤波器，从图中可以看出通带波动和阻带衰减都满足设计要求。 (4)分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字带通滤波器，已知fs = 30kHz ,其等效的模拟滤波器指标为Sv3dB, 2kHzf3kHz, At5dB, f 6kHz, At20dB, f.5kHzo+O.0665z2-O.0363z3+O.0389z