拉普拉斯变换:为什么它能把微分方程变简单?
第一次看拉普拉斯变换,很多人都会被一堆公式劝退。其实刚开始不用急着把公式表全背下来,先理解它到底在做什么。
简单说,拉普拉斯变换就是把原本在时间域里的函数 f(t)f(t)f(t),转换到另一个更方便计算的 sss 域里:
f(t)→F(s)f(t)\\rightarrow F(s)f(t)→F(s)
这样做最大的好处是:原来方程里麻烦的求导、积分,到了 sss 域里会变成更容易处理的代数运算。
所以这篇文章先不追求把所有公式都讲完,而是先抓住几个最核心的点:
f(t)→F(s)f(t)\\rightarrow F(s)f(t)→F(s)
F(s)=L[f(t)]=∫0+∞f(t)e−stdtF(s)=\\mathcal{L}[f(t)]=\\int_0^{+\\infty}f(t)e^{-st}dtF(s)=L[f(t)]=∫0+∞f(t)e−stdt
在零初始条件下,可以简单理解为:
微分→乘以 s\\text{微分}\\rightarrow \\text{乘以 }s微分→乘以 s
积分→除以 s\\text{积分}\\rightarrow \\text{除以 }s积分→除以 s
u(t),t,tn,e±at,sinωt,cosωtu(t),\\quad t,\\quad t^n,\\quad e^{\\pm at},\\quad \\sin\\omega t,\\quad \\cos\\omega tu(t),t,tn,e±at,sinωt,cosωt
其中正弦、余弦很重要,后面控制系统、电路分析、信号分析都会大量用到。
一、拉普拉斯变换是什么?
在时间域中,我们研究的是 f(t)f(t)f(t)。
经过拉普拉斯变换后,它会变成 F(s)F(s)F(s),记作:
F(s)=L[f(t)]F(s)=\\mathcal{L}[f(t)]F(s)=L[f(t)]
其中:
- ttt:时间变量;
- sss:复频域变量;
- f(t)f(t)f(t):原函数;
- F(s)F(s)F(s):拉氏变换后的函数,也叫像函数。
简单说就是:
时间域问题→s 域问题\\text{时间域问题}\\rightarrow \\text{s 域问题}时间域问题→s 域问题
这样做的好处是:很多复杂的微分、积分运算,在 sss 域中可以变成普通代数运算。
二、常用公式是怎么来的?
公式表看起来很多,但它不是凭空背出来的。
所有拉氏变换公式,本质上都来自定义:
F(s)=∫0+∞f(t)e−stdtF(s)=\\int_0^{+\\infty}f(t)e^{-st}dtF(s)=∫0+∞f(t)e−stdt
也就是把不同的 f(t)f(t)f(t) 代进去积分。
新手不需要每个公式都完整推导,但要理解几个典型公式。
1. 单位阶跃函数 u(t)u(t)u(t)
当 t≥0t\\ge0t≥0 时:
u(t)=1u(t)=1u(t)=1
所以:
L[u(t)]=∫0+∞e−stdt\\mathcal{L}[u(t)]=\\int_0^{+\\infty}e^{-st}dtL[u(t)]=∫0+∞e−stdt
积分得:
L[u(t)]=[−1se−st]0+∞\\mathcal{L}[u(t)]=\\left[-\\frac{1}{s}e^{-st}\\right]_0^{+\\infty}L[u(t)]=[−s1e−st]0+∞
当 t→+∞t\\to+\\inftyt→+∞ 时,e−st→0e^{-st}\\to0e−st→0,所以:
L[u(t)]=0−(−1s)=1s\\mathcal{L}[u(t)]=0-\\left(-\\frac{1}{s}\\right)=\\frac{1}{s}L[u(t)]=0−(−s1)=s1
因此:
L[u(t)]=1s\\mathcal{L}[u(t)]=\\frac{1}{s}L[u(t)]=s1
2. 指数函数 e−ate^{-at}e−at
代入定义:
L[e−at]=∫0+∞e−ate−stdt\\mathcal{L}[e^{-at}]=\\int_0^{+\\infty}e^{-at}e^{-st}dtL[e−at]=∫0+∞e−ate−stdt
指数相乘:
e−ate−st=e−(s+a)te^{-at}e^{-st}=e^{-(s+a)t}e−ate−st=e−(s+a)t
所以:
L[e−at]=∫0+∞e−(s+a)tdt=1s+a\\mathcal{L}[e^{-at}]=\\int_0^{+\\infty}e^{-(s+a)t}dt=\\frac{1}{s+a}L[e−at]=∫0+∞e−(s+a)tdt=s+a1
因此:
L[e−at]=1s+a\\mathcal{L}[e^{-at}]=\\frac{1}{s+a}L[e−at]=s+a1
同理:
L[eat]=1s−a\\mathcal{L}[e^{at}]=\\frac{1}{s-a}L[eat]=s−a1
这个地方容易记反:
e−at→1s+ae^{-at}\\rightarrow \\frac{1}{s+a}e−at→s+a1
eat→1s−ae^{at}\\rightarrow \\frac{1}{s-a}eat→s−a1
可以这样记:
指数里是负号,分母里变加号;指数里是正号,分母里变减号。
3. 一次函数 ttt
由定义:
L[t]=∫0+∞te−stdt\\mathcal{L}[t]=\\int_0^{+\\infty}te^{-st}dtL[t]=∫0+∞te−stdt
用分部积分可以得到:
L[t]=1s2\\mathcal{L}[t]=\\frac{1}{s^2}L[t]=s21
进一步推广:
L[tn]=n!sn+1\\mathcal{L}[t^n]=\\frac{n!}{s^{n+1}}L[tn]=sn+1n!
例如:
L[t2]=2s3\\mathcal{L}[t^2]=\\frac{2}{s^3}L[t2]=s32
所以:
L[12t2]=1s3\\mathcal{L}\\left[\\frac{1}{2}t^2\\right]=\\frac{1}{s^3}L[21t2]=s31
4. 正弦和余弦函数
正弦和余弦非常重要,可以用欧拉公式理解:
ejωt=cosωt+jsinωte^{j\\omega t}=\\cos\\omega t+j\\sin\\omega tejωt=cosωt+jsinωt
先求:
L[ejωt]\\mathcal{L}[e^{j\\omega t}]L[ejωt]
根据指数函数的结果:
L[ejωt]=1s−jω\\mathcal{L}[e^{j\\omega t}]=\\frac{1}{s-j\\omega}L[ejωt]=s−jω1
分母有理化:
1s−jω=s+jωs2+ω2\\frac{1}{s-j\\omega}=\\frac{s+j\\omega}{s^2+\\omega^2}s−jω1=s2+ω2s+jω
也就是:
1s−jω=ss2+ω2+jωs2+ω2\\frac{1}{s-j\\omega}=\\frac{s}{s^2+\\omega^2}+j\\frac{\\omega}{s^2+\\omega^2}s−jω1=s2+ω2s+js2+ω2ω
又因为:
ejωt=cosωt+jsinωte^{j\\omega t}=\\cos\\omega t+j\\sin\\omega tejωt=cosωt+jsinωt
所以实部对应余弦,虚部对应正弦:
L[cosωt]=ss2+ω2\\mathcal{L}[\\cos\\omega t]=\\frac{s}{s^2+\\omega^2}L[cosωt]=s2+ω2s
L[sinωt]=ωs2+ω2\\mathcal{L}[\\sin\\omega t]=\\frac{\\omega}{s^2+\\omega^2}L[sinωt]=s2+ω2ω
这两个公式一定要熟:
sinωt→ωs2+ω2\\sin\\omega t\\rightarrow \\frac{\\omega}{s^2+\\omega^2}sinωt→s2+ω2ω
cosωt→ss2+ω2\\cos\\omega t\\rightarrow \\frac{s}{s^2+\\omega^2}cosωt→s2+ω2s
三、新手必须掌握的核心表
刚开始不用背太多,先掌握这张表即可。
| δ(t)\\delta(t)δ(t) | 111 |
| u(t)u(t)u(t) 或 1(t)1(t)1(t) | 1s\\dfrac{1}{s}s1 |
| ttt | 1s2\\dfrac{1}{s^2}s21 |
| tnt^ntn | n!sn+1\\dfrac{n!}{s^{n+1}}sn+1n! |
| e−ate^{-at}e−at | 1s+a\\dfrac{1}{s+a}s+a1 |
| eate^{at}eat | 1s−a\\dfrac{1}{s-a}s−a1 |
| sinωt\\sin\\omega tsinωt | ωs2+ω2\\dfrac{\\omega}{s^2+\\omega^2}s2+ω2ω |
| cosωt\\cos\\omega tcosωt | ss2+ω2\\dfrac{s}{s^2+\\omega^2}s2+ω2s |
四、指数乘正弦、余弦:用位移性质推
后面经常会遇到:
e−atsinωte^{-at}\\sin\\omega te−atsinωt
e−atcosωte^{-at}\\cos\\omega te−atcosωt
这类公式不建议死背,用位移性质更清楚。
如果:
L[f(t)]=F(s)\\mathcal{L}[f(t)]=F(s)L[f(t)]=F(s)
那么:
L[e−atf(t)]=F(s+a)\\mathcal{L}[e^{-at}f(t)]=F(s+a)L[e−atf(t)]=F(s+a)
也就是说,乘上 e−ate^{-at}e−at,就把原公式里的 sss 换成 s+as+as+a。
例如:
L[sinωt]=ωs2+ω2\\mathcal{L}[\\sin\\omega t]=\\frac{\\omega}{s^2+\\omega^2}L[sinωt]=s2+ω2ω
所以:
L[e−atsinωt]=ω(s+a)2+ω2\\mathcal{L}[e^{-at}\\sin\\omega t]=\\frac{\\omega}{(s+a)^2+\\omega^2}L[e−atsinωt]=(s+a)2+ω2ω
同理:
L[e−atcosωt]=s+a(s+a)2+ω2\\mathcal{L}[e^{-at}\\cos\\omega t]=\\frac{s+a}{(s+a)^2+\\omega^2}L[e−atcosωt]=(s+a)2+ω2s+a
如果乘的是 eate^{at}eat,则把 sss 换成 s−as-as−a。
所以记住:
e−at⇒s 换成 s+ae^{-at} \\Rightarrow s\\text{ 换成 }s+ae−at⇒s 换成 s+a
eat⇒s 换成 s−ae^{at} \\Rightarrow s\\text{ 换成 }s-aeat⇒s 换成 s−a
五、拉氏变换最重要的性质
1. 线性性质
线性性质就是:
加减可以拆开,常数可以提出来\\text{加减可以拆开,常数可以提出来}加减可以拆开,常数可以提出来
如果:
L[f1(t)]=F1(s)\\mathcal{L}[f_1(t)]=F_1(s)L[f1(t)]=F1(s)
L[f2(t)]=F2(s)\\mathcal{L}[f_2(t)]=F_2(s)L[f2(t)]=F2(s)
则:
L[af1(t)+bf2(t)]=aF1(s)+bF2(s)\\mathcal{L}[af_1(t)+bf_2(t)]=aF_1(s)+bF_2(s)L[af1(t)+bf2(t)]=aF1(s)+bF2(s)
例如:
L[3t+2e−t]=3L[t]+2L[e−t]\\mathcal{L}[3t+2e^{-t}]=3\\mathcal{L}[t]+2\\mathcal{L}[e^{-t}]L[3t+2e−t]=3L[t]+2L[e−t]
所以:
L[3t+2e−t]=3s2+2s+1\\mathcal{L}[3t+2e^{-t}]=\\frac{3}{s^2}+\\frac{2}{s+1}L[3t+2e−t]=s23+s+12
2. 微分性质
如果这篇文章只能记住一个性质,那优先记住微分性质。
如果:
L[f(t)]=F(s)\\mathcal{L}[f(t)]=F(s)L[f(t)]=F(s)
则:
L[f′(t)]=sF(s)−f(0)\\mathcal{L}[f'(t)]=sF(s)-f(0)L[f′(t)]=sF(s)−f(0)
二阶导数为:
L[f′′(t)]=s2F(s)−sf(0)−f′(0)\\mathcal{L}[f''(t)]=s^2F(s)-sf(0)-f'(0)L[f′′(t)]=s2F(s)−sf(0)−f′(0)
如果初始条件全为 0,则:
L[f′(t)]=sF(s)\\mathcal{L}[f'(t)]=sF(s)L[f′(t)]=sF(s)
L[f′′(t)]=s2F(s)\\mathcal{L}[f''(t)]=s^2F(s)L[f′′(t)]=s2F(s)
这就是拉普拉斯变换最关键的用途:
微分运算→乘以 s 的代数运算\\text{微分运算}\\rightarrow\\text{乘以 }s\\text{ 的代数运算}微分运算→乘以 s 的代数运算
3. 积分性质
如果:
L[f(t)]=F(s)\\mathcal{L}[f(t)]=F(s)L[f(t)]=F(s)
则:
L[∫0tf(τ)dτ]=F(s)s\\mathcal{L}\\left[\\int_0^t f(\\tau)d\\tau\\right]=\\frac{F(s)}{s}L[∫0tf(τ)dτ]=sF(s)
简单记:
微分对应乘 s\\text{微分对应乘 }s微分对应乘 s
积分对应除以 s\\text{积分对应除以 }s积分对应除以 s
4. 延迟性质
如果:
L[f(t)]=F(s)\\mathcal{L}[f(t)]=F(s)L[f(t)]=F(s)
则:
L[f(t−τ)u(t−τ)]=e−τsF(s)\\mathcal{L}[f(t-\\tau)u(t-\\tau)]=e^{-\\tau s}F(s)L[f(t−τ)u(t−τ)]=e−τsF(s)
含义是:
函数延迟 τ\\tauτ 秒出现,在 sss 域中会多乘:
e−τse^{-\\tau s}e−τs
六、初值定理和终值定理
1. 初值定理
用于求 t→0+t\\to0^+t→0+ 时的值:
f(0+)=lims→+∞sF(s)f(0^+)=\\lim_{s\\to+\\infty}sF(s)f(0+)=s→+∞limsF(s)
2. 终值定理
用于求 t→∞t\\to\\inftyt→∞ 时的最终值:
limt→∞f(t)=lims→0sF(s)\\lim_{t\\to\\infty}f(t)=\\lim_{s\\to0}sF(s)t→∞limf(t)=s→0limsF(s)
但终值定理不能乱用。
如果函数最终不稳定,或者一直振荡,比如:
sint,cost\\sin t,\\quad \\cos tsint,cost
就没有最终值,不能用终值定理。
七、拉氏变换如何把微分方程变成代数方程?
以二阶微分方程为例:
y′′(t)+3y′(t)+2y(t)=0y''(t)+3y'(t)+2y(t)=0y′′(t)+3y′(t)+2y(t)=0
设:
Y(s)=L[y(t)]Y(s)=\\mathcal{L}[y(t)]Y(s)=L[y(t)]
对方程两边取拉氏变换:
L[y′′(t)]⏟s2Y(s)−sy(0)−y′(0)+3L[y′(t)]⏟sY(s)−y(0)+2L[y(t)]⏟Y(s)=0\\underbrace{\\mathcal{L}[y''(t)]}_{s^2Y(s)-sy(0)-y'(0)}+3\\underbrace{\\mathcal{L}[y'(t)]}_{sY(s)-y(0)}+2\\underbrace{\\mathcal{L}[y(t)]}_{Y(s)}=\\boxed{0}s2Y(s)−sy(0)−y′(0)L[y′′(t)]+3sY(s)−y(0)L[y′(t)]+2Y(s)L[y(t)]=0
如果是零初始条件:
y(0)=0,y′(0)=0y(0)=0,\\qquad y'(0)=0y(0)=0,y′(0)=0
那么上式变为:
s2Y(s)−s⋅0−0⏟s2Y(s)+3(sY(s)−0)⏟sY(s)+2Y(s)=0\\underbrace{s^2Y(s)-s\\cdot 0-0}_{s^2Y(s)}+3\\underbrace{(sY(s)-0)}_{sY(s)}+2Y(s)=\\boxed{0}s2Y(s)s2Y(s)−s⋅0−0+3sY(s)(sY(s)−0)+2Y(s)=0
整理得:
s2Y(s)+3sY(s)+2Y(s)=0s^2Y(s)+3sY(s)+2Y(s)=\\boxed{0}s2Y(s)+3sY(s)+2Y(s)=0
提取 Y(s)Y(s)Y(s):
(s2+3s+2)Y(s)=0(s^2+3s+2)Y(s)=\\boxed{0}(s2+3s+2)Y(s)=0
这样就把原来的微分方程:
y′′(t)+3y′(t)+2y(t)=0y''(t)+3y'(t)+2y(t)=0y′′(t)+3y′(t)+2y(t)=0
变成了 sss 域中的代数方程:
(s2+3s+2)Y(s)=0(s^2+3s+2)Y(s)=0(s2+3s+2)Y(s)=0
重点是这三个对应关系:
y′′(t)→s2Y(s)−sy(0)−y′(0)y''(t)\\rightarrow s^2Y(s)-sy(0)-y'(0)y′′(t)→s2Y(s)−sy(0)−y′(0)
y′(t)→sY(s)−y(0)y'(t)\\rightarrow sY(s)-y(0)y′(t)→sY(s)−y(0)
y(t)→Y(s)y(t)\\rightarrow Y(s)y(t)→Y(s)
在零初始条件下,它们会简化成:
y′′(t)→s2Y(s)y''(t)\\rightarrow s^2Y(s)y′′(t)→s2Y(s)
y′(t)→sY(s)y'(t)\\rightarrow sY(s)y′(t)→sY(s)
y(t)→Y(s)y(t)\\rightarrow Y(s)y(t)→Y(s)
八、总结
这一篇抓住下面几件事即可:
f(t)→F(s)f(t)\\rightarrow F(s)f(t)→F(s)
F(s)=∫0+∞f(t)e−stdtF(s)=\\int_0^{+\\infty}f(t)e^{-st}dtF(s)=∫0+∞f(t)e−stdt
把微分方程变成代数方程\\text{把微分方程变成代数方程}把微分方程变成代数方程
L[sinωt]=ωs2+ω2\\mathcal{L}[\\sin\\omega t]=\\frac{\\omega}{s^2+\\omega^2}L[sinωt]=s2+ω2ω
L[cosωt]=ss2+ω2\\mathcal{L}[\\cos\\omega t]=\\frac{s}{s^2+\\omega^2}L[cosωt]=s2+ω2s
L[f′(t)]=sF(s)−f(0)\\mathcal{L}[f'(t)]=sF(s)-f(0)L[f′(t)]=sF(s)−f(0)
L[f′′(t)]=s2F(s)−sf(0)−f′(0)\\mathcal{L}[f''(t)]=s^2F(s)-sf(0)-f'(0)L[f′′(t)]=s2F(s)−sf(0)−f′(0)
学完拉氏变换后,下一步就是学习 拉氏逆变换:
F(s)→f(t)F(s)\\rightarrow f(t)F(s)→f(t)
也就是如何从 sss 域回到时间域。



