- ベストアンサー
対称行列の固有値
次の対称行列 1 -3 2 -3 2 1 2 1 3 の固有値xを求めたいのですが、 固有方程式が x^3-6x^2-3x+42=0 となり因数分解ができません。 対称行列の固有値は必ず実数となる らしいので私の計算間違いかと思いますが、 何度計算しても同じ固有方程式になってしまいます(T T) どなた様か御指導願います。
- 数学・算数
- 回答数5
- ありがとう数8
- みんなの回答 (5)
- 専門家の回答
質問者が選んだベストアンサー
これは実対称行列ですから、固有値は3つの実数です。固有方程式の判別式は正ですので、3つの異なる実根になります。Mathematicaはカルダノの公式を使っています。これは不還元の場合になりますので、外見は複素数ですが、実際は実数です。ですから、Mathematicaの解もNo2さんの解も正しい(等しい)のです。このことを確かめるには、簡単な3実根をもつ3次方程式をカルダノの公式にあてはめてみれば良いでしょう。
その他の回答 (4)
- Tacosan
- ベストアンサー率23% (3656/15482)
#3 の回答に対する補足についてですが, 「代数方程式の解を求めるため, その代数方程式を固有方程式として持つ行列の固有値を求める」ということがあったそうです.
- ojisan7
- ベストアンサー率47% (489/1029)
5次以上の代数方程式は数値計算で近似値を求めます。4次以下の代数方程式も根の公式を使うより、数値計算した方が実用的ですよね。数値計算の方法はいくつかありますので、数値解析の教科書等で調べてください。サイトにも紹介されていると思います。「高次方程式の数値計算」をキーワードにして検索してください。
- colder
- ベストアンサー率43% (30/69)
固有方程式はあっています。 実数≠有理数 ということは理解していますか 因数分解は常にできるとは限りません。 固有値は、 4.80560、3.61323、-2.41883 付近です。
お礼
ご回答ありがとうございます。
補足
上の問題はある大学院の過去問題です。 自力で解くことができず、さらに解答がないため、 Mathematica5.1に Solve[x^3 - 6 x^2 - 3 x + 42 == 0, x] と入力し、解かせてみました。 すると複素数の解を出力しました。 3つの複素解は複雑ですが規則性があるため、 人手で解く方法があるのではないかと 調べていました。 colder様が示された数値計算による 解を固有値に採用すると、確かに 元の固有方程式と似た係数をもつ方程式 になるのですが、 Mathematicaが出力した複素解と 対称行列の固有値は必ず実数となる という性質が矛盾しているように 思えたので OKWaveに投稿しました。 Mathematicaの扱いもよく知らないので、 間違っているかもしれません。 どなたか人手で解く方法を教えてくださりませんか。
- Willyt
- ベストアンサー率25% (2858/11131)
固有方程式は下記です。 (1-λ)(2-λ)(3-λ)-4(2-λ)-9(3-λ)-(1-λ)=0 となります。後は腕力の問題です(^_-) 定数項は30になりますね。行列式の計算で右上がりのときにはマイナスをつけなければいけないのを間違えていませんか?
補足
行列式の性質を用いたり、 サラスの方法を用いても 固有方程式は、 x^3-6x^2-3x+42=0 となってしまいました。
関連するQ&A
- 対称行列の固有ベクトル
対称行列の固有ベクトルは互いに垂直という性質がありますが、 固有ベクトル AX1=λ1 X1、 AX2=λ2 X2 の式から n次の対称行列Aは次のように書き表すことができます A= λ1 X1 X1^t +λ2 X2 X2^t+ ・・・ +λn Xn Xn^t なぜ固有ベクトルの式から対称行列の式が表すことができるのでしょうか? 証明を教えてください。よろしくお願いします。
- 締切済み
- 数学・算数
- 対称行列 固有値 実根
n次実対称行列の固有方程式の解(固有値)がn個の実数となるのは、なぜでしょうか。 参考書など見てみましたが、記載されていないようで困っています。 簡単な証明や理論など、どなたかご存知の方がいらっしゃれば教えていただきたいです。
- 締切済み
- 数学・算数
- 3×3行列の固有値と固有ベクトル
以下の行列Aの固有ベクトルを求めようとしているのですが,解を見つけられないでいます. 2 1 0 1 2 0 0 0 -2 計算を進めた結果,固有値λは3,1,-2となり,λ=3,1に対応する固有ベクトルはそれぞれ[1,1,0]t,[1,-1,0]tとなったのですが,λ=-2の場合で求めた固有ベクトル[1,1,k]t(kは任意の実数)がAx=λxに対応しない値になってしまいます.私の計算に何か問題があるのでしょうか? また,行列Aは対称行列なのでそれぞれの固有ベクトルの内積は0になると思うのですが,固有ベクトルの値が得られないことと何か関係があるのでしょうか? 回答よろしくお願いします.
- ベストアンサー
- 数学・算数
- 行列の対角化 固有値を求める
次の行列の固有値、固有ベクトルの作る行列Pを求めて、対角行列に変換せよ。 A= 7 4 -16 -6 1 12 2 2 -5 と言う問題で、 固有値を求めるとき、|A-λE|より (7-λ) 4 -16 -6 (1-λ) 12 2 2 (-5-λ) となって =(7-λ)(1-λ)(-5-λ)+(-6)*2*(-16)+2*4*12-・・・・ としてから展開すると、計算も大変で、そのあとの 因数分解もわかりません;; どうすれば、もっと簡単に固有値を求められるでしょうか? お願いします。
- 締切済み
- 数学・算数
- 行列の固有ベクトルの問題
3X3行列T (0 1 0) (0 0 1) (-2 1 2) がある。 行列Tの固有ベクトルを求めよ。ただし、各固有ベクトルは、第1成分を1とするものを求めよ。 次に 3X1 行列(3) (2) (6) を求めたTの固有ベクトルの線形和の形で表せ。 という問題です。 私は固有ベクトルが (1) (1) (1) (1) (-1) (0.5) (1) (1) (1) と求まりましたが、これは (3) (2) (6) を表せないのではないかと思いますが。 どこか計算ミスをしたのですか? ご指導を待っております!お願いします!
- ベストアンサー
- 数学・算数
- 行列、固有値
2×2の行列Z=aI+bJを考える。 ただし、I=|1 0|および、J=|0 1|である。 |0 1| |-1 0| また、aとbは実数であり、aとbが同時に0になることはないとする。 (1)Z[t]を計算し、a,b,I,Jで表せ。ここで添え字[t]は行列の転置を表す。 これはZ=(a b)Z[t]=(a -b)であるので (-b a) (b a) Z[t]=aI-bJでよいでしょうか?? (2)Zの逆行列を求め、a,b,I,Jで表しなさい。 Z^-1=1/(a^2+b^2)(a -b) なので (b a) 1/(a^2+b^2)*(aI-bJ)でよいでしょうか?? (3)JとZの固有値を求めなさい。 J=(0 1)なので|A-xE|=0より|-x 1| =0 (-1 0) |-1 -x| 従ってx^2+1=0 x^2=-1 x=±i なので Jの固有値は i と -i でよいでしょうか?? Z=(a b) |a-x b |=0 (a-x)^2+b^2=0 ここから、因数分解できないためもとめることができません。 (-b a) |-b a-x| 長々と書いてしまいすみません。よろしくお願いします。
- 締切済み
- 数学・算数
お礼
ご回答、誠にありがとうございます。
補足
>Mathematicaはカルダノの公式を使っています それは存じませんでした。確かにカルダノの公式だと 3次正方行列の固有方程式を解くことができますね。 忘れていました。ありがとうございます。 ojisan7様のご回答を受けて思ったのですが、 5次以上になってくると方程式を 代数的に解くことができないと聞いたことが あります。そのような場合はどのように 求めるのでしょうか。求められる場合の 解き方をご存じの方いますか? 厚かましい追加の質問をお許しください。