njtynhj

拉普拉斯变换（英语：Laplace transform）是应用数学中常用的一种积分变换，又名拉氏轉換，其符號為 $displaystyle mathcal Lleftf(t)right$ 。拉氏變換是一個線性變換，可將一個有實數变量 $displaystyle t(tgeq 0)$ 的函數轉換為一個变量為複數 $s$ 的函數：

displaystyle F(s)=int _0^infty f(t)e^-st,mathrm d t.

拉氏變換在大部份的應用中都是對射的，最常見的 $displaystyle f(t)$ 和 $displaystyle F(s)$ 組合常印製成表，方便查閱。拉普拉斯变换得名自法國天文學家暨數學家皮埃尔-西蒙·拉普拉斯（Pierre-Simon marquis de Laplace），他在機率論的研究中首先引入了拉氏變換。

拉氏變換和傅里叶变换有關，不過傅里叶变换將一個函數或是信號表示為許多弦波的疊加，而拉氏變換則是將一個函數表示為許多矩的疊加。拉氏變換常用來求解微分方程及積分方程。在物理及工程上常用來分析線性非時變系統，可用來分析電子電路、諧振子、光学仪器及機械設備。在這些分析中，拉氏變換可以作時域和頻域之間的轉換，在時域中輸入和輸出都是時間的函數，在頻域中輸入和輸出則是複變角頻率的函數，單位是弧度每秒。

對於一個簡單的系統，拉氏變換提供另一種系統的描述方程，可以簡化分析系統行為的時間^[1]。像時域下的線性非時變系統，在頻域下會轉換為代數方程，在時域下的捲積會變成頻域下的乘法。

形式定义

对于所有实数 $displaystyle tgeq 0$ ，函数 $displaystyle f(t)$ 的拉普拉斯变换是函数 $displaystyle F(s)$ ，定义为：

displaystyle F(s)=int _0^infty e^-stf(t),mathrm d t

参数 $s$ 是一个复数：

s=sigma +iomega ,,

sigma

和

omega

为实数。

拉普拉斯变换的其他表示法中使用 $displaystyle mathcal Lf$ 或 $displaystyle mathcal L_tleftf(t)right$ 而非 $F$ 。 $mathcal L$ 是一个运算符号。

积分的含义取决于函数的类型。该积分存在的一个必要条件是在 $f$ 必须在 $displaystyle [0,infty )$ 上局部可积。对在无穷大处衰减的局部可积函数或指数式（英语：exponential type），该积分可以理解为（恰当）勒贝格积分。然而，在很多应用中有必要将其视作在 $displaystyle infty$ 处条件收敛的反常积分。更一般的，积分可以在较弱的意义上理解，在下面会去处理。

可以用勒贝格积分定义拉普拉斯变换为有限博雷尔测度 $mu$ ^[2]

displaystyle mathcal Lmu (s)=int _[0,infty )e^-st,mathrm d mu (t).

一种特殊情况是当 $mu$ 为概率测度，或者更具体地说，是[[狄拉克 $delta$ 函数]]时。在运算微积（英语：operational calculus）中，拉普拉斯变换的测度常常被视作由分布函数 $f$ 带来的测度。在这种情况下，为了避免混淆，一般写作

displaystyle mathcal Lf(s)=int _0^-^infty e^-stf(t),mathrm d t,

其中是 0⁻ 的下限的简化符号

lim _varepsilon downarrow 0int _-varepsilon ^infty .

这个极限强调任何位于 0 的质点都被拉普拉斯变换完全捕获。虽然使用勒贝格积分，没有必要取这个极限，但它可以更自然地与拉普拉斯–斯蒂尔吉斯变换（英语：Laplace–Stieltjes transform）建立联系。

拉普拉斯逆变换

兩個相異的可積函數，只有在其差的勒貝格測度為零時，才會有相同的拉普拉斯变换。因此以轉換的角度而言，存在其反轉換。包括可積分函數在內，拉普拉斯变换是单射映射，將一個函數空間映射到其他的函數空間。典型的函數空間包括有界連續函數、函數空間 L^∞(0, ∞)、或是更廣義，在 (0, ∞) 區間內的缓增广义函数（函數的最壞情形是多項式成長）。

拉普拉斯逆变换（英语：Inverse Laplace transform）有许多不同的名称，如维奇积分、傅立叶-梅林积分、梅林逆公式，是一个复积分：

