基于Galerkin方法的Helmholtz方程有限差分格式-计算数学专业论文.docxVIP

独创性声明 本人 声 明 所呈 交 的 学位 论 文 是本 人 在 导师 指 导 下进 行 的 研 究工作 及 取得的 研 究成果 。 据我所 知 ，除了 文 中特别 加 以标注 和致谢 的 地方外 ， 论文中 不 包含其 他 人已经 发 表或撰 写 过的研 究成果 ， 也不包 含 为获得 电 子科技 大 学或其 它 教育机 构 的学位 或证书 而 使用过 的 材料。 与 我一同 工 作的同 志 对本研 究 所做的 任何贡献均已在论文中作了明确的说明并表示谢意。 签名： 日期： 年 月 日 论文使用授权 本学 位 论 文作 者 完 全了 解 电 子科 技 大 学有 关 保 留、 使 用 学 位论文 的 规定， 有 权保留 并 向国家 有 关部门 或 机构送 交 论文的 复印件 和 磁盘， 允 许论文 被 查阅和 借 阅。本 人 授权电 子 科技大 学 可 以 将 学 位 论 文 的 全 部 或 部 分 内 容 编 入 有 关 数 据 库 进 行 检 索，可 以 采用影 印 、缩印 或 扫描等 复 制手段 保 存、汇 编 学位论 文。 （保密的学位论文在解密后应遵守此规定） 签名： 导师签名： 日期： 年 月 日 摘 要 Helmholtz 方程是一类非常重要的椭圆型偏微分方程，在物理学电磁辐射、地 震学和声学等相关研究领域里有着十分广泛应用。大量文献给出了关于方程不同 边界的数值解方法，本文在之前研究的基础上，利用 Galerkin 变分原理和有限差 分结合的方法，推导出新的基于 Galerkin 的 Helmholtz 方程有限差分格式。 本文首先给出基于变分原理的差分格式的基本理论，为之后推导新的差分格 式准备充分的基础。然后针对一维的 Helmholtz 方程的三种不同形式，在离散区间 上的构造合适的基函数，导出其二阶精度的差分格式。之后在其基础上，推到四 阶精度的差分格式，并在第一类边界条件下，对一维情况的差分格式做收敛性分 析，并根据特征值分析了一维情况下的条件数。 对于二维 Helmholtz 方程，利用一维导出的结论，在离散后的区域上构造基函 数，利用拉格朗日插值方法逼近函数，通过数值积分计算导出相应的二维情况的 四阶精度的差分格式。讨论了 Dirichlet 边界条件下的差分格式，并分析其收敛性 和特征值。 根据边界条件选取数值例子进行试验，利用正弦变换算法，通过与之前的差 分格式比较，分别得到差分格式的误差和精度，分析导出线性方程组矩阵的条件 数，分析系统的稳定性，实验数据显示新的差分格式具有较好的稳定性，验证了 此差分格式的优越性。 关键词：Helmholtz 方程，Galerkin 原理，有限差分格式，插值方法，数值积分， 正弦变换算法 ABSTRACT Helmholtz equation is an important elliptic partial differential equation. It has the extremely widespread application in the fields of physics, electromagnetic radiation, seismology and acoustic. A large number of researches about this equation are given with different boundary. Under the previous studies, we derive a finite difference scheme based on cosine function of Helmholtz equation which combines Galerkin finite element with finite difference method in this paper. At first, we show the basic theory of the difference scheme based on variational principle as the basis for a new difference scheme. Then we construct an appropriate primary function in discrete intervals of one dimensional Helmholtz equation and then derive the difference scheme with numerical integration method and linear interpolation method. At the same ti