《信号与系统实验之连续线性时不变系统的分析.docx》由会员分享，可在线阅读，更多相关《信号与系统实验之连续线性时不变系统的分析.docx（16页珍藏版）》请在三一文库上搜索。
1、信号与系统实验报告连续线性时不变系统的分析专 业：电子信息工程（实验班）姓 名：曾雄学 号： 14122222203班 级：电实 12-1BF目录验原理与目3二、实 验 过 程 及 结 果 测试 3三、思考题 10四、实验总结 10五、参考文献 11一、实验原理与目的深刻理解连续时间系统的系统函数在分析连续系统的时域特性、 频域特性及稳定性中的重要作用及意义。掌握利用MATLA分析连续系统的时域响应、频响特性和零极点的基本方法。二、实验过程及结果测试1描述某线性时不变系统的微分方程为：y(t) 3y(t) 2y(t)f(t) 2f(t)且 f(t)=t 2， y(0-)=1 ， y (0-)=
2、1 ；试求系统的单位冲激响应、单位阶跃响应、全响应、零状态响应、零输入响应、自由响应和强迫响应。编写相应MATLAB?序，画出各波形图。( 1)单位冲激响应：程序如下：%求单位冲激响应a=1,3,2;b=1,2;sys=tf(b,a);t=0：0.01：10；h=impulse(sys,t);卿画图函数plot()画单位冲激响应的波形plot(h);%单位冲激响应曲线xlabel(t);ylabel(h);title( 单位冲激响应h(t)程序运行所得波形如图一：1 单位冲激响应h( t )IULo. 9_p. 8_q. 7.%r. 6h o.5_0.4.0.3-0.20.1_0 LEEEE0
3、20040060080010001200t图一单位冲激响应的波形(2)单位阶跃响应：程序如下：%t单位阶跃响应a=1,3,2;b=1,2;sys=tf(b,a);t=0:0.01:10;G=step(sys,t);卿画图函数plot()画单位阶跃响应的波形plot(G);%单位阶跃响应曲线xlabel(t);ylabel(g);title(单位阶跃响应g(t) 程序运行所得波形如图二:0.90.80.70.6g 0.50.40.30.20.10g( t )单位阶跃响应020040060080010001200t图二单位阶跃响应的波形(3)零状态响应：程序如下：%t零状态响应yzs=dsolve
4、( D2y+3*Dy+2*y=2*t+2*tV, y(0)=0,Dy(0)=0)卿符号画图函数ezplot() 画各种响应的波形t=0:0.01:3;ezplot(yzs,t);%?状态响应曲线axis(0,3,-1 5);title(零状态响应曲线yzs);ylabel( yzs);程序运行所得波形如图三：t图三零状态响应的波形(4)零输入响应：程序如下：%t零输入响应yzi=dsolve( D2y+3*Dy+2*y=0 , y(0)=1,Dy(0)=1)卿符号画图函数ezplot() 画零输入响应的波形t=0:0.01:3;ezplot(yzi,t);%B输入响应曲线axis(0,3,-1
5、,2);title(零输入响应yzi);ylabel( yzi);程序运行所得波形如图四:21.510.50-0.5-1t图四零输入响应的波形(5)全响应:程序如下:%t全响应y=dsolve( D2y+3*Dy+2*y=2*t+2*tA2,y(0)=1,Dy(0)=1)卿符号画图函数ezplot() 画全响应响应的波形 t=0:0.01:3;ezplot(y,t); 麻响应曲线axis(0,3,-1,5);title(全响应 y);ylabel( y);程序运行所得波形如图五:54310t-1 0y 2图五全响应的波形(6)自由响应:程序如下: 痂由响应y=dsolve( D2y+3*Dy+
6、2*y=2*t+2*P2,y(0)=1,Dy(0)=1); 必响应yht=dsolve( D2y+3*Dy+2*y=0 ,y(0)=1,Dy(0)=1);% 求齐次通解yt=dsolve( D2y+3*Dy+2*y=2*t+2*tA2, y(0)=0,Dy(0)=0); % 求非齐次通解yp=yt-yht;yh=y-yp; %求齐次解，即自由响应t=0:0.01:3;ezplot(yh,t);title(自由响应yh);ylabel( yh);程序运行所得波形如图六：t图六自由响应的波形(7)强迫响应:程序如下: %虽迫响应yht=dsolve( D2y+3*Dy+2*y=0 ,y(0)=1,
7、Dy(0)=1);% 求齐次通解yt=dsolve( D2y+3*Dy+2*y=2*t+2*tA2, y(0)=0,Dy(0)=0); % 求非齐次通解yp=yt-yht; %求特解，即强迫响应 t=0:0.01:3; ezplot(yp,t);title(强迫响应yp); ylabel( yp);程序运行所得波形如图七：强 迫 响 应 yp5EccF43.p2-1 _:0_-1 _ r 一_ ._ - 一00.511.522.53t图七强迫响应的波形2 .给定一个连续线性时不变系统，描述其输入输出之间关系的微分方程为:编写MATLAB?序，绘制系统的幅频响应、相频响应、频率响应的实部和频率响
8、 应的虚部的波形，确定滤波器的类型。(1)幅频响应：程序如下：%昌频响应b=1；a=1,3,2;H,w=freqs(b,a);plot(w,abs(H);title(幅频响应);xlabel(频率(rad/s);y1abel(幅度)；y(t) 3y(t) 2y(t) x(t)幅 频 响 应程序运行所得波形如图八：0.50.450.40.350.30.20.150.10.050012345678910频 率 (rad/s)图八幅频响应的波形(2)相频响应：程序如下：%1频响应b=1；a=1,3,2;H,w=freqs(b,a);phai=angle(H);plot(w,phai);%1出相频响应
9、title(相频响应);xlabel( 3 (rad/s);ylabel( ();程序运行所得波形如图九：相频响应0r:c1r匚Cr-0.5_-1-4-1.5-2-2.5_一-J-3EEELEECCE:0123456789103 (rad/s)图九相频响应的波形(3)频率响应的实部：程序如下：哪率响应的实部的波形b=1;a=1,3,2;H,w=freqs(b,a);Hr=real(H);plot(w,Hr);%画出频率响应的实部的波形title(频率响应的实部)；xlabel( a (rad/s);ylabel( r);程序运行所得波形如图十：0.6频 率 响 应 的 实 部0.50.40.3
10、0.20.10-0.1s (rad/s)图十频率响应的实部的波形(4)频率响应的虚部：程序如下：哪率响应的虚部的波形b=1；a=1,3,2;H,w=freqs(b,a);Hi=imag(H);plot(w,Hi); 喇出频率响应的虚部的波形title(频率响应的虚部)；xlabel( a (rad/s);ylabel( i);程序运行所得波形如图十一：-0.05-0.1-0.15-0.2-0.25-0.3-0.35频 率 响 应 的 虚 部s (rad/s)图H一 频率响应的虚部的波形(5)答：由系统的幅频响应的波形可知该滤波器的类型为低通滤波器3.已知系统函数为H(s)Y 0.8 ,试用MA
11、TLAIB出系统的零极点s3 2s2 2s 1分布图，冲激响应波形、阶跃响应波形、幅频响应曲线和相频响应曲线，并判断 系统的稳定性。(1)冲激响应：程序如下：%中激响应b=1,-2,0.8;a=1,2,2,1;sys=tf(b,a);t=0:0.01:10;H=impulse(sys,t);plot(H);xlabel( t);ylabel( H);title(冲激响应 H(t);程序运行所得波形如图十二：图十二冲激响应的波形(2)阶跃响应:程序如下：t喻跃响应b=1,-2,0.8； a=122j;sys=tf(b,a)6;0.4t=0:0LO1:10;H 0.2G=step(sys,t)0-
12、0.2t4)-0G)0title(阶跃响应 G(t);plot(G); xlabel( ylabel(程序运行所得波形如图十三:1阶 跃 响 应 G(t)Gt图十三阶跃响应的波形(3)零极点分布图：程序如下：极点分布图b=1,-2,0.8;a=1,2,2,1;z=roots(b);点p=roots(a);抽点sys=tf(b,a);pzmap(sys);%出零极点分布图title(零极点分布图)；xlabel( Real Axis);ylabel( Imaginary Axis );程序运行所得波形如图十四:-1.5零 极点 分 布 图0.80.60.40.20-0.2-0.4-0.6-0.8
13、-1-1-0.50Real Axis (seconds0.5T )1.5图十四零极点分布图的波形(4)幅频响应：程序如下：b=1,-2,0.8;a=1,2,2,1;H,w=freqs(b,a);plot(w,abs(H);title(幅频响应);xlabel(频率(rad/s);ylabel(幅度);程序运行所得波形如图十五：图十五幅频响应的波形(5)相频响应: 程序如下：%f频响应b=1,-2,0.8;a=1,2,2,1;H,w=freqs(b,a);phai=angle(H);plot(w,phai);%1出相频响应title(相频响应);xlabel( 3 (rad/s);ylabel(
14、 ();程序运行所得波形如图十六：相 频响 应43 - -2.1_40 _- 1_- 2- 3- 4EcccEE012345678910s (rad/s)图十六相频响应的波形(6)答：由零极点分布图的波形可知，该系统函数的极点全部位于虚轴的左边，所以该系统是稳定的。三、思考题1 .系统函数的零极点对系统频率特性有何影响？答：(1)极点主要影响系统频率特性的峰值，如果极点越靠近单位圆，那么其峰 值就越尖锐。(2)零点主要影响系统频率特性的谷值， 如果零点越靠近单位圆，那么谷值 就越深，特别是当零点在单位圆上时，其系统频率特性会变为零。2 .系统函数的零极点对系统冲激响应有何影响？答：(1)如果极
15、点位于s右半平面或在虚轴上具有二阶以上极点， 那么冲激响应 是发散的。(2)如果极点全部位于s左半平面，那么冲激响应是收敛的。(3)如果极点位于虚轴上且只具有一阶时， 那么冲激响应是非零数值或等幅 振荡。3 .如何根据系统的幅频特性计算出系统的截止频率?答： 截止频率的含义是系统对频率高于某一频率的信号分量具有过滤的功能，而频率低于该频率的信号分量则可直接通过或略有衰减。 在系统的幅频特性中， 我们规定幅频特性达到最大幅度的 0.707 倍并继续下降时的临界频率为系统的截止频率， 所以只要在系统的幅频特性曲线中找到最大幅度的 0.707 倍所对应的那个频率即可得系统的截止频率。四、实验总结回顾这次实验，我们主要是用 MATLABE分析连续系统的时域响应、频响特性和零极点分布。 在实验过程中， 由于本人可以清晰的看到该实验目的是让我们在深刻理解连续时间系统的系统函数在分析连续系统的时域特性、 频域特性及稳定性中的重要作用及意义下。通过利用MATLA来分析连续系统的时域响应、频响特性和零极点分布。（注：专业文档是经验性极强的领域，无法思考和涵盖全面，素材和资料部分 来自网络，供参考。可复制、编制，期待你的好评与关注）