27、边界值问题的有限差分法求解

边界值问题的有限差分法求解

1. 有限差分近似

在直角坐标系的网格中，我们可以通过组合之前的方程来得到一些有限差分近似。
- 拉普拉斯算子的近似 ：对于(\nabla^{2}z=\frac{\partial^{2}z}{\partial x^{2}}+\frac{\partial^{2}z}{\partial y^{2}})，其有限差分近似为：
- (\nabla^{2}z_{i}\approx\frac{z_{l}-2z_{i}+z_{r}}{h^{2}}+\frac{z_{a}-2z_{i}+z_{b}}{k^{2}}\approx\frac{r^{2}z_{l}+r^{2}z_{r}+z_{a}+z_{b}-2(1 + r^{2})z_{i}}{r^{2}h^{2}})，其中(r = \frac{k}{h})。
- 当(r = 1)时，(\nabla^{2}z_{i}\approx\frac{z_{l}+z_{r}+z_{a}+z_{b}-4z_{i}}{h^{2}})。这些中心差分近似的误差为(O(h^{2}))。
- 二阶混合导数的近似 ：对于(\frac{\partial^{2}z}{\partial x\partial y})，其有限差分近似为(D_{xy}z_{i}\approx\frac{(z_{r}-z_{l}) {a}/(2h)-(z {r}-z_{l}) {b}/(2h)}{2k}\approx\frac{z {ar}-z_{al}-z_{br}+z_{bl}}{4hk})。

我们还可以在其他坐标系（如斜坐标系和极坐标