Skip to content

Laplace 变换 ​

拉普拉斯变换(Laplace Transform)是算子振幅分析与线性系统理论中的强有力工具。它能够将时域中的微分、积分运算转化为复频域中的代数运算,从而极大地简化了解微分方程(组)的步骤。本章将详细讨论拉普拉斯变换的定义、核心性质、利用复变函数中的留数理论计算逆变换的方法,以及其在求解常微分方程初值问题中的典型应用。

—

Laplace 变换的定义与存在定理 ​

Laplace 变换的定义 ​

定义 Laplace 变换 设函数 f(t)f(t) 在 t≥0t \ge 0 时有定义, s=σ+iωs = \sigma + i\omega 是一个复变量。若广义积分

F(s)=∫0+∞f(t)e−stdtF(s) = \int_{0}^{+\infty} f(t) e^{-st} dt

在 ss 的某一区域内收敛,则称此积分确定了一个复变量 ss 的函数 F(s)F(s),并称 F(s)F(s) 为函数 f(t)f(t) 的拉普拉斯变换(或像函数),记为:

F(s)=L[f(t)]F(s) = \mathcal{L}[f(t)]

相应地,称 f(t)f(t) 为 F(s)F(s) 的拉普拉斯逆变换(或原函数),记为:

f(t)=L−1[F(s)]f(t) = \mathcal{L}^{-1}[F(s)]

Laplace 变换的存在定理 ​

并非所有的函数都存在拉普拉斯变换。以下定理给出了积分收敛的充分条件。

定理 存在定理 若函数 f(t)f(t) 满足下列条件:

  1. 在 t≥0t \ge 0 的任何有限区间上分段连续;
  1. 当 t→+∞t \to +\infty 时,其增长速度是有限的,即存在常数 M>0M > 0 和 c≥0c \ge 0,使得:

∣f(t)∣≤Mect(称 f(t) 为指数级增长函数) |f(t)| \le M e^{ct} \quad (\text{称 } f(t) \text{ 为指数级增长函数})

则当 Re(s)>c\text{Re}(s) > c 时,拉普拉斯变换积分必绝对收敛,且像函数 F(s)F(s) 在右半平面 Re(s)>c\text{Re}(s) > c 内是解析的。

—

Laplace 变换的基本性质 ​

线性性质 ​

定理 线性性质 设 L[f1(t)]=F1(s)\mathcal{L}[f_1(t)] = F_1(s),L[f2(t)]=F2(s)\mathcal{L}[f_2(t)] = F_2(s),α,β\alpha, \beta 为任意常数,则:

L[αf1(t)+βf2(t)]=αF1(s)+βF2(s)\mathcal{L}[\alpha f_1(t) + \beta f_2(t)] = \alpha F_1(s) + \beta F_2(s)

微分性质(时域微分) ​

定理 微分定理 设 f(t)f(t) 连续,且其导数 f′(t)f'(t) 满足存在定理的条件,若 L[f(t)]=F(s)\mathcal{L}[f(t)] = F(s),则:

