アフィリエイト広告を利用しています

2015年12月04日

Excel VBA ルジャンドル多項式 Legendre polynomial

【Excelショートカット豆知識】
 [Shift] + [Space] 行を選択する
 [Ctrl] + [Space] 列を選択する

@ルジャンドル多項式 Legendre polynomial


 ルジャンドル微分方程式

ルジャンドル微分方程式

 の解がルジャンドル多項式であり、 [−1, 1] の範囲で

ルジャンドル多項式

と表されます。ルジャンドル多項式は次の漸化式を満たします:

    (2n + 1)xPn(x) = (n + 1)Pn + 1(x) − nPn - 1(x)

 具体的には次のような形になります:

    P0(x) = 1
    P1(x) = x
    P2(x) = (3x2 − 1) / 2
    P3(x) = (5x3 − 3x) / 2
    ・・・・・・・・・

AExcel でルジャンドル多項式のグラフを描きます


 Excel でグラフを描いてみると・・・・・・

 ルジャンドル多項式.png

 この図からわかるように、一般に

    Pn(1) = 1

を満たします。P0(x) において n を大きくしてゆくと、[−1, 1] の区間にたくさんの波が閉じ込められている様子がよくわかります。たとえば P50(x) のグラフを描いてみると・・・・・・

 P50(x).png

 このように区間を限定して 各項の係数をうまく調整すれば、その内部ではあたかも cosx や sinx のように振る舞う疑似的な周期関数を作ることができます。仮にルジャンドル多項式の定義域を拡大したとしても、[−1, 1] を超えた場所に波は存在しないので、単調増加あるいは単調減少となり、そのような関数はほとんど用途がありません。任意の関数を展開するためにも波形の存在する区間に限っておく必要があります。別記事にある正規化されたエルミート多項式 u(x) は [−∞, ∞] で定義されていましたが、非常に収束の早い関数ですから、グラフを見ても実質上は −4 〜 4 のあたりに存在する波であると考えることができます。

Bルジャンドル多項式の直交性


 ルジャンドル多項式の直交性は、

ルジャンドル多項式の直交性

によって表されます。また [−1, 1] で値をもつ任意の関数はその範囲内で

ルジャンドル多項式の展開

というように展開でき、その係数 cn

ルジャンドル多項式の係数

で与えられます。

ルジャンドル多項式の応用

 変数 x を cosx 或いは sinx に置き換えると、−1 ≦ cosx, sinx ≦ 1 ですから、実数全域で定義される使いやすい関数となります。Pn(cosx) のグラフを描いてみると:

 P(cosx).png

 小さな波と大きな波が交互に現れる関数です。電磁気学で多重極展開(電荷の作る電位を求める表式)と呼ばれる手法にはこの形が用いられます。
⇒ ルジャンドル多項式を計算する Function マクロ(ユーザー定義関数)  
posted by Blog Cat at 01:29 | Comment(0) | TrackBack(0) | n次関数
検索
Excel VBA 数学教室
数学問題集(解答付き)
下剋上算数
ベクトル解析
サッカーマティクス
Excelで学ぶ統計解析
和算的思考力
学び直し
整数論の理論と演習
大人が手こずる算数
東大生の知恵袋
フーリエ変換
インド式秒算術
Excelで学ぶ微分積分
Excel 数学シミュレーション
オイラーの贈物


ファン
最新記事
カテゴリーアーカイブ