最近复习到自动控制理论,对各种信号和其拉氏变换还不太熟悉,推导一下加深印象。

 在描述一个系统或环节运动规律时,微分方程是其数学模型的最基本形式,然而对于一个二阶以上的系统,其微分方程难以得到解,为了便于研究,引入拉普拉斯变换将时域内的微分方程变成S域中的传递函数进行求解。

1.拉普拉斯变换

 首先给出拉普拉斯变换定义式L[r(t)]=∫0∞r(t)e−stdt=R(s)L{[r(t)}]=\int_{0}^{\infty}r(t)e^{-st}dt=R(s)L[r(t)]=0r(t)estdt=R(s)
其中r(t)r(t)r(t)是关于时间的函数,可以是输入信号。

性质

拉氏变换相对于微分方程的优越性在于其对微分方程中导数和积分量的变换

(1).微分性质:

∫0∞r′(t)e−stdt=∫0∞e−std[r(t)]=e−str(t)∣0∞+s∫0∞r(t)e−stdt=sR(s)−r(0)\int_{0}^{\infty}r^{'}(t)e^{-st}dt=\int_{0}^{\infty}e^{-st}d[r(t)]=e^{-st}r(t)|_{0}^{\infty}+s\int_{0}^{\infty}r(t)e^{-st}dt=sR(s)-r(0)0r(t)estdt=0estd[r(t)]=estr(t)0+s0r(t)estdt=sR(s)r(0)
由于在自动控制原理中传递函数均为0初始条件,因此最终得到L[r′(t)]=sR(s)L[r^{'}(t)]=sR(s)L[r(t)]=sR(s)

(2).积分性质:

∫0∞(∫r(t))e−stdt=−1s∫0∞(∫r(t))d(e−st)=−1sr(t)e−st∣0∞+∫0∞r(t)e−stdt=R(s)s+r(0)s\int_{0}^{\infty}(\int{r(t)})e^{-st}dt=-\frac{1}{s}\int_{0}^{\infty}(\int{r(t)})d(e^{-st})=-\frac{1}{s}r(t)e^{-st}|_{0}^{\infty}+\int_{0}^{\infty}r(t)e^{-st}dt=\frac{R(s)}{s}+\frac{r(0)}{s}0(r(t))estdt=s10(r(t))d(est)=s1r(t)est0+0r(t)estdt=sR(s)+sr(0)
由于在自动控制原理中传递函数均为0初始条件,因此最终得到L[∫r(t)]=R(s)sL[\int{r(t)}]=\frac{R(s)}{s}L[r(t)]=sR(s)
  有了上述两个定理后,就可以对微分方程进行零初始条件下的拉氏变换,进行S域上的研究。

(3).位移性质:

L[r(t−τ)]=∫0∞r(t−τ)e−stdt=∫τ∞r(t−τ)e−stdtL[r(t-\tau)]=\int_{0}^{\infty}r(t-\tau)e^{-st}dt=\int_{\tau}^{\infty}r(t-\tau)e^{-st}dtL[r(tτ)]=0r(tτ)estdt=τr(tτ)estdt
L[r(t)]=∫0∞r(t)e−stdt=R(s)L[r(t)]=\int_{0}^{\infty}r(t)e^{-st}dt=R(s)L[r(t)]=0r(t)estdt=R(s)
∫τ∞r(t−τ)e−stdt=e−sτ∫τ∞r(t−τ)e−s(t−τ)d(t−τ)=e−sτ∫0∞r(u)e−sud(u)=e−sτR(s)\int_{\tau}^{\infty}r(t-\tau)e^{-st}dt=e^{-s\tau}\int_{\tau}^{\infty}r(t-\tau)e^{-s(t-\tau)}d(t-\tau)=e^{-s\tau}\int_{0}^{\infty}r(u)e^{-su}d(u)=e^{-s\tau}R(s)τr(tτ)estdt=esττr(tτ)es(tτ)d(tτ)=esτ0r(u)esud(u)=esτR(s)
注意性质的正反使用。

2.常用输入信号

首先进行一个归纳:

输入信号 拉氏变换
单位脉冲信号δ(t)\delta(t)δ(t) 1
单位阶跃信号1(t)1(t)1(t) 1s\frac{1}{s}s1
单位速度信号 1s2\frac{1}{s^2}s21
单位加速度信号 1s3\frac{1}{s^3}s31
正弦信号 ws2+w2\frac{w}{s^2+w^2}s2+w2w
(1)单位脉冲信号δ(t)\delta(t)δ(t)

δ(t)={∞t=00t≠0\delta(t)=\begin{cases} \infty & t=0 \\ 0 & t\neq0 \\ \end{cases}δ(t)={0t=0t=0
为了便于积分,也可写作下式:
δ(t)={lim⁡ε→01εt=00t≠0\delta(t)=\begin{cases} \lim\limits_{\varepsilon\rightarrow0}\frac{1}{\varepsilon} & t=0 \\ 0 & t\neq0 \\ \end{cases}δ(t)={ε0limε10t=0t=0
对其进行拉氏变换:
L[δ(t)]=∫0∞δ(t)e−stdt=lim⁡ε→0∫0ε1εe−stdt=lim⁡ε→0∫0ε−1εsd(e−st)=lim⁡ε→0−1εse−st∣0ε=lim⁡ε→01εs(1−e−εs)L[\delta(t)]=\int_{0}^{\infty}\delta(t)e^{-st}dt= \lim\limits_{\varepsilon\rightarrow0}\int_{0}^{\varepsilon}\frac{1}{\varepsilon}e^{-st}dt=\lim\limits_{\varepsilon\rightarrow0}\int_{0}^{\varepsilon}-\frac{1}{\varepsilon{}s}d(e^{-st})=\lim\limits_{\varepsilon\rightarrow0}{-\frac{1}{\varepsilon{s}}e^{-st}|_{0}^{\varepsilon}}=\lim\limits_{\varepsilon\rightarrow0}{\frac{1}{\varepsilon{s}}(1-e^{-\varepsilon{s}})}L[δ(t)]=0δ(t)estdt=ε0lim0εε1estdt=ε0lim0εεs1d(est)=ε0limεs1est0ε=ε0limεs1(1eεs)
e−εse^{-\varepsilon{s}}eεs一阶展开展开得e−εs=1−εs+o(εs)e^{-\varepsilon{s}}=1-{\varepsilon{s}}+o({\varepsilon{s}})eεs=1εs+o(εs)
代入上式得L[δ(t)]=1L[\delta(t)]=1L[δ(t)]=1

(2)单位阶跃信号1(t)1(t)1(t)

1(t)={1t>00t≤01(t)=\begin{cases} 1 & t>0 \\ 0 & t\le0 \end{cases}1(t)={10t>0t0
进行拉氏变换:
L[1(t)]=∫0∞e−stdt=−1se−st∣0∞=1sL[1(t)]=\int_{0}^{\infty}e^{-st}dt=-{\frac{1}{s}}e^{-st}|_{0}^{\infty}=\frac{1}{s}L[1(t)]=0estdt=s1est0=s1

(3)单位斜坡(速度)信号

r(t)=t        t≥0r(t)=t~~~~~~~~t\ge0r(t)=t        t0
进行拉氏变换:
L[r(t)]=∫0∞te−stdt=−∫0∞1std(e−st)=−1ste−st∣0∞+1s∫0∞e−stdt=1s2L[r(t)]=\int_{0}^{\infty}te^{-st}dt=-\int_{0}^{\infty}\frac{1}{s}td(e^{-st})=-\frac{1}{s}te^{-st}|_{0}^{\infty}+\frac{1}{s}\int_{0}^{\infty}e^{-st}dt=\frac{1}{s^2}L[r(t)]=0testdt=0s1td(est)=s1test0+s10estdt=s21

(4)单位加速度信号

r(t)=t22        t≥0r(t)=\frac{t^2}{2}~~~~~~~~t\ge0r(t)=2t2        t0
进行拉氏变换:
L[r(t)]=12∫0∞t2e−stdt=−12∫0∞1st2d(e−st)=−12st2e−st∣0∞+1s∫0∞te−stdt=1s3L[r(t)]=\frac{1}{2}\int_{0}^{\infty}t^2e^{-st}dt=-\frac{1}{2}\int_{0}^{\infty}\frac{1}{s}t^2d(e^{-st})=-\frac{1}{2s}t^2e^{-st}|_{0}^{\infty}+\frac{1}{s}\int_{0}^{\infty}te^{-st}dt=\frac{1}{s^3}L[r(t)]=210t2estdt=210s1t2d(est)=2s1t2est0+s10testdt=s31

(5)正弦信号

r(t)=Asinwtr(t)=Asinwtr(t)=Asinwt
根据欧拉公式ejwt=coswt+jsinwte^{jwt}=coswt+jsinwtejwt=coswt+jsinwt
得出sinwt=ejwt−e−jwt2jsinwt=\frac{e^{jwt}-e^{-jwt}}{2j}sinwt=2jejwtejwt
进行拉氏变换:
L[r(t)]=12j∫0∞(ejwt−e−jwt)e−stdt=12j∫0∞e(jw−s)t−e−(jw+s)tdtL[r(t)]=\frac{1}{2j}\int_{0}^{\infty}(e^{jwt}-e^{-jwt})e^{-st}dt=\frac{1}{2j}\int_{0}^{\infty}e^{(jw-s)t}-e^{-(jw+s)t}dtL[r(t)]=2j10(ejwtejwt)estdt=2j10e(jws)te(jw+s)tdt
=12j[1jw−se(jw−s)t∣0∞+1jw+se−(jw+s)t∣0∞]=12j(1s−jw−1s+jw)=ws2+w2=\frac{1}{2j}[\frac{1}{jw-s}e^{(jw-s)t}|_{0}^{\infty}+\frac{1}{jw+s}e^{-(jw+s)t}|_{0}^{\infty}]=\frac{1}{2j}(\frac{1}{s-jw}-\frac{1}{s+jw})=\frac{w}{s^2+w^2}=2j1[jws1e(jws)t0+jw+s1e(jw+s)t0]=2j1(sjw1s+jw1)=s2+w2w

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