信号处理实验资料报告材料.doc

《信号处理实验资料报告材料.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

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

2、)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-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([a<>0,

3、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)结果4、b=[0.2,0.1,0.3,0.1,0.2]a=[1,-1.1,1.5,-0.7,0.3

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)')2.编写程序实现1.求z变换：（1）（2）代码：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)1.求z反变换：（1）（2）

5、symsnzaX1=z/(z-a);x1=iztrans(X1)X2=z/(z-a)^2;x2=iztrans(X2)x1=piecewise([a<>0,a*(a^n/a-kroneckerDelta(n,0)/a)+kroneckerDelta(n,0)])x2=piecewise([a<>0,a*(kroneckerDelta(n,0)/a^2+(a^n*(n-1))/a^2)+a^n/a-kroneckerDelta(n,0)/a])3.用部分分式法求解下列系统函数的z反变换，写出x（n）的表示式，用图形表示出来，取前10个点作图。代码：c=

6、[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)结果：4.已知离散时间系统函数分为：求该系统的零极点分布图，并判断系统的因果稳定性。代码：b=[4,-1.6,-1.6,4]a=[1,0.4,0.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('ϵͳ

7、µÄ³å¼¤ÏìÓ¦');xlabel('n'),ylabel('h(n)')