有限差分方法matlab,有限差分方法的MATLAB编程

第 卷 第 期年 月 《新疆师范大学学报 》 自然科学版 , , 有限差分方法的 编程 王 飞 , 裴永祥 新班师范大学数理信息学院 乌鲁木齐 , 摘 要 实践中应用的电磁场 , 其场域的边 界大多数 比较复杂 , 使解析法难 以应用 虽然有些 电磁 问题经简化后 能解析 求解 , 但解的精度往往难 以达 到实践应用 的要求 。 因此 , 数值计算 电磁场问题成 了一项重要 的工作 。 本文用 材月 软件编程计算二维场 , 方便简便并可 画 出图形 , 更直观地反映二维场的性质 。 关键词 二维场 中图分类号 团 、 〕 文献标识码 文章编号 一 。。 于 卜 一 一 对于大多数科研 工作者来说 , 意味着兴趣 。 软件将人们从乏味的 、 编程 中解放出来 , 使他们真正地把精力放在科研和设计 的核心 问题上 , 大大提高了科研产 出和发现 的步伐 。 在 问世 近 年的时间里 , 是 国际公认的最优秀的科技应用软件之一 , 其系列产品风靡全球 , 成为众多科研工作者和工程设计人员 的首选工作平 台 。 大学生学习 软件 , 一方面使教学 中对难度大而且很抽象的数学 问题形象化 , 另一方面为大学生今后工作打下 良好的基础 。 有 限差分法是 以差分原理为基础的数值计算方法 。 它具有简单 、 直观的特点 , 最早在 电磁场数值场数值分析领域得到 了广泛 的应用 。 即使在 目前有限元法 日益得到广泛应用的时代 , 有限差分法 以其 固有特点仍然是一种不可忽视 的方法 。 对于包括 电磁场在 内的各种物理问题 , 应用有限差分法进行计算的步骤通常为采用一定的网格划分格式离散化场域 , 把实际连续的场离散为有 限 的多个点 , 用这些离散点上 的参数近似描述实际上 连续 的场 , 也就是所谓场 的“ 离散 ” 。 基 于差分原理 的应用 , 对场域内偏微分方程 以及场域边界上的边界条件 包括场域内不同媒质分界面上 的边界条件 进行差分离散化处理 , 即用差商代替偏导数 , 给出相应的差分计算格式 。 结合一定 的代数方程组 的解法 , 编制计算程序 , 求解 由上述步骤所得到的对应 于待求边值问题的差分方程 的解 , 即求 出节点的位 函数值 , 据此进一步求出场强分布 。 差分 设一 函数 , 当其独立变量 有很小的 △ 一 , 所 以相应地函数 的增量为 △ 一 一 △ 被称为函数 的差分 一阶 。 它与微分不同 , 因为差分是有限量的差 , 故通常也被称为有限差分 。 但只要增量 很小 , 差分 △ 与微分 之间的差异将很小根据差分 的定义 , 在差分运算 中还 常用一阶中心差分 , 其公式为 〔收稿 日期 」 〔作者简介 〕王 飞 一 , 男 , 汉族 , 新疆师范大学物理系在读本科生 。 第 期 裴永祥等 有 限差分方法 的 编程 、 , 、 自少又 一 又 十 二 少 一 一 下 , 乙 乙 一 阶差分仍是独立变量 的函数 。 用式 计算一阶差分 的类似差分方法求二阶差分 , 得到 △ , 称为原始 函数 的二 阶差分 。 同样 , 当 很小时 , 二阶差分 △, 很接近 于二阶微分 少 。 平行平面场 中的差分格式 在偏微分方程 中常常提 出如下 的第一边值 问题 即在平面 区域 内 厂 上分别 给 出 了微分方程和边 值条件 挤 铲凸 一 二二 十 二下 万 一 中 , 一 , , 当 , 任 当 , 〔 其 中 中 是在边界 上 的连续 函数 。 差分格式 的构造 我们在 区域 上讨论最简单的泊松方程 挤 又, 及 又 , 〔理