displaystyle f(t)=mathcal L^-1F(t)=frac 12pi ilim _Tto infty int _gamma -iT^gamma +iTe^stF(s),mathrm d s

其中 $gamma$ 是一个使 $displaystyle F(s)$ 的积分路径在收敛域内的实数。另一個拉普拉斯逆变换的公式是由Post反演公式（英语：Post's inversion formula）而來。

在實務上一般會配合查表，將函數的拉普拉斯变换分換為許多已知函數的拉普拉斯变换，再利用觀察的方式產生其拉普拉斯逆变换。在微分方程中會用到拉普拉斯逆变换，會比用傅利葉轉換的處理方式要簡單。

性质和定理

函数 $displaystyle f(t)$ 和 $displaystyle g(t)$ 的拉普拉斯变换分别为 $displaystyle F(s)$ 和 $displaystyle G(s)$ ：

beginalignedf(t)&=mathcal L^-1F(s)\g(t)&=mathcal L^-1G(s)endaligned

下面的表格是一系列单边拉普拉斯变换的性质：^[3]

单边拉普拉斯变换的性质
	时域	s域	注释
线性叠加	$af(t)+bg(t)$	$aF(s)+bG(s)$	可以用积分的基本规则证明。
s域一阶微分	$tf(t)$	$-F'(s)$	$displaystyle F'$ 是 $F$ 的一阶导数。
s域一般微分	$t^nf(t)$	$(-1)^nF^(n)(s)$	更一般的形式是 $displaystyle F(s)$ 的 $n$ 阶导数。
时域一阶微分	$f'(t)$	$sF(s)-f(0)$	$f$ 是一个可微函数，并且其导数为指数类型。这条性质可以通过分部积分得到。
时域二阶微分	$f''(t)$	$s^2F(s)-sf(0)-f'(0)$	$f$ 为二阶可微且二阶导数是指数型的。通过对 $displaystyle f(t)$ 应用微分性质可得。
时域一般微分	$f^(n)(t)$	$s^nF(s)-sum _k=1^ns^k-1f^(n-k)(0)$	$f$ 为 $n$ 阶可微，其 $n$ 阶导数是指数型的。通过数学归纳法证明。
s域积分	$frac 1tf(t)$	$displaystyle int _s^infty F(sigma ),mathrm d sigma$	这是由s域微分和条件收敛推导出来的。
时域积分	$displaystyle int _0^tf(tau ),mathrm d tau =(u*f)(t)$	$1 over sF(s)$	$displaystyle u(t)$ 是阶跃函数，注意到 $displaystyle (u*f)(t)$ 是 $displaystyle u(t)$ 和 $displaystyle f(t)$ 的卷积。
时间标度	$f(at)$	$frac 1aFleft(s over aright)$	$a>0$
s域平移	$e^atf(t)$	$F(s-a)$
时域平移	$f(t-a)u(t-a)$	$e^-asF(s)$	$displaystyle u(t)$ 表示阶跃函数
乘法	$f(t)g(t)$	$displaystyle frac 12pi ilim _Tto infty int _c-iT^c+iTF(sigma )G(s-sigma ),mathrm d sigma$	积分沿完全处在 $F$ 收敛域内的竖直线 $displaystyle Re(sigma )=c$ 。^[4]
卷积	$displaystyle (f*g)(t)=int _0^tf(tau )g(t-tau ),mathrm d tau$	$F(s)cdot G(s)$
复共轭	$f^*(t)$	$F^(s^)$
互相关	$f(t)star g(t)$	$F^(-s^)cdot G(s)$
周期函数	$f(t)$	$displaystyle 1 over 1-e^-Tsint _0^Te^-stf(t),mathrm d t$	$displaystyle f(t)$ 是一个周期为 $T$ 的周期函数，于是对所有 $displaystyle tgeq 0$ ，有 $displaystyle f(t)=f(t+T)$ 。这条性质是时域平移和几何级数的结果。

初值定理：

f(0^+)=lim _sto infty sF(s).

，要求

F(s)

为真分式，即分子的最高次小于分母的最高次，否则使用多项式除法将

F(s)

分解

终值定理：

f(infty )=lim _sto 0sF(s)

，要求sF(s)的所有极点都在左半复平面或原点为单极点。

由于终值定理无需经过部分分式分解或其他困难的代数就能给出长期的行为，它就很有用。如果

F(s)

在右侧面或虚轴上有极点，（例如

f(t)=e^t

或

displaystyle f(t)=sin(t)

