拉普拉斯变换
第一部分:拉普拉斯变换的基本定义与物理意义
1.1 积分变换的通用公式
积分变换(Integral Transform)的本质是通过一个选定的核函数(Kernel) K(s,t),将时域中的连续函数 f(t) 映射到另一个更容易求解的复频域空间(如 s 域或 ω 域):
F(s)=∫abf(t)K(s,t)dt
- 在傅里叶变换(Fourier Transform)中:核函数 K(s,t)=e−jωt,积分区间为实数全域 (−∞,∞)。
- 在拉普拉斯变换(Laplace Transform)中:核函数 K(s,t)=e−st。
1.2 单边拉普拉斯变换的严格定义
在实际工程分析(如电路分析、信号与系统、控制理论)中,由于我们往往关注系统自 t=0 之后加入输入信号的响应,因此主要采用单边拉普拉斯变换(One-sided Laplace Transform):
L[f(t)]=F(s)=∫0−∞f(t)e−stdt
💡 订正说明(关于 0− 下限): 你的手写笔记中写作 ∫0∞。在工程上,为了能够完整、无遗漏地覆盖到可能在 t=0 时刻发生突变或作用的单位冲激函数 δ(t),通常将单边积分的下限严格定义为 0−(即 0 的左侧极限处)。
1.3 算子 s 的物理意义
在变换中,自变量 s 是一个复变量(Complex Variable),称为复频率(Complex Frequency):
s=σ+jω
- σ(实部): 代表衰减因子(Neper frequency),单位为奈培/秒(Np/s),反映信号的指数衰减或增长速度。
- ω(虚部): 代表角频率(Angular frequency),单位为弧度/秒(rad/s),反映信号的等幅振荡速度。
复频域(s-Domain)的量纲分析
因为公式中指数项 s⋅t 必须是一个无量纲(dimensionless)的纯数:
- 时间 t 的单位为秒(s),
- 因此复频率 s 的单位必须是时间秒的倒数,即 s−1(频率单位)。
- 因此,s-domain 也被称为复频域。
将复频域信号 F(s) 还原为时域信号 f(t) 的数学公式为:
L−1[F(s)]=f(t)=2πj1∫σ1−j∞σ1+j∞F(s)estds(t≥0)
该公式称为 Bromwich 积分。由于该积分涉及复变函数中的复轮廓积分(围道积分),计算极其繁琐。在实际工程与做题中,我们几乎从不直接使用此积分,而是通过 部分分式展开法(Partial Fraction Expansion) 结合 常用变换对 进行逆变换求解。
第二部分:为什么要引入拉普拉斯变换?
在电路暂态分析和控制系统分析中,拉普拉斯变换具有压倒性的优势:
- 化微积分为代数:它能够将时域中复杂的微积分运算(微分方程),转化为复频域中极其简单的代数乘除运算(代数方程)。
- 暂态/瞬态分析的绝对利器:相比于经典时域求解微分方程法,拉普拉斯变换在求导性质中,可以天然、自动地将电路元件(如电容、电感)的初始状态(Initial Conditions)一步带入计算,无需在求出通解后再去繁琐地套用边界条件求特解。
- 对相量法(Phasor Method)的终极推广:
- 相量法(Phasor Domain):仅适用于求解正弦稳态电路,它只能处理单一频率的正弦信号,相当于将时域直接映射到了特定频率的相量域。
- 拉普拉斯变换(s-Domain):不仅可以分析稳态,更能分析各种输入源(突变步骤信号、冲激信号、指数信号等)作用下的暂态响应(Transient Response)。当复频率 s=σ+jω 中的实部 σ=0 且系统稳定时,拉普拉斯变换就会平滑地退化为稳态下的傅里叶变换/相量分析法。
下面是高频常用的拉普拉斯变换对照表,要求熟练记忆(默认时域函数均在 t≥0 下作用,即隐含乘上 u(t)):
| 信号名称 | 时域函数 f(t) (t≥0) | 复频域函数 F(s) | 备注 |
|---|
| 单位冲激函数 | δ(t) | 1 | 频谱在全频段等幅分布 |
| 单位阶跃函数 | u(t) | s1 | 最基础的控制输入信号 |
| 斜坡函数(幂函数) | t | s21 | 一般通式:L[tn]=sn+1n! |
| 二次幂函数 | t2 | s32 | 代入上述通式 n=2 的结果 |
| 指数衰减函数 | e−at | s+a1 | 极点位于 s=−a |
| 正弦信号 | sin(ωt) | s2+ω2ω | 极点位于纯虚轴 ±jω |
| 余弦信号 | cos(ωt) | s2+ω2s | 分子为 s |
第四部分:拉普拉斯变换的核心性质
性质一:时域微分性质(Time Differentiation)
在一阶微分情况下:
L[dtdf(t)]=sF(s)−f(0−)
推广到更高阶微分(如二阶微分):
L[dt2d2f(t)]=s2F(s)−sf(0−)−f′(0−)
🔑 物理诠释: 在时域中进行一次求导,相当于在复频域直接乘以 s,但需要减去与初值相关的修正项。若系统的初始状态为零(零状态,即 f(0−)=0),则时域的求导直接等效于 s-Domain 里的乘法:dtd⟺×s。
性质二:时域积分性质(Time Integration)
在时域上从 0− 到 t 进行积分:
L[∫0−tf(τ)dτ]=sF(s)
🔑 物理诠释: 时域积分为累积效应,映射到复频域表现为直接除以 s。即:∫dt⟺÷s。
性质三:复平移性质(频移性质,S-Domain Translation)
若 L[f(t)]=F(s),则时域函数乘以指数衰减因子 e−at,对应在复频域中将自变量 s 平移 a:
L[f(t)e−at]=F(s+a)
💡 解决你的笔记疑惑:“我总感觉和 Fourier transform 转向相反”
这是一个非常经典的直觉,但实际上两者的数学转向在本质上是完美一致的!
-
傅里叶变换的频移性质:
f(t)ejω0t⟺F(ω−ω0)
(时域乘以正指数 ejω0t,频域轴发生右移(减去 ω0))
-
拉普拉斯变换的频移性质:
f(t)e−at⟺F(s+a)
(时域乘以负指数 e−at,复频域发生左移(加上 a),极点向左侧半平面移动)
如果我们令复指数中的实部为正,将拉普拉斯变换平移改写为:L[f(t)es0t]=F(s−s0),你会发现它的符号和傅里叶变换是一模一样的!之所以手写笔记中是加号 s+a,是因为时域相乘的是负指数衰减项 e−at(即这里的 s0=−a)。
第五部分:典型例题(性质的应用)
【例题 6.1 (1)】
题目:已知时域信号 v(t)=t2e−at,求其拉普拉斯变换 V(s)。
解析:
-
首先找出基础时域信号的变换。我们知道幂函数:
L[t2]=s32!=s32
-
应用复平移性质:当 t2 乘上指数项 e−at 时,只需将复频域中所有的 s 直接替换为 (s+a)。
-
由此得出:
V(s)=(s+a)32
【例题 6.1 (2)】
题目:已知时域信号 v(t)=sin(ωt+θ),求其拉普拉斯变换 V(s)。
解析:
-
由于初相位 θ 导致其不是标准的正弦信号,因此我们可以使用三角恒等式将其在时域中展开:
v(t)=sin(ωt)cosθ+cos(ωt)sinθ
-
注意到 cosθ 和 sinθ 都是常系数。根据拉普拉斯变换的线性性质(Linearity),常数可以直接提到变换算子外部:
V(s)=L[sin(ωt)cosθ]+L[cos(ωt)sinθ]
V(s)=cosθ⋅L[sin(ωt)]+sinθ⋅L[cos(ωt)]
-
代入标准正弦、余弦变换对照:
V(s)=cosθ⋅s2+ω2ω+sinθ⋅s2+ω2s
-
合并分子,得到最终结果:
V(s)=s2+ω2ssinθ+ωcosθ
【例题 6.1 (3)】
题目:已知时域信号 v(t)=dtd(e−atsinωt),求其拉普拉斯变换 V(s)。
解析(提供两种完全等价的推导方法,帮助加深理解):
方法一:直接套用一阶微分性质(最简洁、推荐)
-
设内部子信号为 g(t)=e−atsinωt。
-
首先应用复平移性质求出 G(s)。因为已知 L[sinωt]=s2+ω2ω,乘以 e−at 后,自变量 s→s+a:
G(s)=L[e−atsinωt]=(s+a)2+ω2ω
-
计算子信号在 t=0− 时刻的初始值:
g(0−)=e−atsinωt∣t=0−=e0⋅sin(0)=0
-
根据一阶微分性质,微分对应的变换为:
V(s)=sG(s)−g(0−)=s⋅(s+a)2+ω2ω−0
-
得出最终结果:
V(s)=(s+a)2+ω2sω
方法二:先在时域求导,再进行变换(物理推导法)
-
在时域中对 v(t) 使用乘法求导法则(Leibniz rule):
v(t)=dtd(e−atsinωt)=−ae−atsinωt+ωe−atcosωt
-
由于线性性质,我们分别求两项的拉普拉斯变换:
V(s)=−aL[e−atsinωt]+ωL[e−atcosωt]
-
根据复平移性质代入变换对:
- L[e−atsinωt]=(s+a)2+ω2ω
- L[e−atcosωt]=(s+a)2+ω2s+a
-
代入并通分合并:
V(s)=(s+a)2+ω2−aω+ω(s+a)=(s+a)2+ω2−aω+sω+aω=(s+a)2+ω2sω
(两种求解方法的答案高度一致,证明了微积分化为代数性质的正确性)
第六部分:部分分式展开法(Partial Fraction Expansion)
对于实际电路系统函数,最常见的形式是两个多项式的商 F(s)=D(s)N(s)。若该分式为真分式(分母 D(s) 的最高次数 n 大于分子 N(s) 的最高次数 m),我们根据分母极点(分母的根)的不同,分为三种展开模式进行逆变换。
情况 ①:单极点情况(Distinct Real Roots / Simple Poles)
如果分母 D(s)=0 拥有 n 个各不相同的实根(简单极点) p1,p2,…,pn,则 F(s) 可以被拆解为:
F(s)=s−p1K1+s−p2K2+⋯+s−pnKn
其中,每一个拆解系数 Ki(称为该极点处的留数 Residue)可以使用掩盖法(Heaviside Method)计算:
Ki=(s−pi)F(s)∣s=pi
时域逆变换为指数信号的和:
f(t)=(K1ep1t+K2ep2t+⋯+Knepnt)u(t)
【例题 1】
题目:已知复频域表达式如下,求其时域逆变换 v(t):
V(s)=(s+3)(s+4)(s+5)7s2+35s+144
解析:
-
极点为三个不同的简单极点:s=−3,−4,−5。
-
设展开形式为:
V(s)=s+3K1+s+4K2+s+5K3
-
利用公式计算留数 K1,K2,K3:
-
计算 K1:
K1=(s+3)V(s)∣s=−3=(s+4)(s+5)7s2+35s+144s=−3
K1=(−3+4)(−3+5)7(−3)2+35(−3)+144=(1)(2)63−105+144=2102=51
-
计算 K2:
K2=(s+4)V(s)∣s=−4=(s+3)(s+5)7s2+35s+144s=−4
K2=(−4+3)(−4+5)7(−4)2+35(−4)+144=(−1)(1)112−140+144=−1116=−116
-
计算 K3:
K3=(s+5)V(s)∣s=−5=(s+3)(s+4)7s2+35s+144s=−5
K3=(−5+3)(−5+4)7(−5)2+35(−5)+144=(−2)(−1)175−175+144=2144=72
-
代回展开式:
V(s)=s+351−s+4116+s+572
-
进行逆变换(由于是单边变换,时域需乘上阶跃函数 u(t)):
v(t)=(51e−3t−116e−4t+72e−5t)u(t)
【例题 2】
题目:已知系统函数 F(s)=s(s+8)(s+6)96(s+5)(s+12),求其拉普拉斯逆变换 f(t)。
解析:
-
分母根为 s=0,−8,−6。写出展开形式:
F(s)=sK1+s+8K2+s+6K3
-
计算留数:
-
计算 K1:
K1=sF(s)∣s=0=(s+8)(s+6)96(s+5)(s+12)s=0=(8)(6)96(5)(12)=485760=120
-
计算 K2:
K2=(s+8)F(s)∣s=−8=s(s+6)96(s+5)(s+12)s=−8=−8(−8+6)96(−8+5)(−8+12)
K2=1696(−3)(4)=6⋅(−12)=−72
-
计算 K3:
K3=(s+6)F(s)∣s=−6=s(s+8)96(s+5)(s+12)s=−6=−6(−6+8)96(−6+5)(−6+12)
K3=−1296(−1)(6)=8⋅6=48
-
代回展开式:
F(s)=s120−s+872+s+648
-
时域逆变换:
f(t)=(120−72e−8t+48e−6t)u(t)
情况 ②:共轭复根情况(Complex Conjugate Poles)
当分母 D(s) 中存在无法在实域分解的二次项时,系统极点必为一对共轭复数:
p1,2=α±jβ
这种情况下,展开式对应的系数也是共轭复数(K2 和 K2∗):
V(s)=s−p1K1+s−(α+jβ)K2+s−(α−jβ)K2∗
其中,通过极坐标系表示复留数:K2=∣K2∣ejθ。 此时,后面的两项共轭复极点在时域合并后,会直接组合为一个衰减余弦信号:
vconjugate(t)=2∣K2∣eαtcos(βt+θ)u(t)
【例题 3】
题目:已知系统电压的 s 域函数如下,求时域响应 v(t):
V(s)=(s+6)(s2+6s+25)100(s+3)
解析:
-
极点分析:
- 实极点:s=−6。
- 共轭复极点:令 s2+6s+25=0⟹(s+3)2+42=0⟹s=−3±j4(其中 α=−3,β=4)。
-
写出部分分式展开式:
V(s)=s+6K1+s−(−3+j4)K2+s−(−3−j4)K2∗
-
计算实极点留数 K1:
K1=(s+6)V(s)∣s=−6=s2+6s+25100(s+3)s=−6
K1=(−6)2+6(−6)+25100(−6+3)=36−36+25−300=25−300=−12
-
计算复极点留数 K2:
K2=[s−(−3+j4)]V(s)∣s=−3+j4=(s+6)(s+3+j4)100(s+3)s=−3+j4
代入 s=−3+j4(注意:时域代入后分子变为 100(j4)):
K2=(−3+j4+6)(j8)100(j4)=(3+j4)(j8)j400=−32+j24j400
分子和分母同除以 8 以化简:
K2=−4+j3j50
分子分母同乘以分母共轭项 −4−j3 展开:
K2=(−4)2+32j50(−4−j3)=25−j200−j2(150)=25150−j200=6−j8
-
转化为复数的指数极坐标形式:
- 幅值 ∣K2∣=62+(−8)2=100=10
- 相角 θ=arctan(6−8)≈−53.13∘
- 所以复留数为:K2=10e−j53.13∘
-
时域组合还原:
v(t)=K1e−6t+K2e(−3+j4)t+K2∗e(−3−j4)t
v(t)=−12e−6t+10e−j53.13∘e−3tej4t+10ej53.13∘e−3te−j4t
v(t)=−12e−6t+10e−3t[ej(4t−53.13∘)+e−j(4t−53.13∘)]
利用欧拉公式 ejϕ+e−jϕ=2cosϕ 合并共轭对:
v(t)=[−12e−6t+20e−3tcos(4t−53.13∘)]u(t) V
情况 ③:重极点情况(Repeated Poles) —— 补充内容
若分母中含有重数为 r 的重复实根(重极点),即 D(s) 包含 (s−p1)r 项,那么该部分的展开式必须展开为降幂分布的多项式之和:
F(s)=(s−p1)rDr+(s−p1)r−1Dr−1+⋯+s−p1D1+其他极点项
其中,各系数的通用求导公式为:
-
最高次数系数 Dr:
Dr=(s−p1)rF(s)∣s=p1
-
第 k 级降幂系数 Dr−k:
Dr−k=k!1dskdk[(s−p1)rF(s)]s=p1
【新增例题 4(重极点补充案例)】
题目:已知复频域信号 F(s)=(s+1)2(s+2)2s+5,求其拉普拉斯逆变换 f(t)。
解析:
-
极点分析:
- 单极点:s=−2。
- 重极点:s=−1,重数 r=2。
-
设展开分式:
F(s)=(s+1)2D2+s+1D1+s+2K
-
计算简单极点留数 K:
K=(s+2)F(s)∣s=−2=(s+1)22s+5s=−2=(−2+1)22(−2)+5=11=1
-
计算重极点最高阶系数 D2:
D2=(s+1)2F(s)s=−1=s+22s+5s=−1=−1+22(−1)+5=13=3
-
计算重极点一阶导数系数 D1:
D1=dsd[(s+1)2F(s)]s=−1=dsd[s+22s+5]s=−1
利用分式求导法则 (vu)′=v2u′v−uv′ 求解导数部分:
dsd(s+22s+5)=(s+2)22(s+2)−(2s+5)(1)=(s+2)22s+4−2s−5=(s+2)2−1
代入 s=−1:
D1=(s+2)2−1s=−1=(−1+2)2−1=−1
-
代回部分分式表达式:
F(s)=(s+1)23−s+11+s+21
-
求时域逆变换: 利用变换对照表,因为已知时域斜坡乘以指数:L[te−t]=(s+1)21。 所以时域表达式为:
f(t)=(3te−t−e−t+e−2t)u(t)