PDE数值解（二） 截断误差

既然是PDE的数值解，就避免不了一个问题——截断误差。如何理解PDE中的截断误差呢？

对于一个差分格式，事实上只需要将$u$本身直接代入就能够得到截断误差。举例来说，某扩散方程的差分格式为
$$\frac{u_j^{n+1}-u_j^n}{\tau }-a\frac{u_{j+1}^n-u_j^n}{h}=0$$只需要将$u_j^n$替换为$u(x_j,t_n)$，上式即变为了
$$\frac{u(x_j,t_{n+1})-u(x_j,t_n)}{\tau }-a\frac{u(x_{j+1},t_n)-u(x_j,t_n)}{h}=\xi$$接下来就是对上式的化简，显然可以利用刻在骨子里面的Taylor展开。
$$\frac{1}{\tau }\left[\sum _{i=0}^{+\infty }\frac{\frac{{\partial }^{i}u}{\partial t^i}}{i!}{\tau }^i-u\right]-\frac{a}{h}\left[\sum _{i=0}^{+\infty }\frac{\frac{{\partial }^{i}u}{\partial x^i}}{i!}{h}^i-u\right]=\xi$$结合原PDE方程化简一下就可以得到如下结果
$$\xi =\sum _{i=2}^{+\infty }\frac{\frac{{\partial }^{i}u}{\partial t^i}}{i!}{\tau }^{i-1}-a\sum _{i=2}^{+\infty }\frac{\frac{{\partial }^{i}u}{\partial x^i}}{i!}{h}^{i-1}=o(\tau +h）$$这就是此PDE数值解的截断误差。事实上是套路性非常强的一个东西。
（未完待续）