基于Matlab的零输入和零状态响应.docVIP

课程设计题目： 信号系统的时域分析 二、课程设计目的： 1、学习MATLAB 软件的使用。 2、使学生掌握利用工具软件来实现信号系统基本概念、基本原理的方法。 3、通过编程对matlab软件的具体应用有了更好的了解，进一步加强了对函数卷积，零输入，零响应状态这三种函数状态的理解。 二、 信号的卷积运算在系统分析中主要用于求解系统的零状态响应。一般情况，卷积积分的运算比较困难，但在MATLAB中则变得十分简单，MATLAB中是利用conv函数来实现卷积的。 Conv函数功能：实现二个函数 和 的卷积。 格式： 说明： ??? ?表示二个函数， ????? 表示两个函数的卷积结果。 例：已知两信号 f1(t)= ε(t+1)+ε(t-2) f2(t)=ε(t-3)+ε(t-8)??????? 求卷积 ???? ?? MATLAB程序如下： t1=-1:0.01:2; f1=ones(size(t1));??%高度为一的门函数，时间从t=-1到t=2 t2=3:0.01:8; f2=ones(size(t2));??%高度为一的门函数，时间从t=3到t=8 g=conv(f1,f2);??????%对f1和f2进行卷积 t3=2:0.01:10; subplot(3,1,1),plot(t1,f1),title(‘f1’);grid 　%画f1的波形 subplot(3,1,2),plot(t2,f2) ,title(‘f2’);grid?%画f2的波形 subplot(3,1,3),plot(t3,g) ,title(‘f1*f2’);grid?? %画g的波形 2、零输入，零状态响应 一般的连续时间系统分析有以下几个步骤: ①求解系统的零输入响应; ②求解系统的零状态响应; ③求解系统的全响应; ④分析系统的卷积；⑤画出它们的图形. 下面以具体的微分方程为例说明利用MATLAB 软件分析系统的具体方法. 2.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。 以下面式子为例： y″(t)+3y′(t)+6y(t)=6f′(t)-8f′(t) 初始条件为y(0_)=0,y′(0_)=10; MATLAB程序： a=[1,3,6]; n=length(a)-1;Y0=[0,10]; p=roots(a);V=rot90(vander(p));c=V\Y0; dt=0.002;te=9; 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(零输入响应); 2.2连续时间系统零状态响应的数值计算 我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述， 例如，对于以下方程: 可用输入函数，得出它的冲激响应h ，再根据LTI系统的零状态响应y（t）是激励u（t）a或b中的对应元素应为零，不能省略不写，否则出错。 求函数的零状态响应 y″(t)+3y′(t)+6y(t)=6f′(t)-8f′(t) 及初始状态。输入函数。 先求出系统的冲激响应，写出其特征方程 求出其特征根为p和p，及相应的留数r，r；则冲激响应为 输入y（t）可用输入u（t）与冲激响应h（t）的卷积求得。 MATLAB程序： a=[1,3,6]; n=length(a)-1; Y0=[0,10]; b=[1,6,8]; dt=0.001;te=9; t=0:dt:te; u=sin(2*t)+cos(5*t); p=roots(a); %求出齐次解 V=rot90(vander(p)); %生成2×2矩阵并旋转90度