连续时间系统的lti系统的时域仿真零输入响应与零状态响应毕业论文

《连续时间系统的lti系统的时域仿真零输入响应与零状态响应毕业论文》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、※※※※※※※※※※※※※※※※※※※※※※※※2009级信号系统课程设计信号与系统课程设计报告书课题名称连续时间系统的LTI系统的时域仿真零输入响应与零状态响应姓名学号院、系、部电气工程系专业电子信息工程指导教师2011年12月27日连续时间系统的LTI系统的时域仿真零输入响应与零状态响应　　　　　　　　　　　　　　　　　　一、设计目的１、掌握信号经过LTI系统的时域分析方法。２、掌握连续时不变信号处理的基本概念、基本理论和基本方法；３、学会MATLAB的使用，掌握MATLAB的程序设计方法；４、学会用MATLA

2、B对信号进行分析和处理；５、编程实现卷积积分或卷积和，零输入响应，零状态响应； ６、、巩固已经学过的知识，加深对知识的理解和应用，加强学科间的横向联系，７、学会应用MATLAB对实际问题进行仿真。学会对带有非零起始状态的LTI系统进行仿真。二、设计内容及要求（1）根据实际问题建立系统的数学模型，找一个实际的电路，建立如课本第二章例2-8所示的数学模型，并计算其完全响应；（2）用MATLAB描述此系统；（3）仿真实现并绘制输出信号的波形。要求用两种方法仿真实现完全响应。对仿真结果进行比较，并与理论值比较。三、设计方法

3、与步骤：一般的连续时间系统分析有以下几个步骤:①求解系统的零输入响应;②求解系统的零状态响应;③求解系统的全响应;④分析系统的卷积；⑤画出它们的图形。1．连续时间系统的零输入响应描述n阶线性时不变（LTI）连续系统的微分方程为：已知y及各阶导数的初始值为y(0),y(1)(0)，…y(n-1)(0)，求系统的零输入响应。建模当LIT系统的输入为零时，其零输入响应为微分方程的其次解（即令微分方程的等号右端为零），其形式为（设特征根均为单根）其中p1,p2,…pn是特征方程a1λn+a2λn-1+…+an-1λ+an=

4、0的根，它们可以用root(a)语句求得。各系数由y及其各阶导数的初始值来确定。对此有…………………………………………写成矩阵形式为:P1n-1C1+P2n-1C2+…+Pnn-1Cn=Dn-1y0即V•C=Y0其解为：C=V/Y0式中V为范德蒙矩阵，在matlab的特殊矩阵库中有vander。2．卷积的计算连续时间信号和的卷积运算可用信号的分段求和来实现，即：如果只求当t=nD（n为整数）时f(t)的值f(nD)，则上式可得：式中的实际上就是连续时间信号和经等时间间隔均匀抽样的离散序列和的卷积和。当D足够小时，就

5、是卷积积分的结果——连续时间信号f(t)的较好数值近似。建模利用MATLAB实现连续信号卷积的通用程序conv()，程序在计算出卷积积分的数值近似的同时，还绘制出f(t)的时域波形图。3．连续时间系统零状态响应的数值计算我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述，例如，对于以下方程:可用输入函数，得出它的冲击响应h，再根据LTI系统的零状态响应y（t）是激励u（t）与冲击响应h（t）的卷积积分。注意，如果微分方程的左端或右端表达式中有缺项，则其向量a或b中的对应元素应为零，不能省略不写，否则出错

6、。4．连续时间系统的全响应计算上面通过对LTI系统函数的描述，我们可以得知：如果在系统的初始状态不为零，在激励f（t）的作用下，LTI系统的响应称为全响应，它是零输入响应和零状态响应之和，即故可先求出零输入响应和零状态响应，再把两者相加，得到全响应。但简单的相加可能由于零输入与零状态的矩阵不同而不能的出正确的结果，这就需要对矩阵进行截取，使它们的阶数相同。四、程序课本例2-8，已知系统方程式r＇（ｔ）+3ｒ（ｔ）＝3ｅ（ｔ），初状态为ｒ（０-）＝1.5,激励信号为e（t）=u（t），求系统的完全响应。程序如下：a=

7、[1,3];n=length(a)-1;Y0=[3/2];b=[3];dt=0.01;te=6;t=0:dt:te;jieyue=zeros(size(t));jieyue(t>0)=1;jieyue(t==0)=1/2;u=jieyue;p=roots(a);V=rot90(vander(p));c=Y0/V;y1=zeros(1,length(t));fork=1:ny1=y1+c(k)*exp(p(k)*t);endte=t(end);dt=te/(length(t)-1);[r,p,k]=residue(b

8、,a);h=r(1)*exp(p(1)*t);y2=conv(u,h)*dt;y=y1(1:length(t))+y2(1:length(t));figure(1);subplot(3,1,1),plot(t,y1),gridxlabel('t');ylabel('y1');title('零输入响应');subplot(3,1,2),plot(t,y2(1: