拉普拉斯变换与卷积

前面描述 卷积,本文由卷积引入拉普拉斯变换。

拉普拉斯变换就是给傅里叶变换的 iωt 加了个实部,也可以反着理解,原函数乘以 e − β t e^{-\beta t} e−βt 再做傅里叶变换,本质上都是傅里叶变换的扩展。

加入实部的拉普拉斯变换在复平面展开,因此拉普拉斯变换除了有虚部的频率特性,还有实部的衰减特性。但不管怎么说,和傅里叶变换一样,它们均将时域分析转换为频域分析,这就是我经常采用的策略,横竖一颠倒,晴空万里。

看一个频域的运算特性,为什么时域的卷积运算可以转换为频域的乘法运算?看下图:

时域中,信号持续时间和信号生成时间是相对的,必须卷起来运算,而频域是个对称的空间维度,不同频率信号以及针对每个频率信号的处理全都 "摆在那里",因此可通过精巧阵法一下子算出来(如上图矩阵只是一个说明,本质上还是复数的乘法)。

展示一张实操图:

以 g/G 往 f/F 涂颜料为例,上图的卷积操作需要一边卷紧纸一边在 f 曲线上色,这也体现了时间域的特征,而下图的乘积操作则是一个空间域并行操作,颜料涂在 G 半边的曲线,对折一次就能将 F 半边染色。

再看拉普拉斯变换对运算的实际简化效果。

设 u 为输入,g 为系统作用函数,时域上 u * g 的卷积获得输出 x,由上面的时域,频域转换分析获知,拉普拉斯变换将卷积化为了乘法,即 F*G,再经拉普拉斯逆变换就能获得 x(t) 的以下形式:

  • 若特征根没有虚部,要么是 x ( t ) = e a t ± e b t ± e c t . . . x(t)=e^{at}\pm e^{bt}\pm e^{ct}... x(t)=eat±ebt±ect...,要么收敛,要么发散,无震荡;
    • 若有虚部,欧拉公式,要么是 x ( t ) = e a t ( sin ⁡ ( b t ) . . . ) x(t)=e^{at}(\sin(bt)...) x(t)=eat(sin(bt)...),a > 0,震荡发散,a < 0,震荡收敛,a = 0,持续震荡。

发散的都不稳定,收敛的都稳定,区别在于是否震荡。

无需太过于关注极点的物理意义,拉普拉斯变换作为工具就只当它为工具,至于物理意义,一百个物理场景就能有一百个物理意义。

将极点看作逆变换后 x(t) 收敛,发散,震荡特征的总结是高尚的,而不是反过来。从 x(t) 的表达式可明确归纳出 a,b,c 等系数和函数收敛,发散,震荡的关系,而 a,b,c 等系数恰恰就是逆变换前 G(s) 分母为 0 的特征方程解的组成,将这些解在复平面的位置称作 "极点"。

因此根据 x(t) 的收敛,发散,震荡特性与极点实部虚部一一对应关系,通过分析 G(s) 的极点位置就能直接获得系统特性,无需再进行逆变换:

  • 有极点实部 > 0,发散不稳定,因为 e 的正次方总趋向无穷大;
  • 所有极点实部 < 0,收敛稳定,收敛速度取决于绝对值;
  • 极点虚部 ≠ 0,系统震荡,因为 x(t) 包含一个 sin/cos。

这就是拉普拉斯变换 G(s) 作为工具的意义,称 G 为传递函数。

浙江温州皮鞋湿,下雨进水不会胖。

相关推荐
00000cj1 天前
CoAtNet(NeurIPS 2023, Google)论文解读
人工智能·深度学习·计算机视觉·transformer·卷积
西西弗Sisyphus3 天前
卷积层的输出
深度学习·神经网络·计算机视觉·卷积
逼子格25 天前
31、matlab卷积运算:卷积运算、二维卷积、N维卷积
深度学习·神经网络·matlab·信号处理·卷积
微小冷4 个月前
scipy一维卷积函数convolve1d
数值计算·scipy·卷积·科学计算·ndimage·convolve1d·一维卷积
QomolangmaH4 个月前
【深度学习】Pytorch教程(八):PyTorch数据结构:2、张量的数学运算(6):高维张量:乘法、卷积(conv2d~四维张量;conv3d~五维张量)
数据结构·pytorch·python·深度学习·卷积
chairon5 个月前
basic CNN
pytorch·python·深度学习·cnn·卷积
暗然而日章6 个月前
信号与线性系统翻转课堂笔记6——卷积
笔记·信号与系统·卷积
智慧医疗探索者7 个月前
pytorch中Conv1d、Conv2d与Conv3d详解
人工智能·pytorch·卷积
此星光明8 个月前
GEE中核函数在不同缩放级别下的区别
开发语言·前端·javascript·gee·kernel·核函数·卷积