连续时间系统的冲激响应和零状态响应

2024-05-10 03:14:53

实验原理

线性时不变系统用微分方程描述：
y(n)(t)+an-1y(n-1)(t)+…+a1y(1)(t)+a0y(t)=bmf(m)(t)+bm-1f(m-1)(t)+…+b1f(1)(t)+b0f(t) (2-1)

（1）系统的零状态响应可通过求解初始状态为零的微分方程得到。

在MATLAB中，用lsim函数求解LTI的零状态响应，其调用形式为： y=lsim(sys,f,t) (2-2)
其中，t表示计算系统响应的抽样点向量，f是系统输入信号向量，sys是LTI系统模型，用来表示微分方程，差分方程，状态方程。在求解微分方程时，微分方程的LTI系统模型sys要借助MATLAB中的tf函数来获得，其调用形式为： sys=tf(b,a) (2-3)
式中，向量a表示响应及其各阶导数的系数：a=[ an an-1 … a1 a0],一般an=1；
向量b表示激励及其各阶导数的系数： b=[ bm bm-1 … b1 b0]
t表示计算响应所需的时间向量。例如t＝c∶d∶e,其中c是起始时间,e是结束时间,d是时间增量。为了做出平滑的曲线,向量t的增量需要取得足够小。需要注意的是,这里向量a和b的样本值都是按照导数阶数递减的次序排列的。
也可以把式（3-2）、（3-3）合在一起用，即调用格式为： y=lsim(b,a,f,t) (2-4)

（2）系统的单位冲激响应是单位冲激信号作用在系统产生的零状态响应。

在MATLAB中，用impulse函数求解LTI的冲激响应，调用形式为： y=impulse (b,a,t)

例题

微分方程为是
，输入激励为
，求其零状态响应和单位冲激响应,并与用解析式表示的结果进行比较。
用MATLAB函数求解冲激响应的程序如下：
clc;clear
t=[0:0.1:10];
b=[1];
a=[1 3];
x=exp(-1t);
yf=lsim(b,a,x,t);
yf1=(1/2)(exp(-1t)-exp(-3t));
plot(t,yf,‘o’,t,yf1,‘r’);
title(‘零状态响应’);
xlabel(‘t’);ylabel(‘yf(t)’);
legend(‘lsim函数计算结果’,‘解析式计算结果’);
h=impulse(b,a,t);
h1=exp(-3*t);
figure;
plot(t,h,‘o’,t,h1,‘r’);
title(‘单位冲激响应’);
xlabel(‘t’);ylabel(‘h(t)’);
legend(‘impulse函数计算结果’,‘解析式计算结果’);

连续时间系统的冲激响应和零状态响应相关推荐

Matlab计算信号和的卷积，求解微分方程的阶跃响应、冲激响应和零状态响应
参考资料: 连续时间系统的时域分析涉及到的Matlab函数: step:用于计算连续时间系统的单位阶跃响应. impulse:用于计算连续时间系统的单位冲激响应. lsim:用于计算连续系统在任意输入 ...
《信号与系统》连续时间系统零状态响应的 MATLAB 实现
3.5.1 连续时间系统零状态响应的 MATLAB 实现参考书籍 <信号与系统> Matlab 库函数中的 **lsim()**能对微分方程描述的 LTI 连续时间系统的响应进行仿真. ...
使用Matlab求解连续时间系统的零状态响应
连续时间LTI系统以常系数微分方程描述, 系统的零状态响应可通过求解初始状态为零的微分方程得到. 在Matlab中, 控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数lsim. 其调用方 ...
MATLAB：零状态响应(lsim(连续)；filter(离散))、冲激响应(impulse或impz)和阶跃响应(step)、卷积(conv)
例1:求系统y''(t)+2y'(t)+100y(t)=10x(t)的零状态响应,已知x(t)=sin(2pt)u(t).[连续时间系统零状态响应] 例2:求系统y''(t)+2y'(t)+100y( ...
卷积法求解系统的零状态响应_因果系统的零状态响应的一种简易计算方法
年月宜春师专学报年第期 . 因果系统的零状态响应的一种简易计算方法王卫林连续线性系统的研究方法 , 可以按系统的数学模型的求解方式分为时域法和变换域法两大类时域法是直接处理系统的数学模型 ...
matlab求解零输入响应,基于matlab的零输入和零状态响应
基于matlab的零输入和零状态响应一. 课程设计题目: 信号系统的时域分析二.课程设计目的: 1.学习MATLAB 软件的使用. 2.使学生掌握利用工具软件来实现信号系统基本概念.基本原理的方 ...
MATLAB——求系统的零状态响应
题目: 已知系统的微分方程为: d 2 y ( t ) d t 2 + 2 d y ( t ) d t + 100 y ( t ) = f ( t ) \dfrac{d^2y(t)}{dt^2}+2\ ...
卷积法求解系统的零状态响应_求解LTI系统零状态响应过程中的解法辨析
求解 LTI 系统零状态响应过程中的解法辨析马丁 [摘要] 摘要:利用卷积计算 LTI 系统的零状态响应时,单位冲激响应 h ( t ) 的求解是其中至关重要的一步.但用冲激平衡法及傅里叶反变换 ...
卷积法求解系统的零状态响应_连续LTI系统零状态响应求解方法的分析
连续 LTI 系统零状态响应求解方法的分析张淑敏 [摘要] [摘要]零状态响应是电子技术相关课程中的一个重要概念,本文将通过对时域分析法和(复)频域分析法求解连续 LTI 系统零状态响应的分析 ...

最新文章

热门文章