python实现拉普拉斯函数 怎么进行拉普拉斯变换 如何进行拉普拉斯变换

由于自动控制研究的是动态过程，因此常常用微分方程来描述。非线性微分方程的求解十分困难，因此常常在正常工作点附近将非线性微分方程线性化来得到线性微分方程，再通过拉普拉斯变换将微分方程转换为代数方程进行求解。

在自控原理基础中，拉普拉斯变换（拉氏变换）作为一个工具，在大多数情况下不需要考虑数学意义上严谨的使用条件。在此不介绍具体知识，仅给出常用的拉氏变换对及性质。

拉普拉斯变换：

函数$f(t)$的拉氏变换：$F(s) = \int_{0_-}^\infty f(t)e^{-st}dt = \mathscr{L}\{f(t)\}$

函数$F(s)$的拉氏逆变换： $f(t)=\frac{1}{2\pi j} \int_{c -j\infty}^{c +j\infty}F(s)e^{st}ds=\mathscr{L}^{-1}\{F(s)\}\ \ \ \ \ \ t>0$

注意此处有隐含要求：当 $t<0$ 时，$f(t)=0$.控制理论中涉及的时间函数 $f(t)$ 在时间为负时总是等于 $0$，因此采用的是单边拉普拉斯变换。有时积分下限设为 $0_-$ ，这是是为了处理有些时间函数在 $t=0$ 的点上不连续的情况。
在本篇中为描述方便，也记
$\mathscr{L}_+\{f(t)\}=\int_{0_+}^{\infty}f(t)e^{-st}dt \ \ \ \ \ , \ \ \ \ \ \mathscr{L}_-\{f(t)\}=\int_{0_-}^{\infty}f(t)e^{-st}dt$

对于自控原理基础的内容，在一般的情况下很少需要通过定义来计算时间函数的拉氏变换，大多数情况下只需要根据常用的拉氏变换对以及性质计算即可。

常用拉氏变换对

$f(t)$	$F(s)$
脉冲函数 $\delta(t)$	1
阶跃函数 $u(t)$	$\frac{1}{s}$
$e^{-at}·u(t)$	$\frac{1}{s+a}$
$sin(\omega{t})·u(t)$	$\frac{\omega}{s^2+\omega^2}$
$cos(\omega{t})·u(t)$	$\frac{s}{s^2+\omega^2}$
$e^{-at}sin(\omega{t})·u(t)$	$\frac{\omega}{(s+a)^2+\omega^2}$
$e^{-at}cos(\omega{t})·u(t)$	$\frac{s+a}{(s+a)^2+\omega^2}$
$t·u(t)$	$\frac{1}{s^2}$
$\frac{1}{2}t^2·u(t)$	$\frac{1}{s^3}$
$t^n·u(t)$	$\frac{n!}{s^{n+1}}$

常用性质

1. 线性定理
   若 $\mathscr{L}\{f_1(t)\}=F_1(s)\ ,\ \mathscr{L}\{f_2(t)\}=F_2(s)\ ,\ \alpha、\beta$ 为常数，则有：$\mathscr{L}\{\alpha f_1(t) \pm \beta f_2(t)\}=\alpha F_1(s) \pm \beta F_2(s)$
2. 比例定理
   若$\mathscr{L}\{f(t)\}=F(s)\ ,a$为常数，则有：$\mathscr{L}\{f\left(\frac{t}{a}\right)\}=aF(as)\ \ ,\ \ \mathscr{L}\{f(at)\}=\frac{1}{a}F\left(\frac{s}{a}\right)$
