和PDE数值解相关的线性代数方程组求解.pdf

摘 要 在 偏 微 分 方 程 数值 解 中，常采 用 有 限差分 法 、有 限元法 以及 有 限体积 法 等 方 法 ，其 实质 都 是将连 续 的 问题 离散化 ，最后 问题 常 归结为方程 组 的求解 ，比较 常 见 的是大 型线 性 方程 组 ，非 线性 方程 组也往 往要 转 换成 线性 代 数 方程 组 的求 解 . 因此 对 于数 值 求解 偏 微 分方程 而 启 ，‘线代 数 方程 组求解 是 非 常重 要 的一环 . 方程 组 求解 方法 主要 是数值 代 数 研 究 的领 域 ，对 于偏微 分方 程 数 值解 研 究方 向的科 技 人 员而 言 ，往 往 并不 亲 自去研 究方程 组 的求解方法 以及 理 论 分析 ，而 注 重 的是算 法 的使 用 、算法 的结论 以及 算法 的计算 效 果 (计 算精 度 、快 慢 等) . 本 文 写作 目的 即是对 大 型 线代 数 方 程 组 求 解 的各 种 方 法 进 行讨 论 ，理 论 与实 践 相 结 合 ，对 一些 尚不 完善 的理 论 问题进 行 深 入 分析 ，并选 择 Po isson 方 程 差分格 式对 各种 方 法进 行 综合 实验 比对 ，提 供 求 解 方程 组 时算法 选 择 的依 据 . 在 直接 法 部 分 ，首 先应 该指 出在 条件 数大 小 与病 态/ 良态 问题 的关系 、判 断 一个 具体 问题 是 否病 态 的理 论分析 方 法 等还有 许 多值 得 细致研 究 的 问题 ，本 文 只 就通 常 的数 值 分析 教 材 中不 明确 或 易混淆 的一些地 方进 行 了说 明 ，并对 一些 理 论 问题 进 行 了深 入探 讨 ，提 出一些疑 问. 其 次针对 高斯 列 主元 消 去法进 行 了具 体 的 理 论 与 实践相 结合 的分 析 : 对 稠密矩 阵进 行 了数值 实验 ，对 该算 法 适 应 方 程 组类 型 、可求解 方程 组 的规模 、计 算速度 、计算 结果精度 等 问题给 出 了答 案 ; 分析 了 列 主元 消去法 的舍 入 误 差 理论 ，给 出 了因舍入 误 差 影 响导致 计 算 失效 的一个 例 子 ，指 出了舍 入误 差理论 不完善 的地 方 ; 针对 典型 的偏微分方程数值解 中的 Poi sson 方程 问题进 行 了进 一步 实验 ，了解 将连 续 问题 经 由差分格 式离散 化后 产 生 的大 型稀 疏 型方程 组用直接法来编 程求解 的过程 以及 最后 的计 算 效果 .最后 结 合 一维两 点边值 问题 数值 实验讨论 了三对 角方程 组求解 的追赶 法 . 在迭 代 法 一节 ，分 析 了各种基 本 迭代 算法 的收敛 条件 、收敛 速 度 ，为 比较 算 法 之 间计算 效 果 ，仍 采 用 PoiSSon 方 程 问题对 各种 算 法进 行 了数 值 实验 综 合 比 对 . 进 一步地 ，针对 Jac ob i 方法 收敛速度慢 的 问题 ，本文给 出了 Jacobi 方法 结 合 Chebyshev 半迭 代 算 法 ，数值 实验 验证 了它 的计 算速度 和精 度 的提 高 ; 针 对 逐 次超 松 弛 (SOR )迭代 法 最优 松 弛 因子 在 理 论上 很难 找 到 的的 问题 ，本 文用 非线 性 最 小 二 乘 法 进 行矩 形 区域 上 不 同剖 分 下 的数据 拟 合 ，得 到最优 松 弛 因子 。(n) 的 有 理 回归公式 ，并通 过 数值 实验验证 了它 的精确性 . 最后又对共扼梯度法 (CG 法)进行 了实验 ，将其 与前面 的直接法和各种迭代 法进 行 了数值 比较 ，表 明该 算 法 的高 效性 . 关键 词 : 偏 微 分 方 程 数值 解 ; 线 性 方 程 组 ; 直 接 法 ; 迭 代 法 ; 二 维 Poi SSon 方程 ; A b str a e t W e o ften u se th e fin ite d iffe ren ee m eth o d ，fin ite elem en t m eth o d an d th e fin ite v o lu m e m eth o d an d o th er m eth o d s in tl