拉普拉斯变换
物理问题引入:为什么初值动力学问题需要拉普拉斯变换?
在处理力学振动和电网络时,我们常常面临两类实际困难:
- 真实物理实验通常有明确的「启动时刻」:在 𝑡 <0t < 0 时系统处于静止,在 𝑡 =0t = 0 时闭合开关或受到外力撞击,系统带有初始位移 𝑥(0) =𝑥0x(0) = x_0 和初始速度 ˙𝑥(0) =𝑣0x(0) = v_0;
- 傅里叶变换收敛门槛高:如果信号包含阶跃函数、常数甚至线性发散趋势(如直流电、失稳发散振动),普通傅里叶积分 ∫∞−∞𝑓(𝑡)𝑒−𝑖𝜔𝑡d𝑡int_-inf^inf f(t)e^-iomega tdt 往往直接发散无法计算.
皮埃尔 - 西蒙·拉普拉斯(Pierre-Simon Laplace) 引入了带指数衰减抑制因子的复积分:通过令复频率 𝑠 =𝜎 +𝑖𝜔s = sigma + iomega,将 微分方程的初始条件在变换过程中直接自动代入,成为求解常微分方程初值问题的最强代数利器.
1. 单边拉普拉斯变换定义
设因果函数 𝑓(𝑡)f(t) 在 𝑡 ≥0t >= 0 有定义且在 𝑡 <0t < 0 时为零,其 单边拉普拉斯变换 定义为:
𝐹(𝑠)=L[𝑓(𝑡)]=∫∞0𝑓(𝑡)𝑒−𝑠𝑡d𝑡(𝑠=𝜎+𝑖𝜔∈ℂ)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(𝑠) =𝜎Re(s) = sigma 足够大,衰减因子 𝑒−𝜎𝑡e^-sigma t 就能压制绝大多数增长函数(如 𝑡𝑛,𝑒𝑎𝑡t^n, e^at)的发散,使积分在复半平面 Re(𝑠) >𝜎0Re(s) > sigma_0 内严格收敛.
2. 皇冠性质:初始条件自动代入
拉普拉斯变换最令人惊叹的特性,在于对导数进行变换时,初始状态会被天然纳入代数方程:
一阶导数变换
利用分部积分法:
L[𝑓′(𝑡)]=∫∞0𝑓′(𝑡)𝑒−𝑠𝑡d𝑡=𝑓(𝑡)𝑒−𝑠𝑡∣∞0−∫∞0𝑓(𝑡)(−𝑠𝑒−𝑠𝑡)d𝑡=𝑠𝐹(𝑠)−𝑓(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[𝑓″(𝑡)]=𝑠2𝐹(𝑠)−𝑠𝑓(0−)−𝑓′(0−)L[f''(t)] = s^2 F(s) - s f(0^-) - f'(0^-)对比传统 ODE 解法
在传统微积分中,求解二阶非齐次常微分方程需要经历:求齐次通解 →-> 猜特解 →-> 代入求和 →-> 解关于待定常数的线性方程组. 而使用拉普拉斯变换,初始条件直接在第一步就变成了代数方程中的已知数!求解只需解一个一次代数方程,再通过查表或留数求逆变换即可一步到位.
3. 常见物理函数的拉普拉斯变换对
| 时域函数 𝑓(𝑡) (𝑡 ≥0)f(t) , (t >= 0) | 象函数 𝐹(𝑠)F(s) | 收敛域 | 对应物理模型 |
|---|
| 单位脉冲 𝛿(𝑡)delta(t) | 11 | 全复平面 | 瞬时瞬间冲量碰撞 |
| 单位阶跃 Θ(𝑡)Theta(t) | 1𝑠1s | Re(𝑠) >0Re(s) > 0 | 𝑡 =0t=0 时合上直流电源开关 |
| 指数衰减 𝑒−𝑎𝑡e^-at | 1𝑠+𝑎1s + a | Re(𝑠) > −𝑎Re(s) > -a | 放射性衰变、RC 放电响应 |
| 正弦振荡 sin(𝜔𝑡)sin(omega t) | 𝜔𝑠2+𝜔2s^2 + omega^2 | Re(𝑠) >0Re(s) > 0 | 无阻尼简谐振荡 |
| 余弦振荡 cos(𝜔𝑡)cos(omega t) | 𝑠𝑠2+𝜔2ss^2 + omega^2 | Re(𝑠) >0Re(s) > 0 | 简谐初位移自由释放 |
| 阻尼正弦 𝑒−𝛾𝑡sin(𝜔𝑡)e^-gamma tsin(omega t) | 𝜔(𝑠+𝛾)2+𝜔2(s+gamma)^2 + omega^2 | Re(𝑠) > −𝛾Re(s) > -gamma | 阻尼受迫振子与欠阻尼振荡 |
4. 物理实战:RLC 暂态放电响应求解
考虑由电感 𝐿L、电阻 𝑅R、电容 𝐶C 构成的回路,在 𝑡 =0t=0 时电容带电量为 𝑄0Q_0,初始电流为 00. 根据基尔霍夫电压定律(KVL):
𝐿d2𝑞d𝑡2+𝑅d𝑞d𝑡+1𝐶𝑞=0Ld^2 qdt^2 + Rdqdt + 1Cq = 0两边同时做拉普拉斯变换,代入初始条件 𝑞(0) =𝑄0,˙𝑞(0) =0q(0) = Q_0, q(0) = 0:
𝐿[𝑠2𝑄(𝑠)−𝑠𝑄0]+𝑅[𝑠𝑄(𝑠)−𝑄0]+1𝐶𝑄(𝑠)=0L[s^2 Q(s) - s Q_0] + R[s Q(s) - Q_0] + 1CQ(s) = 0直接提取 𝑄(𝑠)Q(s):
𝑄(𝑠)(𝐿𝑠2+𝑅𝑠+1𝐶)=𝑄0(𝐿𝑠+𝑅)⟹𝑄(𝑠)=𝑄0𝑠+𝑅/𝐿𝑠2+𝑅𝐿𝑠+1𝐿𝐶Q(s) left( L s^2 + R s + 1C right) = Q_0 (Ls + R) implies Q(s) = Q_0 s + R/Ls^2 + RLs + 1LC分母二次式的零点(极点)直接决定了电路的三种物理机制:
- 两个互异实根:过阻尼放电(无振荡指数衰减);
- 一对共轭虚根:欠阻尼放电(衰减电磁振荡);
- 二重实根:临界阻尼放电(能量最快泄放).
5. 传递函数与复平面极点稳定性
定义系统的输入信号为 𝑋(𝑠)X(s),响应输出为 𝑌(𝑠)Y(s),两者的比值称为 传递函数(Transfer Function):
𝐻(𝑠)=𝑌(𝑠)𝑋(𝑠)H(s) = Y(s)X(s)- 𝐻(𝑠)H(s) 在复平面上的极点 𝑠𝑝 =𝜎𝑝 +𝑖𝜔𝑝s_p = sigma_p + iomega_p 决定了固有模态:𝑒𝑠𝑝𝑡 =𝑒𝜎𝑝𝑡𝑒𝑖𝜔𝑝𝑡e^s_p t = e^sigma_p t e^iomega_p t;
- 系统稳定的充要条件:𝐻(𝑠)H(s) 的所有极点必须严格位于 复平面的左半平面(即所有 Re(𝑠𝑝) <0Re(s_p) < 0),否则扰动将随时间指数爆炸失稳.
6. 学习衔接
本页面最近更新:2026/10/1 10:49:01,更新历史
发现错误?想一起完善? 在 GitHub 上编辑此页!
本页面贡献者:Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用