![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?e8efa67)
No.1ベストアンサー
- 回答日時:
ううむ。
これだけじゃ回答しようがないと思うなあ。ガウス・ルジャンドルの数値積分というのは、f(x)を-1~1の区間で積分するときに、n次ルジャンドル関数の零点にあたるxでf(x)をサンプリングして重み付きの和を取るんでした。無論、積分区間内に特異点があったりしたら使えません。一般に積分範囲が x=a~b である場合には
x=((b-a)t+a+b)/2
と変数変換すれば、t=-1~1のtに関する積分になる。そしてdx/dt = (b-a)/2という因子を掛け算しておけば良いですね。n次のガウス・ルジャンドル法は、高々n次の多項式で近似できるf(x)を扱う場合に旨く行きます。
さて、ご質問は、おそらく積分範囲 x=-1~1に対してガウス・ルジャンドルの数値積分を使いたいけれど、次数を2にして、分点、すなわちサンプリングする点を±0.5だけにしたい、という注文です。たぶん、±0.5における被積分関数f(x)の値なら簡単に求められる、というのでしょう。
もちろん、適当な一次式ではない関数g(たとえば3次関数)を用いて
x=g(t)
という変数変換でx=±0.5をt=±0.57.... に移し同時にx=±1をt=±1に移す、ということ自体は簡単です。するとf(g(t))と
dx/dt = g'(t)
の積を被積分関数としてt=-1~1について積分することになります。この場合、被積分関数 f(g(t)) g'(t) がtの2次多項式で近似できるんでないと、2次のガウス・ルジャンドル法を使って精度が出るという保証はありません。
高精度の数値積分をやりたいと仰っている割に、f(x)が高々低次の多項式で近似してしまえるんだったら、何もガウス・ルジャンドル法に拘る必要はないんで、例えばニュートン・コーツ型の数値積分、すなわち分点を等間隔に取る方法でも十分じゃないの?と思うんですが、どうなんでしょうね。
或いは分点の数をもっと増やして良い、というのだったら、代わりに例えば-1~-0.5, -0.5~0.5, 0.5~1の3つの区間に分けてそれぞれ積分するのでも良い。被積分関数の傾きが急な部分でサンプリングを細かくしてやるというのも精度が出ますし、その代わりに適当な変数変換をして等間隔サンプリングしたり、ガウス・ルジャンドル法を使ったり…いろんな処方が考えられます。
ですから、「±0.5」と限定なさる理由をもう少し明確に補足して戴くか、具体的に被積分関数をupして戴かないと、ろくな回答にならないと思います。
回答が遅れました。謝罪します。
幾つかの提案を頂きましたが、
その中の変数変換を用いる方法が私が望んでいたものです。また、実際に適用し以下の結果を得ました。
変換関数は、こうなりました。
g(t)=at^3+bt,(a=0.200..., b=0.799...)
積分は、変換関数とGauss-Legendre公式(2次)により、こうなりました。
∫f(x)dx≒cf(0.5)+cf(-0.5), (c=3a(0.577...)^2+b)
これを用いてやってみます。
PS.
ニュートン・コーツと分点数増の方法でなく、変数変換の方法を採用した理由は、1区間の分割数を2以上にできないため、2分割でより良い精度の積分公式が必要なためです。詳しく説明できなくてごめんなさい。
非常に有益な回答ありがとうございます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・「I love you」 をかっこよく翻訳してみてください
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・昔のあなたへのアドバイス
- ・かっこよく答えてください!!
- ・あなたが好きな本屋さんを教えてください
- ・スタッフと宿泊客が全員斜め上を行くホテルのレビュー
- ・【大喜利】【投稿~8/27】 こんなガソリンスタンド二度と来るか!なぜそう思った?
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・【お題】動物のキャッチフレーズ
- ・【お題】甲子園での思い出の残し方
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・「それ、メッセージ花火でわざわざ伝えること?」
- ・自分用のお土産
- ・人生で一番お金がなかったとき
- ・一番好きなみそ汁の具材は?
- ・泣きながら食べたご飯の思い出
- ・ちょっと先の未来クイズ第1問
- ・ゴリラ向け動画サイト「ウホウホ動画」にありがちなこと
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
dx/dy や∂x/∂y の読み方について
-
偏微分について(熱力学)
-
シュレディンガー方程式は暗記...
-
ポアソン括弧式
-
分点座標が±0.5のGauss-Legendr...
-
電気磁気学の∇(ナブラ)の問題で...
-
二変数関数u(x,t)を偏微分した...
-
下の問題を教えてください! 振...
-
2s軌道の極大値についてまた...
-
関数 y=-2x^6logexの導関数が...
-
量子力学:不確定性原理
-
物理の直列回路の問題について...
-
水素原子の軌道角運動量につい...
-
ポテンシャルエネルギーと力
-
x^(2)-x+a(1-a)<0 因数分解の途...
-
高2の数学の対数関数です。 真...
-
年代と年台・・・どちらが正し...
-
「無限の一つ前の数字は何?」...
-
この極限を求める問題で対数を...
-
【数学】 lim x→a ↑これってど...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報