信号与系统matlab学习（二）

 连续系统的响应求解

由连续时间系统的时域和频域分析方法可知，线性时不变系统的微分方程式为：

    对于复杂信号激励下的线性系统，可以将激励信号在时域中分解为单位脉冲信号或单位阶跃信号，把这些单位激励信号分别加在系统中求其响应，然后把这些响应叠加，即可得到复杂信号加在系统中的零状态响应。

一、实验用函数

1、impulse

功能：求解连续系统的冲激响应。

调用格式：

impulse(b,a)：计算并显示出连续系统的冲激响应h(t)的波形。

impulse(b,a,t)：当用户指定t值时，t为实数，将显示连续时间系统在0~t秒间的冲激响应波形，t为向量(t1:dt:t2)，则显示指定时间范围内的冲激响应波形。

2、step

功能：求解连续系统的阶跃响应。

调用格式：

step(b,a)：计算并显示出连续系统的阶跃响应g(t)的波形。

step(b,a,t)：当用户指定t值时，t为实数，将显示连续时间系统在0~t秒间的阶跃响应波形，t为向量(t1:dt:t2)，则显示指定时间范围内的阶跃响应波形。

3、lsim

功能：求解连续系统的零状态响应。

调用格式：

lsim(b,a,x,t)：计算并显示出连续系统的零状态响应的波形。

二、实例

1.已知一个RLC串联振荡电路系统函数，其中L=22mH，C=2000pF，R=100欧姆，求其时域的冲激响应和阶跃响应。

%impulse和step的应用
L = 22*10^-3;
C = 2000*10^-12;
R = 100;
a = [L*C,R*C,1];
b = [1];
t = 0:0.001:10;
subplot(121);impulse(b,a);
subplot(122);step(b,a);

2.用仿真函数lsim计算连续系统的零状态响应。已知系统的微分方程和系统的输入信号，求其零状态响应。

clear all;
a = [1,2,100];
b = 1;
t = 0:0.01:10;
ft = 10*sin(2*pi*t);
y = lsim(b,a,ft,t);
plot(t,y);

3.已知某连续时间系统的微分方程和输入信号，求该系统的零状态响应r(t)

a = [1,2,1];
b = [1,2];
t = 0:0.01:10;
et = exp(-2*t);
y = lsim(b,a,et,t);
plot(y);