您好,欢迎访问三七文档
1、塔里木大学信息工程学院1数字信号处理课程实验报告课程名称:任课教师:机房:计算机编号:实验日期:实验成绩:实验班级:学生姓名:实验名称:z变换和零极点的求法实验目的:1、掌握z变换和反z变换的应用2、掌握用部分分式法求z变换3、掌握零极点的求法实验步骤与内容:1、求Z变换。2.求z反变换。塔里木大学信息工程学院23.塔里木大学信息工程学院34.塔里木大学信息工程学院4塔里木大学信息工程学院5塔里木大学信息工程学院61.将上述的程序进行验证,并进行思考,掌握其编写的方法,将程序和结果复制到下面:1、symszaw0nx1=a^n;X1=ztrans(x1)x2=n;X2=ztrans(x2)x3=(n*(n-1))/2;X3=ztrans(x3)x4=exp(j*w0*n);X4=ztrans(x4)x5=1/n*(n-1);X5=ztrans(x5)X1=-z/(a-z)X2=z/(z-1)^2X3=(z^2+z)/(2*(z-1)^3)-z/(2*(z-1)^2)X4=z/(z-exp(i*w0))X5=z/(z-1)-ztrans(1/n,n,z)2、symsnzaX1=z/(z。
2、-1);x1=iztrans(X1)X2=a*z/(a-z)^2;x2=iztrans(X2)X3=z/(z-1)^3;x3=iztrans(X3)X4=(1-z^-n)/(1-z^-1);x4=iztrans(X4)x1=1x2=piecewise([a0,a*(a^n/a-kroneckerDelta(n,0)/a)+a^2*(kroneckerDelta(n,0)/a^2+(a^n*(n-1))/a^2)])x3=n+binomial(n-1,2)-1x4=iztrans(1/(z^n*(1/z-1)),z,n)+13、b=[1,0,0];a=[1,-1.5,0.5];[rpc]=residuez(b,a)N=20;n=0:N-1;x=r(1)*p(1).^n+r(2)*p(2).^nstem(n,x)塔里木大学信息工程学院7结果0246810121416182000.20.40.60.811.21.41.61.824、b=[0.2,0.1,0.3,0.1,0.2]a=[1,-1.1,1.5,-0.7,0.3]rz=roots(b)rp=roots(a)subplot(1,2,1。
3、),zplane(b,a);subplot(1,2,2),impz(b,a,20);title('ϵͳµÄ³å¼¤ÏìÓ¦');xlabel('n'),ylabel('h(n)')塔里木大学信息工程学院8-1-0.500.51-2-1.5-1-0.500.511.52RealPartImaginaryPart051015-0.1-0.0500.050.10.150.20.250.30.350.4nh(n)系统的冲激响应2.编写程序实现1.求z变换:(1)nnanx)(1(2))sin()(02nnx代码:symszaw0nx1=n*a^n;X1=ztrans(x1)x2=sin(w0*n);X2=ztrans(x2)结果:X1=z/(a*(z/a-1)^2)X2=(z*sin(w0))/(z^2-2*cos(w0)*z+1)2.求z反变换:(1)azzzX)(1(2)22)()(azzzXsymsnzaX1=z/(z-a);x1=iztrans(X1)X2=z/(z-a)^2;x2=iztrans(X2)塔里木大学信息工程学院9x1=piecewise([a0,a*(。
4、a^n/a-kroneckerDelta(n,0)/a)+kroneckerDelta(n,0)])x2=piecewise([a0,a*(kroneckerDelta(n,0)/a^2+(a^n*(n-1))/a^2)+a^n/a-kroneckerDelta(n,0)/a])3.用部分分式法求解下列系统函数的z反变换,写出x(n)的表示式,用图形表示出来,取前10个点作图。32111219812010)(zzzzzX代码:c=[10,20,0,0]d=[1,8,19,12][r,p,c]=residuez(c,d)N1=10;n=0:N1-1;x=r(1)*p(1).^n+r(2)*p(2).^n+r(3)*p(3).^nstem(n,x)结果:0123456789-7-6-5-4-3-2-1012x1064.已知离散时间系统函数分为:3213214.035.04.0146.16.14)(zzzzzzzH求该系统的零极点分布图,并判断系统的因果稳定性。塔里木大学信息工程学院10代码:b=[4,-1.6,-1.6,4]a=[1,0.4,0.。
5、35,-0.4]rz=roots(b)rp=roots(a)subplot(1,2,1),zplane(b,a);subplot(1,2,2),impz(b,a,20);title('ϵͳµÄ³å¼¤ÏìÓ¦');xlabel('n'),ylabel('h(n)')-1-0.500.51-2-1.5-1-0.500.511.52RealPartImaginaryPart051015-4-202468nh(n)系统的冲激响应。
本文标题:信号处理实验报告
链接地址:https://www.777doc.com/doc-6990110 .html