拉普拉斯变换

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

F(s)=\int _{0}^{\infty }f(t)e^{-st}\,\mathrm {d} t.

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

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

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

正式定义

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

F(s)=\int _{0}^{\infty }e^{-st}f(t)\,\mathrm {d} t

其中频率参数 $s$ 是一个复数：

s=\sigma +i\omega ,\,

\sigma

和

\omega

为实数。

除了 $F$ ，有时我们也使用 $\displaystyle {\mathcal {L}}f$ 或 $\displaystyle {\mathcal {L}}_{t}\left\{f(t)\right\}$ 来表示拉普拉斯变换。 ${\mathcal {L}}$ 是一个运算符号。

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

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

{\mathcal {L}}\{\mu \}(s)=\int _{[0,\infty )}e^{-st}\,\mathrm {d} \mu (t).

以上定义的一个特殊情况是 $\mu$ 为概率测度，或者更具体地说，是狄拉克δ函数时。在运算微积分中，拉普拉斯变换的测度常常被视作由分布函数 $f$ 带来的测度。在这种情况下，为了避免混淆，一般写作

{\mathcal {L}}\{f\}(s)=\int _{0^{-}}^{\infty }e^{-st}f(t)\,\mathrm {d} t,

其中积分下限 $0^{-}$ 是

\lim _{\varepsilon \rightarrow 0}\int _{-\varepsilon }^{\infty }

的简化符号。

这个极限强调任何位于 0 的质点都被拉普拉斯变换完全捕获。虽然在使用勒贝格积分时，我们没有必要取这个极限，但它让我们更自然地与拉普拉斯–斯蒂尔吉斯变换建立联系。

拉普拉斯逆变换

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

拉普拉斯逆变换有许多不同的名称，如维奇积分、傅立叶-梅林积分、梅林逆公式，是一个複數积分：

f(t)={\mathcal {L}}^{-1}\{F\}(t)={\frac {1}{2\pi i}}\lim _{T\to \infty }\int _{\gamma -iT}^{\gamma +iT}e^{st}F(s)\,\mathrm {d} s

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

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

性质和定理

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

{\begin{aligned}f(t)&={\mathcal {L}}^{-1}\{F(s)\}\\g(t)&={\mathcal {L}}^{-1}\{G(s)\}\end{aligned}}

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

单边拉普拉斯变换的性质
	时域	s域	注释
线性叠加	$af(t)+bg(t)\$	$aF(s)+bG(s)\$	可以用积分的基本规则证明。
s域一阶微分	$tf(t)\$	$-F'(s)\$	$F'$ 是 $F$ 的一阶导数。
s域一般微分	$t^{n}f(t)\$	$(-1)^{n}F^{(n)}(s)\$	更一般的形式是 $F(s)$ 的 $n$ 阶导数。
时域一阶微分	$f'(t)\$	$sF(s)-f(0)\$	$f$ 是一个可微函数，并且其导数为指数类型。这条性质可以通过分部积分得到。
时域二阶微分	$f''(t)\$	$s^{2}F(s)-sf(0)-f'(0)\$	$f$ 为二阶可微且二阶导数是指数型的。通过对 $f(t)$ 应用微分性质可得。
时域一般微分	$f^{(n)}(t)\$	$s^{n}F(s)-\sum _{k=1}^{n}s^{k-1}f^{(n-k)}(0)\$	$f$ 为 $n$ 阶可微，其 $n$ 阶导数是指数型的。通过数学归纳法证明。
s域积分	${\frac {1}{t}}f(t)\$	$\int _{s}^{\infty }F(\sigma )\,\mathrm {d} \sigma \$	这是由s域微分和条件收敛推导出来的。
时域积分	$\int _{0}^{t}f(\tau )\,\mathrm {d} \tau =(u*f)(t)$	${1 \over s}F(s)$	$u(t)$ 是阶跃函数，注意到 $(u*f)(t)$ 是 $u(t)$ 和 $f(t)$ 的卷积。
时间标度	$f(at)$	${\frac {1}{a}}F\left({s \over a}\right)$	$a>0\$
s域平移	$e^{at}f(t)\$	$F(s-a)\$
时域平移	$f(t-a)u(t-a)\$	$e^{-as}F(s)\$	$u(t)$ 表示阶跃函数
乘法	$f(t)g(t)$	${\frac {1}{2\pi i}}\lim _{T\to \infty }\int _{c-iT}^{c+iT}F(\sigma )G(s-\sigma )\,\mathrm {d} \sigma \$	积分沿完全处在 $F$ 收敛域内的竖直线 $Re(\sigma )=c$ 。[4]
卷积	$(f*g)(t)=\int _{0}^{t}f(\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)$	${1 \over 1-e^{-Ts}}\int _{0}^{T}e^{-st}f(t)\,\mathrm {d} t$	$f(t)$ 是一个周期为 $T$ 的周期函数，于是对所有 $t\geq 0$ ，有 $f(t)=f(t+T)$ 。这条性质是时域平移和几何级数的结果。

