美文网首页
初边值问题的分离变量法

初边值问题的分离变量法

作者: 洛玖言 | 来源:发表于2020-05-15 23:10 被阅读0次

    初边值问题的分离变量法

    考察波动方程的初边值问题
    \begin{aligned} \begin{cases} u_{tt}-a^2u_{xx}=0\cdots①\\ t=0:\;u=\varphi(x),u_{t}=\psi(x)\\ u(0,t)=u(l,t)=0 \end{cases} \end{aligned}
    先求方程①的可以分离变量的非平凡(即不恒等于零)的特解:
    u(x,t)=X(x)T(t)
    带入方程①可得
    X(x)T''(t)-a^2X''(x)T(t)=0
    上式分离变量得
    \dfrac{T''(t)}{a^2T(t)}=\dfrac{X''(x)}{X(x)}
    左边仅是 t 得函数,右边仅是 x 的函数,需要等于同一个常数,记为 -\lambda,于是
    T''(t)+\lambda a^2T(t)=0
    X''(x)+\lambda X(x)=0\cdots②

    又有边界条件知 X(0)=X(l)=0
    求非平凡解:

    • \lambda<0
      ② 的通解可以写成
      X(x)=C_1e^{\sqrt{-\lambda}x}+C_2e^{-\sqrt{-\lambda}x}
      X(0)=C_1+C_2=0
      X(l)=C_1e^{\sqrt{-\lambda}l}+C_2e^{-\sqrt{-\lambda}l}=0
      由于 \begin{vmatrix} 1&1\\ e^{\sqrt{-\lambda}l}&e^{-\sqrt{-\lambda}l} \end{vmatrix}\not=0
      \Rightarrow C_1=C_2=0\Rightarrow\lambda<0 时得不到非平凡解.
    • \lambda=0
      ② 的通解可以写成
      X(x)=C_1+C_2x
      可知 X(x) 也恒为零.
    • \lambda>0
      ② 的通解可以写成
      X(x)=C_1\cos\sqrt{\lambda}x+C_2\sin\sqrt{\lambda}x
      由边界条件知
      X(0)=C_1=0
      X(l)=C_2\sin\sqrt{\lambda}l=0\Rightarrow\lambda=\lambda_{k}=\dfrac{k^2\pi^2}{l^2}\;(k=1,2,\cdots)
      于是得到一族非零解
      X_{k}(x)=C_{k}\sin\dfrac{k\pi}{l}x\quad(k=1,2,\cdots)
      由此可得 T_{k} 的通解为
      T_{k}=A_{k}\cos\dfrac{k\pi a}{l}t+B_{k}\sin\dfrac{k\pi a}{l}t\quad(k=1,2,\cdots)
      其中 A_k,B_k 为任意常数
      由此就得到了满足奇次边界下分离变量形式的特解:
      U_{k}(x,t)=X_{k}(x)T_{k}(t)=\left(A_{k}\cos\dfrac{k\pi a}{l}t+B_{k}\sin\dfrac{k\pi a}{l}t\right)\sin\dfrac{k\pi}{l}x
      这边 C_{k} 和其他两个常数结合成了常数.
      由此得到
      \displaystyle u(x,t)=\sum_{k=1}^{\infty}\left(A_{k}\cos\dfrac{k\pi a}{l}t+B_{k}\sin\dfrac{k\pi a}{l}t\right)\sin\dfrac{k\pi}{l}x\cdots③
      由初值问题得
      \displaystyle u_{t}=\sum_{k=1}^{\infty}\dfrac{k\pi a}{l}\left(-A_{k}\sin\dfrac{k\pi a}{l}t+B_{k}\cos\dfrac{k\pi a}{l}t\right)\sin\dfrac{k\pi}{l}x
      由初始条件
      u(x,0)=\varphi(x)\\u_{t}(x,0)=\psi(x)

      \begin{aligned} \begin{cases} \displaystyle\varphi(x)=\sum_{k=1}^{\infty}A_{k}\sin\dfrac{k\pi}{l}x\\ \displaystyle\psi(x)=\sum_{k=1}^{\infty}\dfrac{B_{k}k\pi a}{l}\sin\dfrac{k\pi}{l}x \end{cases} \end{aligned}
      因此,A_{k}\dfrac{B_{k}k\pi a}{l} 应分别是 \varphi(x)\psi(x)[0,l] 区间中正弦展开得傅里叶级数的系数,即
      \begin{aligned} \begin{cases} \displaystyle A_{k}=\dfrac{2}{l}\int_{0}^l\varphi(\xi)\sin\dfrac{k\pi}{l}\xi\text{d}\xi\\ \displaystyle B_{k}=\dfrac{2}{k\pi a}\int_{0}^l\psi(\xi)\sin\dfrac{k\pi}{l}\xi\text{d}\xi \end{cases} \end{aligned}
      求得 A_{k},B_{k},带入③就得到了初边值问题的解.

    初边值问题
    \begin{aligned} \begin{cases} u_{tt}-a^2u_{xx}=f(x,t)\\ t=0:u=0,\;u_{t}=0\\ u(0,,t)=u(l,t)=0 \end{cases} \end{aligned}

    \displaystyle u(x,t)=\sum_{k=1}^{\infty}\int_{0}^{t}B_{k}(\tau)\sin\dfrac{k\pi a}{l}(t-\tau)\text{d}\tau\cdot\sin\dfrac{k\pi}{l}x
    其中 \displaystyle B_{k}(\tau)=\dfrac{2}{k\pi a}\int_{0}^lf(\xi,\tau)\sin\dfrac{k\pi}{l}\xi\text{d}\xi


    初边值问题:
    \begin{aligned} \begin{cases} u_{tt}-a^2u_{xx}=f(x,t)\\ t=0:\;u=\varphi(x),u_{t}=\psi(x)\\ u(0,t)=u(l,t)=0 \end{cases} \end{aligned}

    Sol:
    \begin{aligned} \begin{cases} u_{tt}-a^2u_{xx}=0\\ t=0:\;u=\varphi(x),u_{t}=\psi(x)\\ u(0,t)=u(l,t)=0 \end{cases} \quad \begin{cases} u_{tt}-a^2u_{xx}=f(x,t)\\ t=0:\;u=0,u_{t}=0\\ u(0,t)=u(l,t)=0 \end{cases} \end{aligned}
    的解为 u_1,u_2,知原初边值问题的解为 u=u_1+u_2.

    相关文章

      网友评论

          本文标题:初边值问题的分离变量法

          本文链接:https://www.haomeiwen.com/subject/hyjkohtx.html