Loading [MathJax]/jax/output/HTML-CSS/jax.js

信号与系统 – 使用拉普拉斯变换求解微分方程


拉普拉斯变换

拉普拉斯变换是一种数学工具,用于将时域中的微分方程转换为频域或s域中的代数方程。

数学上,如果x(t)是时域函数,则其拉普拉斯变换定义为:

L[x(t)]=X(s)=x(t)estdt

使用拉普拉斯变换求解微分方程

线性时不变 (LTI) 系统由常系数微分方程描述,这些方程将系统的输入和输出联系起来。LTI 系统的响应是通过求解这些微分方程得到的。

拉普拉斯变换技术可用于求解描述 LTI 系统的微分方程。使用拉普拉斯变换,时域中的微分方程被转换为s域中的代数方程。通过操作代数方程,在s域中获得解。然后通过对响应进行拉普拉斯反变换,将其转换回时域。

Explore our latest online courses and learn new skills at your own pace. Enroll and become a certified expert to boost your career.

解释

考虑一个输入为x(t),输出为y(t)的 LTI 系统,并由以下微分方程描述:

andndtny(t)+an1dn1dtn1y(t)+...+a0y(t)=bmdmdtmx(t)+bm1dm1dtm1x(t)+...+b0x(t)...(1)

当系统的初始条件为零时,即

y(0)=dy(0)dt=d2y(0)dt2...=d(n1)y(0)dt(n1)=0

以及

x(0)=dx(0)dt=d2x(0)dt2...=d(m1)x(0)dt(m1)=0

现在对式 (1) 中给出的微分方程的两边应用拉普拉斯变换,并忽略初始条件,我们得到:

L[andndtny(t)]+L[an1dn1dtn1y(t)]+...L[a0y(t)]=L[bmdmdtmx(t)]+L[bm1dm1dtm1x(t)]+...L[b0x(t)]

ansnY(s)+an1sn1Y(s)+...+a0Y(s)=bmsmX(s)+bm1sm1X(s)+...+b0X(s)

(ansn+an1sn1+...+a0)Y(s)=(bmsm+bm1sm1+...+b0)X(s)

Y(s)X(s)=(bmsm+bm1sm1+...+b0)(ansn+an1sn1+...+a0)...(2)

方程 (2) 称为 LTI 系统的**传递函数**。因此,LTI 系统的传递函数定义为输出 [Y(s)] 的拉普拉斯变换与输入 [X(s)] 的拉普拉斯变换之比,且初始条件为零。

LTI 系统的传递函数也可以定义为系统脉冲响应的拉普拉斯变换,即

H(s)=Y(s)X(s)...(3)

函数 H(s) 仅取决于描述系统的微分方程的系数,它不取决于输入信号或系统存储的初始能量。

由于系统在s域中的响应由下式给出:

Response,Y(s)=H(s)X(s)...(4)

对式 (4) 进行拉普拉斯反变换,我们得到:

y(t)=L1[H(s)X(s)]...(5)

方程 (5) 称为系统的**零状态响应**。

现在,如果脉冲函数是系统的输入,即 x(t)=δ(t),则 X(s)=1。因此,式 (5) 可以写成

y(t)=L1[H(s)]=h(t)...(6)

方程 (6) 称为系统的**脉冲响应**。

数值示例

求由以下微分方程描述的系统的传递函数:

d2dt2y(t)+5ddty(t)+10y(t)=3ddtx(t)+8x(t)

解答

给定的微分方程为:

d2dt2y(t)+5ddty(t)+10y(t)=3ddtx(t)+8x(t)

忽略初始条件,并在方程两侧进行拉普拉斯反变换,我们得到:

L[d2dt2y(t)+5ddty(t)+10y(t)]=L[3ddtx(t)+8x(t)]

s2Y(s)+5sY(s)+10Y(s)=3sX(s)+8X(s)

(s2+5s+10)Y(s)=(3s+8)X(s)

因此,系统的传递函数为:

H(s)=Y(s)X(s)=(3s+8s2+5s+10)

更新于: 2022年1月11日

2K+ 阅读量

开启你的 职业生涯

通过完成课程获得认证

开始学习
广告