第四章边值问题3.pptVIP

第四章 静态场边值问题的解法 * * 当场域边界的几何形状比较复杂时，很难用解析法进行分析。应采用数值计算法。 4.5 有限差分法 设有一函数f(x)，当独立变量x有一微小增量?x＝h，相应f(x)的增量为：。 通过离散化模型上各离散点的数值解来逼近连续场域内的真实解。 有限差分法将连续场域内的问题转化为离散系统的问题， 1．差分原理 不同于增量为无限小的微分，差分被称为有限差分。 当h很小时，?f(x)?df(x). ?f(x)=f(x+h)-f(x),称为函数f(x)的差分 中心差分?f(x)=f(x+h/2)-f(x-h/2) 一阶差商： 二阶差商 偏导数也可用差商近似表示。因而偏微分方程可表示为差分方程（代数方程）。 2．以二维场为例，将边值问题转化为一组差分方程组（代数方程组）。 (1)决定离散点的分布方式。 按正方网格划分，网格边长（步长）h，网格线的交点称结点。 设结点O上的电位为?(xo,yo)= ?o, 结点1，2，3，4上的电位为?1，?2，?3，?4。 设边值问题是 任一点x的电位 考虑1，3两点x1=xo+h, x3=xo-h 边界条件也可进行离散化处理，对第一类边值，可直接把点函数f(s)的值赋予各边界结点。 3．差分方程的解法 设将场域划分如图. 边界上的值分别为f1,………f16。 在各内点上作出差分，泊松方程变成下列差分方程组 解出关于?