信号处理实验报告

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)()(azzzXsymsnzaX1=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)系统的冲激响应。