）这个公式的行为就是未定义的。

与幂级数的关系

拉普拉斯变换可以看成是幂级数的一个连续模拟。如果 a(n) 是正整数 n 的一个离散函数，那么与 a(n) 相关的幂级数为

sum _n=0^infty a(n)x^n

其中 x 是实变量（参见Z变换）。将对 n 的加和替换成对 t 的积分，则此幂级数的连续形式为

displaystyle int _0^infty f(t)x^t,mathrm d t

其中离散型函数 a(n) 被替换成连续型的 f(t)。（参见下文梅林变换。）改变幂的基底 x 为 e 得

displaystyle int _0^infty f(t)left(e^log xright)^t,mathrm d t

要使这个积分对任何有界函数 f 都收敛，就需要满足 $log x<0$ 。使用 $- s = log x$ 代换就能得到拉普拉斯变换：

displaystyle int _0^infty f(t)e^-st,mathrm d t

换句话说，拉普拉斯变换是幂级数的一个连续模拟，只是把离散参数 n 换成了连续变量 t ， x 换成了 $e - s$ 。

与矩的关系

函数 f 的矩为

displaystyle mu _n=int _0^infty t^nf(t),mathrm d t

如果 f 的前 n 阶矩绝对收敛，则通过反复在积分符号内取微分，就得到 $(-1)^n(mathcal Lf)^(n)(0)=mu _n$ 。这在概率论里是有特别重要的意义的，其中随机变量 X 的矩是 $mu _n=E[X^n]$ 。下面的关系成立：

displaystyle mu _n=(-1)^nfrac mathrm d ^nmathrm d s^nEleft[e^-sXright].

证明函数导数的拉普拉斯变换

很方便用拉普拉斯变换的微分性质来求函数导数的变换。从拉普拉斯变换的基本表达式就可以推导如下：

displaystyle beginalignedmathcal Lleftf(t)right&=int _0^-^infty e^-stf(t),mathrm d t\&=left[frac f(t)e^-st-sright]_0^-^infty -int _0^-^infty frac e^-st-sf'(t),mathrm d tquad text(by parts)\&=left[-frac f(0^-)-sright]+frac 1smathcal Lleftf'(t)right,endaligned

导出

mathcal Lleftf'(t)right=scdot mathcal Lleftf(t)right-f(0^-),

而在双边的情形下，

displaystyle mathcal Lleftf'(t)right=sint _-infty ^infty e^-stf(t),mathrm d t=scdot mathcal Lf(t).

一般化的结果是

mathcal Lleftf^(n)(t)right=s^ncdot mathcal Lleftf(t)right-s^n-1f(0^-)-cdots -f^(n-1)(0^-),

其中 f⁽ⁿ⁾ 表示 f 的 n阶导数，可以由归纳假设得出。

计算广义积分

令 $mathcal Lleftf(t)right=F(s)$ ，则（参见上面的表格）

displaystyle mathcal Lleftfrac f(t)tright=int _s^infty F(p),mathrm d p,

或

displaystyle int _0^infty frac f(t)te^-st,mathrm d t=int _s^infty F(p),mathrm d p.

令 s → 0，假定可以改变取极限顺序，就得到性质

displaystyle int _0^infty frac f(t)t,mathrm d t=int _0^infty F(p),mathrm d p.

即便在不可以交换，此计算依然有暗示性。例如，形式上按此计算得到

displaystyle int _0^infty frac cos at-cos btt,mathrm d t=int _0^infty left(frac pp^2+a^2-frac pp^2+b^2right),mathrm d p=frac 12left.ln frac p^2+a^2p^2+b^2right

这个性质的正确性可以用其他方法证明。它是傅汝兰尼积分（Frullani integral）的一个例子。

例子还有狄利克雷积分。

与其他变换的联系

与傅里叶变换关系

连续傅里叶变换相当于计算令 s = iω 或 s = 2πfi 的双边拉普拉斯变换：

beginalignedhat f(omega )&=mathcal Fleftf(t)right\[1em]&=mathcal Lleftf(t)right

与z变换的联系

z 变换表达式为：

X(z)=sum _n=0^infty x[n]z^-n

其中 $zleftarrow e^sT$ .
比较两者表达式有：

X_q(s)=X(z)Big _z=e^sT.

拉普拉斯变换简表

下表提供了许多常用单变量函数的拉普拉斯变换。 ^[5]^[6] 对于定义和解释，请参见表末的注释。

