当前位置: 首页 > news >正文

从控制系统角度理解拉普拉斯卷积定理

一、什么是拉普拉斯卷积定理?

F_{1}(s)F_{2}(s)分别是函数f_{1}(t)f_{2}(t)的拉普拉斯变换形式,且f_{1}(t)f_{2}(t)满足当t< 0f_{1}(t)=0f_{2}(t)=0。那么有拉普拉斯卷积定理(Laplace convolution theorem):

直接用数学证明该公式往往比较抽象,本文将着力从控制系统的角度来解释该公式的合理性。

二、控制系统和传递函数

控制系统一般用微分方程(differential function)表示,如下所示:

但这个方程往往难以解析,所以往往对该微分方程两边同时作拉普拉斯变换(Laplace intergration transformation)。得到的式子如下所示:

其中,Y(s)是输出量(响应)的拉普拉斯变换形式, X(s)是输入量(激励)的拉普拉斯变换形式。G(s)为两者的比例,被称为传递函数(transfer function)。于是,当Y(s)求出后,输出量函数(响应函数)可以由拉普拉斯逆变换求出。

三、单位冲激响应

当输入函数为单位冲激函数\delta (t)(unit impulse function)时,由于其拉普拉斯变换形式为1,此时,Y(s)满足:

Y(s)=G(s)

也就是说单位冲激函数作为输入量(激励)时,此时系统的输出量(也被称为单位冲击响应,unit impulse response)的拉普拉斯变换形式和传递函数拉普拉斯变换形式相等。由于拉普拉斯变换与原函数存在唯一映射关系,可以说传递函数的原函数等于系统单位冲击响应

规定单位冲激响应为K(s),普通形式为k(t)

四、近似的单位冲激响应

单位冲激函数的有如下性质:

\int_{-\infty }^{\infty}\delta (t)dt=1

单位冲激函数在t=0处接近于无限大,而在别的地方都等于0,且函数在实数轴上积分为1。试想,我们在数轴上某一个比较小区间上,如[t_{1},t_{1}+\Delta t]使得函数f(x)等于某一个非零正常数,而在其他地方均等于0,且在实数轴上积分为1。那么,如果把这个函数作为系统的输入量,系统的响应会如何呢?

此时,可以如下表示输入函数:

x(t)=\frac{1}{\Delta t}[u(t_{1})-u(t_{1}+\Delta t)]

两边作拉普拉斯变换:

X(s)=\frac{1}{\Delta t}(\frac{e^{-st_{1}}}{s}-\frac{e^{-s(t_{1}+\Delta t)}}{s})

其实,当我们让\Delta t无限小时,我们可以把这个函数归结为指数函数e^{-st}t=t_{1}处的导数。X(s)就变成了冲激函数\delta (t-t_{1})的拉普拉斯变换形式,如下所示:

X(s)=e^{-st_{1}}

这样,如果\Delta t比较小,我们可以近似认为输出量和冲激函数作为激励的输出量一样。即X(s)造成的响应为K(s)e^{-st_{1}},普通形式为k(t-t_{1}),即单位冲激响应的变种。

我现在把这种函数f(x)称为近似的单位冲激函数这不是严谨的定义,而是我为了阐述作的规定

五、如何通过求一个不规则激励的响应?

如上图所示,可以把不规则的输入信号拆分为很多很多段持续时间很短的稳定输入量的叠加。具体的说,在0到\Delta t时,可以认为输入量为r(0);在\Delta t2\Delta t时,输入量为r(\Delta t)

我们刚才说,当开始时间为t_{1},持续时间\Delta t很小时,数轴上积分值为1的近似单位冲激函数造成的响应为K(s)e^{-st_{1}},普通形式为k(t-t_{1})

那么,数值上积分为m的近似单位冲激函数的响应为mk(t-t_{1})。现在,从t_{1}t_{1}+\Delta t的积分值为\Delta t·r(t_{1})。所以这段时间的响应表示为r(t_{1})k(t-t_{1})\Delta t

这样,对于这一不规则激励,我们把它进行分段,每一段的响应是可以近似求出的。

六、现在可以开始叠加了

对于一个控制系统而言,任何响应都是由激励造成了。如果没有激励,响应应该严格为0。所以,如果我想求t=t_{1}时的响应,那么我只需要关注[0,t_{1}]范围内的激励即可。将这一段时间平均分为n份(规定n很大),那么t=t_{1}时的响应应当是这n段时间激励造成响应的叠加,如下所示:

y(t_{1})=\sum_{i=1}^{n}r(\tau _{i})k(t_{1}-\tau_{i})\Delta \tau

\Delta \tau无限小时,可以改写为:

y(t_{1})=\int_{0}^{t_{1}} r(\tau)k(t_{1}-\tau)d\tau

也就是:

y(t_{1})=r(t_{1})*k(t_{1})

七、得到等式两边相等

如果激励为r(t),响应为y(t),且两者对应的拉普拉斯变换形式分别是R(s)Y(s),那么其满足

Y(s)=G(s)R(s)

而由第六节我们知道:

Y(s)=L[r(t)*k(t)]

这就有如下式子成立:

L[r(t)*k(t)]=G(s)R(s)=K(s)R(s)

注意,这里用到单位冲激响应的拉普拉斯变换等于传递函数


http://www.mrgr.cn/news/26317.html

相关文章:

  • 2024.9.13 系统运维
  • Java铸基之路:运算符的深入学习!(上)
  • SQL Server 语句日期格式查找方法
  • axure循环介绍
  • 五种IO模型与非阻塞IO
  • 代码随想录Day 42|leetcode题目:188.买卖股票的最佳时机IV、309.最佳买卖股票时机含冷冻期、714.买卖股票的最佳时机含手续费
  • 谷歌创始人谢尔盖·布林回归一线:承认错失先机,每天都在写代码
  • 【25保研经验帖】-更新中
  • C++笔记---继承(上)
  • Spring web mvc入门练习
  • Python 入门教程(3)基础知识 | 3.2、数据类型、变量与常量
  • 天下苦英伟达久矣!PyTorch官方免CUDA加速推理,Triton时代要来?
  • 基于Python实现一个庆祝中秋节的小程序
  • 密码学---黄道十二宫
  • MySQL 日志篇:Redo 文件和自适应检查点
  • 速通GPT-2:Language Models are Unsupervised Multitask Learners全文解读
  • mysql笔记8(多表查询)
  • 【【通信协议之ICMP协议】】
  • 质量体系和质量过程管理及SCIOT平台质量管理功能简介
  • Java项目: 基于SpringBoot+mybatis+maven宠物咖啡馆平台(含源码+数据库+毕业论文)