MATLAB仿真之_连续时间LTI系统仿真和时域解释.docVIP

郑州航空工业管理学院 《电子信息系统仿真》课程设计 级电子信息工程专业班级 题 目 连续时间LTI系统仿真和时域分析 姓 名学号 指导教师 二О一年月日 MATLAB软件简介 MATLAB 是MathWork 公司于1984 年推出的一套面向工程和科学运算的高性能软件，它具有强大的图形处理功能及符号运算功能，为我们实现信号的可视化及系统分析提供了强有力的工具。MATLAB 强大的工具箱函数可以分析连续信号、连续系统，同样也可以分析离散信号、离散系统，并可以对信号进行各种分析域计算，如相加、相乘、移位、反折、傅里叶变换、拉氏变换、Z 变换等等多种计算MATLAB用于算法开发、数据可视化、数据分析以数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 Matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。Simulink是MATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。 设计目的 掌握信号经过LTI系统的时域分析方法。根据连续时不变信号处理的基本概念、理论和方法对信号进行分析和处理，实现卷积积分或卷积和，零输入响应和零状态响应，学会应用MATLAB对实际问题进行仿真原理原理 连续时间零状态响应原理 零状态响应就是在电路初始状态下（动态储能元件储能为零）由外施激励引起的反应，最终状态的为一确定的实数。 LTI连续系统可用如下所示的线性常系数微分方程来描述， Matlab实现及波形输出 卷积分程序及仿真 利用conv函数实现连续时间函数卷积分运算实际上就是先将时间离散化，在计算卷积分，最后根据函数画出输出函数波形，如求f1（t），f2（t）的卷积 ，f (t)=f1（t）*f2（t）。 MATLAB程序： t=0:0.1:2*pi; f1=input(输入函数f1=); f2=input(输入函数f2=); dt=input(dt=); y=conv(f1,f2); plot(dt*([1:length(y)]-1),y); grid on; title(卷积); xlabel(t); ylabel(f1*f) 程序运行结果： 输入以下数据： f1=sin(3*t) f2=cos(3*t+2) dt=0.01 得出图形如图4.1所示： 图4.1 卷积和输出波形图 卷积分运用：利用卷积求零状态 系统零状态响应，其中激励信号e(t),系统单位冲激响应h(t)。 已知系统单位冲激响应，系统激励，使用卷积法求系统零状态响应。 MATLAB程序为： a=1000; t=-5:1/a:5; h=0.5*t.* (stepfun(t,0)-stepfun(t,2) ); e=stepfun(t,-0.5)-stepfun(t,1); r=conv(h,e); t=-10:1/a:10; plot(t,r); title(零状态响应r(t)); xlabel(t); ylabel(r); 零输入程序及仿真 建模 当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(0_)=1,y(0_)=5； MATL