浙江理工大学信号与系统实验报告
实验二 连续时间系统的时域分析
一、实验目的
-
运用 MATLAB 数值求解连续时间系统的冲激响应和阶跃响应;
-
运用 MATLAB 数值求解连续时间系统的零输入响应、零状态响应及全响应;
-
运用 MATLAB 卷积积分法求解连续时间系统的零状态响应。
二、实验仪器
装有 MATLAB 软件的微型计算机 1 台
三、实验原理
对于连续的 LTI 系统,当系统输入为 e(t),输出为 r(t),则输入与输出之间满足如下的线性常系数微分方程:
则可以用向量 a 和 b 来表示该系统,即:
注意,在用向量来表示微分方程描述的连续系统时,向量 a 和 b 的元素一定要以微分方程时间求导的降幂次序来排列,且缺项用 0 来补齐。如对微分方程 2y’‘(t)+2y’(t)+6y(t)=f (t),则表示该系统的对应向量为:
a=[2 2 6]; b=[1]
而对微分方程:y’‘(t)+3y’(t)+2y(t)=f ‘’(t)+f (t),则表示该系统的对应向量应为:
a=[1 3 2]; b=[1 0 1]
当系统输入为单位冲激信号 δ(t)时产生的零状态响应称为系统的单位冲激响应,用 h(t) 表示。若输入为单位阶跃信号 ε(t)时,系统产生的零状态响应则称为系统的单位阶跃响应, 记为 g(t)。系统的单位冲激响应 h(t)包含了系统的固有特性,它是由系统本身的结构及参数所决定的,与系统的输入无关。我们只要知道了系统的冲激响应,即可求得系统在不同激励下产生的响应。因此,求解系统的冲激响应 h(t)对我们进行连续系统的分析具有非常重要的意义。
在 MATLAB 中有专门用于求解连续系统冲激响应和阶跃响应,并绘制其时域波形的函数 impulse( ) 和 step( )。
如果系统输入为 e(t),冲激响应为 h(t),系统的零状态响应为 r(t),则有:
若已知系统的输入信号及初始状态,我们便可以用微分方程的经典时域求解方法,求出系统的响应。但是对于高阶系统,手工计算这一问题的过程非常困难和繁琐。
在 MATLAB 中,应用 lsim( )函数很容易就能对上述微分方程所描述的系统的响应进行仿真,求出系统在任意激励信号作用下的响