1、大连理工大学实验报告学院（系）： 电子信息与电气工程学部 专业： 信息与通信工程 班级： 姓 名： 学号： 组： _ 实验时间： 实验室： 实验台： 指导教师签字： 成绩： 实验程序和结果：【课程练习 1】1、用Matlab求解一下离散时间系统的单位冲激响应。（1）3 +4 + =y(n) y(n-1) y(n-2) x(n)+x(n-1)clccleara=3 4 1;b=1 1;a2=2.5 6 10;n=0:10;impz(b,a,10),grid on;title(系统单位冲击响应h1(n)（2）2.5 y(n)+6y(n-1)+10y(n-2)=x(n)clccleara=2.5 6
2、 10;b=1;n=0:30;impz(b,a,30),grid on;title(系统单位冲击响应h2(n)2、已知某系统的单位冲激响应为 ，试用matlabh(n)=(78)nu(n)u(n1)求当激励信号为 时，系统的零状态响应。x(n)=u(n)u(n5)clcclearnx=-1:6;nh=-1:11;x=uDT(nx)-uDT(nx-5);h=(7.0/8).nh.*(uDT(nh)-uDT(nh-10);y=conv(x,h);subplot(311)stem(x,fill),grid on;xlabel(n),ylabel(x(n);subplot(312)stem(h,fil
3、l),grid on;xlabel(n),ylabel(h(n);subplot(313)stem(y,fill),grid on;xlabel(n),ylabel(y(n)=h(n)*x(n);【课程练习2】1、试用matlab画出下列因果系统的系统函数零极点分布图，并判断系统的稳定性。（1） H(z)=2z2-1.6z-0.9z3-2.5z2+1.96z-0.48（2） H(z)=z-1z4-0.9z3-0.65z2+0.873zclcclearB=0 2 -1.6 -0.9;A=1 -0.25 1.96 -0.48;subplot(211)zplane(B,A),grid on;lege
4、nd(零点 ,极点)title(零极点分布图(1)C=0 0 0 1 -1;D=1 -0.9 -0.65 0.873 0;subplot(212)zplane(C,D),grid on;legend(零点 ,极点)title(零极点分布图(2)由零极点分布图可看出，（1）中有两个极点在单位元外，所以系统不稳定。（2）中有三个极点，且有两个单阶的极点在单位元上，所以系统零界稳定。3、试用matlab绘制系统 的频率响应曲线。H()= 2234+18clcclearb=1 0 0;a=1 -3/4 1/8;H,w=freqz(b,a,400,whole);Hm=abs(H);Hp=angle(H)
5、;subplot(211)plot(w,Hm),grid on;xlabel(w),ylabel(Magnitude)title(离散系统幅频特性曲线 )subplot(212)plot(w,Hp),grid on;xlabel(w),ylabel(Phase)title(离散系统向频特性曲线 )3、自行设计系统函数，验证系统函数零极点分布与其时域特性的关系。（1） ，（2） ，H()=z21.2+0.72 H()= z21.5+1（3） H()=z22+1.5clcclearb1=0 1 0;a1=1 -1.2 0.72;subplot(121)zplane(b1,a1)title()sub
6、plot(122)impz(b1,a1,30);grid on;figureb2=0 1 0;a2=1 -1.5 1;subplot(121)zplane(b2,a2)title()subplot(122)impz(b2,a2,30);grid on;figureb3=0 1 0;a3=1 -2 1.5;subplot(121)zplane(b3,a3)title()subplot(122)impz(b3,a3,30);grid on;figure由图形可看出当极点位于单位圆内时，h(n)为衰减序列，系统是稳定的；当极点位于单位圆上时，h(n)为等幅序列，系统处于零界稳定状态；当极点位于单位圆外时，h(n)为发散序列，此时的系统不稳定。实验体会：通过此次实验，我对离散时间系统的性质有了更好的了解，对Matlab也有了进一步的了解。编程中用到的都是函数，所以对z变换的具体做法不是很了解。
