当前位置: 首页 > news >正文

拉普拉斯变换求解线性常系数微分方程

绝大多数参考来源:通义千问

拉普拉斯定义与性质:

一、什么是拉普拉斯变换?

拉普拉斯变换(Laplace Transform)是一种积分变换,由法国数学家皮埃尔·西蒙·拉普拉斯(Pierre-Simon Laplace)在约1812年提出。它的核心思想是:

将一个关于时间 ttt 的函数 f(t)f(t)f(t) 转换为一个关于复变量 sss 的函数 F(s)F(s)F(s)

这个转换使得原本在时域中难以处理的微分、积分运算,变为在复频域中的代数运算。


二、拉普拉斯变换的定义

对于一个定义在 t≥0t \geq 0t0 上的实值函数 f(t)f(t)f(t),其单边拉氏变换定义为:

F(s)=L{f(t)}=∫0∞f(t)e−stdt F(s) = \mathcal{L}\{f(t)\} = \int_0^\infty f(t) e^{-st} dt F(s)=L{f(t)}=0f(t)estdt

其中:

  • f(t)f(t)f(t):原函数(时域函数)
  • F(s)F(s)F(s):象函数(复频域函数)
  • s=σ+jωs = \sigma + j\omegas=σ+:复变量(σ\sigmaσ 是实部,ω\omegaω 是虚部,j=−1j = \sqrt{-1}j=1
  • e−ste^{-st}est:核函数(变换的核心)

注意:我们通常使用的是单边拉氏变换(从0到∞积分),适用于因果系统(即系统在 t<0t < 0t<0 时无响应)。


三、拉氏反变换

F(s)F(s)F(s) 恢复 f(t)f(t)f(t) 的过程称为拉氏反变换,公式为:

f(t)=L−1{F(s)}=12πj∫c−j∞c+j∞F(s)estds f(t) = \mathcal{L}^{-1}\{F(s)\} = \frac{1}{2\pi j} \int_{c - j\infty}^{c + j\infty} F(s) e^{st} ds f(t)=L1{F(s)}=2πj1cjc+jF(s)estds


四、线性算子

1、什么是“线性算子”?

一个变换(或算子)T\mathcal{T}T线性的,当且仅当它满足两个性质:

  1. 可加性(Additivity)
    T{f(t)+g(t)}=T{f(t)}+T{g(t)} \mathcal{T}\{f(t) + g(t)\} = \mathcal{T}\{f(t)\} + \mathcal{T}\{g(t)\} T{f(t)+g(t)}=T{f(t)}+T{g(t)}

  2. 齐次性(Homogeneity)
    T{af(t)}=aT{f(t)} \mathcal{T}\{a f(t)\} = a \mathcal{T}\{f(t)\} T{af(t)}=aT{f(t)}

如果一个变换同时满足这两个性质,我们就说它是线性算子


2、证明:拉氏变换是线性算子

我们从定义出发,证明它满足可加性齐次性

🔷 1. 可加性(Additivity)

设两个函数 f(t)f(t)f(t)g(t)g(t)g(t),它们的拉氏变换分别为:

  • L{f(t)}=∫0∞f(t)e−stdt\mathcal{L}\{f(t)\} = \int_0^\infty f(t) e^{-st} dtL{f(t)}=0f(t)estdt
  • L{g(t)}=∫0∞g(t)e−stdt\mathcal{L}\{g(t)\} = \int_0^\infty g(t) e^{-st} dtL{g(t)}=0g(t)estdt

考虑它们的和 $ f(t) + g(t) $ 的拉氏变换:
L{f(t)+g(t)}=∫0∞[f(t)+g(t)]e−stdt \mathcal{L}\{f(t) + g(t)\} = \int_0^\infty [f(t) + g(t)] e^{-st} dt L{f(t)+g(t)}=0[f(t)+g(t)]estdt

利用积分的线性性(即 $ \int (a + b) = \int a + \int b $):
=∫0∞f(t)e−stdt+∫0∞g(t)e−stdt=L{f(t)}+L{g(t)} = \int_0^\infty f(t) e^{-st} dt + \int_0^\infty g(t) e^{-st} dt = \mathcal{L}\{f(t)\} + \mathcal{L}\{g(t)\} =0f(t)estdt+0g(t)estdt=L{f(t)}+L{g(t)}

✅ 所以:可加性成立


🔷 2. 齐次性(Homogeneity)

aaa 是任意常数(实数或复数),考虑 af(t)a f(t)af(t) 的拉氏变换:
L{af(t)}=∫0∞[af(t)]e−stdt \mathcal{L}\{a f(t)\} = \int_0^\infty [a f(t)] e^{-st} dt L{af(t)}=0[af(t)]estdt

把常数 aaa 提到积分外面(积分的齐次性):
=a∫0∞f(t)e−stdt=aL{f(t)} = a \int_0^\infty f(t) e^{-st} dt = a \mathcal{L}\{f(t)\} =a0f(t)estdt=aL{f(t)}

✅ 所以:齐次性成立


总结

拉氏变换是线性算子,因为它本质上是一个积分变换,而积分具有天然的线性性质:

∫[af(t)+bg(t)]e−stdt=a∫f(t)e−stdt+b∫g(t)e−stdt \int [a f(t) + b g(t)] e^{-st} dt = a \int f(t) e^{-st} dt + b \int g(t) e^{-st} dt [af(t)+bg(t)]estdt=af(t)estdt+bg(t)estdt

这直接保证了:
L{af(t)+bg(t)}=aF(s)+bG(s) \mathcal{L}\{a f(t) + b g(t)\} = a F(s) + b G(s) L{af(t)+bg(t)}=aF(s)+bG(s)


微分方程求解实例:

题目:

求解以下二阶常系数线性微分方程:

y′′(t)+3y′(t)+2y(t)=e−t,y(0)=1,y′(0)=0 y''(t) + 3y'(t) + 2y(t) = e^{-t}, \quad y(0) = 1, \quad y'(0) = 0 y′′(t)+3y(t)+2y(t)=et,y(0)=1,y(0)=0


第一步:对方程两边取拉氏变换

L{y′′(t)+3y′(t)+2y(t)}=L{e−t} \mathcal{L}\{y''(t) + 3y'(t) + 2y(t)\} = \mathcal{L}\{e^{-t}\} L{y′′(t)+3y(t)+2y(t)}=L{et}

利用拉氏变换的线性性质,展开得:
L{y′′(t)}+3L{y′(t)}+2L{y(t)}=L{e−t} \mathcal{L}\{y''(t) \} + 3\mathcal{L}\{y'(t) \} + 2\mathcal{L}\{y(t)\} = \mathcal{L}\{e^{-t}\} L{y′′(t)}+3L{y(t)}+2L{y(t)}=L{et}

L{y(t)}=Y(s)\mathcal{L}\{y(t)\} = Y(s)L{y(t)}=Y(s)

我们逐项变换:

1. L{y′′(t)}\mathcal{L}\{y''(t)\}L{y′′(t)}

根据微分定理(后续有推导):
L{y′′(t)}=s2Y(s)−sy(0)−y′(0)=s2Y(s)−s⋅1−0=s2Y(s)−s \mathcal{L}\{y''(t)\} = s^2 Y(s) - s y(0) - y'(0) = s^2 Y(s) - s \cdot 1 - 0 = s^2 Y(s) - s L{y′′(t)}=s2Y(s)sy(0)y(0)=s2Y(s)s10=s2Y(s)s

2.L{y′(t)}\mathcal{L}\{y'(t)\}L{y(t)}

L{y′(t)}=sY(s)−y(0)=sY(s)−1 \mathcal{L}\{y'(t)\} = s Y(s) - y(0) = s Y(s) - 1 L{y(t)}=sY(s)y(0)=sY(s)1

3. L{2y(t)}=2Y(s)\mathcal{L}\{2y(t)\} = 2Y(s)L{2y(t)}=2Y(s)
4. L{e−t}\mathcal{L}\{e^{-t}\}L{et}

查表:L{e−at}=1s+a\mathcal{L}\{e^{-at}\} = \frac{1}{s+a}L{eat}=s+a1,所以
L{e−t}=1s+1 \mathcal{L}\{e^{-t}\} = \frac{1}{s+1} L{et}=s+11


第二步:将所有项代入原方程的拉氏变换形式

原方程变为:

[s2Y(s)−s]+3[sY(s)−1]+2Y(s)=1s+1 [s^2 Y(s) - s] + 3[s Y(s) - 1] + 2Y(s) = \frac{1}{s+1} [s2Y(s)s]+3[sY(s)1]+2Y(s)=s+11

展开整理:

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

合并同类项:

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

移项:

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


第三步:解出 Y(s)Y(s)Y(s)

先化简右边:

1s+1+s+3=1+(s+3)(s+1)s+1 \frac{1}{s+1} + s + 3 = \frac{1 + (s+3)(s+1)}{s+1} s+11+s+3=s+11+(s+3)(s+1)

计算分子:
(s+3)(s+1)=s2+4s+3⇒1+s2+4s+3=s2+4s+4 (s+3)(s+1) = s^2 + 4s + 3 \Rightarrow 1 + s^2 + 4s + 3 = s^2 + 4s + 4 (s+3)(s+1)=s2+4s+31+s2+4s+3=s2+4s+4

所以右边为:
s2+4s+4s+1 \frac{s^2 + 4s + 4}{s+1} s+1s2+4s+4

左边系数:
s2+3s+2=(s+1)(s+2) s^2 + 3s + 2 = (s+1)(s+2) s2+3s+2=(s+1)(s+2)

因此:
Y(s)=s2+4s+4(s+1)(s+1)(s+2)=s2+4s+4(s+1)2(s+2) Y(s) = \frac{s^2 + 4s + 4}{(s+1)(s+1)(s+2)} = \frac{s^2 + 4s + 4}{(s+1)^2 (s+2)} Y(s)=(s+1)(s+1)(s+2)s2+4s+4=(s+1)2(s+2)s2+4s+4

注意:分子 s2+4s+4=(s+2)2s^2 + 4s + 4 = (s+2)^2s2+4s+4=(s+2)2

所以:
Y(s)=(s+2)2(s+1)2(s+2)=s+2(s+1)2(约去一个 s+2,前提是 s≠−2) Y(s) = \frac{(s+2)^2}{(s+1)^2 (s+2)} = \frac{s+2}{(s+1)^2} \quad \text{(约去一个 } s+2 \text{,前提是 } s \ne -2\text{)} Y(s)=(s+1)2(s+2)(s+2)2=(s+1)2s+2(约去一个 s+2,前提是 s=2

最终得到:
Y(s)=s+2(s+1)2 Y(s) = \frac{s+2}{(s+1)^2} Y(s)=(s+1)2s+2


第四步:对 Y(s)Y(s)Y(s) 进行部分分式分解(可选)

我们可以尝试分解以便查表反变换:

s+2(s+1)2=As+1+B(s+1)2 \frac{s+2}{(s+1)^2} = \frac{A}{s+1} + \frac{B}{(s+1)^2} (s+1)2s+2=s+1A+(s+1)2B

两边同乘 (s+1)2(s+1)^2(s+1)2
s+2=A(s+1)+B s + 2 = A(s+1) + B s+2=A(s+1)+B

s=−1s = -1s=1
−1+2=A(0)+B⇒B=1 -1 + 2 = A(0) + B \Rightarrow B = 1 1+2=A(0)+BB=1

再令 s=0s = 0s=0
0+2=A(1)+1⇒A=1 0 + 2 = A(1) + 1 \Rightarrow A = 1 0+2=A(1)+1A=1

所以:
Y(s)=1s+1+1(s+1)2 Y(s) = \frac{1}{s+1} + \frac{1}{(s+1)^2} Y(s)=s+11+(s+1)21


第五步:取拉氏反变换

查表:

  • L−1{1s+1}=e−t\mathcal{L}^{-1}\left\{\frac{1}{s+1}\right\} = e^{-t}L1{s+11}=et
  • L−1{1(s+1)2}=te−t\mathcal{L}^{-1}\left\{\frac{1}{(s+1)^2}\right\} = t e^{-t}L1{(s+1)21}=tet

所以:

y(t)=e−t+te−t=e−t(1+t),t≥0 y(t) = e^{-t} + t e^{-t} = e^{-t}(1 + t), \quad t \geq 0 y(t)=et+tet=et(1+t),t0


微分定理

拉氏变换的微分定理是:

如果 f(t)f(t)f(t) 的拉氏变换为 F(s)=L{f(t)}F(s) = \mathcal{L}\{f(t)\}F(s)=L{f(t)},且 f(t)f(t)f(t)t≥0t \geq 0t0 上连续可导,那么:

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

其中 f(0−)f(0^-)f(0)t→0−t \to 0^-t0 时的初值(通常记为 f(0)f(0)f(0),在因果系统中理解为 t=0t=0t=0 时刻的初始状态)。

这个公式是整个拉氏变换应用于微分方程的核心


1、推导过程:从定义出发

我们从拉氏变换的定义开始:

L{f′(t)}=∫0∞f′(t)e−stdt \mathcal{L}\{f'(t)\} = \int_0^\infty f'(t) e^{-st} dt L{f(t)}=0f(t)estdt

这是一个积分,被积函数是两个函数的乘积:f′(t)f'(t)f(t)e−ste^{-st}est
这提示我们可以使用分部积分法(Integration by Parts)。

分部积分公式:

∫u dv=uv−∫v du \int u \, dv = uv - \int v \, du udv=uvvdu

我们设:

  • u=e−st⇒du=−se−stdtu = e^{-st} \quad \Rightarrow \quad du = -s e^{-st} dtu=estdu=sestdt
  • dv=f′(t)dt⇒v=f(t)dv = f'(t) dt \quad \Rightarrow \quad v = f(t)dv=f(t)dtv=f(t)

✅ 注意:∫f′(t)dt=f(t)\int f'(t) dt = f(t)f(t)dt=f(t),这是基本积分性质。

代入分部积分公式:

∫0∞f′(t)e−stdt=[f(t)e−st]t=0t=∞−∫0∞f(t)(−se−st)dt \int_0^\infty f'(t) e^{-st} dt = \left[ f(t) e^{-st} \right]_{t=0}^{t=\infty} - \int_0^\infty f(t) (-s e^{-st}) dt 0f(t)estdt=[f(t)est]t=0t=0f(t)(sest)dt

化简:

=[lim⁡t→∞f(t)e−st−f(0)e0]+s∫0∞f(t)e−stdt = \left[ \lim_{t \to \infty} f(t) e^{-st} - f(0) e^{0} \right] + s \int_0^\infty f(t) e^{-st} dt =[tlimf(t)estf(0)e0]+s0f(t)estdt

即:

=[lim⁡t→∞f(t)e−st−f(0)]+sF(s) = \left[ \lim_{t \to \infty} f(t) e^{-st} - f(0) \right] + s F(s) =[tlimf(t)estf(0)]+sF(s)


2、关键一步:极限项的处理

我们来看这一项:
lim⁡t→∞f(t)e−st \lim_{t \to \infty} f(t) e^{-st} tlimf(t)est

  • 假设 f(t)f(t)f(t) 的增长速度不超过某个指数函数(这是拉氏变换存在的条件之一);
  • e−st=e−σte−jωte^{-st} = e^{-\sigma t} e^{-j\omega t}est=eσtet,其中 s=σ+jωs = \sigma + j\omegas=σ+
  • Re(s)=σ>σ0\text{Re}(s) = \sigma > \sigma_0Re(s)=σ>σ0(某个收敛横坐标)时,e−σte^{-\sigma t}eσt 衰减得比 f(t)f(t)f(t) 增长得快;
  • 所以只要 sss收敛域内,就有:
    lim⁡t→∞f(t)e−st=0 \lim_{t \to \infty} f(t) e^{-st} = 0 tlimf(t)est=0

✅ 这就是为什么拉氏变换要求“在某个右半平面收敛”。

因此,极限项为 0。


3、最终结果

代入上面的结果:

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

这就是我们熟悉的微分定理!

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


4、高阶导数的推广

二阶导数:

L{f′′(t)}=L{ddt(f′(t))} \mathcal{L}\{f''(t)\} = \mathcal{L}\left\{ \frac{d}{dt}(f'(t)) \right\} L{f′′(t)}=L{dtd(f(t))}

套用一次微分定理:

=s⋅L{f′(t)}−f′(0)=s[sF(s)−f(0)]−f′(0)=s2F(s)−sf(0)−f′(0) = s \cdot \mathcal{L}\{f'(t)\} - f'(0) = s[sF(s) - f(0)] - f'(0) = s^2 F(s) - s f(0) - f'(0) =sL{f(t)}f(0)=s[sF(s)f(0)]f(0)=s2F(s)sf(0)f(0)

三阶导数:

L{f′′′(t)}=s3F(s)−s2f(0)−sf′(0)−f′′(0) \mathcal{L}\{f'''(t)\} = s^3 F(s) - s^2 f(0) - s f'(0) - f''(0) L{f′′′(t)}=s3F(s)s2f(0)sf(0)f′′(0)

规律:每求一次导,就多一个初始条件项,且 s 的幂次降低一次


5、直观理解:为什么会有 −f(0)-f(0)f(0) 这一项?

我们可以这样理解:

  • 拉氏变换是从 t=0t=0t=0 开始积分的,它“看不到” t<0t<0t<0 的信息;
  • 但函数的导数包含了从 t=0−t=0^-t=0 到$ t=0^+$ 的跳变信息;
  • 所以在变换中,必须显式地补上初始状态的影响
  • −f(0)-f(0)f(0) 正是对“历史状态”的修正项。

🔍 类比:就像你解微分方程时需要初始条件一样,拉氏变换在转换过程中就把这些条件“打包”进去了。


6、特殊情况:零初始条件

如果系统初始静止,即 f(0)=0f(0) = 0f(0)=0,那么:

L{f′(t)}=sF(s) \mathcal{L}\{f'(t)\} = sF(s) L{f(t)}=sF(s)

这意味着:

时域中的微分操作 ⇔ 复频域中乘以 sss

这正是拉氏变换强大的地方——微分变成了乘法


补充

在我们的实例计算中,可以看到有两次查表,一次是等式右边的e−te^{-t}et做拉氏变换,不是所有的函数都可以做拉氏变换,所以等式右边是有一定的要求(后面会给出),另一次是求到Y(s)Y(s)Y(s)后做拉氏逆变换。其实查表就是提前求出来的结果,以便之后直接使用,就跟求导一样,提前记好各个函数求导的结果是什么。

一、可进行拉氏变换的函数

1、具体条件详解

条件说明
必须是因果信号t<0t < 0t<0 时为 0,通常乘以 u(t)u(t)u(t)
分段连续允许有限个跳跃间断点(如方波)
指数阶增长存在 M,σ0M, \sigma_0M,σ0,使得 r(t)≤Meσ0tr(t) \leq M e^{\sigma_0 t}r(t)Meσ0t
常见工程信号都支持阶跃、冲激、正弦、指数、多项式等

2、 常见支持的输入类型

类型示例拉氏变换 F(s)F(s)F(s)是否支持
常数5555s\frac{5}{s}s5✅ 支持
单位阶跃函数$u(t $1s\frac{1}{s}s1✅ 支持
指数函数e−ate^{-at}eat, e2te^{2t}e2t1s+a\frac{1}{s+a}s+a1, 1s−2\frac{1}{s-2}s21✅ 支持(需注意收敛域)
幂函数tnt^ntnn!sn+1\frac{n!}{s^{n+1}}sn+1n!✅ 支持
正弦/余弦函数sin⁡(ωt)\sin(\omega t)sin(ωt), cos⁡(ωt)\cos(\omega t)cos(ωt)ωs2+ω2\frac{\omega}{s^2 + \omega^2}s2+ω2ω, ss2+ω2\frac{s}{s^2 + \omega^2}s2+ω2s✅ 支持
斜坡函数t⋅u(t)t \cdot u(t)tu(t)1s2\frac{1}{s^2}s21✅ 支持
冲激函数δ(t)\delta(t)δ(t)111✅ 支持(理想激励)
延迟信号f(t−a)u(t−a)f(t - a)u(t - a)f(ta)u(ta)e−asF(s)e^{-as}F(s)easF(s)✅ 支持
分段连续函数方波、三角波(周期信号)可变换(如用 1−e−sTs(1+e−sT)\frac{1 - e^{-sT}}{s(1 + e^{-sT})}s(1+esT)1esT 等)✅ 支持
多项式 × 指数t2e−3tt^2 e^{-3t}t2e3t2(s+3)3\frac{2}{(s+3)^3}(s+3)32✅ 支持

3、不支持或需特殊处理的情况

类型示例问题处理方式
不满足指数阶增长et2e^{t^2}et2增长太快,积分发散❌ 不支持
非因果信号(t<0t<0t<0 非零)sin⁡(t)\sin(t)sin(t)(定义在 −∞<t<∞-\infty < t < \infty<t<拉氏变换默认 t≥0t \geq 0t0需截断为 sin⁡(t)u(t)\sin(t)u(t)sin(t)u(t)
无界振荡sin⁡(t2)\sin(t^2)sin(t2)非指数阶,难以变换❌ 一般不支持
随机信号(白噪声)w(t)w(t)w(t)不是确定函数需用随机过程方法(如功率谱)
狄拉克函数的导数δ′(t)\delta'(t)δ(t)数学上存在,但物理不可实现可形式变换为 sss,但需谨慎

二、可查表的拉氏变换

一、回顾:拉氏变换的定义

拉氏变换的定义是:

L{f(t)}=F(s)=∫0∞f(t)e−stdt \mathcal{L}\{f(t)\} = F(s) = \int_0^\infty f(t) e^{-st} dt L{f(t)}=F(s)=0f(t)estdt

我们要做的,就是对不同的 f(t)f(t)f(t),计算这个积分,得到对应的 F(s)F(s)F(s)。下面我们就逐个推导几个最常用的函数。


二、逐个推导常用函数的拉氏变换

1. 单位阶跃函数 f(t)=u(t)f(t) = u(t)f(t)=u(t) 或简写为 111

定义u(t)={0,t<01,t≥0u(t) = \begin{cases} 0, & t < 0 \\ 1, & t \geq 0 \end{cases}u(t)={0,1,t<0t0

由于我们只关心 t≥0t \geq 0t0,所以 f(t)=1f(t) = 1f(t)=1

L{1}=∫0∞1⋅e−stdt=∫0∞e−stdt \mathcal{L}\{1\} = \int_0^\infty 1 \cdot e^{-st} dt = \int_0^\infty e^{-st} dt L{1}=01estdt=0estdt

这是一个基本积分:

=[e−st−s]0∞=lim⁡t→∞e−st−s−e0−s = \left[ \frac{e^{-st}}{-s} \right]_0^\infty = \lim_{t \to \infty} \frac{e^{-st}}{-s} - \frac{e^{0}}{-s} =[sest]0=tlimsestse0

  • Re(s)>0\text{Re}(s) > 0Re(s)>0(复数 s 的实部) 时,e−st→0e^{-st} \to 0est0
  • 所以第一项为 0

=0−(−1s)=1s = 0 - \left( -\frac{1}{s} \right) = \frac{1}{s} =0(s1)=s1

结论
L{1}=1s,Re(s)>0 \boxed{\mathcal{L}\{1\} = \frac{1}{s}, \quad \text{Re}(s) > 0} L{1}=s1,Re(s)>0


2. 指数函数 f(t)=e−atf(t) = e^{-at}f(t)=eataaa 为常数)

L{e−at}=∫0∞e−ate−stdt=∫0∞e−(s+a)tdt \mathcal{L}\{e^{-at}\} = \int_0^\infty e^{-at} e^{-st} dt = \int_0^\infty e^{-(s+a)t} dt L{eat}=0eatestdt=0e(s+a)tdt

这和上面一样:

=[e−(s+a)t−(s+a)]0∞=0−(−1s+a)=1s+a = \left[ \frac{e^{-(s+a)t}}{-(s+a)} \right]_0^\infty = 0 - \left( -\frac{1}{s+a} \right) = \frac{1}{s+a} =[(s+a)e(s+a)t]0=0(s+a1)=s+a1

收敛条件:Re(s+a)>0⇒Re(s)>−a\text{Re}(s + a) > 0 \Rightarrow \text{Re}(s) > -aRe(s+a)>0Re(s)>a

结论
L{e−at}=1s+a,Re(s)>−a \boxed{\mathcal{L}\{e^{-at}\} = \frac{1}{s+a}, \quad \text{Re}(s) > -a} L{eat}=s+a1,Re(s)>a

📌 特例:当 $ a = 0 $,就是前面的 $ \mathcal{L}{1} = \frac{1}{s} $


3. 线性函数 $ f(t) = t $

L{t}=∫0∞te−stdt \mathcal{L}\{t\} = \int_0^\infty t e^{-st} dt L{t}=0testdt

使用分部积分法

设:

  • $ u = t \Rightarrow du = dt $
  • $ dv = e^{-st} dt \Rightarrow v = -\frac{1}{s} e^{-st} $

则:
∫te−stdt=−tse−st+∫1se−stdt=−tse−st−1s2e−st \int t e^{-st} dt = -\frac{t}{s} e^{-st} + \int \frac{1}{s} e^{-st} dt = -\frac{t}{s} e^{-st} - \frac{1}{s^2} e^{-st} testdt=stest+s1estdt=stests21est

代入上下限 $ 0 \to \infty $:

  • 当 $ t \to \infty :: t e^{-st} \to 0 (因为指数衰减快于线性增长),(因为指数衰减快于线性增长),(因为指数衰减快于线性增长), e^{-st} \to 0 $
  • 当 $ t = 0 :: -\frac{0}{s} \cdot 1 - \frac{1}{s^2} \cdot 1 = -\frac{1}{s^2} $

所以:
L{t}=[−tse−st−1s2e−st]0∞=0−(−1s2)=1s2 \mathcal{L}\{t\} = \left[ -\frac{t}{s} e^{-st} - \frac{1}{s^2} e^{-st} \right]_0^\infty = 0 - \left( -\frac{1}{s^2} \right) = \frac{1}{s^2} L{t}=[stests21est]0=0(s21)=s21

结论
L{t}=1s2,Re(s)>0 \boxed{\mathcal{L}\{t\} = \frac{1}{s^2}, \quad \text{Re}(s) > 0} L{t}=s21,Re(s)>0


4. 幂函数 f(t)=tnf(t) = t^nf(t)=tn(一般形式)

可以证明(数学归纳法或伽马函数):

L{tn}=n!sn+1,n=0,1,2,… \mathcal{L}\{t^n\} = \frac{n!}{s^{n+1}}, \quad n = 0,1,2,\dots L{tn}=sn+1n!,n=0,1,2,

  • n=0n=0n=0L{1}=1s\mathcal{L}\{1\} = \frac{1}{s}L{1}=s1
  • n=1n=1n=1L{t}=1s2\mathcal{L}\{t\} = \frac{1}{s^2}L{t}=s21
  • n=2n=2n=2L{t2}=2s3\mathcal{L}\{t^2\} = \frac{2}{s^3}L{t2}=s32
  • 依此类推

📌 这个公式在部分分式反变换中非常有用。


5. 单位脉冲函数(Dirac Delta)f(t)=δ(t)f(t) = \delta(t)f(t)=δ(t)

定义δ(t)\delta(t)δ(t) 是一个“在 t=0t=0t=0 处无限高、无限窄,但面积为 1”的理想化函数,满足:
∫0∞δ(t)ϕ(t)dt=ϕ(0) \int_0^\infty \delta(t) \phi(t) dt = \phi(0) 0δ(t)ϕ(t)dt=ϕ(0)
对任意连续函数 ϕ(t)\phi(t)ϕ(t) 成立。

所以:
L{δ(t)}=∫0∞δ(t)e−stdt=e−s⋅0=1 \mathcal{L}\{\delta(t)\} = \int_0^\infty \delta(t) e^{-st} dt = e^{-s \cdot 0} = 1 L{δ(t)}=0δ(t)estdt=es0=1

结论
L{δ(t)}=1 \boxed{\mathcal{L}\{\delta(t)\} = 1} L{δ(t)}=1


6. 正弦函数 f(t)=sin⁡(ωt)f(t) = \sin(\omega t)f(t)=sin(ωt)

我们用欧拉公式来推导:

sin⁡(ωt)=ejωt−e−jωt2j \sin(\omega t) = \frac{e^{j\omega t} - e^{-j\omega t}}{2j} sin(ωt)=2jetet

利用线性性

L{sin⁡(ωt)}=12j[L{ejωt}−L{e−jωt}] \mathcal{L}\{\sin(\omega t)\} = \frac{1}{2j} \left[ \mathcal{L}\{e^{j\omega t}\} - \mathcal{L}\{e^{-j\omega t}\} \right] L{sin(ωt)}=2j1[L{et}L{et}]

查指数函数公式:

  • L{ejωt}=1s−jω\mathcal{L}\{e^{j\omega t}\} = \frac{1}{s - j\omega}L{et}=s1
  • L{e−jωt}=1s+jω\mathcal{L}\{e^{-j\omega t}\} = \frac{1}{s + j\omega}L{et}=s+1

所以:
L{sin⁡(ωt)}=12j(1s−jω−1s+jω) \mathcal{L}\{\sin(\omega t)\} = \frac{1}{2j} \left( \frac{1}{s - j\omega} - \frac{1}{s + j\omega} \right) L{sin(ωt)}=2j1(s1s+1)

通分:
=12j⋅(s+jω)−(s−jω)(s−jω)(s+jω)=12j⋅2jωs2+ω2=ωs2+ω2 = \frac{1}{2j} \cdot \frac{(s + j\omega) - (s - j\omega)}{(s - j\omega)(s + j\omega)} = \frac{1}{2j} \cdot \frac{2j\omega}{s^2 + \omega^2} = \frac{\omega}{s^2 + \omega^2} =2j1(s)(s+)(s+)(s)=2j1s2+ω22=s2+ω2ω

结论
L{sin⁡(ωt)}=ωs2+ω2,Re(s)>0 \boxed{\mathcal{L}\{\sin(\omega t)\} = \frac{\omega}{s^2 + \omega^2}, \quad \text{Re}(s) > 0} L{sin(ωt)}=s2+ω2ω,Re(s)>0


7. 余弦函数 f(t)=cos⁡(ωt)f(t) = \cos(\omega t)f(t)=cos(ωt)

同样用欧拉公式:
cos⁡(ωt)=ejωt+e−jωt2 \cos(\omega t) = \frac{e^{j\omega t} + e^{-j\omega t}}{2} cos(ωt)=2et+et

L{cos⁡(ωt)}=12(1s−jω+1s+jω)=12⋅2ss2+ω2=ss2+ω2 \mathcal{L}\{\cos(\omega t)\} = \frac{1}{2} \left( \frac{1}{s - j\omega} + \frac{1}{s + j\omega} \right) = \frac{1}{2} \cdot \frac{2s}{s^2 + \omega^2} = \frac{s}{s^2 + \omega^2} L{cos(ωt)}=21(s1+s+1)=21s2+ω22s=s2+ω2s

结论
L{cos⁡(ωt)}=ss2+ω2,Re(s)>0 \boxed{\mathcal{L}\{\cos(\omega t)\} = \frac{s}{s^2 + \omega^2}, \quad \text{Re}(s) > 0} L{cos(ωt)}=s2+ω2s,Re(s)>0


三、总结:查表是怎么来的?

函数 f(t)f(t)f(t)推导方法
1,e−at1, e^{-at}1,eat直接积分
t,tnt, t^nt,tn分部积分 或 递推
δ(t)\delta(t)δ(t)冲激函数的筛选性质
sin⁡(ωt),cos⁡(ωt)\sin(\omega t), \cos(\omega t)sin(ωt),cos(ωt)欧拉公式 + 指数函数变换
te−att e^{-at}teat分部积分 或 利用 S域微分性质

✅ 所有这些结果都被“打包”进了一张拉氏变换表,我们不需要每次都重新积分,只需查表即可。


四、补充:S域微分性质(快速推导高阶函数)

还有一个高级技巧:S域微分性质

L{tf(t)}=−ddsF(s) \mathcal{L}\{t f(t)\} = -\frac{d}{ds} F(s) L{tf(t)}=dsdF(s)

例如:

  • L{e−at}=1s+a\mathcal{L}\{e^{-at}\} = \frac{1}{s+a}L{eat}=s+a1
  • 所以 L{te−at}=−dds(1s+a)=1(s+a)2\mathcal{L}\{t e^{-at}\} = -\frac{d}{ds} \left( \frac{1}{s+a} \right) = \frac{1}{(s+a)^2}L{teat}=dsd(s+a1)=(s+a)21

这比直接积分更快!


三、可查表的拉氏逆变换

逆变换的理论公式(反演积分)是:
f(t)=L−1{F(s)}=12πj∫c−j∞c+j∞F(s)estds f(t) = \mathcal{L}^{-1}\{F(s)\} = \frac{1}{2\pi j} \int_{c - j\infty}^{c + j\infty} F(s) e^{st} ds f(t)=L1{F(s)}=2πj1cjc+jF(s)estds
这是一个复变函数中的围线积分,在工程中直接计算较难。

但我们在实际应用中“查表”得到的那些逆变换结果,并不是每次都用这个复积分算出来的,而是通过以下几种方法“反向推导”或“构建”出来的:


一、逆变换查表公式的来源

逆变换主要是通过以下 四种方法 得到的:

✅ 方法 1:利用正变换的已知结果(对称性)

这是最直接的方法。

既然我们知道:
L{e−at}=1s+a \mathcal{L}\{e^{-at}\} = \frac{1}{s+a} L{eat}=s+a1
那么自然就有:
L−1{1s+a}=e−at \mathcal{L}^{-1}\left\{\frac{1}{s+a}\right\} = e^{-at} L1{s+a1}=eat

查表的本质是“双向记忆”:正变换和逆变换是一一对应的。

所以,所有正变换推导出的结果,都可以“倒过来”作为逆变换查表的依据。


✅ 方法 2:部分分式分解 + 查基本表

这是工程中最常用的方法。适用于有理函数(多项式之比)。

例子:求

F(s)=1(s+1)(s+2) F(s) = \frac{1}{(s+1)(s+2)} F(s)=(s+1)(s+2)1
的逆变换。

步骤 1:部分分式分解

设:
1(s+1)(s+2)=As+1+Bs+2 \frac{1}{(s+1)(s+2)} = \frac{A}{s+1} + \frac{B}{s+2} (s+1)(s+2)1=s+1A+s+2B

两边乘 (s+1)(s+2)(s+1)(s+2)(s+1)(s+2)
1=A(s+2)+B(s+1) 1 = A(s+2) + B(s+1) 1=A(s+2)+B(s+1)

s=−1s = -1s=11=A(1)⇒A=11 = A(1) \Rightarrow A = 11=A(1)A=1
s=−2s = -2s=21=B(−1)⇒B=−11 = B(-1) \Rightarrow B = -11=B(1)B=1

所以:
F(s)=1s+1−1s+2 F(s) = \frac{1}{s+1} - \frac{1}{s+2} F(s)=s+11s+21

步骤 2:查表

  • L−1{1s+1}=e−t\mathcal{L}^{-1}\left\{\frac{1}{s+1}\right\} = e^{-t}L1{s+11}=et
  • L−1{1s+2}=e−2t\mathcal{L}^{-1}\left\{\frac{1}{s+2}\right\} = e^{-2t}L1{s+21}=e2t

步骤 3:线性叠加

f(t)=e−t−e−2t,t≥0 f(t) = e^{-t} - e^{-2t}, \quad t \geq 0 f(t)=ete2t,t0

✅ 这就是我们查表能直接用的原因:复杂函数被拆成基本函数的和


✅ 方法 3:配方法 + 查三角函数表

适用于分母为二次项的情况。

例子:求

F(s)=1s2+2s+5 F(s) = \frac{1}{s^2 + 2s + 5} F(s)=s2+2s+51
的逆变换。

步骤 1:配方

s2+2s+5=(s+1)2+4=(s+1)2+22 s^2 + 2s + 5 = (s+1)^2 + 4 = (s+1)^2 + 2^2 s2+2s+5=(s+1)2+4=(s+1)2+22

所以:
F(s)=1(s+1)2+22 F(s) = \frac{1}{(s+1)^2 + 2^2} F(s)=(s+1)2+221

步骤 2:利用 S域平移性质

我们知道:
L{sin⁡(2t)}=2s2+4⇒L−1{1s2+4}=12sin⁡(2t) \mathcal{L}\{\sin(2t)\} = \frac{2}{s^2 + 4} \Rightarrow \mathcal{L}^{-1}\left\{\frac{1}{s^2 + 4}\right\} = \frac{1}{2} \sin(2t) L{sin(2t)}=s2+42L1{s2+41}=21sin(2t)

而:
L−1{F(s−a)}=eatf(t)(S域平移) \mathcal{L}^{-1}\left\{F(s - a)\right\} = e^{at} f(t) \quad \text{(S域平移)} L1{F(sa)}=eatf(t)(S域平移)

这里 F(s+1)=1(s+1)2+4F(s+1) = \frac{1}{(s+1)^2 + 4}F(s+1)=(s+1)2+41,所以:

L−1{1(s+1)2+22}=e−t⋅12sin⁡(2t) \mathcal{L}^{-1}\left\{\frac{1}{(s+1)^2 + 2^2}\right\} = e^{-t} \cdot \frac{1}{2} \sin(2t) L1{(s+1)2+221}=et21sin(2t)

即:
f(t)=12e−tsin⁡(2t) f(t) = \frac{1}{2} e^{-t} \sin(2t) f(t)=21etsin(2t)

✅ 所以查表时,我们不仅记“形式”,还记“如何变形”。


✅ 方法 4:利用性质推导新公式(如 S域微分、时域微分)

我们可以用已知变换 + 性质,推导出新的逆变换。

例子:求

L−1{1(s+a)2} \mathcal{L}^{-1}\left\{\frac{1}{(s+a)^2}\right\} L1{(s+a)21}

我们知道:
L{e−at}=1s+a \mathcal{L}\{e^{-at}\} = \frac{1}{s+a} L{eat}=s+a1

利用 S域微分性质
L{tf(t)}=−ddsF(s)⇒L−1{−ddsF(s)}=tf(t) \mathcal{L}\{t f(t)\} = -\frac{d}{ds} F(s) \Rightarrow \mathcal{L}^{-1}\left\{-\frac{d}{ds} F(s)\right\} = t f(t) L{tf(t)}=dsdF(s)L1{dsdF(s)}=tf(t)

所以:
−dds(1s+a)=1(s+a)2⇒L−1{1(s+a)2}=te−at -\frac{d}{ds} \left( \frac{1}{s+a} \right) = \frac{1}{(s+a)^2} \Rightarrow \mathcal{L}^{-1}\left\{\frac{1}{(s+a)^2}\right\} = t e^{-at} dsd(s+a1)=(s+a)21L1{(s+a)21}=teat

✅ 所以我们就能在表中写下:
1(s+a)2↔te−at \frac{1}{(s+a)^2} \leftrightarrow t e^{-at} (s+a)21teat


三、复变函数视角:反演积分(高级,了解即可)

理论上,逆变换是通过复平面上的围线积分计算的:
f(t)=12πj∫c−j∞c+j∞F(s)estds f(t) = \frac{1}{2\pi j} \int_{c - j\infty}^{c + j\infty} F(s) e^{st} ds f(t)=2πj1cjc+jF(s)estds

对于有理函数 $ F(s) $,可以用**留数定理(Residue Theorem)**计算:

f(t)=∑Res[F(s)est,极点] f(t) = \sum \text{Res}[F(s) e^{st}, \text{极点}] f(t)=Res[F(s)est,极点]

例子:F(s)=1(s+1)(s+2)F(s) = \frac{1}{(s+1)(s+2)}F(s)=(s+1)(s+2)1

极点:s=−1,−2s = -1, -2s=1,2

  • s=−1s = -1s=1 的留数:
    lim⁡s→−1(s+1)⋅est(s+1)(s+2)=e−t1=e−t \lim_{s \to -1} (s+1) \cdot \frac{e^{st}}{(s+1)(s+2)} = \frac{e^{-t}}{1} = e^{-t} s1lim(s+1)(s+1)(s+2)est=1et=et
  • s=−2s = -2s=2 的留数:
    lim⁡s→−2(s+2)⋅est(s+1)(s+2)=e−2t−1=−e−2t \lim_{s \to -2} (s+2) \cdot \frac{e^{st}}{(s+1)(s+2)} = \frac{e^{-2t}}{-1} = -e^{-2t} s2lim(s+2)(s+1)(s+2)est=1e2t=e2t

所以:
f(t)=e−t−e−2t f(t) = e^{-t} - e^{-2t} f(t)=ete2t

✅ 这和部分分式结果一致!

🔍 这说明:部分分式本质上是留数法的代数实现


http://www.dtcms.com/a/355648.html

相关文章:

  • 数字接龙(dfs)(蓝桥杯)
  • npm install 安装离线包的方法
  • 【论文阅读】健全个体无辅助运动期间可穿戴传感器双侧下肢神经机械信号的基准数据集
  • 如何打造品牌信任护城河?
  • Spark入门:从零到能跑的实战教程
  • 腾讯云重保流程详解:从预案到复盘的全周期安全防护
  • ♻️旧衣回收小程序|线上模式新升级
  • 网页爬虫的实现
  • 苹果ImageIO零日漏洞分析:攻击背景与iOS零点击漏洞历史对比
  • 2025 深度洞察!晶圆背面保护膜市场全景调研与投资机遇解析
  • 推荐一款JTools插件Crypto
  • 基于Spring Session + Redis + JWT的单点登录实现
  • Redis使用简明教程
  • SQL 查询优化全指南:从语句到架构的系统性优化策略
  • 初识分布式事务
  • week5-[一维数组]归并
  • 数据结构与算法-算法-42. 接雨水
  • AI 如何 “看见” 世界?计算机视觉(CV)的核心技术:图像识别、目标检测与语义分割
  • Scrapy 框架实战:构建高效的快看漫画分布式爬虫
  • 试验铁地板在现代工业中的应用与特性
  • AI医疗影像诊断新突破:从肺部CT结节识别到眼底病变筛查,提升疾病早诊效率
  • MTK Linux DRM分析(十四)- Mediatek KMS实现mtk_drm_drv.c(Part.2)
  • 工业机器人如何通过Modbus TCP转CanOpen网关高效通信!
  • 机器学习基本介绍
  • 【练习九】Java实现加油站支付小程序:存款与消费
  • 健永科技RFID技术在高压电厂机器狗巡检中的应用
  • Access token(访问令牌:以JWT格式无状态存储)和Refresh token(刷新令牌:有状态存储于Redis/DB)区别与联系、Redis黑名单
  • C#-mqtt通讯,服务端和客户端,以及esp32-mqtt
  • 第二十节:3D文本渲染 - 字体几何体生成与特效
  • 神经网络 | 基于matlab的LSTM详解