11. 可分离的偏微分方程 PDE

  • 受张力的绳 suspended chain

    • 首先位移函数是一个双变量函数 u(x,t)u(x,t)

      • 考虑二维平面,具有时间属性
    • 我们假定 utu_t 表示对 tt 求导结果,同理 uttu_{tt} 表示对时间求二阶导的结果

    • 微分方程等式 utt=c2uxxu_{tt} = c^2\cdot u_{xx}, c2=Ftenρc^2 = \frac{F_{ten}}{\rho}

      • FtenF_{ten} 表示的是绳子内部的张力,假设绳子静止 (time-independent)且绳子重力效果不大
      • ρ\rho 表示绳子密度
      • 为了细化公式,定义 Fgen=ρgxF_{gen} = \rho gx
    • 由于绳子的两端固定,这就是经典的 BVP 边界问题

  • 分离思路

    • u(x,t)=X(x)T(t)u(x,t) = X(x)\cdot T(t) 来获得两个单变量方程的积

    • 我们要证明这个分离是存在且唯一的

    • 回到上面的 suspended chain 例子

      • u(x,t)=X(x)T(t)u(x,t) = X(x)T(t)

      • 利用基本等式 utt=c2uxxXT=c2XTu_{tt} = c^2 u_{xx}\Rightarrow XT'' = c^2X''T

      • 这下就是可分离变量了,我们写成等式 1c2TT=1XX\frac{1}{c^2T}T'' = \frac 1 X X''

        • 由于两边变量不同,那么两边就是常数,令之为 λ-\lambda
        • 这下就有两个特征方程
      • 同时由于两端固定,我们有等式 {X(0)T(0)=0X(l)T(0)=0\begin{cases}X(0)T(0)=0\\X(l)T(0)=0\end{cases}

        • 从这个可以看出要么 T(t)=0T(t) =0 恒定,或者看出 XX 在两端为 0
        • 于是就有 X+λX=0,X(l)=X(0)=0X'' + \lambda X = 0, X(l) = X(0) = 0
      • 特征值解为 λn=(nπl)2\lambda_n = (\frac{n\pi}{l})^2

      • 采用上一讲的 SL 试根得到解集 Xn=Cnsin(λnx)=Cnsin(nπxl)X_n = C_n\sin(\sqrt{\lambda_n}x) = C_n\sin(\frac{n\pi x}{l})

        • 一般这个时候我们有边界问题进行约束,没有的话解集就是上式
        • 同理 Tn=Dncoscnπtl+EnsincnπtlT _n = D_n\cos\frac{cn\pi t}l +E_n\sin\frac{cn\pi t}l
      • 总的结果是 u=Σn=1(Fncoscnπtl+Gnsincnπtl)sinnπxlu= \Sigma_{n=1}^\infty (F_n\cos\frac{cn\pi t}{l}+G_n\sin\frac{cn\pi t}l)\sin\frac{n\pi x}{l}

        • 边界问题 u(x,0)=ΣFnsinnπxlu(x,0) = \Sigma F_n\sin\frac{n\pi x}{l} , ut(x,0)=ΣGncnπlsinnπxlu_t (x,0)=\Sigma G_n\frac{cn\pi}{l}\sin\frac{n\pi x}l

        • 这里需要将这个 Fn,GnF_n,G_n 变成傅里叶级数会更好分析一点

        • f(x)=u(x,0)=2l0lf(x)sin(nπxl)dxsin(nπxl)f(x) = u(x,0) = \frac 2 l \int_0^l f(x)\sin(\frac{n\pi x}l)dx\cdot \sin(\frac{n\pi x}l)

          • 这下中间部分就是 FnF_n 了,后面不要忘记 × sin\sin
          • 注意最终结果是回到 uu 本身,也就是上面加粗处
    • 唯一解定理:对于形如 uttc2uxx=F(x,t)u_{tt} -c^2u_{xx} = F(x,t),拥有两个边界问题 u(0,t)=a(t),u(l.t)=b(t)u(0,t) = a(t),u(l.t)=b(t) 和初始条件 u(x,0)=f(x),ut(x,0)=g(x)u(x,0)=f(x),u_t(x,0)= g(x)

  • 第一类贝塞尔方程的解系