ガウス求積法(求積法選択)

[1-4] /4件

表示件数

ご意見・ご感想: ガウス求積法に関し、n次式が 1/(2n-1) の積分点で済む理由を明快に説明しているのを、これまで見たことがありません。ウイキペディアでも ” f(x) が 2n − 1 次以下の多項式であれば上記の式が厳密に成立する ”と解説されていますが、理由は明快には説明していないようです。私なら、1次元の局所原点を中点に置けば、奇関数項の積分は左右の象限で相殺してゼロになるから、偶関数項のみの積分で済む、がその理由である。よって、2Dでは約1/4に、3Dでは約1/8の計算量で済む、と説明します。なお、積分はべき級数のみ可能で、他に方法はない。よって、指数関数でも何でも、すべてべき級数展開して積分する。Taylor展開して積分すると理解すればよい。非線形の代表はステップ関数ですが、近似的に高次積分するのが現実で、それより良い方法は無いようだ。（あればご教示願いたい。）厳密解は手計算（あるいは領域区分してガウス求積）できるので、精度検証も可能。なお、X-FEMと呼ぶ技法が在るが、その真髄には上述で辿れる、筈。以上

使用目的: 数値積分プログラム作成の為
バグの報告: 以下の部分が間違っています。
(4) Laguerre [0,∞)
正しくは、
(4) Laguerre [0,∞]
です。
keisanより: Laguerre-Gauss求積法の積分区間は、[0,∞)となります。
参考リンク
Weisstein, Eric W. "Laguerre-Gauss Quadrature." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/Laguerre-GaussQuadrature.html

使用目的: 　日頃、数値処理をして開発の検討を行う業務についているので、数値処理法の勉強に役立てております。
ご意見・ご感想: 　ガウス・クロンロッド求積法は、関数電卓に採用されていると存じますが、是非、その内容を理解したく、評価点Xiと重みWiの求め方をご開示頂けないでしょうか。　ガウス・ルジャンドルの評価点ｎにｎ＋１個の評価点を追加されていることまでは分かったのですが、ｎ＋１追加の評価点の求め方やそれぞれの重みの求め方が分かりませんでした。　調査の仕方が悪いのでしょうか、ガウス・クロンロッド法の解説文を見つけることができませんでした。　以前の関数電卓では、シンプソンの積分公式を採用されていたと存じますが、御社関数電卓（プログラム電卓）愛用者として、この求積法にかえられてから速く・精度よくなったと感謝しております。

積分区間（a,b）の積分を様々なガウス求積法で計算します。
(ガウス‐ルジャンドル積分、第1種チェビシェフ積分、第2種チェビシェフ積分、ラゲール積分、エルミート積分、ヤコビ積分、ロバート積分、クロンロッド積分)