实验三 连续时间LTI系统的时域分析
连续时间LTI系统的时域分析
![连续时间LTI系统的时域分析](https://img.taocdn.com/s3/m/03f483f276a20029bd642d90.png)
0.6
1
2 Time(sec)
3
4
第六章 周期信号的傅里叶级数及频谱
6.1 周期信号的傅里叶级数
设周期信号f (t ),其周期是T,角频率为0 2f 0 2 , T 则该信号可展开为三角 式的傅里叶级数,即 f (t ) a0 a1 cos0t a2 cos0t ... b1 sin 0t b2 cos0t ... a0 (an cos n0 bn sin n0t )
解:求的零输入响应和零状态响应后,完全响应则为二者之和
>> eq='D2y+3*Dy+2*y=0'; %齐次解求零输入相应 >> cond='y(0)=1,Dy(0)=2'; >> yzi=dsolve(eq,cond);yzi=simplify(yzi) yzi = -3*exp(-2*t)+4*exp(-t) >> eq1='D2y+3*Dy+2*y=Dx+3*x'; %零状态响应求解 >> eq2='x=exp(-3*t)*heaviside(t)'; >> cond='y(-0.001)=0,Dy(-0.001)=0'; %起始条件 >> yzs=dsolve(eq1,eq2,cond);yzs=simplify(yzs.y) yzi = heaviside(t)*(-exp(-2*t)+exp(-t)) >> yt=simplify(yzi+yzs) %完全响应 yt = -3*exp(-2*t)+4*exp(-t)-exp(-2*t)*heaviside(t)+exp(-t)*heaviside(t)
连续时间LTI系统的时域分析-副本解读
![连续时间LTI系统的时域分析-副本解读](https://img.taocdn.com/s3/m/2ad8d3435901020207409c78.png)
1、已知系统的微分方程和激励信号,使用MATLAB的符号运算命令画出系统的零状态、零输入响应和完全响应:y’’+4y’+4y=x’(t)+3x(t);clear alleq1='D2y+4*Dy+4*y=0';eq2='D2y+4*Dy+4*y=Dx+3*x';eq3='x=exp(-t)*heaviside(t)';cond1='y(0)=1,Dy(0)=1';cond2='y(-0.01)=0,Dy(-0.01)=0';answer1=dsolve(eq1,cond1);answer2=dsolve(eq2,eq3,cond2);answer3.y=answer1+answer2.y;subplot(3,1,1);ezplot(answer2.y);title('零状态响应');subplot(3,1,2);ezplot(answer1);title('零输入响应');subplot(3,1,3);ezplot(answer3.y);title('完全响应');00.51 1.522.533.544.550.10.20.3t零状态响应-6-5-4-3-2-10123x 105t零输入响应-6-4-224-4-205t 完全响应2、某连续系统的输入为e(t),输出为y(t),系统的微分方程为:y ’’(t)+5y ’(6)+6y(t)=3f ’(t)+2f(t) ❶绘出该系统的单位冲激响应和单位阶跃响应时域图形 ❷若,绘出系统的零状态响应时域图形clear allsys=tf([3 2],[1 5 6]); t=0:0.01:3; h=impulse(sys,t); g=step(sys,t); subplot(2, 1 ,1); plot(t,h);title('单位冲激响应'); subplot(2 ,1 ,2); plot(t,g);title('单位阶跃响应');-10123单位冲激响应00.20.40.60.8单位阶跃响应(2)clear allsys=tf([3 2],[1 5 6]); t=0:0.01:3;f=exp(-2*t).*heaviside(t); h=lsim(sys,f,t); plot(t,h);title('零状态响应');零状态响应3、已知描述系统的微分方程和激励信号f(t )分别如下,试用数值方法求系统的单位冲激响应h(t)和零状态响应y(t ),并根据y(t)=f(t)*h(t)验证系统零状态响应的波形。
信号与系统实验报告实验三连续时间LTI系统的频域分析报告
![信号与系统实验报告实验三连续时间LTI系统的频域分析报告](https://img.taocdn.com/s3/m/4678a485ce2f0066f53322d2.png)
. .实验三 连续时间LTI 系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义;2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法,理解具有不同频率响应特性的滤波器对信号的滤波作用;3、学习和掌握幅度特性、相位特性以及群延时的物理意义;4、掌握用MATLAB 语言进行系统频响特性分析的方法。
基本要求:掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法,深刻理解LTI 系统的频率响应特性的物理意义,理解滤波和滤波器的概念,掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。
二、实验原理及方法1 连续时间LTI 系统的频率响应所谓频率特性,也称为频率响应特性,简称频率响应(Frequency response ),是指系统在正弦信号激励下的稳态响应随频率变化的情况,包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。
上图中x(t)、y(t)分别为系统的时域激励信号和响应信号,h(t)是系统的单位冲激响应,它们三者之间的关系为:)(*)()(t h t x t y =,由傅里叶变换的时域卷积定理可得到:)()()(ωωωj H j X j Y =3.1或者: )()()(ωωωj X j Y j H =3.2)(ωj H 为系统的频域数学模型,它实际上就是系统的单位冲激响应h(t)的傅里叶变换。
即⎰∞∞--=dt e t h j H tj ωω)()( 3.3 由于H(j )实际上是系统单位冲激响应h(t)的傅里叶变换,如果h(t)是收敛的,或者说是绝对可积(Absolutly integrabel )的话,那么H(j )一定存在,而且H(j )通常是复数,连续时间LTI 系统的时域及频域分析图系统LTI )(t h )(ωj H )(t y )(ωj X )(ωj Y )(t x坐标形式:)()()(ωϕωωj ej H j H = 3.4上式中,)j (ωH 称为幅度频率相应(Magnitude response ),反映信号经过系统之后,信号各频率分量的幅度发生变化的情况,)(ωϕ称为相位特性(Phase response ),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。
实验三连续时间LTI系统的时域分析实验报告
![实验三连续时间LTI系统的时域分析实验报告](https://img.taocdn.com/s3/m/5d4872b003d276a20029bd64783e0912a2167cfd.png)
实验三连续时间LTI系统的时域分析实验报告一、实验目的通过实验三的设计和实现,达到如下目的:1、了解连续时间LTI(线性时不变)系统的性质和概念;2、在时域内对连续时间LTI系统进行分析和研究;3、通过实验的设计和实现,了解连续时间LTI系统的传递函数、共轭-对称性质、单位冲激响应等重要性质。
二、实验原理在常见的线性连续时间系统中,我们知道采用差分方程的形式可以很好地表示出该系统的性质和特点。
但是,在本实验中,我们可以采用微分方程的形式来进行相关的研究。
设系统的输入为 x(t),输出为 y(t),系统的微分方程为:其中,a0、a1、…、an、b0、b1、…、bm为系统的系数,diff^n(x(t))和diff^m(y(t))分别是输入信号和输出信号对时间t的n阶和m阶导数,也可以记为x^(n)(t)和y^(m)(t)。
系统的单位冲激响应函数 h(t)=dy/dx| x(t)=δ(t),则有:其中,h^(i)(t)表示h(t)的第i阶导数定义系统的传递函数为:H(s)=Y(s)/X(s)在时域内,系统的输出y(t)可以表示为:其中,Laplace^-1[·]函数表示Laplace逆变换,即进行s域到t域的转化。
三、实验步骤1、在Simulink中,构建连续时间LTI系统模型,其中系统的微分方程为:y(t)=0.1*x(t)-y(t)+10*dx/dt2、对系统进行单位冲激响应测试,绘制出系统的单位冲激响应函数h(t);4、在S函数中实现系统单位冲激响应函数h(t)的微分方程,并使用ODE45框图绘制出系统单位冲激响应函数h(t)在t=0~10s之间的图像;6、利用数据记录栏,记录系统在不同的参数下的变化曲线、阶跃响应函数u(t)和单位冲激响应函数h(t)的变化规律。
四、实验数据分析1、单位冲激响应测试那么,当输入信号为单位冲激函数δ(t)时,根据系统的微分方程,可以得知输出信号的形式为:即单位冲激响应函数h(t)为一个包含了单位冲激函数δ(t)在内的导数项序列。
连续LTI系统的时域分析
![连续LTI系统的时域分析](https://img.taocdn.com/s3/m/6d5ccd26192e45361066f585.png)
连续LTI 系统的时域分析 1、 实验目的(1) 熟悉连续LTI 系统在典型激励信号下的响应及其特征(2) 掌握连续LTI 系统单位冲激响应的求解方法 (3) 掌握用卷积法计算连续时间系统的零状态响应(4) 掌握matlab 相关函数的调用格式和作用 (5) 能够用matlab 对系统进行时域分析 2、 实验原理对于一般的n 阶LTI 连续系统,其微分方程的形式可以写成()(1)(1)()(1)(1)110110()()...()()()()...()()n n m m n m m y t a y t a y t a y t b f t b f t b f t b f t ----++++=++++式中系数均为常数,y(t)为系统的相应变量,f(t)为系统的激励信号 对于一个动态系统而言,其响应不仅与激励有关,还与系统的初始状态有关,对于线性系统通常分为零输入响应和零状态响应两部分。
对于低阶系统一般可以通过解析法得到响应,但是对于高阶系统,手工计算比较困难,利用matlab 强大的计算功能就能方便的得到系统的冲激响应、阶跃响应、零输入响应、零状态响应和全响应。
(1) 直接求解法 涉及的函数有impulse(冲激响应),step (阶跃响应),roots (零输入响应),lsim(零状态响应)等,在matlab 中,要求以系数向量的形式输入系统的微分方程,因此在使用前必须对系统的微分方程进行变换,得到其传递函数,分别用向量a 和b 表示分母多项式和分子多项式的系统(按S 的降幂排列) (2) 卷积计算法 用卷积积分极端线性时不变系统的零状态响应。
若对于LTI 系统,则系统的冲激响应为h(t),当输入为f(t)时系统的零状态响应()()()()*()zs y t f h t d f t h t τττ+∞-∞=-=⎰若系统为因果系统,则输入信号从t=0时加入,上式变为 0()()()yzs t f h t d τττ∞=-⎰3、 涉及的matlab 函数 (1) lsim 函数 功能:计算并画出系统在任意输入下的零状态响应 调用格式:lsim(b,a,x,t)其中b 和a 是由描述系统的微分方程的系数决定的表示该系统的两个行向量,x 和t 是表示输入信号的行向量(t 表示输入信号时间范围的向量,x 表示输入信号在向量t 定义的时间点上的取样值)。
实验三 连续时间LTI系统的时域分析实验报告
![实验三 连续时间LTI系统的时域分析实验报告](https://img.taocdn.com/s3/m/3b91d92f4028915f814dc29c.png)
实验三连续时间L TI系统的时域分析实验报告实验三连续时间LTI系统的时域分析一、实验目的1、学会使用符号法求解连续系统的零输入响应和零状态响应2、学会使用数值法求解连续系统的零状态响应3、学会求解连续系统的冲激响应和阶跃响应二、实验原理及实例分析1、连续时间系统零输入响应和零状态响应的符号求解连续时间系统可以使用常系数微分方程来描述,其完全响应由零输入响应和零状态响应组成。
MATLAB符号工具箱提供了dsolve函数,可以实现对常系数微分方程的符号求解,其调用格式为:dsolve(‘eq1,eq2…’,’cond1,cond2,…’,’v’)其中参数eq表示各个微分方程,它与MATLAB符号表达式的输入基本相同,微分和导数的输入是使用Dy,D2y,D3y来表示y的一价导数,二阶导数,三阶导数;参数cond表示初始条件或者起始条件;参数v表示自变量,默认是变量t。
通过使用dsolve函数可以求出系统微分方程的零输入响应和零状态响应,进而求出完全响应。
2、连续时间系统零状态响应的数值求解在实际工程中使用较多的是数值求解微分方程。
对于零输入响应来说,其数值解可以通过函数initial来实现,而该函数中的参量必须是状态变量所描述的系统模型,由于现在还没有学习状态变量相关内容,所以此处不做说明。
对于零状态响应,MATLAB控制系统工具箱提供了对LTI系统的零状态响应进行数值仿真的函数lsim,利用该函数可以求解零初始条件下的微分方程的数值解。
其调用格式为:y=lsim(sys,f,t),其中t表示系统响应的时间抽样点向量,f是系统的输入向量;sys表示LTI系统模型,用来表示微分方程、差分方程或状态方程。
在求解微分方程时,sys是有tf函数根据微分方程系数生成的系统函数对象,其语句格式为:sys=tf(a,b)。
其中,a和b分别为微分方程右端和左端的系数向量。
例如,对于微分方程a3y'''(t)?a2y''(t)?a1y'(t)?a0y(t)?b3f'''(f)?b2f''(t)?b1f'(t)?b0f(t) 可以使用a?[a3,a2,a1,a0];b?[b3,b2,b1,b0];sys?tf(b,a)获得其LTI模型。
3_2连续时间LTI系统响应的时域分析
![3_2连续时间LTI系统响应的时域分析](https://img.taocdn.com/s3/m/5651c7c676eeaeaad1f330c4.png)
y (0) A B 1 / 3 1
y '(0) 2 A 4 B 1 / 3 1 2 4 t 1 t 2 t y ( t ) 2e e e , t 0 3 3
解: (3) 求方程的全解
y (t ) y h (t ) y p (t ) Ae
2 t
Be
4 t
A= 11/4,B=-7/4 1 y '(0) 2 A 4 B 2 2 11 7 1 y ( t ) e 2 t e 4 t te 2 t , t 0 4 4 2
[例] 已知描述某连续时间LTI系统的微分方程 y" (t ) 6 y' (t ) 8 y(t ) x(t ), t 0 初始条件y(0)=1, y' (0)=2, 输入信号x(t)=e-2t u(t),求全解y(t)。
解: (2) 求方程 y''(t)+6y'(t)+8y(t) = x(t)的特解yp(t)
2.基于零输入响应和零状态响应的方法
系统响应 y(t ) = 零输入响应 yzi (t )+ 零状态响应yzs (t ) 求解齐次微分方程
y zs ( t ) x ( t ) h ( t )
※ 零输入响应求解
※ 冲激响应的求解 ※ 零状态响应求解
当输入x(t)=0
y( n ) (t ) an1 y( n1) (t ) a1 y '(t ) a0 y(t ) 0
上式为齐次方程,因此零输入响应具有齐次解的形式
MATLAB与信号实验——连续LTI系统的时域分析
![MATLAB与信号实验——连续LTI系统的时域分析](https://img.taocdn.com/s3/m/04f1c686ba4cf7ec4afe04a1b0717fd5370cb254.png)
MATLAB与信号实验-——-连续LTI系统的时域分析在信号处理中,MATLAB是一个强大的工具,它提供了许多功能,使我们能够模拟和分析各种信号系统。
对于连续LTI系统,时域分析是一个重要的方法,它允许我们直接观察系统的输入和输出信号之间的关系。
下面是一个关于连续LTI系统的时域分析的实验。
一、实验目的本实验的目的是验证连续LTI系统的时域响应,通过使用MATLAB模拟系统,我们可以观察到不同的输入信号产生的输出信号,从而了解系统的特性。
二、实验步骤1.定义系统:首先,我们需要定义我们的连续LTI系统。
这可以通过使用MATLAB中的lti函数来完成。
我们需要提供系统的传递函数,它描述了系统的输入和输出之间的关系。
2.设置输入信号:为了观察系统的行为,我们需要设置一个合适的输入信号。
在MATLAB中,我们可以使用square函数来生成一个方波信号,该信号具有固定的频率和幅度。
3.模拟系统:使用MATLAB的lsim函数,我们可以模拟我们的连续LTI系统。
这个函数将输入信号和系统的传递函数作为参数,然后计算出系统的输出信号。
4.分析结果:我们可以使用MATLAB的图形功能来观察输入和输出信号。
这可以帮助我们理解系统的行为,并验证我们的模型是否正确。
三、实验结果与分析在实验中,我们使用了不同的输入信号(如方波、正弦波等)来测试我们的连续LTI系统。
对于每种输入信号,我们都观察了系统的输出信号,并记录了结果。
通过对比不同的输入和输出信号,我们可以得出以下结论:1.对于方波输入,系统的输出信号是带有延迟的方波,这表明系统对突变信号的响应是瞬时的。
2.对于正弦波输入,系统的输出信号是与输入信号同频同相位的正弦波,这表明系统对正弦波的响应是具有稳定性的。
这些结果验证了连续LTI系统的基本特性:即对于单位阶跃函数(突变信号)的输入,系统的响应是瞬时的;而对于周期性输入(如正弦波),系统的响应具有稳定性。
这些结果与我们在理论上学到的知识相符,从而验证了我们的模型是正确的。
实验三 线性时不变(LTI)连续系统的时域分析
![实验三 线性时不变(LTI)连续系统的时域分析](https://img.taocdn.com/s3/m/8a95bcd4240c844769eaeefd.png)
执行结果
实验任务 1:LTI系统的微分方程y''(t)+2y'(t)+y(t)=f'(t)+2f(t),激励f (t)=e-2tε(t),
(1) 利用 impulse 函数获得冲激响应; (2) 利用 lsim 函数求取零状态响应; (3) 用卷积分析法计算其零状态响应; 要求:在一个图形窗口里以 3 个子图形式绘制冲激响应和两种方法得到的零状 态响应的波形。 (4)改变系统的 a 系数矩阵,观察冲激响应和零状态响应时域波形的变化情况。 建议 a 系数向量分别如下取值讨论。
用。
一、实验目的
1. 掌握系统时域分析常用函数的使用方法; 2. 理解系统特征根对系统时域特性的影响;
二、实验原理及内容 2.1 连续系统的时域分析 2.1.1 连续系统时域分析的几个常用函数
设 LTI 连续时间系统的微分方程为
a 2 y''(t)+a 1 y'(t)+a 0 y(t)=b 2 f''(t)+b 1 f'(t)+b 0 f(t)
将积分变量离散化即将用n替代d用替代只要时域取样间隔足够小上式可近似为再把观察响应时刻离散化即将t用k替换只要足够小通常将fn简记为fnhkn简记为hknyk简记为yk这样上式便可表示为因此两个连续信号ft和ht的卷积yt可用ft和ht的取样信号f均取整数
实验三 线性时不变(LTI)连续系统的时域分析
2.1.2 连续信号卷积的近似计算
连续系统的零状态响应 y(t)可通过输入信号 f(t)与系统冲激响应 h(t)的卷积求 得;但是计算机只能处理数字信号,不能直接处理模拟信号,因此,卷积积分不 能直接用计算机计算。为了解决这个问题,可以将连续信号用取样信号来近似表 示,利用卷积和近似求得卷积积分。下面就连续信号卷积积分的近似计算进行简 单推导。
实验三 LTI连续系统的时域分析
![实验三 LTI连续系统的时域分析](https://img.taocdn.com/s3/m/bf1e2d6e7e21af45b307a85a.png)
实验三线性时不变连续系统的时域分析一、实验目的1.加深对线性时不变系统中零状态响应概念的理解,掌握其求解方法;2.掌握求解给定连续系统的冲激响应和阶跃响应的方法。
二、实验原理1. 线性时不变(LTI) 连续时间系统用常系数线性微分方程进行描述,系统的零状态响应就是在系统初始状态为零条件下微分方程的解。
MATLAB控制系统工具箱提供了一个lsim函数来求解连续时间系统的零状态响应。
设系统方程为:,该方程左边、右边的系数向量分别为,,所对应的系统模型sys可借助MATLAB中的tf函数得到:sys = tf(b, a) .这样,系统的零状态响应为:y = lsim(sys, f, t) ,其中f是输入信号向量,t是与f对应的时间变量。
2. 连续系统的冲激响应、阶跃响应分别是输入信号为和所对应的零状态响应。
MATLAB控制系统工具箱专门提供了两个函数求解连续系统的冲激响应和阶跃响应。
冲激响应:y = impulse(sys, t) ;阶跃响应:y = step(sys, t) .其中sys, t的含义同上。
三、程序示例示例1:已知系统的微分方程为,。
求零状态响应。
a = [1 4 4]; %将y(t)各阶导数的系数放在向量a中b = [1 3]; %将f(t)各阶导数的系数放在向量b中sys = tf(b, a); %求系统函数systd = 0.01; %定义时间间隔t = 0 : td : 10; %定义时间向量f = exp(-t); %将f(t)表示出来y = lsim(sys, f, t); %求系统的零状态响应yplot(t, y); %绘出零状态响应的波形xlabel('t(sec)'); % 这行代码是给出x坐标的标签ylabel('y(t)'); % 这行代码是给出y坐标的标签grid on % 这行代码是在图上显示方格程序运行结果见下图。
示例2:利用impulse 函数和step 函数,求连续系统冲激响应和阶跃响应。
信号与系统仿真—连续时间LTI系统的时域分析
![信号与系统仿真—连续时间LTI系统的时域分析](https://img.taocdn.com/s3/m/cc7f68de31126edb6e1a1066.png)
电子信息工程系实验报告成绩:课程名称:计算机仿真技术指导教师(签名):实验项目名称:实验4 信号与系统仿真—连续时间LTI系统的时域分析实验时间:2011-11-22 班级:电信092 姓名:XXX 学号:910706201一、实验目的:学会运用MATLAB符号求解连续系统的零输入响应和零状态响应;学会运用MATLAB数值求解连续系统的零状态响应;学会运用MATLAB求解连续系统的冲激响应和阶跃响应;思考运用MATLAB卷积积分法求解系统的零状态响应。
二、实验环境:硬件:PC机,PII以上 CPU,内存1G;软件:Matlab7.1三、实验原理:连续时间系统零状态响应的数值求解:对于零状态响应,MATLAB控制系统工具箱提供了对LTI系统的零状态响应进行数值仿真的函数lsim,该函数可求解零初始条件下微分方程的数值解,语句格式为:y=lsim(sys,f,t)其中,t表示计算系统响应的时间抽样点向量;f是系统的输入信号向量;sys表示LTI系统模型,用来表示微分方程、差分方程或者状态方程。
在求微分方程时,sys是由MATLAB的tf函数根据微分方程系数生成的系统函数对象,其语句格式为:sys = tf(b,a)其中,b和a分别为微分方程右端和左端的系数向量。
例如,对于微分方程可用;;获得其LTI模型。
注意,如果微分方程的左端或右端表达式中有缺项,则其向量a 或者b中对应的元素应为0,不能省略不写,否则会出错。
连续时间系统冲激响应和阶跃响应的求解:在MATLAB中,对于连续LTI系统的冲激响应和阶跃响应的数值解,可分别用控制系统工具箱提供的函数impulse和step来求解。
其语句格式分别为:y = impulse(sys,t);y = step(sys,t);其中,t表示计算机系统响应的时间抽样点向量,sys表示LTI系统模型。
四、实验内容及过程:1.1、已知系统的微分方程和激励信号如下,试用MATLAB命令绘出系统零状态响应的时域仿真波形图。
连续时间系统的时域分析实验报告
![连续时间系统的时域分析实验报告](https://img.taocdn.com/s3/m/035172735f0e7cd1842536a1.png)
实验二连续时间系统的时域分析一、实验目的通过使用MATLAB 软件对连续时间线性非时变系统的时域特性进行仿真分析,熟悉IT 系统在典型激励下的响应及特征,熟悉相应MATLAB 函数的调用格式和作用,熟悉井掌握用MATLAB 函数求解冲激响应、阶跃响应、零输入响应、零状态响应及全响应的方法。
二、实验原理(一)连续时间系统的时域分析方法 连续时间线性非时变系统(LTI )的输入()t f 与输出()t y 可以用线性常系数微分方程来描述:()()()()()()()()()()t f b t f b t fb t y a t y a t y a t y a m m n n n n 0'10'111++=++++--如果已知系统的输入信号()t f 及系统的初始条件为()()()()()-----0,,0,0,01'''n y y y y ,就可以利用解析方法求出系统的响应。
线性系统的全响应由零输入响应分量和零状态响应分量组成。
零输入响应是指当输入为零时仅由t=0的初始条件产生的系统响应,零状态响应是当初始条件(在t=0)假定为零时仅由0≥t 时的输入产生的系统响应分量。
零输入响应(单极点时)为:()∑==+++=nk t k tn ttx k n e c ec ec ec t y 12121λλλλ f式中,n c c c 、、、 21为任意待定常数,由初始条件确定。
零状态响应为:()()()τττd t h f t y f -=⎰∞∞-此式是对任意输入()t f ,用单位冲激响应()t h 形式表示的零状态响应()t y f 的公式。
已知()t h 就可确定任意输入()t f 的零状态响应()t y f ,即系统对任意输入的响应都可以用单位冲激响应确定。
系统总响应为:()()()()()τττλd t h f ec t y t y t y tnj j f x j -+=+=⎰∑∞∞-=1对于高阶系统,手工计算非常繁琐。
连续时间LTI系统的时域分析报告
![连续时间LTI系统的时域分析报告](https://img.taocdn.com/s3/m/9dfb3d73192e45361166f52b.png)
一、课程设计题目:基于 MATLAB 的连续时间LTI 系统的时域分析二、基本要求:① 掌握连续时不变信号处理的基本概念、基本理论和基本方法; ② 学会 MATLAB 的使用,掌握 MATLAB 的程序设计方法; ③ 学会用 MATLAB 对信号进行分析和处理;④ 编程实现卷积积分或卷积和,零输入响应,零状态响应; ⑤ 撰写课程设计论文,用信号处理基本理论分析结果。
三、设计方法与步骤:一般的连续时间系统分析有以下几个步骤: ①求解系统的零输入响应; ②求解系统的零状态响应; ③求解系统的全响应; ④分析系统的卷积;⑤画出它们的图形. 下面以具体的微分方程为例说明利用MATLAB 软件分析系统的具体方法.1.连续时间系统的零输入响应描述n 阶线性时不变(LTI )连续系统的微分方程为:已知y 及各阶导数的初始值为y(0),y (1)(0),… y (n-1)(0), 求系统的零输入响应。
建模当LIT 系统的输入为零时,其零输入响应为微分方程的其次解(即令微分方程的等号右端为零),其形式为(设特征根均为单根)其中p 1,p 2,…,p n 是特征方程a 1λn +a 2λn-1+…+a n λ+a n =0的根,它们可以用root(a)语句求得。
各系数 由y 及其各阶导数的初始值来确定。
对此有………………………………………………………………………………………写成矩阵形式为: P 1n-1C 1+ P 2n-1C 2+…+ P n n-1C n =D n-1y 01121111n 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 dtdt dt dt dt -++-++⋅⋅⋅⋅⋅++=+⋅⋅⋅⋅++即V•C=Y0其解为:C=V\Y0式中V为德蒙矩阵,在matlab的特殊矩阵库中有vander。
以下面式子为例:y(0_)=1,y(0_)=5;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);endplot(t,y);gridxlabel('t') ;ylabel('y');title('零输入响应');程序运行结果:用这个通用程序来解一个三阶系统,运行此程序并输入a=[1,5,4] Y0=[1,5] dt=0.01 te=6结果如下图:根据图可以分析零输入响应,它的起始值与输入函数无关,只与它的初始状态值有关,其起始值等于y(0_)的值。
MATLAB 3.连续时间LTI系统的时域分析
![MATLAB 3.连续时间LTI系统的时域分析](https://img.taocdn.com/s3/m/d1663687d0d233d4b04e6907.png)
1
3.连续时间 LTI 系统的时域分析
作业题
2、某连续系统的输入为 e(t),输出为 y(t),系统的微分方程为:y ’’(t)+5y ’(6)+6y(t)=3f ’(t)+2f(t) 绘出该系统的单位冲激响应和单位阶跃响应时域图形 若f t = ������ −2������ ������(������),绘出系统的零状态响应时域图形
将激励信号和单位冲击响应进行卷积,由下图发现结果一致。
6
3.连续时间 LTI 系统的时域分析
作业题
4、如下图所示的电路中,已知 电压 y(t ) 作为系统输出:
R1 R2 R3 4() , L1 L2 1( H ) ,如果以电阻 R3 上
请画出系统的冲击响应 h(t) 请画出系统在激励e(������) = 12������(������)(v)作用下的零状态响应,要求写出系统微分方程的求 解步骤。 对所求结果进行验证。
零输入响应: eq1='D2y+4*Dy+4*y=0'; cond='Dy(-0.01)=1,y(-0.01)=1'; answer=dsolve(eq1,cond) answer = 103/(100*exp(1/50)*exp(2*t))+ (3*t)/(exp(1/50)*exp(2*t)) >> 全响应: eq1='D2y+4*Dy+4*y=Df+3*f'; eq2='f=exp(-t)*heaviside(t)'; cond='Dy(-0.01)=1,y(-0.01)=1'; answer=dsolve(eq1,eq2,cond); simplify(answer.y) ans = (3*t+2*exp(t+1/50)*heaviside(t)-2*exp(1/50)*heaviside(t)-t*exp(1/50)* heaviside(t) + 103/100)/exp(2*t + 1/50)
MATLAB仿真之连续时间LTI系统仿真和时域分析
![MATLAB仿真之连续时间LTI系统仿真和时域分析](https://img.taocdn.com/s3/m/0d8fa674182e453610661ed9ad51f01dc28157af.png)
r=conv(h,e);t=-10:l/a:10;PlOt(I.r);title('零状态响应r(t)'); xlabel('t');ylabel('r');零输入程序及仿真建模当UT系统的输入为零时,其零输入响应为微分方程的其次解(即令微分方程的等号右端为零),其形式为(设特征根均为单根)Mf)=GeR÷Ge网+••••+”'其中PbP2,,∙∙,Pn是特征方程alλn+a2λn-l+∙∙∙+anλ+an=O的根,它们可以用rool(八)语句求得。
各系数IIIy及其各阶导数的初始值来确定。
对此有G+G+•…+G=NOp l C l+p2C2+--+P ll C n=Dy0PFG+〃2”工+•…3Y写成矩阵形式为:PJC+IY1C J+∙∙∙+PJC=D"*PlPi - P nC*-∣C∣t-∣JU-IP∣Pi…P…V为范德蒙矩阵,在matlab的特殊矩阵库中有Vandero以下面式子为例:√(r)+5y(0+4y(r)=2∕(∕)-4∕(r)y(OJ=l,y(OJ=5:MAT1.AB程序:a=input(,输入分母系数a=[al,a2,...]=');n=length(八)-l;YO=inputC输入初始条件向量YO=[yO,DyO,D2yO,.p=roots(八);V=rot90(vander(p));c=V∖Y0';dt=inρut('dt=');te=inpιιt('te-);t=O:dt:te;y=zeros(1,length(t));fork=kny=y÷c(k)*exp(p(k)*t);endplot(t,y);gridon:xlabel(,t');ylabel('y');litle('零输入响应');程序运行结果:用这个通用程序来解一个三阶系统,运行此程序并输入a=[l,5,4]Y0=[l,5]dt=O.Olte=6结果如下列图:依据图可以分析零输入响应,它的起始值与输入函数无关,只与它的初始状态值有关,其起始值等于y(0_)的值。
实验三 连续时间LTI系统的时域分析
![实验三 连续时间LTI系统的时域分析](https://img.taocdn.com/s3/m/2e62a826647d27284b7351fd.png)
实验三 连续时间LTI 系统的时域分析 一、实验目的1.学会用MA TLAB 求解连续系统的零状态响应; 2. 学会用MATLAB 求解冲激响应及阶跃响应; 3.学会用MA TLAB 实现连续信号卷积的方法; 二、实验原理1.连续时间系统零状态响应的数值计算我们知道,LTI 连续系统可用如下所示的线性常系数微分方程来描述,()()0()()NMi j i j i j a yt b f t ===∑∑在MA TLAB 中,控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数lsim 。
其调用格式y=lsim(sys,f,t)式中,t 表示计算系统响应的抽样点向量,f 是系统输入信号向量,sys 是LTI 系统模型,用来表示微分方程,差分方程或状态方程。
其调用格式sys=tf(b,a)式中,b 和a 分别是微分方程的右端和左端系数向量。
例如,对于以下方程:''''''''''''32103210()()()()()()()()a y t a y t a y t a y t b f t b f t b f t b f t +++=+++可用32103210[,,,];[,,,];a a a a a b b b b b == (,)sys tf b a = 获得其LTI 模型。
注意,如果微分方程的左端或右端表达式中有缺项,则其向量a 或b 中的对应元素应为零,不能省略不写,否则出错。
例3-1 已知某LTI 系统的微分方程为 y’’(t)+ 2y’(t)+100y(t)=f(t)其中,'(0)(0)0,()10sin(2)y y f t t π===,求系统的输出y(t).解:显然,这是一个求系统零状态响应的问题。
其MATLAB 计算程序如下: ts=0;te=5;dt=0.01; sys=tf([1],[1,2,100]); t=ts:dt:te;f=10*sin(2*pi*t); y=lsim(sys,f,t); plot(t,y);xlabel('Time(sec)'); ylabel('y(t)');2.连续时间系统冲激响应和阶跃响应的求解在MATLAB 中,对于连续LTI 系统的冲激响应和阶跃响应,可分别用控制系统工具箱提供的函数impluse 和step 来求解。
信号与系统实验报告实验三连续时间LTI系统的频域分析
![信号与系统实验报告实验三连续时间LTI系统的频域分析](https://img.taocdn.com/s3/m/5cab651265ce050877321372.png)
实验三 连续时间LTI 系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义;2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法,理解具有不同频率响应特性的滤波器对信号的滤波作用;3、学习和掌握幅度特性、相位特性以及群延时的物理意义;4、掌握用MATLAB 语言进行系统频响特性分析的方法。
基本要求:掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法,深刻理解LTI 系统的频率响应特性的物理意义,理解滤波和滤波器的概念,掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。
二、实验原理及方法1 连续时间LTI 系统的频率响应所谓频率特性,也称为频率响应特性,简称频率响应(Frequency response ),是指系统在正弦信号激励下的稳态响应随频率变化的情况,包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。
上图中x(t)、y(t)分别为系统的时域激励信号和响应信号,h(t)是系统的单位冲激响应,它们三者之间的关系为:)(*)()(t h t x t y =,由傅里叶变换的时域卷积定理可得到:)()()(ωωωj H j X j Y =3.1或者: )()()(ωωωj X j Y j H =3.2)(ωj H 为系统的频域数学模型,它实际上就是系统的单位冲激响应h(t)的傅里叶变换。
即⎰∞∞--=dt et h j H tj ωω)()( 3.3由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换,如果h(t)是收敛的,或者说是绝对可积(Absolutly integrabel )的话,那么H(j ω)一定存在,而且H(j ω)通常是复数,因此,也可以表示成复数的不同表达形式。
在研究系统的频率响应时,更多的是把它表示成极坐标形式:)()()(ωϕωωj ej H j H = 3.4上式中,)j (ωH 称为幅度频率相应(Magnitude response ),反映信号经过系统之后,信号各频率分量的幅度发生变化的情况,)(ωϕ称为相位特性(Phase response ),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验三 连续时间LTI 系统的时域分析
一、实验目的
1、学会使用符号法求解连续系统的零输入响应和零状态响应
2、学会使用数值法求解连续系统的零状态响应
3、学会求解连续系统的冲激响应和阶跃响应
二、实验原理及实例分析
1、连续时间系统零输入响应和零状态响应的符号求解
连续时间系统可以使用常系数微分方程来描述,其完全响应由零输入响应和零状态响应组成。
MATLAB 符号工具箱提供了dsolve 函数,可以实现对常系数微分方程的符号求解,其调用格式为:
dsolve(‘eq1,eq2…’,’cond1,cond2,…’,’v’)
其中参数eq 表示各个微分方程,它与MATLAB 符号表达式的输入基本相同,微分和导数的输入是使用Dy ,D2y ,D3y 来表示y 的一价导数,二阶导数,三阶导数;参数cond 表示初始条件或者起始条件;参数v 表示自变量,默认是变量t 。
通过使用dsolve 函数可以求出系统微分方程的零输入响应和零状态响应,进而求出完全响应。
[实例1]试用Matlab 命令求齐次微分方程0)()(2)(='+''+'''t y t y t y 的零输入响应,已知起始条件为2)0(,1)0(,1)0(=''='=---y y y 。
3、连续时间系统冲激响应和阶跃响应的求解
在连续时间LTI系统中,冲激响应和阶跃响应是系统特性的描述。
在MATLAB中,对于冲激响应和阶跃响应的数值求解,可以使用控制工具箱中提供的函数impulse和step来求解。
)
, (
) ,
(
t
sys
step
y t
sys
impulse
y
=
=
其中t表示系统响应的时间抽样点向量,sys表示LTI系统模型。
三、实验内容
1、已知系统的微分方程和激励信号,使用MATLAB命令画出系统的零状态响应和零输入响应(零状态响应分别使用符号法和数值法求解,零输入响应只使用符号法求解)。
要求题目2必做,题目1选做。
(1)
2、已知系统的微分方程,使用MATLAB命令画出系统的冲激响应和阶跃响应(数值法)。
要求题目2必做,题目1选做。
(1)
四、实验报告要求
1、格式:实验名称、实验目的、实验原理、实验环境、实验内容(上述几部分代码及结果图形)、实验思考等。