信号与系统(Python) 学习笔记 (6) 拉普拉斯变换 Laplace Transform

最新推荐文章于 2024-01-31 21:22:03 发布

Varalpha

最新推荐文章于 2024-01-31 21:22:03 发布

阅读量7.5k

点赞数 10

分类专栏： # 信号与系统文章标签：信号处理傅立叶分析

本文链接：https://blog.csdn.net/varalpha/article/details/105096437

版权

信号与系统专栏收录该内容

13 篇文章 85 订阅

订阅专栏

6. 拉普拉斯变换

6.1. 拉普拉斯变换 Laplace Transform

6.1.1 双边拉普拉斯变换的定义

有些函数不满足绝对可积条件，求解傅里叶变换困难。为此，可用一衰减因 ${\color{blue}e^{-\sigma t}}$ ( ${\color{red}\sigma}$ 为实常数）乘信号 $f (t)$ , 适当选取 $\sigma$ 的值，使乘积信号 $e^{-\sigma t}$ 当 $t\to \infty$ 时信号幅度趋近于 $0$ , 从而使 $e^{-\sigma t}$ 的傅里叶变换存在。
$\begin{aligned}F_b ({\color{red}\sigma} + j \omega) & =\mathfrak{F}\big[ f(t) {\color{red}e^{-\sigma t} }\big] \\ & = \int^{\infty}_{-\infty}f(t) {\color{red}e^{-\sigma t}} e^{-j\omega t}dt \\ & = \int^{\infty}_{-\infty}f(t) e^{-({\color{red}\sigma} + j\omega)t} dt \end{aligned}$
相应的傅里叶逆变换为:
$f(t){\color{red}e^{-\sigma t}} = \frac{1}{2\pi} \int^{\infty}_{-\infty} F_b ({\color{red}\sigma} + j\omega) e^{j\omega t} d \omega$
$\frac{1}{2\pi} \int^{\infty}_{-\infty} F_b ({\color{red}\sigma} + j\omega) e^{({\color{red}\sigma} +j\omega) t} d \omega$
令 ${\color{red}s = \sigma + j\omega}, \; d \omega = ds/j$ 有:
$F_b ({\color{red}s}) = \int^{\infty}_{-\infty}f(t) e^{-{\color{red}s}t} dt$
$f({\color{red}t}) = \frac{1}{2\pi {\color{blue}j}} \int^{{\color{blue}\sigma+ j}\infty}_{{\color{blue}\sigma -j}\infty} F_b ({\color{blue}s}) e^{{\color{blue}s}{\color{red} t}} d {\color{blue}s}$
$F_b(s)$ 称为 $f (t)$ 的双边拉氏变换（或象函数），
$f (t)$ 称为 $F_b(s)$ 的双边拉氏逆变换（或原函数）

6.1.2 收敛域

只有选择适当的 $\sigma$ 值才能使积分收敛，信号 $f (t)$ 的双边拉普拉斯变换存在

收敛域：使 $f (t)$ 拉氏变换存在的 $\sigma$ 取值范围。
例1: 因果信号 $f_1(t) = e^{\alpha t} \varepsilon (t)$ , 求其拉普拉斯变换:
$\begin{aligned}F_{1b}(s) & = \int^{\infty}_0 e^{\alpha t} e^{-st}dt\\& = \frac{1}{s-\alpha} \big[ 1-\lim_{t\to\infty} e^{-(\sigma-\alpha)t} e^{-j\omega t}\big] \\ & = \begin{cases}\frac{1}{s-\alpha} ,\; &\mathcal{Re}[s] = \sigma > \alpha \\ 不定, \; &\mathcal{Re}[s] = \sigma = \alpha \\ 无界 ,\; &\mathcal{Re}[s] = \sigma < \alpha \end{cases}\end{aligned}$
- 可见，对于因果信号，仅当 $\mathcal{Re}[s]=\sigma>\alpha$ 时，其拉氏变换存在。收敛域如图所示。
例2: 反因果信号 $f_2(t) = e^{\beta t} \varepsilon (-t)$ , 求其拉普拉斯变换:
$\begin{aligned}F_{2b}(s) & = \int^0_{-\infty} e^{\beta t} e^{-st}dt\\& = \frac{-1}{s-\beta} \big[ 1-\lim_{t\to-\infty} e^{-(\sigma-\beta)t} e^{-j\omega t}\big] \\ & = \begin{cases}\frac{-1}{s-\beta} ,\; &\mathcal{Re}[s] = \sigma < \beta \\ 不定, \; &\mathcal{Re}[s] = \sigma = \beta \\ 无界 ,\; &\mathcal{Re}[s] = \sigma > \beta \end{cases}\end{aligned}$
- 可见，对于反因果信号，仅当 $\mathcal{Re}[s]=\sigma<\beta$ 时，其拉氏变换存在。收敛域如图所示。
例3: 双边信号
$\begin{aligned}f_3(t) = f_1(t) + f_2(t) = \begin{cases}e^{\beta t}, & t<0 \\ e^{\alpha t} &t>0 \end{cases}\end{aligned}$
- 仅当 $\beta > \alpha$ 其收敛域为 $\alpha< \mathcal{Re} [s] <\beta$ 的一个带状区域，如图所示。
双边拉氏变换必须标出收敛域。
对于双边拉普拉斯变换而言， $F_b(s)$ 和收敛域一起，可以唯一地确定 $f (t)$ 。即
$\overset{\text{一一对应}}{\longleftrightarrow} F_b(S) + {\color{blue} \text{收敛域}}$
不同的信号可以有相同的 $F_b(s)$ ，但收敛域不同。

6.1.3 单边拉氏变换的定义

通常遇到的信号都有初始时刻，不妨设其初始时刻为坐标原点。
这样, $t < 0$ 时, $f (t) = 0$ 。从而拉氏变换式写为:
$\int^{\infty}_{0_-} f(t) e^{-st} dt$
- 称为单边拉氏变换。简称拉氏变换。
- 其收敛域一定是 $Re[s]>\alpha$ ，可以省略。
$\mathfrak{L}[f(t)]$
$\mathfrak{L}[f(t)] = F(s) \overset{\text{def}}{=} \int^{\infty}_{0} f(t) e^{-st} dt$
$\mathfrak{L}^{-1}[F(s)]$
$\mathfrak{L}^{-1}[F(s)] =f(t) \overset{\text{def}}{=} \Big[\frac{1}{2\pi j} \int^{\sigma+ j\infty}_{\sigma -j\infty} F_b (s) e^{s t} d s\Big]\varepsilon(t)$
${\color{red}f(t)\overset{\text{1-to-1}}{\longleftrightarrow} F(s)}$

6.1.4 单边拉氏变换与傅里叶变换的关系

$\int^{\infty}_{0} f(t) e^{-st} dt, \; \mathcal{Re}[s] >\sigma_0$
$F(j\omega) = \int^{\infty}_{0} f(t) e^{-j\omega t} dt$

要讨论其关系， $f (t)$ 必须为因果信号:
根据收敛坐标 $\sigma_0<0$ 的值可分为以下三种情况:
1. $\sigma_0<0$ ，即 $F (s)$ 的收敛域包含 $j\omega$ 轴，则 $f (t)$ 的傅里叶变换存在，并且
  $F(j\omega)=F(s)\big\vert_{s=j\omega}$
  - 如 $e^{-2t} \varepsilon(t) \longleftrightarrow F(s) = 1/(s+2), \; \sigma >-2$
  - 则 $F(j\omega) = 1/(j\omega+2)$
2. $\sigma_0=0$ ，即 $F (s)$ 的收敛边界为 $j\omega$ 轴，则
  $F(j\omega)=\lim_{\sigma\to0}F(s)$
  - 如 $\varepsilon(t) \longleftrightarrow F(s) = 1/s$
  - 则
    $\begin{aligned} F(j\omega) &=\lim_{\sigma\to0}\frac{1}{\sigma+j\omega}\\ & = \lim_{\sigma\to0} \frac{\sigma}{\sigma^2 + \omega^2} + \lim_{\sigma\to0}\frac{-j\omega}{\sigma^2+\omega^2} \\&= \pi \delta(\omega) +\frac{1}{j\omega}\end{aligned}$
3. $\sigma_0>0$ ，即 $F(j\omega)$ 不存在。
  - 如 $e^{2t} \varepsilon(t) \longleftrightarrow F(s) = 1/(s-2), \; \sigma >2$
  - 则其傅里叶变换 $F(j\omega)$ 不存在

6.1.5 常见信号的拉普拉斯变换

$\begin{aligned} \displaystyle \delta(t) \longleftarrow & \longrightarrow 1,\; \sigma > -\infty \\ \varepsilon(t)\longleftarrow & \longrightarrow \frac{1}{s},\; \sigma > 0 \\ 1 \longleftarrow & \longrightarrow \frac{1}{s},\; \sigma >0 \\ e^{s_0 t} \longleftarrow & \longrightarrow \frac{1}{s-s_0}, \; \sigma > \mathcal{Re}[s_0] \\ \cos(\omega_0 t) = \frac{e^{j\omega_0 t} + e^{-j\omega_0 t}}{2} \longleftarrow & \longrightarrow \frac{s}{s^2 +\omega_0^2} \\ \sin(\omega_0 t) = \frac{e^{j\omega_0 t} - e^{-j\omega_0 t}}{2j} \longleftarrow & \longrightarrow \frac{\omega_0}{s^2 +\omega_0^2} \\ f_T(t) \longleftarrow & \longrightarrow \frac{1}{1-e^{-sT}} \int^{T}_{0} f_T(t) e^{-st} dt \\ \delta_T(t) \longleftarrow & \longrightarrow \frac{1}{1-e^{-sT}} \end{aligned}$

周期信号 $f_T(t)$ 解释:
$\begin{aligned}F_T(s) & = \int^{\infty}_0 f_T(t) e^{-st} dt \\ &= \int^{T}_0 f_T(t) e^{-st} dt + \int^{2T}_T f_T(t) e^{-st} dt + \cdots \\ & = \sum^{\infty}_{n=0} \int^{(n+1)T}_{nT} f_T(t) e^{-st} dt\end{aligned}$
令 $t = t + n T$ :
$\sum^{\infty}_{n=0} e^{-nsT} \int^{T}_{0} f_T(t) e^{-st} dt = \frac{1}{1-e^{-sT}} \int^{T}_{0} f_T(t) e^{-st} dt$

6.2. 拉普拉斯变换的性质

6.2.1 线性性质

若
$\begin{aligned} \displaystyle f_1(t) \longleftarrow & \longrightarrow F_1(s),\; & \mathcal{Re}[s]>\sigma_1 \\ f_2(t) \longleftarrow & \longrightarrow F_2(s),\; & \mathcal{Re}[s]>\sigma_2 \\ a_1f_1(t) + a_2f_2(t) \longleftarrow & \longrightarrow a_1F_1(s) + a_2F_2(s) ,\;& \mathcal{Re}[s]>\max(\sigma_1,\sigma_2)\\ \end{aligned}$

6.2.2 尺度变换

若
$\begin{aligned} \displaystyle f(t) \longleftarrow & \longrightarrow F(s), \; &{\color{red} 实数 \alpha >0} ,\; \mathcal{Re}[s]>\sigma_0\\ f(\alpha t) \longleftarrow & \longrightarrow \frac{1}{\alpha}F(\frac{s}{\alpha}) ,\; &\mathcal{Re}[s]>\alpha\sigma_0\\ \end{aligned}$

6.2.3 时移性质

若
$\begin{aligned} \displaystyle f(t) \longleftarrow & \longrightarrow F(s),\; \mathcal{Re}[s]>\sigma_0,\; &{\color{red} 实常数 t_0 >0} ,\; \mathcal{Re}[s]>\sigma_0\\ f(t-t_0){\color{red}\varepsilon(t-t_0)} \longleftarrow & \longrightarrow e^{-st_0} F(s) ,\; & \mathcal{Re}[s]>\sigma_0\\ f(\alpha t-t_0){\color{red}\varepsilon(\alpha t-t_0)} \longleftarrow & \longrightarrow \frac{1}{\alpha} e^{-\frac{t_0}{\alpha} s} F(\frac{s}{\alpha}) ,\; &{\color{red} 实数 \alpha >0} ,\; \mathcal{Re}[s]>\sigma_0\\ \end{aligned}$
若 $f (t)$ 为 因果信号,
$f(t-t_0)\longleftarrow \longrightarrow e^{-st_0} F(s)$

6.2.4. 复频移特性

若
$\begin{aligned} \displaystyle f(t) \longleftarrow & \longrightarrow F(s), \; &{\color{red} 复常数 s_\alpha = \sigma_\alpha+ j\omega_\alpha} ,\; \mathcal{Re}[s]>\sigma_0\\ f(t)e^{s_\alpha t} \longleftarrow & \longrightarrow F(s-s_\alpha),\; &\mathcal{Re}[s]>\sigma_0+\sigma_\alpha\\ \end{aligned}$

6.2.5. 时域微分特性

若
$\begin{aligned} \displaystyle f(t) \longleftarrow & \longrightarrow F(s), \; &\mathcal{Re}[s]>\sigma_0\\ f^{\prime}(t) \longleftarrow & \longrightarrow sF(s)-f(0_-)\\ f^{\prime\prime}(t) \longleftarrow & \longrightarrow s^2 F(s)-sf(0_-)-f^{\prime}(0_-)\\\end{aligned}$
若 $f (t)$ 为 因果信号，则
$f^{(n)} (t) \longleftarrow \longrightarrow s^n F(s)$

6.2.6. 时域积分特性

若
$\begin{aligned} \displaystyle f(t) \longleftarrow & \longrightarrow F(s), \; &\mathcal{Re}[s]>\sigma_0\\ \int^{t}_{0_-} f(x)dx \longleftarrow & \longrightarrow \frac{1}{s}F(s)\\ \Big(\int^{t}_{0_-}\Big)^n f(x)dx \longleftarrow & \longrightarrow \frac{1}{s^n}F(s)\\ f^{(-1)}(t) = \int^{t}_{-\infty} f(x)dx \longleftarrow & \longrightarrow s^{-1}F(s)+s^{-1}f^{(-1)}(0_-)\\\end{aligned}$
若 $f (t)$ 为 因果信号，则
$\longleftarrow \longrightarrow \frac{F_n(s)}{s^n}$

6.2.7. 复频域微分和积分

若
$\begin{aligned} \displaystyle f(t) \longleftarrow & \longrightarrow F(s), \; &\mathcal{Re}[s]>\sigma_0\\ (-t) f(t) \longleftarrow & \longrightarrow \frac{d F(s)}{ds}\\ (-t)^n f(t) \longleftarrow & \longrightarrow \frac{d^n F(s)}{d s^n}\\ \frac{f(t)}{t} \longleftarrow & \longrightarrow \int^{\infty}_{s} F(\eta)d\eta\\ \end{aligned}$

6.2.8. 时域卷积定理

若因果函数:
$\begin{aligned} \displaystyle f_1(t) \longleftarrow & \longrightarrow F_1(s),\; & \mathcal{Re}[s]>\sigma_1 \\ f_2(t) \longleftarrow & \longrightarrow F_2(s),\; & \mathcal{Re}[s]>\sigma_2 \\ f_1(t) \star f_2(t) \longleftarrow & \longrightarrow F_1(s) \cdot F_2(s) ,\;& \mathcal{Re}[s]>\max(\sigma_1,\sigma_2)\\ \end{aligned}$

6.2.9. 复频域卷积定理

若因果函数:
$\begin{aligned} \displaystyle f_1(t) \longleftarrow & \longrightarrow F_1(s),\; & \mathcal{Re}[s]>\sigma_1 \\ f_2(t) \longleftarrow & \longrightarrow F_2(s),\; & \mathcal{Re}[s]>\sigma_2 \\ f_1(t) \cdot f_2(t) \longleftarrow & \longrightarrow \frac{1}{2\pi j} \int^{c+j\infty}_{c-j\infty} F_1(\eta) \cdot F_2(s-\eta)d\eta ,\;& \mathcal{Re}[s]>\max(\sigma_1,\sigma_2)\\ \end{aligned}$

6.2.10. 初值终值定理

初值定理和终值定理常用于由 $F (s)$ 直接求 $f (0 +)$ 和 $f(\infty)$ ,而不必求出原函数 $f (t)$ 。

初值定理:
- 设函数 $f (t)$ 不含 $\delta(t)$ 及其各阶导数（即 $F (s)$ 为真分式，若 $F (s)$ 为假分式化为真分式），则:
  $f(0_+) = \lim_{t\to0_+} f(t) = \lim_{s\to\infty} s F(s)$
终值定理:
- 若 $f (t)$ ，当 $t\to \infty$ 时存在, 并且 $\leftrightarrow F(s)$ , $\mathcal{Re}[s]>\sigma_0$ , $\sigma_0<0$ , 则:
  $f(\infty) =\lim_{s\to 0} sF(s)$

6.3. 拉普拉斯反变换

$\mathfrak{L}^{-1}[F(s)] \longleftrightarrow \mathfrak{L} [f(t)] = F(s)$

6.3.1. 拉普拉斯反变换

直接利用定义式求反变换—复变函数积分，比较困难。
通常的方法:
1. 查表；
2. 利用性质；
3. 部分分式展开 --结合。
若象函数 $F (s)$ 是 $s$ 的有理分式, 可写为:
${\color{blue}F(s) = \displaystyle \frac{b_m s^m + b_{m-1}s^{m-1} + \cdots + b_1 s + b_0}{s^n + a_{n-1}s^{n-1} + \cdots +a_1s+a_0}}$
若 $m\geq n$ （假分式），可用多项式除法将象函数 $F (s)$ 分解为
${\color{blue}有理多项式 P(s)+ 有理真分式}$
${\color{blue}F(s) = P(s) + \frac{B_0(s)}{A(s)}}$
$P (s)$ 的拉普拉斯逆变换由冲激函数及其各阶导数构成。
- 例: $P(s)\to a_1 s^2 + a_2 s + a_3 \to a_1\delta^{\prime\prime}(t) + a_2 \delta^{\prime}(t) +a_3\delta(t)$
下面主要讨论有理真分式。

6.3.2. 部分分式展开法

若 $F (s)$ 是 $s$ 的实系数有理真分式 $(m < n)$ ，则
${\color{blue}F(s) = \displaystyle \frac{B(s)}{A(s)} = \displaystyle \frac{b_m s^m + b_{m-1}s^{m-1} + \cdots + b_1 s + b_0}{s^n + a_{n-1}s^{n-1} + \cdots +a_1s+a_0}}$
- 式中 $A (s)$ 称为 $F (s)$ 的特征多项式 (characteristic polynomial)，方程 $A (s) = 0$ 称为特征方程，它的根称为特征根，也称为 $F (s)$ 的固有频率（或自然频率）。 $n$ 个特征根 $p_i$ 称为 $F (s)$ 的极点

$F (s)$ 为单极点（单根）
$\displaystyle \frac{B(s)}{A(s)} = \frac{K_1}{s-p_1}+\frac{K_2}{s-p_2}+ \cdots + \frac{K_i}{s-p_i}+ \cdots + \frac{K_n}{s-p_n}$
$K_i = (s - p_i) F(s) \big\vert _{s=pi}$
$\mathfrak{L}^{-1} \big[ \frac{1}{s-p_i} \big] = e^{p_i t} \varepsilon(t)$
- 特例 $F (s)$ 包含共轭复根时 ( $p_{1,2} = -\alpha\pm j\beta$ ):
  $\begin{aligned} F(s) &= \displaystyle \frac{B(s)}{D(s)[(s+\alpha)^2 + \beta^2]}\\ &= \displaystyle \frac{B(s)}{D(s)(s+\alpha-j\beta)(s+\alpha+j\beta)}\\ &= \displaystyle \frac{K_1}{s+\alpha-j\beta}+\frac{K_2}{s+\alpha+j\beta}+F_2(s)\\ K_1 &= [(s+\alpha - j\beta)F(s)]\big\vert_{s=-\alpha +j\beta} \\ &= \lvert K_1\rvert e^{j\theta} \\ &=A+jB \\ K_2 &= K_1^* = \lvert K_1 \rvert e^{-j\theta} = A -jB\\ F_1(s) &= \displaystyle \frac{K_1}{s+\alpha-j\beta}+\frac{K_2}{s+\alpha+j\beta}\\ & = \displaystyle \frac{\lvert K_1\rvert e^{j\theta}}{s+\alpha-j\beta}+\frac{\lvert K_1 \rvert e^{-j\theta}}{s+\alpha+j\beta}\\ f_1(t) &= 2 \lvert K_1\rvert e^{-\alpha t} \cos(\beta t + \theta) \varepsilon(t)\end{aligned}$
  $K_{1,2} = A \pm jB, \; f_1(t) = 2 e^{-\alpha t} [A \cos(\beta t) - B \sin(\beta t) ] \varepsilon(t)$
$F (s)$ 有重极点（重根）
- 若 $A (s) = 0$ 在 $s=p_1$ 处有 $r$ 重根,
  $\displaystyle \frac{B(s)}{A(s)} = \frac{K_{11}}{(s-p_1)^r}+ \frac{K_{12}}{(s-p_1)^r}+ \cdots + \frac{K_{1r}}{(s-p_1)^r}$
  $K_{11} = \displaystyle [(s - p_1)^r F(s)] \big\vert _{s=p1}$
  $K_{12} = \displaystyle \frac{d[(s - p_1)^r F(s)]}{ds} \big\vert _{s=p1}$
  $K_{1i} = \displaystyle \frac{1}{(i-1)!} \frac{d^{i-1}}{ds^{i-1}}[(s - p_1)^r F(s)] \Big\vert _{s=p_1}$
  $\mathfrak{L}[t^n \varepsilon(t)] = \frac{n!}{s^{n+1}}$
  $\mathfrak{L}^{-1}\big[\frac{1}{(s-p_1)^{n+1}}\big] = \frac{1}{n!} t^n e^{p_1 t} \varepsilon(t)$

To TOP 至目录

Varalpha

关注

10
点赞
踩
35

收藏

觉得还不错? 一键收藏
0
评论
信号与系统(Python) 学习笔记 (6) 拉普拉斯变换 Laplace Transform

6. 拉普拉斯变换6.1. 拉普拉斯变换 Laplace Transform6.1.1 双边拉普拉斯变换的定义有些函数不满足绝对可积条件，求解傅里叶变换困难。为此，可用一衰减因 e−σt{\color{blue}e^{-\sigma t}}e−σt (σ{\color{red}\sigma}σ 为实常数）乘信号 f(t)f(t)f(t), 适当选取 σ\sigmaσ 的值，使乘积信号...
复制链接

扫一扫