質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.35%
Mathematica

Mathematicaは、ウルフラム・リサーチによって開発されている数式処理システムです。

Q&A

0回答

834閲覧

固有値がEigenvalues関数を使用すると求まらないのに、定義に従った求め方だと求めることができる理由について

MF0524

総合スコア50

Mathematica

Mathematicaは、ウルフラム・リサーチによって開発されている数式処理システムです。

0グッド

0クリップ

投稿2021/10/26 08:54

###発生している問題
mathematicaを使用して、次のような変数mを持つ4x4の行列の固有値を求めようとしています。

mathematica

1testmat = {{m, 0, 0, 1}, {1, m, 0, 0}, {0, 1, m, 0}, {0, 0, 1, 1 + m}}

この行列に対して、Eigenvalues関数を用いて固有値を計算した場合、固有値は計算されず、以下のような解けなかった固有方程式が出力されます。

mathematica

1>{Root[-1 + m^3 + m^4 + (-3 m^2 - 4 m^3) #1 + (3 m + 6 m^2) #1^2 + (-1 - 4 m) #1^3 + #1^4 &, 1], 2 Root[-1 + m^3 + m^4 + (-3 m^2 - 4 m^3) #1 + (3 m + 6 m^2) #1^2 + (-1 - 4 m) #1^3 + #1^4 &, 2], 3 Root[-1 + m^3 + m^4 + (-3 m^2 - 4 m^3) #1 + (3 m + 6 m^2) #1^2 + (-1 - 4 m) #1^3 + #1^4 &, 3], 4 Root[-1 + m^3 + m^4 + (-3 m^2 - 4 m^3) #1 + (3 m + 6 m^2) #1^2 + (-1 - 4 m) #1^3 + #1^4 &, 4]} 5

上記の固有方程式は変数について4次の方程式であるため、解が求まると思うのですが、なぜ求まらないのでしょうか?

###試したこと
固有値は以下のようにEigenvalues関数を使用せず固有値を定義に従って求めることもできるはずで、こちらの場合は解が求まりました。

mathematica

1xmat = {{x, 0, 0, 0}, {0, x, 0, 0}, {0, 0, x, 0}, {0, 0, 0, x}}; 2eigentest = testmat - xmat; 3Solve[Det[eigentest] == 0, x] 4

ちなみにこの方法の場合に解く対象の行列式は、以下になり、Eigenvalues関数で出力された式と一致します。

mathematica

1Collect[Expand[Det[eigentest]], x] 2>-1 + m^3 + m^4 + (-3 m^2 - 4 m^3) x + (3 m + 6 m^2) x^2 + (-1 - 4 m) x^3 + x^4

Eigenvalues関数で解が求まらない理由は何かあるのでしょうか?
どなたかわかる方は教えて頂けると幸いです。
どうぞよろしくお願い致します。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだ回答がついていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.35%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問