初值定理：

f(0^{+})=\lim _{s\to \infty }{sF(s)}

，要求

{F(s)}

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

{F(s)}

分解

终值定理：

f(\infty )=\lim _{s\to 0}{sF(s)}

，要求

sF(s)

的所有极点都在左半复平面或原点为单极点。

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

F(s)

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

f(t)=e^{t}

或

f(t)=\sin(t)

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

与幂级数的关系

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

\sum _{n=0}^{\infty }a(n)x^{n}

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

\int _{0}^{\infty }f(t)x^{t}\,\mathrm {d} t

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

\int _{0}^{\infty }f(t)\left(e^{\log {x}}\right)^{t}\,\mathrm {d} t

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

\int _{0}^{\infty }f(t)e^{-st}\,\mathrm {d} t

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

与矩的关系

函数 f 的矩为

\mu _{n}=\int _{0}^{\infty }t^{n}f(t)\,\mathrm {d} t

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

\mu _{n}=(-1)^{n}{\frac {\mathrm {d} ^{n}}{\mathrm {d} s^{n}}}E\left[e^{-sX}\right].

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

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

{\begin{aligned}{\mathcal {L}}\left\{f(t)\right\}&=\int _{0^{-}}^{\infty }e^{-st}f(t)\,\mathrm {d} t\\&=\left[{\frac {f(t)e^{-st}}{-s}}\right]_{0^{-}}^{\infty }-\int _{0^{-}}^{\infty }{\frac {e^{-st}}{-s}}f'(t)\,\mathrm {d} t\quad {\text{(by parts)}}\\&=\left[-{\frac {f(0^{-})}{-s}}\right]+{\frac {1}{s}}{\mathcal {L}}\left\{f'(t)\right\},\end{aligned}}

导出

