「統計学関連なんでもあり」の過去ログ---003

★ General Polynomial Fit　について教えてください。 ★

　214　General Polynomial Fit　について教えてください。　　松井　秀子　　1999/10/27 (水) 19:49
　　224　Re: General Polynomial Fit　について教えてください。　　堀　啓造　　1999/10/28 (木) 18:54
　　215　Re: General Polynomial Fit　について教えてください。　　マンボウ　　1999/10/27 (水) 23:44

214.　General Polynomial Fit　について教えてください。　　松井　秀子　　1999/10/27 (水) 19:49

プログラムを作っているんですが　General Polynomial Fit　を使わなくてはなりません。そこで教えていただきたいのですが，英語で以下のように説明されています。
そのアルゴリズムに定数で:
Mathematics Functions
polynomial order must be greater than or equal to zero: 0 」 m < n-1, where n is the number of sample points, and m is the polynomial order. If polynomial order is less than zero, the VI sets Polynomial Fit Coefficients to an empty array and returns an error. polynomial order defaults to 2.

algorithm has six selections:
0: SVD
1: Givens
2: Givens2
3: Householder
4: LU decomposition
5: CholeskyIt defaults to SVD (0).

ここに書かれている6種類のアルゴリズムの定義についてご存知であれば教えていただけませんでしょうか。

申し訳ありませんが　どなたかご教授していただければ幸いです。
宜しくお願いいたします。

　　　　　[このページのトップへ]

224.　Re: General Polynomial Fit　について教えてください。　　堀　啓造　　1999/10/28 (木) 18:54

> 0: SVD

Singular value decomposition 特異値分解
これは行列が特異の場合（正則でない場合を含む）でも固有値・固有ベクトルに分解します。

これは数値計算法に紹介されていますし，インターネットでも手に入れることができます。
http://www.netlib.org/svdpack/index.html
にたくさんのプログラムがあります。

> 5: CholeskyIt defaults to SVD (0).

コレスキー分解のことかな。

　　　　　[このページのトップへ]

215.　Re: General Polynomial Fit　について教えてください。　　マンボウ　　1999/10/27 (水) 23:44

> ここに書かれている6種類のアルゴリズムの定義についてご存知であれば教えていただけませんでしょうか。

これは統計学というよりも数値計算の分野ですね。
連立方程式の解を求める方法（およびそれに関連する方法;行列の対角化とか三角化）ですね。
奥村晴彦氏の，「C言語によるアルゴリズム事典」等も参考になるでしょう。全部が出ていたかどうかは分かりませんが。
ソースなどはインターネット上に公開されています。以下を参照してみて下さい。
http://www.matsusaka-u.ac.jp/~okumura/algo/

ハードディスクの肥やしになっているファイルを調べたら，
house.c -- Householder (ハウスホルダー) 変換
lu.c -- LU分解
はあるみたいですね。

その他に，行列計算ライブラリみたいな本をさがすと良いかも。

#しかし，なにも，これら6つの手法を用意する必要なんかないと思うんですが。
#お仕事関係で，仕様書に定められているからということなんでしょうか?
#お仕事関係なら，このような基礎的なアルゴリズムは知っておかねばならないだろうし，そうでなくても，少なくとも自分で調べることができなきゃ...

#一番簡単には，掃き出し法とかでいいんじゃないだろうか。

　　　　　[このページのトップへ]