首页 > 其他分享 >拉普拉斯变换10.24

拉普拉斯变换10.24

时间:2024-10-24 17:21:11浏览次数:8  
标签:infty right frac 变换 拉普拉斯 10.24 left

目录

1.拉普拉斯变换

\[\mathscr{L} [f(t)]=F(s)= \int ^ \infty _0 f(t) e^{-st}dt \]

$s=\sigma + j \omega $​​ 是一个复数

简单例子 ,证明 \(\mathscr{L} [e^{-at}]=\frac {1}{s+a}\)

\[\mathscr{L} [e^{-at}] =\int ^ \infty _0 e^{-at} e^{-st} dt =\int ^ \infty _0 e^{-(a+s)t} dt \\ =-\frac{1}{s+a} e^{-(a+s)t} \Big|^ \infty _0 = \frac{-a}{s+a} \\ = - \frac{1}{s+a} \lim_{b \to \infty} e^{-(s+a)t} -(-\frac{1}{s+a}e^0 \\ =\frac {1}{s+a} \]

2. 拉普拉斯收敛域

把$s=\sigma + j \omega $ 代入 \(F(s)=\int ^ \infty _0 e^{-at} e^{-st} dt\)

\[F(s)=\int ^ \infty _0 e^{-at} e^{-st} dt =\int ^ \infty _0 e^{-at} e^{-(\sigma + j \omega)t} dt =\int ^ \infty _0 e^{-(a+\sigma)t} e^{- j \omega t} dt \]

上式积分项中的第二项 $ e^{-j \omega t } $是旋转矢量,可以根据欧拉公式展开,是周期函数,不影响收敛。

积分项中的第一项的指数 $-(a+\sigma)<0 $ ,才能保证积分收敛。所以, $ \sigma >-a$ 是上述变换的收敛域。

3.导数的拉普拉斯变换

函数 $ f(t) $ 的导数的拉普拉斯变换是拉普拉斯变换理论中的一个重要性质。如果 $ f(t) $ 及其导数 $ f'(t) $ 都是分段光滑且指数阶的,那么 $ f'(t) $ 的拉普拉斯变换可以表示为:

\[\mathscr{L}\{f'(t)\} = sF(s) - f(0^+) \]

其中,$ F(s) $ 是 $ f(t) $ 的拉普拉斯变换,$ f(0^+) $ 是 $ f(t) $ 在 $ t = 0 $ 处的值。

这个性质说明,函数导数的拉普拉斯变换等于 $ s $ 乘以该函数的拉普拉斯变换减去该函数在 $ t = 0 $ 时的值。这个公式在求解微分方程的拉普拉斯变换时非常有用。

推导过程

假设 $ f(t) $ 是一个分段光滑且指数阶的函数,那么它的拉普拉斯变换定义为:

\[F(s) = \mathscr{L}\{f(t)\} = \int_0^\infty e^{-st} f(t) \, dt \]

现在,我们要求 $ f'(t) $ 的拉普拉斯变换:

\[\mathscr{L}\{f'(t)\} = \int_0^\infty e^{-st} f'(t) \, dt \]

使用分部积分法,令 $ u = e^{-st} $ 和 $ dv = f'(t) dt $,则 $ du = -se^{-st} dt $ 和 $ v = f(t) $。应用分部积分:

\[\int_0^\infty e^{-st} f'(t) \, dt = \left[ e^{-st} f(t) \right]_0^\infty - \int_0^\infty (-s e^{-st}) f(t) \, dt \]

\[= \left[ e^{-st} f(t) \right]_0^\infty + s \int_0^\infty e^{-st} f(t) \, dt \]

由于 \(f(t)\) 是指数阶的,当 $ t \to \infty $ 时,$ e^{-st} f(t) \to 0 $。因此,第一项计算结果为 $ f(0^+) $:

\[= f(0^+) + sF(s) \]

所以,我们得到:

\[\mathscr{L}\{f'(t)\} = sF(s) - f(0^+) \]

这就是函数导数的拉普拉斯变换的公式。这个公式可以推广到更高阶的导数。对于 $ n $ 阶导数,拉普拉斯变换是:

\[\mathscr{L}\{f^{(n)}(t)\} = s^n F(s) - s^{n-1} f(0^+) - s^{n-2} f'(0^+) - \cdots - f^{(n-1)}(0^+) \]

这个公式在求解微分方程和控制系统分析中非常有用。

5.传递函数

系统输出的L变换与系统输入的L变换之比,定义传递函数

\[G(s)=\frac{X(s)}{U(s)} \]

6.电感电阻电路

动态方程

电流的动态微分方程为

\[e\left(t\right)=L\frac{\mathrm{d}i\left(t\right)}{\mathrm{d}t}+Ri\left(t\right) \]

image

定义 此动态 系统的 输入为 电压其 \(u(t)=e(t)\), 输出为 电流\(x(t)=i(t)\), 则

\[u(t)=L \frac{d x(t)}{dt}+R(t) \]

拉氏变换

对两边进行L变换

\[\mathcal{L}[u\left(t\right)]=\mathcal{L}\left[L \frac{\mathrm{d}x\left(t\right)}{\mathrm{d}t}+Rx\left(t\right)\right]\\\Rightarrow U\left(s\right)=L\left(sX\left(s\right)-x\left(0\right)\right)+RX\left(s\right) \]

考虑零初始状态,\(x(0)=0\) ,则

\[U(s)=LsX(s)+RX(s)=(Ls+R)X(s) \]

其传递函数为

\[G(s)=\frac{X(s)}{U(s)}=\frac{1}{Ls+R} \]

常数输入

\(u(t)=C\)作用系统,L变换为

\[U(s)=\mathcal{L}[u(t)]=\mathcal{L}[C]=\mathcal{L}[C\mathrm{e}^{-0t}]=C \frac{1}{s+0}=C \frac{1}{s} \]

代入传函 \(G(s)=\frac{1}{Ls+R}\)

\[X(s)=U(s)G(s)=C\left(\frac{1}{s}\right)\left(\frac{1}{Ls+R}\right)=\frac{C}{L} \frac{1}{s\left(s+\frac{R}{L}\right)} \]

分式分解,设待定系数A,B

\[X(s)=\frac{C}{L}\frac{1}{s\left(s+\frac{R}{L}\right)}=\frac{C}{L}\left(\frac{A}{s}+\frac{B}{s+\frac{R}{L}}\right) \]

解得\(A=\frac LR,B=-\frac LR\),则

\[X\left(s\right)=\frac{C}{L}\left(\frac{L}{R}-\frac{L}{R}\right)=\frac{C}{R}\left(\frac{1}{s}-\frac{1}{s+\frac{R}{L}}\right) \]

L逆变换

对其进行L逆变换

\[x\left(t\right)=\mathcal{L}^{-1}\left[X\left(s\right)\right]=\frac{C}{R}\left(\mathrm{e}^{0t}-\mathrm{e}^{-\frac{R}{L^{2}}}\right)=\frac{C}{R}-\frac{C}{R}\mathrm{e}^{-\frac{R}{L^{2}}t} \]

image

7. 控制系统传递函数

如图所示的开环控制系统(Open Loop Control System):其中\(R(s)\)是参考值(Reference)或目标值,\(C(s)\)是控制器,原动态系统的传递函数\(G\)(s)被称为控制系统的开环传递函数(Open Loop Transfer Function)。控制量是\(U(s)\),也就是原动态系统的输入。控制系统的输出等于原动态系统的输出\(X(s)\)​。

image

控制系统本质上也是一个动态系统,从参考值\(R(s)\)(它同时也是该控制系统的输入,又称参考输人)到系统输出\(X(s)\)是串联的结构,即

\[X(s)=U(s)G(s)=R(s)C(s)G(s) \]

若将输出\(X(s)\)反馈到输人端,则可以形成一个闭环控制系统,如图 所示。其中, 参考值与输出之间的差称为误差(Error),\(E\left(s\right)=R\left(s\right)-X\left(s\right)\),其对应的时间函数是\(e\left(t\right)=\) \(r(t)-x\left(t\right)\),控制器 \(C(s)\)将根据误差决定控制量\(U(s)\)。

image

根据传递兩数的代数性质,可得

\[X(s)=U(s)G(s)=E(s)C(s)G(s) \]

将\(E(s)=R(s)-X(s)\) 代入上式

\[\begin{aligned}&X\left(s\right)=\left(R\left(s\right)-X\left(s\right)\right)C\left(s\right)G\left(s\right)\\&\Rightarrow(1+C(s)G(s))X(s)=C(s)G(s)R(s)\\&\Rightarrow X(s)=\frac{C(s)G(s)R(s)}{1+C(s)G(s)}\end{aligned} \]

定义闭环传递函数

\[G_{\mathrm{cl}}(s)=\frac{X(s)}{R(s)}=\frac{C(s)G(s)}{1+C(s)G(s)} \]

得到简化的闭环控制系统框图

8.非零初始状态的传递函数

在 传递函 数的定义中宥 一个先 决条件 ,即零初始 条件。 但在实际情 况中, 往往需要处理非零初始 状态的 系统。

考虑一阶微分方程

\[\frac{\mathrm{d}x(t)}{\mathrm{d}t}+ax(t)=u(t) \]

L变换

\[\mathcal{L}\left[\frac{\mathrm{d}x\left(t\right)}{\mathrm{d}t}+ax\left(t\right)\right]=\mathcal{L}[u\left(t\right)]\\sX\left(s\right)-x\left(0\right)+aX\left(s\right)=U\left(s\right) \]

在零初始条件下\((x(0)=0)\), 可写成 s\(X(s)+aX(s)=U(s)\),系统的传递函数是

\[G(s)=\frac{X(s)}{U(s)}=\frac{1}{s+a} \]

此时定义新的系统输入\(:U_1(s)=U(s)+x(0)\),代人式中,得到

\[G(s)=\frac{X(s)}{U_{1}(s)}=\frac{1}{s+a} \]

image

这两个系统的传递函数是相同的,其中非零初始条件系统多出一个输人,而这个输人的拉普拉斯变换等于其初始条件\(x(0)\)。对它进行拉普拉斯逆变换可以得到其原函数,即

\[\mathcal{L}^{-1}[x(0)]=x(0)\delta(t) \]

其中,\(\delta(t)\)是单位冲激函数,可以把它理解为在很短的时间内释放出的一个单位的能量。将它乘以一个系数 \(x\)(0),则相当于在一瞬间对系统施加了\(x\)(0)个单位的能量(系统的输出也将叠加\(x(0)h(t))\)​。因为这个能量是瞬间的,并不持续,所以它不会影响到系统的稳定性分析与特征分析。高阶系统的非零初始条件的分析则比较复杂,但是其理念与一阶系统相同,系统的初始状态可以理解为瞬时间赋予系统的“能量”。

标签:infty,right,frac,变换,拉普拉斯,10.24,left
From: https://www.cnblogs.com/redufa/p/18500016

相关文章

  • 24.10.24
    A大家使用了整体二分+可撤销并查集,倍增等方法...考虑线段树合并。在跑Kruskal时,如果一个询问的两个点在同一个连通块内,那么这个询问就是可回答的,但是可回答不一定要回答,因为如果此后加的边权相同那么其实里面的点还能再往外走。所以在加边时如果新加的边权大于连通块边权,那......
  • 10.24Python_pandas_基础
    一、基础1、概述Pandas是一个开源的第三方Python库,从Numpy和Matplotlib的基础上构建而来Pandas名字衍生自术语“paneldata”(面板数据)和“Pythondataanalysis”(Python数据分析)Pandas已经成为Python数据分析的必备高级工具,它的目标是成为强大、灵活、可以......
  • 10.24
    考前挂分是个好迹象,至少不像啥也不会那么绝望是不是/A.城市间交通第一眼整体二分+可撤销并查集,觉得有点难写,而且两个\(\log\)。再看一眼,发现最小生成树+倍增优秀单\(\log\)做法。B.最小公倍数第一眼这不是我们P3911最小公倍数之和吗?坏消息是忘了怎么莫反了。于是写了......
  • 2024.10.24 1234版
    起于《海奥华预言》的思考◆地球管理结构和参考持续更新中...... 英文地址:https://github.com/zhuyongzhe/Earth/tags中文地址:https://www.cnblogs.com/zhuyongzhe85作者:朱永哲 ---------------------------------------------------------------------------------......
  • 10.24
    实验3:工厂方法模式本次实验属于模仿型实验,通过本次实验学生将掌握以下内容:1、理解工厂方法模式的动机,掌握该模式的结构;2、能够利用工厂方法模式解决实际问题。[实验任务一]:加密算法目前常用的加密算法有DES(DataEncryptionStandard)和IDEA(InternationalDataEncryptionAlgo......
  • matlab求冲击输入拉普拉斯变换,和响应
    在MATLAB中,求解冲击输入的拉普拉斯变换以及系统的冲击响应通常涉及以下几个步骤:定义冲击输入:在MATLAB中,冲击输入通常用狄拉克δ函数(Diracdeltafunction)表示,可以使用dirac函数来创建。计算拉普拉斯变换:使用laplace函数计算冲击输入的拉普拉斯变换。计算系统冲击响应:如果......
  • LTE 基于快速哈达玛hadamard变换SSS辅同步信号SSS检测之hadamard变换公式推导
    LTESSSs序列生成的阶数为31阶,所以hadamard矩阵的阶数为32阶,定义一个32阶的hadamard矩阵H32,下面进行hadamard快速变换公式推导继续分解后面公式的推导小编还在继续......
  • 10.24
    一.填空题(共3题,60分)(填空题)支持向量到超平面的距离之和称之为。我的答案:(1)间隔(填空题)支持向量机的核心思想是(最大化/最小化)间隔。我的答案:(1)最大化(填空题)函数可以作为核函数。我的答案:(1)满足Mercer定理条件......
  • 2022.10.24
    练习情况P8593「KDOI-02」一个弹的投题目拆分为两个问题,一个是求每个炸弹的威力,另一个是求最多减少多少威力。根据物理知识可知,当且仅当\(y_i=y_j\)时,这两枚导弹才有可能相遇。将落地点离散化。使用权值树状数组求逆序对。Code:P8593CF1311FMovingPoints将速度离散......
  • 对图像进行傅立叶变换的目的是什么
    傅立叶变换是一种在图像处理和信号处理中非常重要的技术,它能够将信号从时域(时间域)转换到频域(频率域),从而揭示出信号的频谱信息。对图像进行傅立叶变换的主要目的包括以下几个方面: 1.分析图像的频率成分:通过傅立叶变换,可以将图像从空间域转换到频率域,得到图像的频谱图。这有......