由于拉普拉斯变换是一个线性算子：

和的拉普拉斯变换等于各项的拉普拉斯变换的总和。

mathcal Lleftf(t)+g(t)right=mathcal Lleftf(t)right+mathcal Lleftg(t)right

一个函数的倍数的拉普拉斯变换等于该函数的拉普拉斯变换的倍数。

mathcal Lleftaf(t)right=amathcal Lleftf(t)right

使用这个线性性质，以及各种三角、双曲、和复数（等）的性质，可以从其他拉普拉斯变换得到一些拉普拉斯变换，这会比直接通过使用定义更快。

单边拉普拉斯变换取时域为非负实数的函数作为输入，这就是下表中所有时域函数都乘以单位阶跃函数 u(t) 的原因。表中涉及时间延迟 τ 的条目必须是因果的（即 τ > 0）。因果系统是 t = 0 之前的冲激响应 h(t) 都为零的一个系统。在一般情况下，因果系统的收敛区域和反因果系统是不相同的。

函数

时域

f(t)=mathcal L^-1leftF(s)right

拉普拉斯s域

F(s)=mathcal Lleftf(t)right

收敛区域

参考

单位脉冲函数

delta (t)

1

所有 s

检验

延迟脉冲函数

delta (t-tau )

e^-tau s

单位脉冲函数
的时移

单位阶跃函数

u(t)

1 over s

Re(s) > 0

对单位冲激函数积分

延迟单位阶跃函数

u(t-tau )

frac 1se^-tau s

Re(s) > 0

单位阶跃函数
的时移

斜坡函数

tcdot u(t)

frac 1s^2

Re(s) > 0

两次积分
单位脉冲函数

n次幂
（n 为整数）

t^ncdot u(t)

n! over s^n+1

Re(s) > 0
(n > −1)

n次积分
单位阶跃函数

q次幂
（q 为复数）

t^qcdot u(t)

Gamma (q+1) over s^q+1

Re(s) > 0
Re(q) > −1

^[7]^[8]

n次方根

sqrt[n]tcdot u(t)

1 over s^frac 1n+1Gamma left(frac 1n+1right)

Re(s) > 0

由前一条性质中令 q = 1/n 得到。

频移的 n 次方

t^ne^-alpha tcdot u(t)

frac n!(s+alpha )^n+1

Re(s) > −α

对单位阶跃函数积分，
应用频移

延迟的频移的 n 次方

(t-tau )^ne^-alpha (t-tau )cdot u(t-tau )

frac n!cdot e^-tau s(s+alpha )^n+1

Re(s) > −α

对单位阶跃函数积分，
应用频移，
应用时移

指数衰减

e^-alpha tcdot u(t)

1 over s+alpha

Re(s) > −α

单位阶跃函数
的频移

双侧指数衰减
（仅对于双边变换）

e^

2alpha over alpha ^2-s^2

−α < Re(s) < α

单位阶跃函数
的频移

指数趋近

(1-e^-alpha t)cdot u(t)

frac alpha s(s+alpha )

Re(s) > 0

单位阶跃函数
减去指数衰减

正弦

sin(omega t)cdot u(t)

omega over s^2+omega ^2

Re(s) > 0

Bracewell 1978，第227页

余弦

cos(omega t)cdot u(t)

s over s^2+omega ^2

Re(s) > 0

Bracewell 1978，第227页

双曲正弦

sinh(alpha t)cdot u(t)

alpha over s^2-alpha ^2

Re(s) > |α|

Williams 1973，第88页

双曲余弦

cosh(alpha t)cdot u(t)

s over s^2-alpha ^2

Re(s) > |α|

Williams 1973，第88页

指数衰减
正弦波

e^-alpha tsin(omega t)cdot u(t)

omega over (s+alpha )^2+omega ^2

Re(s) > −α

Bracewell 1978，第227页

指数衰减
余弦波

e^-alpha tcos(omega t)cdot u(t)

s+alpha over (s+alpha )^2+omega ^2

Re(s) > −α

Bracewell 1978，第227页

自然對數

ln(t)cdot u(t)

-1 over s,left[ln(s)+gamma right]

Re(s) > 0

Williams 1973，第88页

n 阶第一类贝塞尔函数

J_n(omega t)cdot u(t)

frac left(sqrt s^2+omega ^2-sright)^nomega ^nsqrt s^2+omega ^2

Re(s) > 0
(n > −1)

Williams 1973，第89页