3. 时域平移定理
   若$\mathscr{L}\{f(t)\}=F(s)\ ,\ \tau$为正常数，将$f(t)$纯滞后$\tau$秒：$f(t-\tau)u(t-\tau)=\left\{ \begin{array}{c} 0 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ t<\tau \\ f(t-\tau) \ \ \ \ t>\tau\end{array}\right.$ ， 则有：$\mathscr{L}\{f(t-\tau)u(t-\tau)\}=e^{-\tau s}F(s)$
4. 复域平移定理
   若$\mathscr{L}\{f(t)\}=F(s)\ ,\ \tau$为实数或负数 ，则有：$\mathscr{L}\{e^{-\tau t}f(t)\}=F(s+\tau)$
5. 实微分定理
   若$\mathscr{L}\{f(t)\}=F(s)$ ，$f(t)$的$1到n$阶导数分别存在，且其拉氏变换也都存在，则有：$\mathscr{L}_\pm\left\{\frac{d^n}{dt^n}f(t)\right\}=s^nF(s)-\sum_{k=1}^ns^{n-k}f^{(k-1)}(0^\pm)$
6. 实积分定理
   若$\mathscr{L}\{f(t)\}=F(s)$ ，则有多重积分的拉氏变换：$\mathscr{L}_\pm\left\{\int…\int f(t)(dt)^n\right\}=\frac{1}{s^n}F(s)+\sum_{k=1}^n\frac{1}{s^{n-k+1}}\left[\int…\int f(t)(dt)^k\right]_{t=0^\pm}$多重定积分的拉氏变换：$\mathscr{L}\left\{\int_{0^\pm}^t…\int_{0^\pm}^tf(t)(dt)^n\right\}=\frac{1}{s^n}F(s)$
7. 复微分定理
   若$\mathscr{L}\{f(t)\}=F(s)$ ，则除$F(s)$的极点外，有：$\mathscr{L}\{t^nf(t)\}=(-1)^n\frac{d^n}{ds^n}F(s)$
8. 复积分定理
   若$\mathscr{L}\{f(t)\}=F(s)$ ，则除$F(s)$的极点外，有：$\mathscr{L}\left\{\frac{1}{t^n}f(t)\right\}=\int_s^\infty…\int_s^\infty F(s) (ds)^n$
9. 卷积定理
   函数$f_1(t)$与函数$f_2(t)$的卷积运算为：$f_1(t)*f_2(t)=\int_0^tf_1(t-\tau)f_2(\tau)d\tau=\int_0^tf_1(\tau)f_2(t-\tau)d\tau=f_2(t)*f_1(t)$若$\mathscr{L}\{f_1(t)\}=F_1(s)$ ，$\mathscr{L}\{f_2(t)\}=F_2(s)$ ，则有：$\mathscr{L}\{f_1(t)*f_2(t)\}=F_1(s)F_2(s)$
10. 函数乘积的拉氏变换
    若$\mathscr{L}\{f(t)\}=F(s)$ ，$\mathscr{L}\{g(t)\}=G(s)$，则有：$\mathscr{L}\{f(t)g(t)\}=\frac{1}{2\pi j}\int_{c-j\infty}^{c+j\infty}F(p)G(s-p)dp$
11. 终值定理
    若$\mathscr{L}\{f(t)\}=F(s)$ ，$sF(s)$的所有极点均位于左半开平面，则$\lim_{t\to\infty}f(t)$存在，成立终值定理：$\lim_{t\to\infty}f(t)=f(\infty)=\lim_{s\to0}[sF(s)]$
12. 初值定理
    若$f(t)、f^{$的拉氏变换存在、$\lim_{s\to\infty}[sF(s)]$存在，则有初值定理：$\lim_{t\to0^+}f(t)=f(0^+)=\lim_{s\to\infty}[sF(s)]$

拉氏反变换

除根据常见拉氏变换对、使用拉氏变换公式和性质进行拉式反变换外，一般用部分分式展开法。
设 $F(s)$ 是关于 $s$ 的有理真分式，$F(s)=\frac{B(s)}{A(s)}=\frac{b_ms^m+b_{m-1}s^{m-1}+···+b_1s+b_0}{a_ns^n+a_{n-1}s^{n-1}+···+a_1s+a_0}\ \ \ \ \ \ m<n$其中 $b_m、b_{m-1}、···、b_1、b_0，a_n、a_{n-1}、···、a_1、a_0$ 是实常数，$m、n$ 是正整数，将 $F(s)$

情况一：特征方程 $A(s)=0$
在此情况下 $F(s)$ 可以展开成：$F(s)=\frac{B(s)}{A(s)}=\frac{c_1}{s-s_1}+···+\frac{c_{n-1}}{s-s_{n-1}}+\frac{c_n}{s-s_n}=\sum_{i=1}^n\frac{c_i}{s-s_i}\ \ \ \ i=1,···,n$其中 $s_i$ 为特征方程 $A(s)=0$ 的第 $i$ 个单根，$c_i$ 为 $F(s)$ 在 $s_i$ 处的留数，其计算方法如下：$c_i=\left\{(s-s_i)F(s)\right\}_{s=s_i}$注意这里的 $c_i$

情况二：特征方程$A(s)=0$有重根
设 $A(s)$ 有 $r$ 重根 $s_1$ ,则在此情况下有：$F(s)=\frac{B(s)}{A(s)}=\frac{B(s)}{(s-s_1)^r(s-s_{r+1})···(s-s_n)}$于是 $F(s)$ 可以展开成：$F(s)=\frac{c_1}{s-s_1}+···+\frac{c_r}{(s-s_1)^r}+\frac{c_{r+1}}{s-s_{r+1}}+··+\frac{c_{n-1}}{s-s_{n-1}}+\frac{c_n}{s-s_n}=\sum_{j=1}^r\frac{c_j}{(s-s_1)^j}+\sum_{i=r+1}^n\frac{c_i}{s-s_i}$其中 $j=1,···,r\ ;\ i=r+1,···,n$

留数 $c_i$ 按照情况一中方法计算，留数 $c_j$ 的计算方法如下：$c_r=\left\{(s-s_1)^rF(s)\right\}_{s=s_1}$$c_{r-1}=\left\{\frac{d}{ds}[(s-s_1)^rF(s)]\right\}_{s=s_1}$$c_{r-j}=\frac{1}{j!}\left\{\frac{d^j}{ds^j}[(s-s_1)^rF(s)]\right\}_{s=s_1}$$c_1=\frac{1}{(r-1)!}\left\{\frac{d^{r-1}}{ds^{r-1}}[(s-s_1)^rF(s)]\right\}_{s=s_1}$注意这里的 $c_i、c_j$

以情况二为一般情况，可以得到 $F(s)$ $f(t)=\mathscr{L}^{-1}\{F(s)\}=\mathscr{L}^{-1}\left\{\sum_{j=1}^r\frac{c_j}{(s-s_1)^j}+\sum_{i=r+1}^n\frac{c_i}{s-s_i}\right\}$$f(t)=\mathscr{L}^{-1}\{F(s)\}=\left\{c_1+c_2t+···+\frac{c_{r-1}}{(r-2)!}t^{r-2}+\frac{c_r}{(r-1)!}t^{r-1}\right\}e^{s_1t}+\sum_{i=r+1}^nc_ie^{s_it}$

结语

最后需要说明的是：
有些函数是无法进行拉普拉斯变换的（在数学上，函数 $f(t)$ 存在拉氏变换的条件是拉普拉斯积分收敛），
但是物理上可以实现的信号总是存在拉氏变换的。