信号与系统课程设计报告书.docVIP

信号与系统课程设计报告书 课题名称 零输入响应与零状态响应 姓 名 学 号 院、系、部 电气工程系 专 业 电子信息工程 指导教师 2013年1月8日 连续时间系统的LTI系统的时域仿真 -------零输入响应与零状态响应 一掌握信号经过LTI系统的时域分析方法。 ?巩固已经学过的知识，加深对知识的理解和应用，加强学科间的横向联系，学会应用MATLAB对实际问题进行仿真。二要求）: ①求解系统的零输入响应; ②求解系统的零状态响应; ③求解系统的全响应; ④分析系统的卷积；⑤画出它们的图形. 下面以具体的微分方程为例说明利用MATLAB 软件分析系统的具体方法. 1．连续时间系统的零输入响应 描述n阶线性时不变（LTI）连续系统的微分方程为： 已知y及各阶导数的初始值为y(0),y(1)(0)，… y(n-1)(0)， 求系统的零输入响应。 建模 当LIT系统的输入为零时，其零输入响应为微分方程的其次解（即令微分方程的等号右端为零），其形式为（设特征根均为单根） 其中p1,p2,…,pn是特征方程a1λn+a2λn-1+…+anλ+an=0的根，它们可以用root(a)语句求得。各系数 由y及其各阶导数的初始值来确定。对此有 写成矩阵形式为: P1n-1C1+ P2n-1C2+…+ Pnn-1Cn=Dn-1y0 即 V?C=Y0 其解为：C=V\Y0 式中 V为范德蒙矩阵，在matlab的特殊矩阵库中有vander。 以下面式子为例： 初始条件为,； MATLAB程序： a=input(输入分母系数a=[a1,a2,...]=); n=length(a)-1; Y0=input(输入初始条件向量 Y0=[y0,Dy0,D2y0,...]=); p=roots(a);V=rot90(vander(p));c=V\Y0; dt=input(dt=);te=input(te=); t=0:dt:te;y=zeros(1,length(t)); for k=1:n y=y+c(k)*exp(p(k)*t);end plot(t,y);grid xlabel(t) ;ylabel(y); title(零输入响应); 程序运行结果： 用这个通用程序来解一个三阶系统，运行此程序并输入 a=[1,2,3] Y0=[2,7] dt=0.002 te=7 结果如下图：根据图可以分析零输入响应，它的起始值与输入函数无关，只与它的初始状态值有关，其起始值等于y(0_)的值。随着时间的推移，最后零输入响应的值无限的趋近于0。 2．连续时间系统零状态响应的数值计算 我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述， 例如，对于以下方程: 可用输入函数，得出它的冲激响应h ，再根据LTI系统的零状态响应y（t）是激励u（t）a或b中的对应元素应为零，不能省略不写，否则出错。 求函数的零状态响应 及初始状态。输入函数。 建模 先求出系统的冲激响应，写出其特征方程 求出其特征根为p和p，及相应的留数r，r；则冲激响应为 输入y（t）可用输入u（t）与冲激响应h（t）的卷积求得。 MATLAB程序： a=input(输入分母系数a=[a1,a2,...]=); b=input(输入输入信号系数b=[b1,b2,...]=); dt=input(dt=);te=input(te=); t=0:dt:te; u=input(输入函数u=); te=t(end); dt=te/(length(t)-1); [r,p,k]=residue(b,a); h=r(1)*exp(p(1)*t)+r(2)*exp(p(1)*t); subplot(2,1,1),plot(t,h);grid title(冲激函数); y=conv(u,h)*dt; subplot(2,1,2), plot(t,y(1:length(t)));grid title(零状态响应); 程序运行结果 执行这个程序，取a=[1,2,3] b=[4,6] dt=0.001 te=7 得出图形如下： 由于初始状态为零，所以零状态的起始值也为零，即h(t)包含了连续系统的固有特性，与系统的输入无关。只要知道了系统的冲激响应，即可求得系统在不同输入时产生的输出。因此，求解系统的冲激响应h对进行连续时间系统的分析具有非常重要的意义 3．连续时间系统的全响应计算 上面通过对LTI系统函数的描述，我们可