信号处理实验报告.docx

1、信号处理实验报告数字信号处理课程实验报告课程名称：任课教师：机 房：计算机编号：实验日期：实验成绩：实验班级：学生XX：实验名称：z变换和零极点的求法实验目的：1、掌握z变换和反z变换的应用2、掌握用局部分式法求z变换3、掌握零极点的求法实验步骤与内容：1、求Z变换。2.求z反变换。3.4.1.将上述的程序进展验证，并进展思考，掌握其编写的方法，将程序和结果复制到下面：1、syms zaw0nx1=an; 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

2、*(n-1); X5=ztrans(x5)X1 =-z/(a - z)X2 =z/(z - 1)2X3 =(z2 + 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、syms nzaX1=z/(z-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 =1 x2 =piecewise(a 0, a*(an/a

3、 - kroneckerDelta(n, 0)/a) + a2*(kroneckerDelta(n, 0)/a2 + (an*(n - 1)/a2) x3 =n + binomial(n - 1, 2) - 1 x4 =iztrans(1/(zn*(1/z - 1), z, n) + 1 3、b=1,0,0;a=1,-1.5,0.5;r p c=residuez(b,a)N=20; n=0: N-1;x=r(1)*p(1).n+r(2)*p(2).nstem(n,x)结果4、b=0.2,0.1,0.3,0.1,0.2a=1,-1.1,1.5,-0.7,0.3rz=roots(b)rp=root

4、s(a)subplot(1,2,1),zplane(b,a);subplot(1,2,2),impz(b,a,20);title();xlabel(n),ylabel(h(n)2.编写程序实现1.求z变换：1 2代码：syms zaw0nx1=n*an; X1=ztrans(x1)x2=sin(w0*n); X2=ztrans(x2)结果：X1 =z/(a*(z/a - 1)2)X2 =(z*sin(w0)/(z2 - 2*cos(w0)*z + 1)2.求z反变换：12syms nzaX1=z/(z-a); x1=iztrans(X1)X2=z/(z-a)2;x2=iztrans(X2)x1

5、 =piecewise(a 0, a*(an/a - kroneckerDelta(n, 0)/a) + kroneckerDelta(n, 0) x2 =piecewise(a 0, a*(kroneckerDelta(n, 0)/a2 + (an*(n - 1)/a2) + an/a - kroneckerDelta(n, 0)/a)3.用局部分式法求解以下系统函数的z反变换，写出xn的表示式，用图形表示出来，取前10个点作图。代码：c=10,20,0,0d=1,8,19,12r,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)结果：4.离散时间系统函数分为：求该系统的零极点分布图，并判断系统的因果稳定性。代码：b=4,-1.6,-1.6,4a=1,0.4,0.35,-0.4rz=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)