偏微分方程（二）

（未完，待补充）

大致思路和一般步骤

起源：机械振动或电磁振动（波动方程的解）总可以分解为具有各种频率和振幅的简谐运动的叠加，而每个简谐运动具有形式：
$$e^{iw(t+sx)}=e^{iwt}{e}^{iswx}$$
这正是物理上的驻波。这就表明波动方程的解可写为驻波的叠加，而每个驻波都只是只含变量 x 的函数和只含变量 t 的函数的乘积，即具有变量分离的形式。

适用范围：有界区域下的混合问题

一般步骤：分离变量得到关于不同变量的常微分方程组 $⟶$ 由边界条件寻找特征值和相应的特征函数 $⟶$ 由叠加原理得到带未知系数的解 $⟶$ 利用初值条件确定未知系数得到最终结果

相容性条件：略，但很重要

双曲型方程的混合问题

（1）一维齐次波动方程：
$$\{\begin{array}{ll}{u}_{tt}-a^2{u}_{xx}=0,& 0<x<l,t>0\\ & \\ u{|}_{t=0}=\phi (x),u_t{|}_{t=0}=\psi (x),& 0\le x\le l\\ & \\ u{|}_{x=0}=0,u{|}_{x=l}=0,& t\ge 0\end{array}$$

解：

1. 变量分离：$u(x,t)=X(x)T(t)$, 代回 $u_{tt}-a^2{u}_{xx}=0$，有 $\frac{X^{\prime \prime }(x)}{X(x)}=\frac{1}{a^2}\frac{T^{\prime \prime }(t)}{T(t)}=-\lambda$
   $$\{\begin{array}{l}{X}^{\prime \prime }(x)+\lambda X(x)=0,\\ \\ T^{\prime \prime }(t)+\lambda a^{2}T(t)=0\end{array}$$

2. 利用边界条件，有
   $$\{\begin{array}{l}{X}^{\prime \prime }(x)+\lambda X(x)=0,\\ \\ X(0)=X(l)=0\end{array}$$
   为了获得平凡解，最终得到
   $$(c_1=0,c_2\ne 0)\sin \sqrt{\lambda }l=0\text{\hspace{1em}(*)}$$
   得到特征值和相应特征函数，
   $$\begin{gathered} {\lambda }_k=(\frac{k\pi }{l}{)}^2,(k=1,2,\cdots ) \\ {X}_k(x)=c_k\sin \frac{k\pi x}{l},(k=1,2,\cdots ) \end{gathered}$$

3. 对于每一个特征值，关于 $T$ 的方程的通解可写为
   $$T_k(t)=a_k\cos \frac{k\pi a}{l}t+b_k\sin \frac{k\pi a}{l}t,k=1,2,\cdots$$
   于是，对任意的 $A_k=c_k{a}_k,B_k=c_k{b}_k$，函数
   $$u_k(x,t)=X_k(x)T_k(t)=(A_k\cos \frac{k\pi a}{l}t+B_k\sin \frac{k\pi a}{l}t)\sin \frac{k\pi x}{l}$$
   满足方程和边界条件，由叠加原理（方程和边界条件都是线性的）
   $$u(x,t)=\sum _{k=1}^{\infty }(A_k\cos \frac{k\pi a}{l}t+B_k\sin \frac{k\pi a}{l}t)\sin \frac{k\pi x}{l}$$

4. 利用初值条件确定未知系数 $A_k,B_k$
   $$\{\begin{array}{l}\phi (x)=u(x,0)={\sum }_{k=1}^{\infty }{A}_k\sin \frac{k\pi }{l}x\\ \\ \psi (x)=\frac{\partial }{\partial t}u(x,0)={\sum }_{k=1}^{\infty }{B}_k\frac{k\pi a}{l}\sin \frac{k\pi }{l}x\end{array}$$
   利用 $\phi ,\psi$ 在 $[0,l]$ 上的傅里叶级数，可得
   $$\{\begin{array}{l}{A}_k=\frac{2}{l}{\int }_0^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}\\ \\ B_k=\frac{2}{k\pi a}{\int }_0^l\psi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}\end{array}$$

• 注1：物理意义——两端固定的弦的自由振动

• 注2：（*）—— 实际上要对 $\lambda$ 的不同取值分情况讨论
  $$\{\begin{array}{ll}X(x)=c_1{e}^{\sqrt{-\lambda x}}+c_2{e}^{-\sqrt{-\lambda x}},& \lambda <0\\ & \\ X(x)=c_1+c_{2}x,& \lambda =0\\ & \\ X(x)=c_1\cos \sqrt{\lambda }x+c_2\sin \sqrt{\lambda }x,& \lambda >0\end{array}$$

• 注3： fourier级数只展开成含有正弦项（余弦项）的级数
  解释为什么 $A_k=\frac{2}{l}{\int }_0^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}$，函数要展开成只有正弦项的级数，说明这是个奇函数（在后面的推理中会解释），这里，我们定义在 $[0,l]$ 上的函数 $\phi$ 需要改造成一个奇函数，因此我们进行了奇延拓，把 $[0,l]$ 扩大到 $[-l,l]$，此时便有
  $$\begin{array}{rl}{A}_k& =\frac{1}{l}{\int }_{-l}^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}\\ & =\frac{1}{l}[{\int }_0^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}+{\int }_{-l}^0\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}]\\ & =\frac{1}{l}[{\int }_0^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}+{\int }_{-l}^0\phi (-x)\sin \frac{-k\pi x}{l}\mathtt{d}\mathtt{x}]\\ & =\frac{1}{l}[{\int }_0^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}+{\int }_0^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}]\\ & =\frac{2}{l}{\int }_0^l\phi (x)\sin \frac{k\pi x}{l}\mathtt{d}\mathtt{x}\end{array}$$
  而
  $$\begin{array}{rl}{B}_k& =\frac{1}{l}{\int }_{-l}^l\phi (x)\cos \frac{k\pi x}{l}\mathtt{d}\mathtt{x}\\ & =\frac{1}{l}[{\int }_0^l\phi (x)\cos \frac{k\pi x}{l}\mathtt{d}\mathtt{x}+{\int }_{-l}^0\phi (x)\cos \frac{k\pi x}{l}\mathtt{d}\mathtt{x}]\\ & =\frac{1}{l}[{\int }_0^l\phi (x)\cos \frac{k\pi x}{l}\mathtt{d}\mathtt{x}-{\int }_{-l}^0\phi (-x)\cos \frac{-k\pi x}{l}\mathtt{d}\mathtt{x}]\\ & =\frac{1}{l}[{\int }_0^l\phi (x)\cos \frac{k\pi x}{l}\mathtt{d}\mathtt{x}-{\int }_0^l\phi (x)\cos \frac{k\pi x}{l}\mathtt{d}\mathtt{x}]\\ & =0\end{array}$$