{\mathcal {L}}\left\{f'(t)\right\}=s\cdot {\mathcal {L}}\left\{f(t)\right\}-f(0^{-}),

而在双边的情形下，

{\mathcal {L}}\left\{{f'(t)}\right\}=s\int _{-\infty }^{\infty }e^{-st}f(t)\,\mathrm {d} t=s\cdot {\mathcal {L}}\{f(t)\}.

一般化的结果是

{\mathcal {L}}\left\{f^{(n)}(t)\right\}=s^{n}\cdot {\mathcal {L}}\left\{f(t)\right\}-s^{n-1}f(0^{-})-\cdots -f^{(n-1)}(0^{-}),

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

计算广义积分

令 ${\mathcal {L}}\left\{f(t)\right\}=F(s)$ ，则（参见上面的表格）

{\mathcal {L}}\left\{{\frac {f(t)}{t}}\right\}=\int _{s}^{\infty }F(p)\,\mathrm {d} p,

或

\int _{0}^{\infty }{\frac {f(t)}{t}}e^{-st}\,\mathrm {d} t=\int _{s}^{\infty }F(p)\,\mathrm {d} p.

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

\int _{0}^{\infty }{\frac {f(t)}{t}}\,\mathrm {d} t=\int _{0}^{\infty }F(p)\,\mathrm {d} p.

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

\int _{0}^{\infty }{\frac {\cos at-\cos bt}{t}}\,\mathrm {d} t=\int _{0}^{\infty }\left({\frac {p}{p^{2}+a^{2}}}-{\frac {p}{p^{2}+b^{2}}}\right)\,\mathrm {d} p={\frac {1}{2}}\left.\ln {\frac {p^{2}+a^{2}}{p^{2}+b^{2}}}\right|_{0}^{\infty }=\ln b-\ln a.

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

例子还有狄利克雷积分。

与傅里叶变换关系

连续傅里叶变换相当于计算令 $s=\imath \omega$ 或 $s=2\pi f\imath$ 的双边拉普拉斯变换：

{\begin{aligned}{\hat {f}}(\omega )&={\mathcal {F}}\left\{f(t)\right\}\\[1em]&={\mathcal {L}}\left\{f(t)\right\}|_{s=i\omega }=F(s)|_{s=i\omega }\\[1em]&=\int _{-\infty }^{\infty }e^{-\imath \omega t}f(t)\,\mathrm {d} t.\\\end{aligned}}

与z变换的联系

z 变换表达式为：

X(z)=\sum _{n=0}^{\infty }x[n]z^{-n}

其中 $z\leftarrow e^{sT}\$ 。比较两者表达式有：

X_{q}(s)=X(z){\Big |}_{z=e^{sT}}.

拉普拉斯变换简表

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

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

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

{\mathcal {L}}\left\{f(t)+g(t)\right\}={\mathcal {L}}\left\{f(t)\right\}+{\mathcal {L}}\left\{g(t)\right\}

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

{\mathcal {L}}\left\{af(t)\right\}=a{\mathcal {L}}\left\{f(t)\right\}

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

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

函数

时域

f(t)={\mathcal {L}}^{-1}\left\{F(s)\right\}

拉普拉斯s域

F(s)={\mathcal {L}}\left\{f(t)\right\}

收敛区域

参考

单位脉冲函数

\delta (t)\

1

所有

s

检验

延迟脉冲函数

\delta (t-\tau )\

e^{-\tau s}\

单位脉冲函数
的时移

单位阶跃函数

u(t)\

{1 \over s}

\Re (s)>0

对单位冲激函数积分

延迟单位阶跃函数

u(t-\tau )\

{\frac {1}{s}}e^{-\tau s}

\Re (s)>0

单位阶跃函数
的时移

斜坡函数

t\cdot u(t)\

{\frac {1}{s^{2}}}

\Re (s)>0

两次积分
单位脉冲函数

n次幂
（n 为整数）

t^{n}\cdot u(t)

{n! \over s^{n+1}}

\Re (s)>0{\text{ for }}n>-1

n次积分
单位阶跃函数

q次幂
（q 为复数）

t^{q}\cdot u(t)

{\Gamma (q+1) \over s^{q+1}}

\Re (s)>0

Re(q) > −1

[7][8]

n次方根

{\sqrt[{n}]{t}}\cdot u(t)

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

\Re (s)>0

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

频移的 n 次方

t^{n}e^{-\alpha t}\cdot u(t)

{\frac {n!}{(s+\alpha )^{n+1}}}

Re(s) > −α

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

延迟的频移的 n 次方

(t-\tau )^{n}e^{-\alpha (t-\tau )}\cdot u(t-\tau )

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

Re(s) > −α

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

指数衰减

e^{-\alpha t}\cdot u(t)

{1 \over s+\alpha }

Re(s) > −α

单位阶跃函数
的频移

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

e^{-\alpha |t|}\

{2\alpha  \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 t}\sin(\omega t)\cdot u(t)\

{\omega  \over (s+\alpha )^{2}+\omega ^{2}}

Re(s) > −α

Bracewell 1978，第227頁

指数衰减
余弦波

e^{-\alpha t}\cos(\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}}}-s\right)^{n}}{\omega ^{n}{\sqrt {s^{2}+\omega ^{2}}}}}

\Re (s)>0

(n > −1)

Williams 1973，第89頁

误差函数

\mathrm {erf} (t)\cdot u(t)

{\frac {1}{s}}e^{{\frac {1}{4}}s^{2}}\left(1-\operatorname {erf} {\frac {s}{2}}\right)

\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, 3rd, Cambridge University Press: 455, 2010, ISBN 978-0-521-86153-3
J.J.Distefano, A.R. Stubberud, I.J. Williams, 2nd, Schaum's outlines: 78, 1995, ISBN 0-07-017052-5
Lipschutz, S.; Spiegel, M. R.; Liu, J., , 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., 2nd, McGraw-Hill Companies, 1967, ISBN 0-07-035370-0.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[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, 3rd, Cambridge University Press: 455, 2010, ISBN 978-0-521-86153-3

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

[7] Lipschutz, S.; Spiegel, M. R.; Liu, J., , 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

拉普拉斯变换

正式定义

拉普拉斯逆变换

性质和定理

与幂级数的关系

与矩的关系

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

计算广义积分

与傅里叶变换关系

与z变换的联系

拉普拉斯变换简表

变换及其性质的应用实例

在工程学上的应用

相關條目

參考書目、資料來源