24、偏微分方程的有限差分法求解

偏微分方程的有限差分法求解

1. 双曲型偏微分方程

双曲型偏微分方程中，我们考虑一维波动方程：
[c^{2}\frac{\partial^{2}u}{\partial x^{2}}=\frac{\partial^{2}u}{\partial t^{2}}, 0 < x < L, t > 0]
此方程描述了拉紧弦中的波，其中 (c) 是弦中波的传播速度。通过与标准二阶偏微分方程对比，可知 (B^{2}-4AC=-4c^{2}(-1)>0)，所以该方程是双曲型的，并且需要 (x = 0) 和 (x = L) 处的边界条件以及 (t = 0) 时的初始条件。

1.1 有限差分近似

将 (L) 分成 (n) 段，使得 (h = L/n)，并考虑时间步长为 (k)。在节点 ((i,j)) 处，基于中心有限差分近似，可得：
[c^{2}\left(\frac{u_{i - 1,j}-2u_{i,j}+u_{i + 1,j}}{h^{2}}\right)=\left(\frac{u_{i,j - 1}-2u_{i,j}+u_{i,j + 1}}{k^{2}}\right)]
经过整理，得到：
[u_{i,j + 1}=\alpha^{2}(u_{i - 1,j}+u_{i + 1,j})+2(1 - \alpha^{2})u_{i,j}-u_{i,j - 1}]
其中 (\alpha^{2}=c^{2}k^{2}/h^{2})。当 (j = 0) 时，有：
[u_{i,1}=\alpha^{2}(u_{i - 1,0}+u_{i + 1,0})+2(1 - \alpha^{2})u_{i,0}-u