- ベストアンサー
ヤコビ法とQR法について
stomachmanの回答
- stomachman
- ベストアンサー率57% (1014/1775)
大きい固有値だけに限って求められれば十分、という応用が沢山あります。また、制御工学では不安定な固有値成分だけ取り出したい。そういう場合どうするのか、ちょっと調べてみたら 「ランチョス法に、等角写像を組み合わせる。」 というテクニックがあるみたいです。 等角写像 y = (x+c) / (x-c) （x,yは複素数、cは適当な正の定数） によって、不安定なやつを単位円の外に、安定な奴を単位円の中に写像しておいて、ランチョス法（絶対値が大きい固有値から順に出てきます）を適用するらしい。 まだ調査中ですが、お急ぎのようなのでとりあえず。
関連するQ&A
- 固有値問題の解法について
固有値問題を解くプログラムを作りたいと考えています。 対象とする行列は「実対称行列（正方行列）」で、251次の 行列を考えており、この行列の固有値、固有ベクトルを 求めたいです。そこで、現在調べたところ、ヤコビ法、 QR法などがあるのですが、どの方法がより高速に求められる でしょうか。ご教示願いたいと思います。また、上記の方法 よりも高速に求められる解法などありましたら、教えて頂き たいです。宜しくお願いします。
- ベストアンサー
- 数学・算数
- QR法を使った行列の固有値問題について
行列の固有値と固有ベクトルを求めるプログラムをJavaで 作っています。そこで、アルゴリズムとしてQR法を用いて いるのですが、この方法では固有値のみを求めており、 固有ベクトルは求めていないことがわかりました。そこで、 QR法を用いつつ、固有ベクトルも同時に求める方法がありましたら、 教えて頂きたいと思います。宜しくお願いします。
- 締切済み
- 数学・算数
- ＱＲ法について
ＱＲ法により，固有値を求めるプログラムを作っています． ＱＲ法でqjを求めるときに，rjjの逆数が必要ですが，rjj＝０になったとき，うまく求めることができません． このような場合に，どうしたら解決できますか？ ちなみに言語はＣ言語ですが，解決法さえいただければ自力で組むつもりです．「教育」の「数学」のカテゴリか迷いましたが・・・・
- 締切済み
- その他（プログラミング・開発）
- scilab ヤコビ法について、errorがでます
ヤコビ法で固有値、固有ベクトルを求めるプログラムを色々調べながら作っていたのですが、errorが出てしまい手詰まりになりました。あまり、scilab自体詳しくなく、多々間違いがあるかもしれませんが、学校でどうしても必要なため、どうか、解答の方をよろしくお願いします。 //初期化 A=[2 1 3;1 0 1;3 1 0] Pinf=eye(3,3) //ヤコビ法：Aの対角成分が固有値、 Pinfの行ベクトルが固有ベクトルになる for i=1:21, i if abs(A(1,2))>abs(A(1,3)) then theta=atan(2*A(1,2))/(A(1,2)-A(1,3))/2 P=[cos(theta),sin(theta);-sin(theta),cos(theta)]; elseif abs(A(1,3))>abs(A(2,3)) then theta=atan(2*A(1,3))/(A(1,3)-A(2,3))/2 P=[cos(theta),sin(theta);-sin(theta),cos(theta)]; else theta=atan(2*A(2,3)/(A(2,3)-A(1,2)))/2 P=[cos(theta),sin(theta);-sin(theta),cos(theta)]; end Pinf=P*Pinf A=P*A*P' ,end 現在、下から4行目のendを入れなければ、end または 「else が抜けています...」とerrorが出てしまい、また、endを入れると、その下の Pinf=P*Pinf の所で、「一貫性がない掛け算です。」とerrorがでます。 以上が、現在困っている状況になります。 よろしくおねがいします。
- ベストアンサー
- その他（プログラミング・開発）
- 大学数学でのヤコビ行列について
大学数学でのヤコビ行列について 現在ヤコビ行列を求める以下の問題をやっています． 問題 次の写像のヤコビ行列を求めなさい． f:R^2->R^3 ※２次元から３次元への写像 f×(x,y) = xe^y + cosy x x + e^y ※(x,y):列ベクトル 自分なりに解いてみたところ，右辺の各成分をxとyで偏微分して， f = e^y ,xe^y-sin y 1 , 0 1 , e^y と求めました． このやり方と答えで合っているのでしょうか？ 参考にしたURLを載せておきます． http://www.ne.jp/asahi/search-center/internationalrelation/mathWeb/Differentiation/PartialDifferential2VarFnctn/GradHessian.htm どなたかご教授お願いします．
- 締切済み
- 数学・算数
- ニュートン法に関して
数値計算初心者です。数値計算で分からないことがあるので質問します。よろしくお願いします。 y=f(x,a)という関数があってパラメータaを非線形最小２乗法のニュートン法やマルカート法を使って求めたいのですが、計算過程でf(x)を各パラメータで偏微分してヤコビ行列を求める必要があると思われます。 例えばf(x)が複雑な関数で偏微分するのに困難な関数であった場合、 偏微分をしなくてΔxを決定するにはどのような方法があるのでしょうか？
- 締切済み
- 数学・算数
- 一般化逆行列と最小二乗法
最小二乗法は割と簡単に理解することができますし、式の誘導も簡単ですが、分数が出てきたら分母がゼロでないとか、逆行列が存在するとか理想的な条件を仮定しているように思います。そこでその理想的な条件が存在しない場合、すなわち逆行列が存在しない場合、”一般化逆行列を用いて計算する”とサラリと書いてある本がありました。データ解析ソフトRなどもそれに対応しているかもしれません。一般化逆行列というのはすんなり受け入れられるものでしょうか。何か別の指標があってそれを最小化するとか何らかのペナルティとか損失を甘受した上で計算していると思うのですが、いきなりピンチヒッターとして出てくることができるみたいに書いてありました。数理統計の本には共線性がある場合とか行列式が極めて小さな値になるとかの場合に出てくるようです。少し読んでみると固有値・固有ベクトル（正規直交行列を構成）で行列を展開したもののような記述もあり、これはこれで普通のことのように思うのですが。一般化逆行列とはどのようなものだと思えばいいでしょうか。
- 締切済み
- 数学・算数