误差函数

mathrm erf (t)cdot u(t)

frac 1se^frac 14s^2left(1-operatorname erf frac s2right)

Re(s) > 0

Williams 1973，第89页

注释：

u(t) 表示单位阶跃函数。

$delta (t),$ 表示狄拉克δ函数。

Γ(z) 表示Γ函数。

γ 是欧拉-马歇罗尼常数。

t 为实数，通常表示时间，
尽管它可以表示任意独立空间。

s 为复角频率，而 Re(s) 是它的实部。

α, β, τ 和 ω 是实数。

n 是整数。

变换及其性质的应用实例

拉普拉斯变换在物理学和工程中是常用的；线性时不变系统的输出可以通过卷积单位脉冲响应与输入信号来计算，而在拉氏空间中执行此计算将卷积通过转换成乘法来计算。后者是更容易解决，由于它的代数形式。

拉普拉斯变换也可以用来解决微分方程，这被广泛应用于电气工程。拉普拉斯变换把线性差分方程化简为代数方程，这样就可以通过代数规则来解决。原来的微分方程可以通过施加逆拉普拉斯变换得到其解。英国电气工程师奧利弗·黑維塞第一次提出了一个类似的计划，虽然没有使用拉普拉斯变换；以及由此产生的演算被誉为黑維塞演算。

在工程学上的应用

应用拉普拉斯变换解常变量齐次微分方程，可以将微分方程化为代数方程，使问题得以解决。在工程学上，拉普拉斯变换的重大意义在于：将一个信号从时域上，转换为复频域（s域）上来表示，对于分析系统特性，系统稳定有着重大意义；在线性系统，控制自动化上都有广泛的应用。

參考書目、資料來源

^ Korn & Korn 1967，§8.1

^ Feller 1971，§XIII.1

^ Korn & Korn 1967，第226–227页

^ Bracewell 2000，Table 14.1, p. 385

^ K.F. Riley, M.P. Hobson, S.J. Bence, Mathematical methods for physics and engineering 3rd, Cambridge University Press: 455, 2010, ISBN 978-0-521-86153-3

^ J.J.Distefano, A.R. Stubberud, I.J. Williams, Feedback systems and control 2nd, Schaum's outlines: 78, 1995, ISBN 0-07-017052-5

^ Lipschutz, S.; Spiegel, M. R.; Liu, J., Mathematical Handbook of Formulas and Tables, Schaum's Outline Series 3rd, McGraw-Hill: 183, 2009, ISBN 978-0-07-154855-7 - provides the case for real q.

^ http://mathworld.wolfram.com/LaplaceTransform.html - Wolfram Mathword provides case for complex q

電機電子類科《工程數學》，ISBN 957-584-377-0，作者陳锡冠、胡曦、周祯晖老師，高立出版社。

Korn, G. A.; Korn, T. M., Mathematical Handbook for Scientists and Engineers 2nd, McGraw-Hill Companies, 1967, ISBN 0-07-035370-0 .

[1] Korn & Korn 1967，§8.1

[2] Feller 1971，§XIII.1

[3] Korn & Korn 1967，第226–227页

[4] Bracewell 2000，Table 14.1, p. 385

[5] K.F. Riley, M.P. Hobson, S.J. Bence, Mathematical methods for physics and engineering 3rd, Cambridge University Press: 455, 2010, ISBN 978-0-521-86153-3

[6] J.J.Distefano, A.R. Stubberud, I.J. Williams, Feedback systems and control 2nd, Schaum's outlines: 78, 1995, ISBN 0-07-017052-5

[7] Lipschutz, S.; Spiegel, M. R.; Liu, J., Mathematical Handbook of Formulas and Tables, Schaum's Outline Series 3rd, McGraw-Hill: 183, 2009, ISBN 978-0-07-154855-7 - provides the case for real q.

[8] ttp://mathworld.wolfram.com/LaplaceTransform.html - Wolfram Mathword provides case for complex q

搜尋此網誌

njtynhj

拉普拉斯变换

目录

形式定义

拉普拉斯逆变换

性质和定理

与幂级数的关系

与矩的关系

证明函数导数的拉普拉斯变换

计算广义积分

与其他变换的联系

与傅里叶变换关系

与z变换的联系

拉普拉斯变换简表

变换及其性质的应用实例

在工程学上的应用

相關條目

參考書目、資料來源

Popular posts from this blog

The Dalles, Oregon

眉山市

清晰法令