抛物型方程的混合问题

$$\{\begin{array}{ll}{u}_t-a^2{u}_{xx}=0,& 0<x<l,0<t<T\\ & \\ u{|}_{t=0}=\phi (x),& 0\le x\le l\\ & \\ u_x{|}_{x=0}=0,u_x{|}_{x=l}=0,& 0\le t\le T\end{array}$$
解：（模仿双曲型）

1. 变量分离：$u(x,t)=X(x)T(t)$, 代回 $u_t-a^2{u}_{xx}=0$，有 $\frac{X^{\prime \prime }(x)}{X(x)}=\frac{1}{a^2}\frac{T^{\prime }(t)}{T(t)}=-\lambda$
   $$\{\begin{array}{l}{X}^{\prime \prime }(x)+\lambda X(x)=0,\\ \\ T^{\prime }(t)+\lambda a^{2}T(t)=0\end{array}$$

2. 利用边界条件，有
   $$\{\begin{array}{l}{X}^{\prime \prime }(x)+\lambda X(x)=0,\\ \\ X(0)=X(l)=0\end{array}$$
   为了获得平凡解，

   当 $\lambda =0$ 时，${\lambda }_0=0,u_0(x,t)=X_0(x)T_0(t)=A_0$，

   当 $\lambda >0$ 时，
   $$(c_2=0,c_1\ne 0)\sin \sqrt{\lambda }l=0$$
   得到特征值和相应特征函数，
   $$\begin{gathered} {\lambda }_k=(\frac{k\pi }{l}{)}^2,(k=1,2,\cdots ) \\ {X}_k(x)=B_k\cos \frac{k\pi x}{l},(k=1,2,\cdots ) \end{gathered}$$

3. 对于每一个特征值，关于 $T$ 的方程的通解可写为
   $$T_k(t)=C_k{e}^{-(\frac{ak\pi }{l}{)}^{2}t},k=1,2,\cdots$$
   于是，对任意的 $A_k=B_k{C}_k$，函数
   $$u_k(x,t)=X_k(x)T_k(t)=A_k{e}^{-(\frac{ak\pi }{l}{)}^{2}t}\cos \frac{k\pi x}{l},k=1,2,\cdots$$
   满足方程和边界条件，由叠加原理（方程和边界条件都是线性的）
   $$u(x,t)=\frac{A_0}{2}+\sum _{k=1}^{\infty }{A}_k{e}^{-(\frac{ak\pi }{l}{)}^{2}t}\cos \frac{k\pi x}{l},$$

4. 利用初值条件确定未知系数 $A_k,k=0,1,2,\cdots$
   $$u(x,0)=\frac{A_0}{2}+\sum _{k=1}^{\infty }{A}_k\cos \frac{k\pi }{l}x=\phi (x)$$
   利用 $\phi$ 在 $[0,l]$ 上的傅里叶级数，可得
   $$A_k=\frac{2}{l}{\int }_0^l\phi (x)\cos \frac{k\pi x}{l}\mathtt{d}\mathtt{x},k=0,1,2,\cdots$$

5. 最终得到形式解
   $$u(x,t)=\sum _{k=0}^{\infty }(\frac{2}{l}{\int }_0^l\phi (\xi )\cos \frac{k\pi \xi }{l}\mathtt{d}\xi )e^{-(\frac{ak\pi }{l}{)}^{2}t}\cos \frac{k\pi x}{l}$$

• 当 $\lambda =0$ 时，方程的通解为 $X(x)=c_1+c_{2}x\Rightarrow c_2=0$，取 $c_1=1$，我们可以得到一个非平凡解 $X_0(x)\equiv 1$. 对于 ${\lambda }_0=\lambda =0$，$T_0(t)=A_0$（$A_0$ 为某个常数），此时就有 $u_0(x,t)=X_0(x)T_0(t)=A_0$