拉普拉斯逆变换
拉普拉斯逆变换是从 s 域函数 F(s) 恢复时域信号 f(t) 的运算。逆变换的一般公式为:
f(t)=2πj1∫σ−j∞σ+j∞F(s)estds,t>0
其中积分路径 Re[s]=σ 位于 F(s) 的收敛域内。直接计算该复变积分通常较为复杂,实际中常用部分分式法和留数法进行求解。
部分分式法
基本思想
对于有理函数 F(s)(分子分母均为 s 的多项式之比):
F(s)=D(s)N(s)=ansn+an−1sn−1+⋯+a0bmsm+bm−1sm−1+⋯+b0
将 F(s) 分解为多个简单分式之和,每个简单分式对应一个已知的拉普拉斯变换对,通过查表即可求得各分量的时域形式,再利用线性性质求和。
处理步骤
- 化真分式:若 m≥n(分子次数不低于分母次数),需先用多项式除法将其化为多项式与真分式之和
- 分母因式分解:将 D(s) 分解为一次因式和二次因式的乘积
- 部分分式展开:根据极点类型写出展开形式
- 求待定系数
- 逐项查表求逆变换
极点分类与展开形式
(1)单实极点
若 D(s) 有 n 个互异的单实根 p1,p2,…,pn,则:
F(s)=s−p1k1+s−p2k2+⋯+s−pnkn
系数 ki 可通过下式计算:
ki=(s−pi)F(s)s=pi
逆变换结果:
f(t)=i=1∑nkiepitu(t)
示例:求 F(s)=s2+3s+23s+5 的逆变换。
分母因式分解:s2+3s+2=(s+1)(s+2)
F(s)=s+1k1+s+2k2
k1=(s+1)F(s)s=−1=s+23s+5s=−1=12=2
k2=(s+2)F(s)s=−2=s+13s+5s=−2=−1−1=1
f(t)=(2e−t+e−2t)u(t)
(2)共轭复极点
若 D(s) 有一对共轭复极点 p1,2=α±jω,通常不直接使用复数系数,而是将对应分式合并为:
F(s)=(s+α)2+ω2As+B+⋯
逆变换结果(查正弦/余弦变换对):
L−1{(s+α)2+ω2ω}=e−αtsin(ωt)u(t)
L−1{(s+α)2+ω2s+α}=e−αtcos(ωt)u(t)
示例:求 F(s)=s2+2s+5s+1 的逆变换。
配方法:s2+2s+5=(s+1)2+4
F(s)=(s+1)2+22s+1
由变换对查表得:
f(t)=e−tcos(2t)u(t)
(3)重极点
若 D(s) 在 p0 处有 r 重极点,则部分分式展开中包含:
F(s)=s−p0k11+(s−p0)2k12+⋯+(s−p0)rk1r+⋯
系数计算公式:
k1j=(r−j)!1dsr−jdr−j[(s−p0)rF(s)]s=p0,j=1,2,…,r
逆变换结果:
L−1{(s−p0)r1}=(r−1)!tr−1ep0tu(t)
示例:求 F(s)=s2(s+1)1 的逆变换。
F(s)=sA+s2B+s+1C
计算系数:
- B=s2F(s)s=0=s+11s=0=1
- C=(s+1)F(s)s=−1=s21s=−1=1
- 求 A:通分后比较分子系数,得 A=−1
F(s)=−s1+s21+s+11
f(t)=(−1+t+e−t)u(t)
留数法
基本原理
留数法直接利用逆变换的复变积分公式,将积分转化为对 F(s)est 所有极点的留数之和:
f(t)=2πj1∫σ−j∞σ+j∞F(s)estds=k=1∑nRes[F(s)est,pk],t>0
其中 Res[F(s)est,pk] 为 F(s)est 在极点 pk 处的留数。
留数计算方法
单极点
若 pk 为 F(s) 的一阶极点,则:
Res[F(s)est,pk]=s→pklim(s−pk)F(s)est
r 重极点
若 pk 为 F(s) 的 r 阶极点,则:
Res[F(s)est,pk]=(r−1)!1s→pklimdsr−1dr−1[(s−pk)rF(s)est]
部分分式法与留数法的比较
| 方法 | 优点 | 缺点 | 适用场景 |
|---|
| 部分分式法 | 直观简便,无需复变函数知识 | 重极点系数计算较繁琐 | 有理函数 F(s),极点较少 |
| 留数法 | 统一公式,适用于各种极点类型 | 需复变函数知识,计算留数 | 极点较多或有重极点 |
注意事项
- 收敛域的重要性:相同的 F(s) 表达式若收敛域不同,对应的时域信号也不同。例如 F(s)=s+11,收敛域 Re[s]>−1 对应 e−tu(t),收敛域 Re[s]<−1 对应 −e−tu(−t)
- 因果性:在单边拉普拉斯变换中,默认 t<0 时 f(t)=0,因此无需关注收敛域对逆变换的影响
- 部分分式法查表:熟练掌握常用变换对是快速求解逆变换的关键
链接到