(Library) Legendre陪関数計算ライブラリ

(C) Takeshi Horinouchi

LICENSE: Ruby's

概要

RubyでLegendre陪関数(球面調和関数の緯度依存部分)を計算するライブラリです. 通常の定義と違いノルムが2に正規化されています. P_n^m について,まず P_m^m, P_m+1^m を漸化式で求め, ついで各mについて,nの漸化式を計算します.

ダウンロード

assoc_legendre.rb ライブラリ. テストプログラムがついています.テストはノルムが2に なる確認(ただしガウス緯度でしてない簡略版)と, 簡単な図示.

ynm_global_plot.rb 描画アプリケーション.経度について三角関数をかけた球面調和関数を図示します.

依存ライブラリ

NArray および RubyDCL (電脳Rubyホームページを参照のこと)

使用法と描画結果

ライブラリファイル assoc_legendre.rb を,ruby コマンドの直接の引数として実行すると,テストが実行されます. 実行は,

% ruby assoc_legendre.rb

テストでは標準出力にいろいろ垂れ流されるほか,次の図が表示されます.

assoc_legendre.png

なお,

% ruby assoc_legendre.rb 4

とすると,PNG画像が出力されます.

また,次を実行すると

% ruby ynm_global_plot.rb

次の図が表示されます.

ynm_global_plot.png

こちらも引数4でPNG出力です.

参考

ノルムを2に正規化したLegendre陪関数について:

石岡圭一,「スペクトル法による数値計算入門」, 東京大学出版会, 2004. 
更新日時:2008/12/02 18:41:38
キーワード:[特殊関数]
参照: