高阶高斯积分节点的高精度数值计算∫b-中国工程科学.PDFVIP

高阶高斯积分节点的高精度数值计算 １ ２ １ ３ 张庆礼 ，王晓梅 ，殷绍唐 ，江海河 （１畅 中国科学院安徽光学精密机械研究所，合肥 ２３００３１ ；２畅 中国科学院环境光学与技术重点实验室，合肥 ２３００３１ ； ３畅 中国科学院合肥物质科学研究院，合肥 ２３００３１） ［摘要］　 在工程数值计算、X射线衍射线形分析、光谱学等领域常使用高斯数值积分，高斯积分的节点及权重 因子是数值积分的必须数据。 研究了高次勒让德、拉盖尔和厄米多项式的零点，即高斯－ 勒让德、高斯－ 拉 盖尔、高斯－ 厄米积分的节点的计算方法，给出了一种有效的高精度数值算法——— 搜索迭代方法（scan唱iteration method，SIM）。 根据勒让德、拉盖尔、厄米多项式的特点，对拉盖尔多项式、厄米多项式的定义稍做变化后，获 得了计算多项式值的稳定递推关系。 求它们的根时，先在一定范围内以一定的步长搜索根所在的区间，获得 所有根的各自区间范围后，再通过常用的迭代方法如割线法、二分法进行求解。 数值实验表明，这种方法是 非常有效的，可获得高次勒让德、拉盖尔、厄米多项式的全部高精度根值。 ［关键词］　 高斯积分；勒让德多项式；拉盖尔多项式；厄米多项式；求根 ［中图分类号］　 O１７４畅６；O２４１　 ［文献标识码］　 A　 ［文章编号］　１００９－１７４２（２００８）０２－００３５－０６ 式（１）对次数≤ ２ n－ １ 的多项式精确成立，即为高斯 1　 前言 积分公式，它的代数精度为 ２ n－ １ 。 在所有的 n 个 　 　 高斯数值积分是一种有效的高精度数值计算方 节点积分公式中，高斯积分有最高的代数精度。 法。 在工程数值计算、固体物理学、光电子学等领域 可以证明，积分公式（１）为高斯积分的充要条件 ［２］ 常用来对难以获得解析解的积分进行数值计算。 例 是积分节点为正交多项式的根 。 不同的高斯积分 ［１］ 公式如高斯－ 勒让德、高斯－ 拉盖尔、高斯－ 厄米、 如，在工程电磁场计算中常用的贝塞尔函数 ，如果 采用其级数解来计算，需考虑很复杂的递推稳定性 高斯－ 切比雪夫积分的积分区间和正权重函数不 问题；反之，若从它的积分形式出发来计算，采用数 同［３］ 。 高斯－ 勒让德、高斯－ 拉盖尔、高斯－ 厄米、 个节点的高斯积分，即可获得满意的数值精度。 在 高斯－ 切比雪夫的正权重函数和积分区间分别为 １ ２ X 射线衍射中的线形分析和光谱学中，Voigt 函数能 和［ － １ ，１］ ，e－ x 和［０ ，∞ ） ，e－ x 和 （ － ∞ ，＋ ∞ ） ，（１ － 较好地描述X 射线衍射峰和光谱峰的线形，但该函 x２ ）－ １桙２ 和［ － １ ，１］ 。 数没有解析的积分解；固体物理学中热容表达式中 很多文献可以查到高斯积分公式的节点 x 和