matlab求解零输入响应,基于matlab的零输入和零状态响应

基于matlab的零输入和零状态响应

一、 课程设计题目： 信号系统的时域分析 二、课程设计目的： 1、学习MATLAB 软件的使用。 2、使学生掌握利用工具软件来实现信号系统基本概念、基本原 理的方法。 3、通过编程对 matlab 软件的具体应用有了更好的了解，进一步 加强了对函数卷积，零输入，零响应状态这三种函数状态的理解。 二、 基本原理 1、卷积 信号的卷积是数学上的一种积分运算，两个信号的卷积定义为： 信号的卷积运算在系统分析中主要用于求解系统的零状态响应。一 般情况，卷积积分的运算比较困难，但在 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系统的输入为零时，其零输入响应为微分方程的其次解 (即令微分方程的等号右端为零) ，其形式为(设特征根均为单根) 1 1 2 1 1 1 1 n n m n n m m n n m d y d y dy d u du a a a a y b b b u dt dt dt dt dt                     其中 p 1 ,p 2 ,…,p n 是特征方程 a 1 λ n +a 2 λ n-1 +…+a n λ+a n =0 的根，它们 可以用root(a)语句求得。各系数 由y及其各阶导数的初始值来 确定。对此有 写成矩阵形式为: P 1 n-1 C 1 + P 2 n-1 C 2 +…+ P n n-1 C n =D n-1 y 0 1 0 1 2 2 0 1 1 1 1 1 2 0 1 1 1 n n n n n n n C y p p p C Dy p p p C D y                                                   A A A A A A 即 V•C=Y 0 其解为：C=V\Y 0式中 V为范德蒙矩阵，在matlab的特殊矩阵库中有vander。 以下面式子为例： y″(t)+3y′(t)+6y(t)=6f′(t)-8f′(t) 初始条件为y(0_)=0,y′(0_)=10; MATLAB程序： a=[1,3,6]; 1 2 1 2 ( ) n p t pt p t n y t Ce C e C e         1 2 0 n C C C y         1 1 2 2 0 n n pC p C p C Dy         1 1 1 1 1 1 2 2 0 n n n n n n p C p C p C D y             1 2 1 1 1 1 2 1 1 1 n n n n n p p p V p p p                          A A A A   1 2 n C C C C     1 0 0 0 n C y Dy D y         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连续系统可用如下所示的线性常系数微分方程来 描述， ( ) ( ) 0 0 ( ) ( ) N M i j i j i j a y t b f t      例如，对于以下方程: 3 2 1 0 3 2 1 0 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) a y t a y t a y t a y t b f t b f t b f t b f t        可用 输入函数 ，得出它的 3 2 1 0 3 2 1 0 [ , , , ], [ , , , ], a a a a a b b b b b   ( ) u f t  冲激响应h ，再根据LTI系统的零状态响应y(t)是激励u(t) 与冲激响应h(t)的卷积积分。 注意，如果微分方程的左端或右端表达式中有缺项，则其向量 a 或 b中的对应元素应为零，不能省略不写，否则出错。 求函数的零状态响应y″(t)+3y′(t)+6y(t)=6f′(t)-8f′(t) 及初始状态 。输入函数 。 (0 ) (0 ) 0 zs zs y y     ) * 5 cos( ) * 2 sin( t t y   先求出系统的冲激响应，写出其特征方程0 3 2 2      求出其特征根为p和p，及相应的留数r，r；则冲激响应为1 2 1 2 ( ) pt p t h t re re   输入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;