matlab偏微分方程数值解误差_偏微分方程数值解不稳定性的本质

应该说，传统偏微分方程数值解的思路就是用微分代替差分，与常微分方程数值解无异。但是真正属于这个领域特色的东西则是格式收敛性、稳定性的分析。初学者通常会惊异于差分格式居然在某些情况下会不稳定，出现 blow up 现象，得到面目全非的解。对此，国内外绝大多数教科书只是用傅里叶分析、能量估计法等给出分析稳定性的判别条件，而未能深入探讨数值不稳定性产生的原因。

笔者自今年入坑PDE数值解以来，对此问题也没有一个清晰的答案，直到接触的问题相对多了一些，最近才给出了一个至少能让自己满意的回答，笔者自己将其归结为“过度线性化”。本文我们通过一些最简单的例子，用最浅显的数学语言来解释数值不稳定性的本质。