信号与系统实训指导(MATLAB)第3章连续时间系统的时域分析.ppt

　　 本章介绍连续时间系统的时域分析以及如何用MATLAB函数计算连续线性时不变系统的响应。3.1节简要介绍连续时间系统时域分析用到的MATLAB函数。3.2节介绍用MATLAB求连续系统的单位阶跃响应和单位冲激响应。3.3节介绍用MATLAB求连续系统对任意输入的零状态响应。3.4节介绍连续信号卷积的计算和用卷积求解连续系统的零状态响应。3.5节为线性时不变系统实例分析。 　　　　连续时间系统的时域分析涉及到的MATLAB函数主要有：　　 step　　 impulse　　 lsim　　 dsolve　　 conv　　下面就简要介绍这些函数。 　　1.step　　函数step用于计算连续时间系统的单位阶跃响应。　　2.impulse　　函数impulse用于计算连续时间系统的单位冲激响应。　　3.lsim　　函数lsim用于计算连续系统在任意输入作用下的响应。　　4.dsolve　　这个函数是SymbolicMathToolbox中的函数，用于求解常微分方程。 　　　　线性时不变连续时间系统可以用微分方程描述： (3.1)式中:x(t)是输入,y(t)是输出,系数ak(k＝0,1,…,N)和bk(k＝0,1,…,M)均为常数。一般情况下,aN＝1,方程(3.1)可以写做  (3.2)　　单位阶跃响应是单位阶跃信号作用在系统产生的零状态响应；单位冲激响应是单位冲激信号作用在系统产生的零状态响应。用step函数或impulse函数求系统相应的响应时,先把式(3.1)中微分方程左右两边的系数用向量a和b表示,向量a表示响应及其各阶导数的系数： 　 a=[ aN aN-1 ... a1 a0]　　向量b表示激励及其各阶导数的系数：　　 b=[ bM bM-1 ... b1 b0]　　然后用命令y＝step(b,a,t)或y＝impulse(b,a,t)求系统的阶跃响应或冲激响应,命令中t表示计算响应所需的时间向量。例如t＝c∶d∶e,其中c是起始时间,e是结束时间,d是时间增量。为了做出平滑的曲线,向量t的增量需要取得足够小。需要注意的是,这里向量a和b的样本值都是按照导数阶数递减的次序排列的。 　　例3-1 求解微分方程 的单位阶跃响应和单位冲激响应,并与用解析式表示的结果进行比较。　　解：用手工计算微分方程得到的阶跃响应的表达式为冲激响应的表达式为 用MATLAB函数求解响应的程序如下： 　　 t＝[0∶0.1∶10］；　　 b＝[1］；　　 a＝[13］；　　 s＝step(b,a,t)；　　 s1＝(1/3)*(1－exp(－3*t；　　 plot(t,s,′o′,t,s1,′r′)；　　 title(′StepResponse′)；ylabel(′y(t)′)；xlabel(′t′)；legend(′step函数′,′解析式′) 　　 figure　　 h＝impulse(b,a,t)；　　 h1＝exp(－3*t)；　　 plot(t,h,′o′,t,h1,′r′)；　　 title(′ImpulseResponse′)；ylabel(′y(t)′)；xlabel(′t′)；legend(′impulse函数′,′解析式′)　　系统的阶跃响应波形如图3-1所示。 　　 图3-1 阶跃响应的波形 　　图3-2 冲激响应的波形 　　　　用MATLAB的lsim函数可以求线性时不变连续系统在任意输入作用下的响应。用lsim函数求系统的非零状态响应时，需要把描述系统的微分方程转化为状态方程，这里主要介绍如何求系统的零状态响应。 系统微分方程如式(3.1)所列,向量a、b和t与step函数中的定义相同,向量x包含在向量t区间内的输入信号x(t)的样本值,利用命令y＝lsim(b,a,x,t)就会得到微分方程表示的系统的零状态响应。由这个命令产生的输出向量y包含向量t中每一时刻上的微分方程输出值。　　例3-2 求解微分方程 的零状态响应,输入信号 并与用解析式表示的结果进行比较。 　　解：用手工计算微分方程得到的零状态响应的表达式为  用MATLAB函数求解响应的程序如下：　　 t＝[0∶0.1∶10］；　　 x＝exp(－1*t)；　　 b＝[1］；　　 a＝[13］；　　 s＝lsim(b,a,x,t)；　　 s1＝(1/2)*(exp(－1*t)－exp(－3*t； 　　 plot(t,s,′o′,t,s1,′r′)；　　 ylabel(′y_{zs}(t)′)；xlabel(′t′)；legend(′lsim函数′,′解析式′)　　系