matlab在有限差分法中的应用,MATLAB在有限差分法数值计算中的应用

第 卷 第 期 四川理工 学院学报 自然科学版仃 闷以 仃 从 年 月 文章编号 一 洲沁 一 在有限差分法数值计算中的应用 赵德奎 ’ , 刘 勇 ’ 四川理工 学院物理 系 , 四川 自贡 四川省人工智能重点 实验室 , 四川 自贡 摘 要 介绍 了应 用有限差分法求电位分布的一般步骤 , 给 出了 仿真的程序设计流程图 , 并针对静 电场 中轴对称情形 下的正六边形 电位分布进行 仿真 。 关键词 有限差分法 正六边形 电位 中图分类号 文献标识码 引 言 关于电磁场的数值计算 , 常用的方法有 有限单元法 , 有限差分法 , 有限体积法等。 由于这些方法 只能获得近似解 , 因此 , 利用软件进行仿真或者求数值解就显得非常必要 。 近年来有大量的文献在这方面进行了探讨 。 例如 文献【 利用 对载流螺线管的磁场分布进行了仿真 文献 利用 对静电场的电位分布进行了仿真 文献 利用 对正六边形二维场域的电位分布给出了数值解 。 本文在综合各文献的基础之上 , 基于 , 对有限差分法求解静电场的电位分布给出一般的程序设计方法 , 并结合具体的例子进行说明 。 有限差分法的基本思想 有限差分法的基本思想是把连续的定解区域用有限个离散点构成的网格来代替 , 这些离散点称作网格的节点 把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似 把原方程和定解条件中的微商用差商来近似 , 积分用积分和来近似 , 于是原微分方程和定解条件就近似地代之以代数方程组 , 即有限差分方程组 , 解此方程组就可 以得到原问题在离散点上的近似解 。 然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解 。 考虑到有限差分法和有限单元法在求解步骤的雷同 , 参照文献【 , 给出有限差分法数值计算的基本步骤 区域的离散或子 区域划分插值函数的选择方程组的建立 力 一 程组的求解 。 一般模型 的建立 与数据流程 图 设电位分布只是 , 的函数 , 沿 方向没有变化 , —十 — 日尹 日 则拉普拉斯方程 击 一 伽“ 其中尹 , 夕 是关于 , 的电位分布函数 。 考虑静电场边值条件 图 有限差分的网络分割 收稿 日期 一 一 作者简介 赵德奎 一 , 男 , 重庆人 , 助教 , 研究方 向 为模式识别与 智能控制 。 四川理工学院学报 自然科学版 年 月 子 ·子 · ‘“ 区域 ”内, 冲 八“ 其中 是 尹 , 在边界 上的分布函数 。 采用 图 的正方形网格进行离散 , 拉普拉斯方程变为 码 仇 仇 势 一 汽 用高斯一塞德而迭代法 , 得到 砚 户、 , 之 , , 图 节点顺序 式 ’一 专广 、 式只 试 试一 , 其中 , 表示节点所在的位置 。 表示迭代次数 。 下同 。 在高斯一塞德尔迭代中 , 网格节点一般按“ 自然顺序 ”排列 , 即先“从左到右 ” , 再“从 下到 匕 , 的顺序排列 , 如图 所示 。 迭代也是按 自然顺序进行 。 若采川逐次超松弛法 , 得到 。户 ’一 。汽, 十号, 九, , , 十 此、 蛤匀 试卫一 此 其中 称为加速收敛因子 , 其值的选取直接影响迭代的收敛速度 。 基 于式 的迭代公式 , 可以给出程序设计流程图 , 如图 所示一 个轴对称情形 下 电位分布的仿真程序 轴对称电磁场问题是电工设备设计分析中常遇到的一大类问题 。 我们以文献 提到的一个例子 , 稍加修改 , 得到下面一个例子 。 边 长为“ 的正六边形 的二维场域内无电荷分布 , 各边的电位依次为 十 , 一 , 一 。