实验三连续时间系统的时域分析

《实验三连续时间系统的时域分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验三连续时间系统的时域分析实验目的：1、熟悉和掌握常用的用于信号与系统时域分析的MATLAB函数；2、掌握如何利用Matlab软件求解一个线性时不变连续时间系统的零状态响应、冲激响应和阶跃响应。实验原理：在信号与线性系统中，LTI(线性时不变)连续时间系统以常系数微分方程描述，系统的零状态响应可以通过求解初始状态为零的微分方程得到。在Matlab中，控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数Isim,其调用形式为：y=lsim{sysjt}式中，t表示计算系统响应的抽样点向量，f是系统输入信号向量(

2、即激励)，是LTI系统模型，用来表示微分方程。在求解微分方程时，微分方程的LTI系统模型sys耍借助Matlab中的tf函数来获得，其调用形式为：sys=tf{b,a)式中，b和a分别为微分方程右端和左端各项的系数向量。例如对于三阶微分方程：%,(，)+“2/(，)+/⑺+aQy(t)=W++/V’⑺+M•⑴可以用以下命令：b=[b3，b2，bl，b0];a=[a3，a2，al，a0];sys=tf(b,a);来获得LTI模型。系统的LTI模型建立后，就可以求出系统的冲激响应和阶跃响应。在连续时间LTI中，冲击响应和阶跃响

3、应是系统特性的描述。输入为单位冲击函数＜5(Z)所引起的零状态响应称为单位冲击响应，简称冲击响应，用/z(Z)表示；输入为单位阶跃函数£«所引起的零状态响应称为单位阶跃响应，简称阶跃响应，用wW表示。求解系统的冲激响应的函数是impulse,求解系统的阶跃响应可以利用函数step,其调用形式分别为：y=bnpulse（sys,0和y=xtep（sys,Z）式中t表示计算系统响应的抽样点向量，sys是LTI系统模型。实验内容：一、已知系统的系统转移算子为//，（妁二，2/?+2，求该系统的零状态响应曲1厂2+3/?+3线。假

4、设系统的激励冲）=sin（z）,t在［0,2川之间，步长0.01。程序代码：波形图：t=0:0.01:2*pi;b=[22j;AP零状态哺应a=[133];uo06ZAsys=tf(b,a);04/•f=sin⑴；02.•y=lsim(sys，f，t);L/plot(t,y)•02/•xlabel(’t’)；-0.4•/ylabel(’y⑴’);-0.6/-title（’零状态响应’）_nftgridon;•08£12345671二、已知系统的系统转移算子为=—，求该系统的零状态响应曲p一+2/?+3线。假设系

5、统的激励冲）=cos⑴，t往［0,2;r］之间，步长0.01。程序代码：波形图：t=0:0.01:2*pi;08零状态杭应2b=[20];a=[123];0.6sys=tf(b，a);04//■f=cos(t);0.2\/-y=lsim(sys,f,t);'/plot(t，y)•02/■xlabelCt’)；•0.4-ylabeKW’)；•0.6title（’零状态响应2’）nq皇觀皇1t1O8()1234567gridon;1三'己知系统的微分方程为:»（輯），求该系统的零状态响应曲线。假设系统的激励冲）=10

6、sin（2;zr）,t在［0,2兀］之间，步长0.01。程序代码：波形图：t=0:0.01:2*pi;nq苓状态析应3b=[6];U.J0.25Aa=[l56];0.2sys=tf（b,a）;015IAAAAAf=10*sin（2*pi*t）;01sAAAA•y=lsim（sys，f，t）;0.05yplot（t,y）0•005V'xlabel（’t’）；ylabelCy（t）1）;titleC零状态响应3’）gridon;•01nicVVvVv■垂■11畫1•U.IO123456712四、己知系统的微分方程为：++=求

7、系统的冲激响drdt应和阶跃响应曲线，将两幅图显示在一个窗口，t在.［0,2;r］之间，步长0.01。程序代码:t=0:0.01:2*pi;b=[10];a=[l2100];sys=tf（b，a）;y1=impulse（sys，t）;y2=step（sys,t）;subplot（2，l，l）;plot（t，yl）;xlabel（’t’）；ylabel（’h⑴’）；title（‘冲激响应•）gridon;subplot（2，l，2）;plot（t，y2）;xlabel（’t’）；ylabel（’u⑴’）；title（‘阶跃响

8、应*）波形图:醐应gridon;ry五、已知系统的微分方程为：^4^+2^4-32r(Z)=-^+16^(Z),假设系统drdtdt的激励e(r)=f，t在［0,2兀］之间,步长0.01。将系统的激励函数、冲激响应和零状态响应显示在一个窗口。程序代码：t=0:0.01:2*pi;b=[l16];a=[l