第二十讲 拉普拉斯变换求解线性ODE

一,拉氏变换公式:

  • \mathcal {L}[f(t)]=\int_{0}^{\infty }f(t)e^{-st}dts> 0

二,保证f(t)的拉氏变换存在的条件:

  • 只有在\int_{0}^{\infty }f(t)e^{-st}dt收敛的情况下,拉氏变换才存在
  • 增长条件:保证拉氏变换存在的前提是f(t)的增长速度不能超过e^{-st}的收敛速度,\lim_{t\rightarrow \infty }\frac{f(t)}{e^{st}}=0
  • 指数型表达法:\left |f(t) \right |\leq Ce^{kt}C> 0k\geq 0t\geq 0
  • 正例:
  1. \left |sin(t) \right |\leq 1\cdot e^{0t}
  2. t^{n}\leq M\cdot e^{t}
  • 反例:
  1. \frac{1}{t},因为\int_{0}^{\infty }\frac{1}{t}e^{-st}dt\approx \int_{0}^{\infty }\frac{1}{t}dt=\left. ln(t) \right |^{\infty }_{0}\rightarrow \infty
  2. e^{t^{2}},因为e^{t^{2}}> e^{kt}t^{2}> kt\lim_{t\rightarrow \infty }\frac{t^{2}}{kt}\rightarrow \infty

三,导数的拉氏变换:

  • 解微分方程:{y}''+A{y}'+By=h(t)
  • 设初始条件:y(0)=y_{0}{y}'(0)={y}'_{0}
  • 一阶导数的拉氏变换:
  1. \mathcal {L}[{y}'(t)]=\int_{0}^{\infty }{y}'(t)e^{-st}dt=\left. y(t)e^{-st} \right |^{\infty }_{0}+s\int_{0}^{\infty }y(t)e^{-st}dt(分部积分)
  2. 要使\lim_{t\rightarrow \infty }\frac{y(t)}{e^{st}}\rightarrow 0,必须满足增长条件:\left |y(t) \right |\leq Ce^{kt}C> 0k\geq 0t\geq 0,或者说s> k
  3. \left. y(t)e^{-st} \right |^{\infty }_{0}=0-y(0)
  4. s\int_{0}^{\infty }y(t)e^{-st}dt=s\cdot \mathcal {L}[y(t)]
  5. \mathcal {L}[{y}'(t)]=s\cdot \mathcal {L}[y(t)]-y(0)
  6. 代入初始条件:y(0)=y_{0}
  • 二阶导数的拉氏变换:
  1. \mathcal {L}[{y}''(t)]=s\cdot \mathcal {L}[{y}'(t)]-{y}'(0)=s[s\cdot \mathcal {L}[y(t)]-y(0)]-{y}'(0)
  2. \mathcal {L}[{y}''(t)]=s^{2}\cdot \mathcal {L}[y(t)]-s\cdot y(0)-{y}'(0)
  3. 代入初始条件:y(0)=y_{0}{y}'(0)={y}'_{0}

四,例题:

  • 求解:{y}''-y=e^{-t}
  • 已知初始条件:y(0)=1{y}'(0)=0
  • 用微分方程法解题步骤:
  1. 解齐次方程{y}''-y=0
  2. 求特征方程r^{2}-r=0
  3. 根据特征根求通解
  4. 解非齐次方程{y}''-y=e^{-t}
  5. 化成特殊形式,根据指数输入定理求特解
  6. 通解=齐次方程通解+非齐次方程特解
  7. 代入初始条件得特解

五,用拉氏变换解如上例题的步骤:

  • 第一步:等式两边拉氏变换:
  1. 左边:\mathcal {L}[{y}''-y]=\mathcal {L}[{y}'']-\mathcal {L}[y]=s^{2}\cdot \mathcal {L}[y]-s\cdot y(0)-{y}'(0)-\mathcal {L}[y]
  2. 代入初始条件:y(0)=1{y}'(0)=0
  3. 左边:\mathcal {L}[{y}''-y]=s^{2}\cdot \mathcal {L}[y]-s-\mathcal {L}[y]
  4. 右边:\mathcal {L}[e^{-t}]=\frac{1}{s+1}(根据公式,查表)
  5. 左边=右边:s^{2}\cdot \mathcal {L}[y]-s-\mathcal {L}[y]=\frac{1}{s+1}
  • 第二步:解出\mathcal {L}[y]
  1. (s^{2}-1)\mathcal {L}[y]=\frac{1}{s+1}+s=\frac{s^{2}+s+1}{s+1}
  2. \mathcal {L}[y]=\frac{s^{2}+s+1}{(s+1)^{2}(s-1)}
  • 第三步:拉式逆变换\mathcal {L}^{-1}[\mathcal {L}[y]]=y(t)
  1. 部分分式:\frac{s^{2}+s+1}{(s+1)^{2}(s-1)}=\frac{A}{(s+1)^{2}}+\frac{B}{(s+1)}+\frac{C}{(s-1)}
  2. 等式两边同乘(s+1)^{2},令s=-1,解得A=-\frac{1}{2}
  3. 等式两边同乘(s-1),令s=1,解得C=\frac{3}{4}
  4. 代入法:代入s=0,得-1=-\frac{1}{2}+B-\frac{3}{4},解得B=\frac{1}{4}
  5. \mathcal {L}^{-1}[-\frac{1}{2}\cdot \frac{1}{(s+1)^{2}}]=-\frac{1}{2}\cdot t\cdot e^{-t}(查表)
  6. \mathcal {L}^{-1}[\frac{1}{4}\cdot \frac{1}{(s+1)}]=\frac{1}{4}\cdot e^{-t}(查表)
  7. \mathcal {L}^{-1}[\frac{3}{4}\cdot \frac{1}{(s-1)}]=\frac{3}{4}\cdot e^{t}(查表)
  8. 结果:y(t)=-\frac{1}{2}\cdot t\cdot e^{-t}+\frac{1}{4}\cdot e^{-t}+\frac{3}{4}\cdot e^{t}

你可能感兴趣的:(微分方程)