资源描述:
《连续时间系统的时分析》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、实验三连续时间系统的时域分析一实验目的:1、熟悉和掌握常用的用于信号与系统吋域分析的MATLAB函数;2、掌握如何利用Matlab软件求解一个线性时不变连续时间系统的零状态响应、冲激响应和阶跃响应。二实验原理:在倍号与线性系统屮,LTI(线性时不变)连续时间系统以常系数微分方程描述,系统的零状态响应可以通过求解初始状态为零的微分方程得到。在Matlab屮,控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数Mm,其调用形式为:y=Isim{sys,f,t)式屮,t表示计算系统响应的抽样点向量,f是系统输入信号向量(即激励),是LTI系统模型,用来表示微分方程。在
2、求解微分方程时,微分方程的LTI系统模型sys要借助Matlab屮的tf函数来获得,其调用形式为:sys=tf{b,a)式屮,b和a分别为微分方程右端和左端各项的系数向量。例如对于三阶微分方程:以3,(,)+%/(,)+“,/(,)+aQy(t)=W)+b2ft)++1)(、f(t)可以用以下命令:b=[b3,b2,bl,b0];a=[a3,a2,al,aO];sys=tf(b,a);来获得LTI模型。系统的LTI模型建立后,就可以求出系统的冲激响应和阶跃响应。在连续吋间LTI屮,冲击响应和阶跃响应是系统特性的描述。输入为单位冲击函数利r)所引起的零状态响应称为单位冲击
3、响应,简称冲击响应,用/咐)表示;输入为单位实验三连续时间系统的时域分析一实验目的:1、熟悉和掌握常用的用于信号与系统吋域分析的MATLAB函数;2、掌握如何利用Matlab软件求解一个线性时不变连续时间系统的零状态响应、冲激响应和阶跃响应。二实验原理:在倍号与线性系统屮,LTI(线性时不变)连续时间系统以常系数微分方程描述,系统的零状态响应可以通过求解初始状态为零的微分方程得到。在Matlab屮,控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数Mm,其调用形式为:y=Isim{sys,f,t)式屮,t表示计算系统响应的抽样点向量,f是系统输入信号向量(即激
4、励),是LTI系统模型,用来表示微分方程。在求解微分方程时,微分方程的LTI系统模型sys要借助Matlab屮的tf函数来获得,其调用形式为:sys=tf{b,a)式屮,b和a分别为微分方程右端和左端各项的系数向量。例如对于三阶微分方程:以3,(,)+%/(,)+“,/(,)+aQy(t)=W)+b2ft)++1)(、f(t)可以用以下命令:b=[b3,b2,bl,b0];a=[a3,a2,al,aO];sys=tf(b,a);来获得LTI模型。系统的LTI模型建立后,就可以求出系统的冲激响应和阶跃响应。在连续吋间LTI屮,冲击响应和阶跃响应是系统特性的描述。输入为单位
5、冲击函数利r)所引起的零状态响应称为单位冲击响应,简称冲击响应,用/咐)表示;输入为单位阶跃函数£(z)所引起的零状态响应称为单位阶跃响应,简称阶跃响应,用《)表示。求解系统的冲激响应的函数是impulse,求解系统的阶跃响应可以利用函数step,其调用形式分别为:y=和=step(sys,t)式屮t表示计算系统响应的杣样点向量,sys是LTI系统模型三实验内容:—、己知系统的系统转移算子为求该系统的零状态响应曲线。假设系统的激励冲)=sin⑴,t在
6、0,2川之间,步长0.01。t=0:0.01:2*pi;、已知系统的系统转移算子为H2(/7)=p2+2/7+3,求该系统
7、的零状态响应曲b=[22];a=[133];sys=tf(b,a};f=sin⑴;y=lsim(sys,f,t);Plot(t,y)xlabelCf);ylabelfy⑴•);gridon;线。假设系统的激励冲)=cos⑺,t在[0,2兀]之间,步长0.01t=0:0.01:2*pi;b=[20];a=[123];sys=tf(b,a);f=cos(t);y=lsim(sys,f,t);plot{t,y)xlabel(Y);ylabelCy⑴gridon;s三'已知系统的微綱为:»仂,,求该系统的零状态响应曲线。假设系统的激励冲)=10sin(2尔),t在[0,2兀]之间
8、,步长0.01。t=0:0.01:2*pi;b=[6];a=[l56];sys=tf(b,a);f=10*sin(2*pi*t);y=lsim(sys,f,t);plot{t,y)xlabel(Y);ylabelfy⑴%gridon;四'已知系统的微分撤芽十呼偷⑺,),求系统的冲激响冲激响应应和阶跃响应曲线,将两幅图显示在一个窗U,t在[0,2;r]之间,步长0.01。t=0:0.01:2*pi;b=[10];a=[l2100];sys=tf(b,a);yl=impulse(sys/t);y2=step(sys/t);subp