波动方程与特征线法

波动方程与特征线法

公式有点多，建议电脑端看。

1、特征线法

（1） 特征线法自我总结

将方程转化为${\partial }_{t}u+p{\partial }_{x}u=q$，则有两个常微分方程：
$$\begin{array}{rl}& \frac{dx}{dt}=p,x(0)=x_0\\ & \frac{dU}{dt}=q,U(0)=含x_0的表达式\end{array}$$
先解出x(t)，再写出U(t)，U(t)中往往含有$x_0$，将$x_0$用含有x的式子表示。

（2）特征线法与初值问题

$$\begin{array}{rl}& {\partial }_{t}u+(x+t){\partial }_{x}u+u=x,u(x,0)=x,x\in \mathbb{R},t>0\\ & \therefore \frac{dx}{dt}=x+t,x(0)=x_0\\ & \frac{dU}{dt}+U=x(t,t_0),U(0)=x_0\end{array}$$

2、齐次化原理与一维初值问题

对于如下一维波动方程的初值问题：
$$\square u={\partial }_{tt}u-a^2\Delta u$$

$$\square u=f,u(x,0)=\phi (x),{\partial }_{t}u(x,0)=\psi (x)$$

根据线性叠加原理，有：
$$\begin{array}{rl}& \square u_1=0,u_1(x,0)=\phi (x),{\partial }_t{u}_1(x,0)=0;\\ & \square u_2=0,u_2(x,0)=0,{\partial }_t{u}_2(x,0)=\psi (x);\\ & \square u_3=f,u_3(x,0)=0,{\partial }_t{u}_3(x,0)=0;\end{array}$$
得到：初值问题的解：$u=u_1+u_2+u_3$
$$\begin{array}{rl}& u_1(x,t)=\frac{1}{2}(\phi (x+at)+\phi (x-at))\\ & u_2(x,t)=\frac{1}{2a}{\int }_{x-at}^{x+at}\psi (\xi )d\xi \\ & u_3(x,t)=\frac{1}{2a}{\int }_0^t{\int }_{x-a(t-\tau )}^{x+a(t-\tau )}d(\xi ,\tau )d\xi d\tau \end{array}$$
所以，一维波动方程的初值问题的解为：
$$u(x,t)=\frac{1}{2}(\phi (x+at)+\phi (x-at))+\frac{1}{2a}{\int }_{x-at}^{x+at}\psi (\xi )d\xi +\frac{1}{2a}{\int }_0^t{\int }_{x-a(t-\tau )}^{x+a(t-\tau )}d(\xi ,\tau )d\xi d\tau$$
当$f\equiv 0$的时候，得到D’Alembert公式：
$$u(x,t)=\frac{1}{2}(\phi (x+at)+\phi (x-at))+\frac{1}{2a}{\int }_{x-at}^{x+at}\psi (\xi )d\xi$$

3、一维半无界定解问题

在半无界问题区域 Q ‾ = { 0 ≤ x < ∞ , 0 ≤ t < ∞ } \overline{Q}=\{0 \leq x< \infty,0\leq t <\infty\} Q​={0≤x<∞,0≤t<∞}上考虑定解问题：

$$\begin{array}{r}\{\begin{array}{rc}\square u& =f(x,t)\\ u(x,0)& =\phi (x)\\ {\partial }_{t}u(x,0)& =\psi (x)\\ u(0,t)& =g(t)\end{array}\end{array}$$
考虑$u=v+g(t)$，所以得到：
$$\begin{array}{r}\{\begin{array}{rc}\square v& =f(x,t)\\ v(x,0)& =\phi (x)\\ {\partial }_{t}v(x,0)& =\psi (x)\\ v(0,t)& =0\end{array}\end{array}$$
延拓法，将$f,\phi ,\psi$延拓成奇函数，转化为初值问题，得到初值问题的解，然后按照$x-at$的正负分类讨论，代入$f,\phi ,\psi$：
$$\begin{array}{rl}& 当x\ge at时，\\ & u(x,t)=\frac{1}{2}(\phi (x+at)+\phi (x-at))+\frac{1}{2a}{\int }_{x-at}^{x+at}\psi (\xi )d\xi +\frac{1}{2a}{\int }_0^t{\int }_{x-(t-\tau )}^{x+a(t-\tau )}f(\xi ,\tau )d\xi d\tau \\ & 当x<at时，\\ & u(x,t)=\frac{1}{2}(\phi (x+at)-\phi (at-x))+\frac{1}{2a}{\int }_{at-x}^{x+at}\psi (\xi )d\xi +\frac{1}{2a}{\int }_{t-\frac{x}{a}}^t{\int }_{x-(t-\tau )}^{x+a(t-\tau )}f(\xi ,\tau )d\xi d\tau \\ & +\frac{1}{2a}{\int }_0^{t-\frac{x}{a}}{\int }_{(t-\tau )-x}^{x+a(t-\tau )}f(\xi ,\tau )d\xi d\tau \end{array}$$

4、高维初值问题

从一维初值问题出发，通过球面平均法推导出三维波动方程初值问题的解的表达式，再通过降维法推导二维波动方程初值问题的解的表达式。
这里只是简单写了球面平均法，降维法以后再敲公式。

高维的初值问题：
$$\begin{array}{r}\{\begin{array}{l}\square u={\partial }_{tt}u-a^2\Delta u=f(x,t)\\ u(x,0)=\phi (x)\\ {\partial }_{t}u(x,0)=\psi (x)\end{array}\end{array}$$

（1）三维波动方程的初值问题——球面平均法

我们同样有齐次化原理，只需要解出$u_2$的方程即可：
$$\square u_2=0,u_2(x,0)=0,{\partial }_t{u}_2(x,0)=\psi (x);$$
我们将直角坐标转换为空间球坐标$(r,\theta ,\phi )$，则有：
$$\frac{1}{a^2}\frac{{\partial }^2{u}_2}{\partial t^2}=\frac{1}{r^2}\frac{\partial }{\partial r}(r^2\frac{\partial u_2}{\partial r})+\frac{1}{r^2\sin \theta }\frac{\partial }{\partial \theta }(\sin \theta \frac{\partial u_2}{\partial \theta })+\frac{1}{r^2{\sin }^2\theta }\frac{{\partial }^2{u}_2}{\partial {\phi }^2}$$

在球对称假设下（即$u_2$不依赖$\theta ,\phi$），可以化简为：
$$\frac{1}{a^2}\frac{{\partial }^2(r{u}_2)}{\partial t^2}=\frac{{\partial }^2(r{u}_2)}{\partial r^2}$$
显然，球对称情况下，$r{u}_2$满足一维波动方程，所以可以求出$u_2$然后再根据齐次化原理求出$u_1,u_3$.