有限差分法的原理与计算步骤.docxVIP

有限差分法的原理与计算步骤 有限差分法的原理与计算步骤 PAGE 有限差分法的原理与计算步骤 有限差分法的原理与计算步骤 原理 基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似， 积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组 ， 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 计算步骤 在采用数值计算方法求解偏微分方程时，若将每一处导数由有限差分近似公式替代，从而把求解偏微分方程的问题转换成求解代数方程的问题，即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下： （1）区域离散化，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格； （2）近似替代，即采用有限差分公式替代每一个格点的导数； （3）逼近求解。换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程 二、有限差分法的程序流程图