信号与系统实验(编程四)new

《信号与系统实验(编程四)new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、大连理工大学实验报告学院（系）：电子信息与电气工程学部专业：信息与通信工程班级：姓名：学号：组：___实验时间：实验室：实验台：指导教师签字：成绩：实验程序和结果：【课程练习1】1、用Matlab求解一下离散时间系统的单位冲激响应。（1）3yn+4yn-1+yn-2=xn+xn-1clccleara=[341];b=[11];a2=[2.5610];n=0:10;impz(b,a,10),gridon;title('系统单位冲击响应h1(n)')（2）2.5yn+6yn-1+10yn-2=xnclccl

2、eara=[2.5610];b=[1];n=0:30;impz(b,a,30),gridon;title('系统单位冲击响应h2(n)'2、已知某系统的单位冲激响应为hn=78nun-un-1，试用matlab求当激励信号为xn=un-un-5时，系统的零状态响应。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'),g

3、ridon;xlabel('n'),ylabel('x(n)');subplot(312)stem(h,'fill'),gridon;xlabel('n'),ylabel('h(n)');subplot(313)stem(y,'fill'),gridon;xlabel('n'),ylabel('y(n)=h(n)*x(n)');【课程练习2】1、试用matlab画出下列因果系统的系统函数零极点分布图，并判断系统的稳定性。（1）Hz=2z2-1.6z-0.9z3-2.5z2+1.96z-0.48（2）Hz

4、=z-1z4-0.9z3-0.65z2+0.873zclcclearB=[02-1.6-0.9];A=[1-0.251.96-0.48];subplot(211)zplane(B,A),gridon;legend('零点','极点')title('零极点分布图(1)')C=[0001-1];D=[1-0.9-0.650.8730];subplot(212)zplane(C,D),gridon;legend('零点','极点')title('零极点分布图(2)')由零极点分布图可看出，（1）中有两个极点在

5、单位元外，所以系统不稳定。（2）中有三个极点，且有两个单阶的极点在单位元上，所以系统零界稳定。3、试用matlab绘制系统Hz=z2z2-34z+18的频率响应曲线。clcclearb=[100];a=[1-3/41/8];[H,w]=freqz(b,a,400,'whole');Hm=abs(H);Hp=angle(H);subplot(211)plot(w,Hm),gridon;xlabel('w'),ylabel('Magnitude')title('离散系统幅频特性曲线')subplot(212

6、)plot(w,Hp),gridon;xlabel('w'),ylabel('Phase')title('离散系统向频特性曲线')3、自行设计系统函数，验证系统函数零极点分布与其时域特性的关系。（1）Hz=zz2-1.2z+0.72，（2）Hz=zz2-1.5z+1，（3）Hz=zz2-2z+1.5clcclearb1=[010];a1=[1-1.20.72];subplot(121)zplane(b1,a1)title('¼«µãÔÚµ¥Î»Ô²ÄÚ')subplot(122)impz(b1,a1,3

7、0);gridon;figureb2=[010];a2=[1-1.51];subplot(121)zplane(b2,a2)title('¼«µãÔÚµ¥Î»Ô²ÉÏ')subplot(122)impz(b2,a2,30);gridon;figureb3=[010];a3=[1-21.5];subplot(121)zplane(b3,a3)title('¼«µãÔÚµ¥Î»Ô²Íâ')subplot(122)impz(b3,a3,30);gridon;figure由图形可看出当极点位于单位圆内时，h(n

8、)为衰减序列，系统是稳定的；当极点位于单位圆上时，h(n)为等幅序列，系统处于零界稳定状态；当极点位于单位圆外时，h(n)为发散序列，此时的系统不稳定。实验体会：通过此次实验，我对离散时间系统的性质有了更好的了解，对Matlab也有了进一步的了解。编程中用到的都是函数，所以对z变换的具体做法不是很了解。