跳转至

拉普拉斯变换

物理问题引入:为什么初值动力学问题需要拉普拉斯变换?

在处理力学振动和电网络时,我们常常面临两类实际困难:

  1. 真实物理实验通常有明确的「启动时刻」:在 t<0t < 0 时系统处于静止,在 t=0t = 0 时闭合开关或受到外力撞击,系统带有初始位移 x(0)=x0x(0) = x_0 和初始速度 x˙(0)=v0x(0) = v_0;
  2. 傅里叶变换收敛门槛高:如果信号包含阶跃函数、常数甚至线性发散趋势(如直流电、失稳发散振动),普通傅里叶积分 ∫−∞∞f(t)e−iωtdtint_-inf^inf f(t)e^-iomega tdt 往往直接发散无法计算.

皮埃尔 - 西蒙·拉普拉斯(Pierre-Simon Laplace) 引入了带指数衰减抑制因子的复积分:通过令复频率 s=σ+iωs = sigma + iomega,将 微分方程的初始条件在变换过程中直接自动代入,成为求解常微分方程初值问题的最强代数利器.


1. 单边拉普拉斯变换定义

设因果函数 f(t)f(t) 在 t≥0t >= 0 有定义且在 t<0t < 0 时为零,其 单边拉普拉斯变换 定义为:

F(s)=L[f(t)]=∫0∞f(t)e−stdt(s=σ+iω∈C)F(s) = L[f(t)] = int_0^inf f(t) e^-s t dt quad (s = sigma + iomega in C)

收敛域(Region of Convergence, ROC)

只要实部 Re(s)=σRe(s) = sigma 足够大,衰减因子 e−σte^-sigma t 就能压制绝大多数增长函数(如 tn,eatt^n, e^at)的发散,使积分在复半平面 Re(s)>σ0Re(s) > sigma_0 内严格收敛.


2. 皇冠性质:初始条件自动代入

拉普拉斯变换最令人惊叹的特性,在于对导数进行变换时,初始状态会被天然纳入代数方程:

一阶导数变换

利用分部积分法:

L[f′(t)]=∫0∞f′(t)e−stdt=f(t)e−st|0∞−∫0∞f(t)(−se−st)dt=sF(s)−f(0−)L[f'(t)] = int_0^inf f'(t) e^-stdt = left. f(t)e^-st right|_0^inf - int_0^inf f(t)(-s e^-st)dt = s F(s) - f(0^-)

二阶导数变换

同理推广到二阶导数:

L[f″(t)]=s2F(s)−sf(0−)−f′(0−)L[f''(t)] = s^2 F(s) - s f(0^-) - f'(0^-)
对比传统 ODE 解法

在传统微积分中,求解二阶非齐次常微分方程需要经历:求齐次通解 →-> 猜特解 →-> 代入求和 →-> 解关于待定常数的线性方程组. 而使用拉普拉斯变换,初始条件直接在第一步就变成了代数方程中的已知数!求解只需解一个一次代数方程,再通过查表或留数求逆变换即可一步到位.


3. 常见物理函数的拉普拉斯变换对

时域函数 f(t)(t≥0)f(t) , (t >= 0)象函数 F(s)F(s)收敛域对应物理模型
单位脉冲 δ(t)delta(t)11全复平面瞬时瞬间冲量碰撞
单位阶跃 Θ(t)Theta(t)1s1sRe(s)>0Re(s) > 0t=0t=0 时合上直流电源开关
指数衰减 e−ate^-at1s+a1s + aRe(s)>−aRe(s) > -a放射性衰变、RC 放电响应
正弦振荡 sin⁡(ωt)sin(omega t)ωs2+ω2s^2 + omega^2Re(s)>0Re(s) > 0无阻尼简谐振荡
余弦振荡 cos⁡(ωt)cos(omega t)ss2+ω2ss^2 + omega^2Re(s)>0Re(s) > 0简谐初位移自由释放
阻尼正弦 e−γtsin⁡(ωt)e^-gamma tsin(omega t)ω(s+γ)2+ω2(s+gamma)^2 + omega^2Re(s)>−γRe(s) > -gamma阻尼受迫振子与欠阻尼振荡

4. 物理实战:RLC 暂态放电响应求解

考虑由电感 LL、电阻 RR、电容 CC 构成的回路,在 t=0t=0 时电容带电量为 Q0Q_0,初始电流为 00. 根据基尔霍夫电压定律(KVL):

Ld2qdt2+Rdqdt+1Cq=0Ld^2 qdt^2 + Rdqdt + 1Cq = 0

两边同时做拉普拉斯变换,代入初始条件 q(0)=Q0,q˙(0)=0q(0) = Q_0, q(0) = 0:

L[s2Q(s)−sQ0]+R[sQ(s)−Q0]+1CQ(s)=0L[s^2 Q(s) - s Q_0] + R[s Q(s) - Q_0] + 1CQ(s) = 0

直接提取 Q(s)Q(s):

Q(s)(Ls2+Rs+1C)=Q0(Ls+R)⟹Q(s)=Q0s+R/Ls2+RLs+1LCQ(s) left( L s^2 + R s + 1C right) = Q_0 (Ls + R) implies Q(s) = Q_0 s + R/Ls^2 + RLs + 1LC

分母二次式的零点(极点)直接决定了电路的三种物理机制:

  1. 两个互异实根:过阻尼放电(无振荡指数衰减);
  2. 一对共轭虚根:欠阻尼放电(衰减电磁振荡);
  3. 二重实根:临界阻尼放电(能量最快泄放).

5. 传递函数与复平面极点稳定性

定义系统的输入信号为 X(s)X(s),响应输出为 Y(s)Y(s),两者的比值称为 传递函数(Transfer Function):

H(s)=Y(s)X(s)H(s) = Y(s)X(s)
  • H(s)H(s) 在复平面上的极点 sp=σp+iωps_p = sigma_p + iomega_p 决定了固有模态:espt=eσpteiωpte^s_p t = e^sigma_p t e^iomega_p t;
  • 系统稳定的充要条件:H(s)H(s) 的所有极点必须严格位于 复平面的左半平面(即所有 Re(sp)<0Re(s_p) < 0),否则扰动将随时间指数爆炸失稳.

6. 学习衔接



评论