L[f′(t)]=sF(s)−f(0)\mathcal{L}[f'(t)] = sF(s) - f(0)

一般地,若 f(t)f(t) 的各阶导数均满足条件,则有 nn 阶微分公式:

L[f(n)(t)]=snF(s)−sn−1f(0)−sn−2f′(0)−⋯−f(n−1)(0)\mathcal{L}[f^{(n)}(t)] = s^n F(s) - s^{n-1}f(0) - s^{n-2}f'(0) - \cdots - f^{(n-1)}(0)

积分性质(时域积分) ​

定理 积分定理 若 L[f(t)]=F(s)\mathcal{L}[f(t)] = F(s),则其变上限积分的拉氏变换为:

L[∫0tf(τ)dτ]=1sF(s)\mathcal{L}\left[ \int_{0}^{t} f(\tau) d\tau \right] = \frac{1}{s} F(s)

相似性质 ​

定理 相似性质 若 aa 为正实数,且 L[f(t)]=F(s)\mathcal{L}[f(t)] = F(s),则:

L[f(at)]=1aF(sa)(Re(sa)>c)\mathcal{L}[f(at)] = \frac{1}{a} F\left(\frac{s}{a}\right) \quad \left(\text{Re}\left(\frac{s}{a}\right) > c\right)

证明 根据定义,有:

L[f(at)]=∫0+∞f(at)e−stdt\mathcal{L}[f(at)] = \int_{0}^{+\infty} f(at)e^{-st} dt

令 u=atu = at,则 t=uat = \frac{u}{a},dt=1adudt = \frac{1}{a}du。代入积分式中得:

L[f(at)]=∫0+∞f(u)e−suad(ua)=1a∫0+∞f(u)e−(sa)udu=1aF(sa)\mathcal{L}[f(at)] = \int_{0}^{+\infty} f(u)e^{-s\frac{u}{a}} d\left(\frac{u}{a}\right) = \frac{1}{a} \int_{0}^{+\infty} f(u)e^{-\left(\frac{s}{a}\right)u} du = \frac{1}{a} F\left(\frac{s}{a}\right)

延伸推论:

  1. 若已知 L[f(t)]=F(s)\mathcal{L}[f(t)] = F(s),则 L[f(ta)]=aF(as)\mathcal{L}\left[f\left(\frac{t}{a}\right)\right] = aF(as)。

  2. 结合衰减性质(位移性质),可得:L[eatf(ta)]=aF(a(s−a))\mathcal{L}\left[e^{at}f\left(\frac{t}{a}\right)\right] = aF(a(s-a))。

卷积性质 ​

定义 卷积 设函数 f(t)f(t) 与 g(t)g(t) 在 [0,+∞)[0, +\infty) 上有定义,称积分

∫0tf(τ)g(t−τ)dτ(t≥0)\int_{0}^{t} f(\tau)g(t-\tau)d\tau \quad (t \ge 0)

为函数 f(t)f(t) 与 g(t)g(t) 的卷积,记为 f(t)∗g(t)f(t) * g(t)。

卷积运算满足交换律,即:f(t)∗g(t)=g(t)∗f(t)f(t) * g(t) = g(t) * f(t)。

定理 卷积定理 若 L[f(t)]=F(s)\mathcal{L}[f(t)] = F(s),L[g(t)]=G(s)\mathcal{L}[g(t)] = G(s),则两个函数时域卷积的拉普拉斯变换等于它们各自像函数的乘积:

L[f(t)∗g(t)]=F(s)⋅G(s)\mathcal{L}[f(t) * g(t)] = F(s) \cdot G(s)

证明 由定义开展二重积分:

L[f(t)∗g(t)]=∫0+∞[∫0tf(τ)g(t−τ)dτ]e−stdt=∫0+∞f(τ)dτ∫τ+∞g(t−τ)e−stdt\mathcal{L}[f(t) * g(t)] = \int_{0}^{+\infty} \left[ \int_{0}^{t} f(\tau)g(t-\tau)d\tau \right] e^{-st} dt = \int_{0}^{+\infty} f(\tau) d\tau \int_{\tau}^{+\infty} g(t-\tau)e^{-st} dt

在内层积分中令 u=t−τu = t - \tau,则 t=u+τt = u + \tau,dt=dudt = du。积分限从 [τ,+∞)[\tau, +\infty) 变为 [0,+∞)[0, +\infty):

L[f(t)∗g(t)]=∫0+∞f(τ)dτ∫0+∞g(u)e−s(u+τ)du=(∫0+∞f(τ)e−sτdτ)⋅(∫0+∞g(u)e−sudu)=F(s)G(s)\mathcal{L}[f(t) * g(t)] = \int_{0}^{+\infty} f(\tau) d\tau \int_{0}^{+\infty} g(u)e^{-s(u+\tau)} du = \left( \int_{0}^{+\infty} f(\tau)e^{-s\tau} d\tau \right) \cdot \left( \int_{0}^{+\infty} g(u)e^{-su} du \right) = F(s)G(s)

逆变换应用:

L−1[F(s)G(s)]=f(t)∗g(t)=∫0tf(τ)g(t−τ)dτ\mathcal{L}^{-1}[F(s)G(s)] = f(t) * g(t) = \int_{0}^{t} f(\tau)g(t-\tau)d\tau

初值定理与终值定理 ​

定理 初值定理 设 L[f(t)]=F(s)\mathcal{L}[f(t)] = F(s),且 lim⁡t→0+f(t)\lim_{t \to 0^+} f(t) 存在,则:

lim⁡t→0+f(t)=lim⁡s→∞sF(s)\lim_{t \to 0^+} f(t) = \lim_{s \to \infty} sF(s)

定理 终值定理 设 L[f(t)]=F(s)\mathcal{L}[f(t)] = F(s),若 lim⁡t→+∞f(t)\lim_{t \to +\infty} f(t) 存在,且像函数 F(s)F(s) 在右半复平面及虚轴上除原点外无奇点(即 F(s)F(s) 的所有奇点均位于左半平面 Re(s)<0\text{Re}(s) < 0 或原点),则:

lim⁡t→+∞f(t)=lim⁡s→0sF(s)\lim_{t \to +\infty} f(t) = \lim_{s \to 0} sF(s)

—

用留数计算 Laplace 逆变换 ​

定理 留数反演定理 设像函数 F(s)F(s) 在整个扩充复平面上除去有限个孤立奇点 s1,s2,⋯ ,sns_1, s_2, \cdots, s_n 外处处解析,且满足条件 lim⁡s→∞F(s)=0\lim_{s \to \infty} F(s) = 0。如果选取实数 β\beta 使得 F(s)F(s) 的所有奇点都位于半平面 Re(s)<β\text{Re}(s) < \beta 内,则当 t>0t > 0 时,F(s)F(s) 的拉普拉斯逆变换 f(t)f(t) 可以表示为 F(s)estF(s)e^{st} 在其所有奇点处的留数之和:

f(t)=L−1[F(s)]=∑k=1nRes[F(s)est,sk](t>0)f(t) = \mathcal{L}^{-1}[F(s)] = \sum_{k=1}^{n} \text{Res}[F(s)e^{st}, s_k] \quad (t > 0)

典型极点处的留数计算法 ​

1. 单极点(一阶极点):

Res[F(s)est,sk]=lim⁡s→sk(s−sk)F(s)est\text{Res}[F(s)e^{st}, s_k] = \lim_{s \to s_k} (s - s_k)F(s)e^{st}

  1. 高阶极点(以 mm 阶极点为例):

Res[F(s)est,sk]=1(m−1)!lim⁡s→skdm−1biographicaldsm−1[(s−sk)mF(s)est]\text{Res}[F(s)e^{st}, s_k] = \frac{1}{(m-1)!} \lim_{s \to s_k} \frac{d^{m-1 biographical}}{ds^{m-1}} \left[ (s - s_k)^m F(s)e^{st} \right]

—

Laplace 变换的应用 ​

拉普拉斯变换的基本应用逻辑是将时域的微分方程,通过导数的拉氏变换公式转化为关于 Y(s)Y(s) 的复频域代数方程,解出 Y(s)Y(s) 后再求其拉氏逆变换得到原微分方程的解。

典型例题 1:求解高阶常微分方程初值问题 ​

例题 求微分方程 y′′−2y′−3y=te−ty'' - 2y' - 3y = te^{-t} 满足初始条件 y∣t=0=0,y′∣t=0=1y|_{t=0} = 0, y'|_{t=0} = 1 的解。

证明 设 L[y(t)]=Y(s)\mathcal{L}[y(t)] = Y(s)。已知 L[te−t]=1(s+1)2\mathcal{L}[te^{-t}] = \frac{1}{(s+1)^2}。 对方程两边同时取拉普拉斯变换:

[s2Y(s)−sy(0)−y′(0)]−2[sY(s)−y(0)]−3Y(s)=1(s+1)2\left[ s^2Y(s) - s y(0) - y'(0) \right] - 2\left[ sY(s) - y(0) \right] - 3Y(s) = \frac{1}{(s+1)^2}

代入初始条件 y(0)=0,y′(0)=1y(0) = 0, y'(0) = 1 得:

(s2−2s−3)Y(s)−1=1(s+1)2(s^2 - 2s - 3)Y(s) - 1 = \frac{1}{(s+1)^2}

由于 s2−2s−3=(s−3)(s+1)s^2 - 2s - 3 = (s-3)(s+1),故:

(s−3)(s+1)Y(s)=1+1(s+1)2=(s+1)2+1(s+1)2(s-3)(s+1)Y(s) = 1 + \frac{1}{(s+1)^2} = \frac{(s+1)^2 + 1}{(s+1)^2}

Y(s)=1(s−3)(s+1)3+1(s−3)(s+1)Y(s) = \frac{1}{(s-3)(s+1)^3} + \frac{1}{(s-3)(s+1)}

为了求 y(t)=L−1[Y(s)]y(t) = \mathcal{L}^{-1}[Y(s)],我们可以分别对两部分用留数反演定理或分式展开进行计算: 对于 Y(s)Y(s),其孤立奇点为 s1=3s_1 = 3(单极点)和 s2=−1s_2 = -1(三阶极点)。 根据留数反演定理:

y(t)=Res[Y(s)est,3]+Res[Y(s)est,−1]y(t) = \text{Res}[Y(s)e^{st}, 3] + \text{Res}[Y(s)e^{st}, -1]

经计算(或结合部分分式展开与时域卷积)可得原方程的唯一时域特解。

典型例题 2:求解常微分方程组初值问题 ​

例题 求解常微分方程组

{x′′−y′′=etx′′−2y′′+y=0\begin{cases} x'' - y'' = e^t \\ x'' - 2y'' + y = 0 \end{cases}

满足初始条件:x(0)=1,x′(0)=0,y(0)=0,y′(0)=0x(0) = 1, x'(0) = 0, y(0) = 0, y'(0) = 0。

证明 设 L[x(t)]=X(s)\mathcal{L}[x(t)] = X(s),L[y(t)]=Y(s)\mathcal{L}[y(t)] = Y(s)。对方程组两边同时进行拉普拉斯变换: 因为 x(0)=1x(0)=1 且其他初值均为 0,所以依据微分性质:

L[x′′]=s2X(s)−sx(0)−x′(0)=s2X(s)−s\mathcal{L}[x''] = s^2X(s) - sx(0) - x'(0) = s^2X(s) - s

L[y′′]=s2Y(s)\mathcal{L}[y''] = s^2Y(s)

代入原方程组得到其复频域代数方程组:

{s2X(s)−s−s2Y(s)=1s−1s2X(s)−s−2s2Y(s)+Y(s)=0\begin{cases} s^2X(s) - s - s^2Y(s) = \frac{1}{s-1} \\ s^2X(s) - s - 2s^2Y(s) + Y(s) = 0 \end{cases}

两式相减(第一式减第二式)可以迅速消去变量 X(s)X(s):

s2Y(s)−Y(s)=1s−1  ⟹  (s2−1)Y(s)=1s−1s^2Y(s) - Y(s) = \frac{1}{s-1} \implies (s^2 - 1)Y(s) = \frac{1}{s-1}

解得像函数 Y(s)Y(s) 的显式表达式:

Y(s)=1(s+1)(s−1)2Y(s) = \frac{1}{(s+1)(s-1)^2}

将其代回第一式,即可解出 X(s)X(s)。最终两变量的像函数形式为:

{X(s)=1s2+1s2(s−1)+1s(s−1)2(s+1)Y(s)=1(s+1)(s−1)2\begin{cases} X(s) = \frac{1}{s^2} + \frac{1}{s^2(s-1)} + \frac{1}{s(s-1)^2(s+1)} \\ Y(s) = \frac{1}{(s+1)(s-1)^2} \end{cases}

最后,利用拉普拉斯逆变换(部分分式展开法或高阶极点留数法)可得时域特解: 对于 Y(s)Y(s):s=1s=1 是二阶极点,s=−1s=-1 是单极点。通过计算留数和反演可得:

y(t)=L−1[1(s+1)(s−1)2]=12tet−34et−14e−t+1y(t) = \mathcal{L}^{-1}\left[\frac{1}{(s+1)(s-1)^2}\right] = \frac{1}{2}te^t - \frac{3}{4}e^t - \frac{1}{4}e^{-t} + 1

同理对 X(s)X(s) 每一项求逆变换即可解出 x(t)x(t)。

用 VitePress 搭建 · 托管在 GitHub Pages