以下の行列Aの固有値に対して、最大実部を与えるkの値を求めたいです。
A=[[1-0.01*k^2,1],[-4,-k^2-2]]※1
ただし仮定としてAのトレースは負であるとします。
固有値をλとした場合、以下の固有方程式が得られます。
λ^2-tr(A)+det(A)=0※2
そのため、固有値は以下のようなkを含んだ式で表されます※3。
λ=0.5(tr(A)+sqrt(tr(A)^2-4*det(A)))
λの実部が最大化する場合は、det(A)が最小の場合であり、
det(A)
=(1-0.01*k^2)*(-k^2-2)+4
=0.01*k^4-0.98*k^2+2
より、k^2=49の場合にdet(A)が最小になり、λが最大化すると思うのですが合っていますか?
pythonで答えを確かめたところ、k=4の場合に固有値の実部が最大化したのですが、どこか間違っているのでしょうか。
どなたかご教授頂けると幸いです。
どうぞよろしくお願い致します。
※1 行列Aはpythonの書き方を参考にしており、例えば最初の[]内には1行目の要素が順番に入っています。
※2 tr: トレース、det: 行列式を表します。
※3 tr(A)が負であり、最大固有値に興味があるので、ルート(sqrt関数で表記)の前の符号は+に限定します。
No.3ベストアンサー
- 回答日時:
> A=[[1-0.01*k^2,1],[-4,-k^2-2]]
x = k^2
と考えたほうが楽だと思う。
> λの実部が最大化する場合は、det(A)が最小の場合
はおかしいな。
λ = (tr(A) + sqrt(D))/2
D = tr(A)^2 - 4det(A)
だから、実部Reλは
Reλ = if D>0 then λ else tr(A)/2
そしてtr(A)は定数じゃないのだから、
dλ/dk = 0 (あるいは dλ/dx = 0)
を真面目に解かなくちゃ。
No.2
- 回答日時:
A の固有方程式 λ^2 - tr(A)λ + det(A) = 0 の判別式は、
D = { tr(A) }^2 - 4 det(A)
= { (1-0.01k^2)+(-k^2-2) }^2 - 4{ (1-0.01k^2)(-k^2-2)-1(-4) }
= 0.9801k^4 + 5.94k^2 - 7
= 0.9801(k^2 + 3.0303)^2 - 16
≧ 0.9801(0 + 3.0303)^2 - 16
= -7.
D ≧ 0 のときは、固有値はふたつとも実数になって
Reλ = λ だが、
大きいほうの固有値は λ = { tr(A) + √(tr(A)^2 - 4det(A)) }/2.
-7 ≦ D < 0 のときは、固有値は共役複素数となって、
その実部は Reλ = tr(A)/2.
質問文中の答案は、この場合分けをしていない。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Java Java モンスターブリーダー 1 2023/02/05 09:44
- 数学 黄色マーカーの意味がよく分かりません。 X軸方向に4倍、Y軸方向に8倍に引き伸ばした らへんは何とな 3 2023/04/10 20:08
- 数学 行列(I-βG)の逆行列が存在することの証明について 1 2023/06/23 01:33
- 数学 連立微分方程式の解き方について 7 2022/12/16 13:39
- 物理学 スピン 行列表示 固有状態 測定値 1 2022/08/16 18:39
- Visual Basic(VBA) Excel(VBA) 特定の条件に該当する行の値、書式を同じセルにコピ&ペーストしたいです 1 2022/05/21 18:18
- Excel(エクセル) 列を挿入しても式の一部を固定するには 6 2023/05/31 21:57
- その他(プログラミング・Web制作) プログラミング pythonの問題について 2 2022/04/19 00:41
- 数学 線形代数 A= 2 -1 1 ( 0 0 2 ) 0 -1 3 の固有値と固有空間を使って、この行列 2 2023/02/03 12:39
- Excel(エクセル) エクセルのオートフィルについて質問があります。 2 2022/04/08 22:31
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
固有値の値について
-
3x3行列の固有値の求め方
-
行列の積の固有値
-
行列の2分の1乗の計算の仕方
-
A^3= E の問題。
-
固有ベクトルが複数の場合
-
3×3行列の固有値が重解をとる...
-
行列・対角化可能の条件は?
-
正定値行列は正則行列
-
正則でない行列の固有値・固有...
-
ジョルダン標準形について
-
主成分分析の途中。固有値の求...
-
固有値と固有ベクトル・重解を...
-
行列A^2からAを求める方法について
-
冪零行列の固有値は 0 のみであ...
-
エルミート行列の対角化の証明
-
行列の問題で、A^(-1)=PDP^(-1)...
-
2次形式 対角化
-
固有値が実数でない計算
-
f(x,y)=x^2-2xy^2+y^4-y^5 この...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報