5.拉普拉斯变换

最新推荐文章于 2025-04-30 18:18:59 发布

恋日

最新推荐文章于 2025-04-30 18:18:59 发布

阅读量4.7k

点赞数 2

CC 4.0 BY-SA版权

分类专栏：给工程师的微分方程课

本文链接：https://blog.youkuaiyun.com/yahhqy/article/details/106816615

给工程师的微分方程课专栏收录该内容

11 篇文章

订阅专栏

0.引言

拉普拉斯变换可以处理常系数ode当非齐次项为非连续或者脉冲函数的情形。

1.定义

函数 $f (x)$ 的拉普拉斯变换 $F(s)=\mathcal{L}\{f(t)\}$ 按照以下积分进行定义：
$F(s)=\int_{0}^{\infty} e^{-s t} f(t) d t\tag{1}$

式中， $s$ 的值需使以上积分收敛。
拉普拉斯变换是一个线性变换(易证)，即 $\mathcal{L}\left\{c_{1} f_{1}(t)+c_{2} f_{2}(t)\right\}=c_{1} \mathcal{L}\left\{f_{1}(t)\right\}+c_{2} \mathcal{L}\left\{f_{2}(t)\right\}\tag{2}$

1.1 例子

$\begin{aligned} \mathcal{L}\left\{e^{a t}\right\} &=\int_{0}^{\infty} e^{-(s-a) t} d t \\ &\left.=-\frac{1}{s-a} e^{-(s-a) t}\right]_{0}^{\infty} \\ &=\frac{1}{s-a} \end{aligned}\tag{3}$

注：要使式(3)的积分收敛，需有 $s > a$

2.常系数微分方程的拉式变换

对于如下非齐次常系数二阶ode
$\ddot{x}+b \dot{x}+c x=g(t), \quad x(0)=x_{0}, \quad \dot{x}(0)=u_{0}\tag{4}$
两边取拉式变换，并考虑拉式变换是一种线性变换，可得
$\mathcal{L}\{\ddot{x}\}+b \mathcal{L}\{\dot{x}\}+c \mathcal{L}\{x\}=\mathcal{L}\{g\}\tag{5}$

函数导数的拉式变换可由分步积分求得：
$\int_{0}^{\infty} e^{-s t} \dot{x} d t=\left.x e^{-s t}\right|_{0} ^{\infty}+s \int_{0}^{\infty} e^{-s t} x d t=s X(s)-x_{0}\tag{6}$

$\int_{0}^{\infty} e^{-s t} \ddot{x} d t=\left.\dot{x} e^{-s t}\right|_{0} ^{\infty}+s \int_{0}^{\infty} e^{-s t} \dot{x} d t=-u_{0}+s\left(s X(s)-x_{0}\right)=s^{2} X(s)-s x_{0}-u_{0}\tag{7}$

注:要使(6)(7)能够收敛，需要 $s > 0$ 。
则式(5)为
$a\left(s^{2} X-s x_{0}-u_{0}\right)+b\left(s X-x_{0}\right)+c X=G\tag{8}$

这是一个关于 $X = X (s)$ 的代数方程，容易求解。之后，对 $X = X (s)$ 进行拉式反变换来求 $x = x (t)$ 。
注：求解示意图

3.单位阶跃函数

单位阶跃函数 $u_{c}(t)$ 定义为
$u_{c}(t)=\left\{\begin{array}{ll} 0, & t<c \\ 1, & t \geq c \end{array}\right.\tag{9}$

要使上式有实际物理意义，需有 $c > 0$
阶跃函数的拉式变换为
$\mathcal{L}\left\{u_{c}(t)\right\}=\int_{0}^{\infty} e^{-s t} u_{c}(t) d t=\frac{e^{-c s}}{s}\tag{9}$

要使上式积分收敛，需要 $s > 0$ 。

3.1 表示滞后函数

阶跃函数可以被用于表示函数 $f (t)$ 向 $t$ 的正方向平移 $c$ 后的函数(该函数和 $f (t)$ 相比滞后了 $c$ ):
$u_{c}(t) f(t-c)=\left\{\begin{aligned} 0, & t<c \\ f(t-c), & t \geq c \end{aligned}\right.\tag{10}$