分离变量法

分离变量法核心思路

分离变量法的核心思路是，将PDE(偏微分方程)变为多个ODE(常微分方程)，那么我们如何进行呢？

PDE和ODE的一个显著不同是：PDE的自变量一定有两个以上，而ODE的自变量只有一个，那么我们能否将多个自变量拆分使求解PDE变为求解ODE呢？如果满足PDE的自变量 $u(x,t)$ 是某种特别的形式，例如
$$u(x,t)=X(x)T(t),$$
可以多个一元函数相乘得到，那么每个一元函数 $X(x)$ 和 $T(t)$ 满足的方程应该都是ODE，只要能求出 $X(x)$ 和 $T(t)$，就可能求出 $u(x,t)$。这样我们就把求解PDE问题转变为求解ODE问题，这就是分离变量法的思路。

ODE求解思路：先求出微分方程的特解，由线性无关的特解叠加出通解，再用定解条件去确定叠加系数。

在力学中，驻波的表达式为
$$u(x,t)=2A\cos \frac{2\pi x}{\lambda }\cos 2\pi \gamma t$$
这使我们自然想到，对于定解问题可设其特解为
$$u(x,t)=X(x)T(t)$$
其中 $X(x)$ 和 $T(t)$ 分别只是变数 $x$ 和 $t$ 的函数。为了弄清楚定解问题究竟有什么样的驻波解，应将 $u(x,t)=X(x)T(t)$ 式分别代人方程和定解条件中。

考虑长为 $l$ 、两端固定的弦的自由振动, 方程及定解条件为
$$\begin{array}{llc}\frac{{\partial }^{2}u}{\partial t^2}-a^2\frac{{\partial }^{2}u}{\partial x^2}=0,& 0<x<l,t>0,& \\ {u|}_{x=0}=0,& {u|}_{x=l}=0,& t⩾0,\\ {u|}_{t=0}=\varphi (x),& {\frac{\partial u}{\partial t}|}_{t=0}=\psi (x),& 0⩽x⩽l.\end{array}$$

分离变量法求解步骤

求解分为四步：分离变量、求解本征值、求特解并叠加出一般解、利用本征函数的正交性叠加系数。

详细求解过程每本数理方法课本上都有详细介绍，证明篇幅也较长，故笔者只对其进行简要说明。

分离变量：我们希望求得的特解具有分离变量的形式，则$u(x,t)=X(x)T(t)$，带回波动方程则有$X(x)T^{\prime \prime }(t)=a^2{X}^{\prime \prime }(x)T(t)$，在波动方程中，左端和右端相等，就必须共同等于一个既与 $x$ 无关、又与 $t$ 无关的常数 $-\lambda$，$\frac{1}{a^2}\frac{T^{\prime \prime }(t)}{T(t)}=\frac{X^{\prime \prime }(x)}{X(x)}=-\lambda$，上面的结果又可以写成
$$\begin{array}{rl}& T^{\prime \prime }(t)+\lambda a^{2}T(t)=0,\\ & X^{\prime \prime }(x)+\lambda X(x)=0.\end{array}$$
这样我们就得到了两个常微分方程，同时引入了一个待定常数 $\lambda$。

求解本征值：

分离变量的结果是得到了多个含有待定常数的齐次ODE和齐次边界条件，它们一同构成常微分方程的本征值问题。当待定参数 $\lambda$ 取某些特定值时，可以有同时满足ODE $X^{\prime \prime }(x)+\lambda X(x)=0$与齐次边界条件${u|}_{x=0}=0,{u|}_{x=l}=0$的非零解 $X(x)$ ，$\lambda$ 的这些特定值称为本征值，相应的非零解$X(x)$称为本征函数。

详细过程略去

最终我们求得了所有的本征值 ${\lambda }_n$ 以及相应的本征函数 $X_n(x)$，有无穷多组，可以用脚标 $n$ 进行标记。

求特值并叠加出一般解：

在求解本征值问题后, 将本征值 ${\lambda }_n$ 代入方程$T^{\prime \prime }(t)+\lambda a^{2}T(t)=0$，求出相应的 $T_n(t)$，$T_n(t)=C_n^{\prime }\cos \frac{n\pi }{l}at+D_n^{\prime }\sin \frac{n\pi }{l}at$。

因此，就得到了同时满足偏微分方程和边界条件的特解
$$u_n(x,t)=\left(C_n\cos \frac{n\pi }{l}at+D_n\sin \frac{n\pi }{l}at\right)\sin \frac{n\pi }{l}x,n=1,2,3,\cdots .$$
这样的特解有无穷多个，每个特解都满足齐次偏微分方程和齐次边界条件，但单独一个特解不足以满足给定的初始条件，我们需要将全部特解叠加起来，得到
$$u(x,t)=\sum _{n=1}^{\infty }(C_n\cos \frac{n\pi }{l}at+$$