推 alamabarry:polyeig 01/03 04:09
→ alamabarry:不然簡單一點eig(a*inv(b)) 01/03 04:11
想請教各位，關於特徵值的算法，以工數上學到的是
det( a-mI )=0
行列式求得m之兩根為特徵值。
在matlab內使用 eig(a) 函數計算。
但如果我今天要求的是
det( a-mb )=0
b=[x,0;0,y] 這種非單位矩陣的矩陣
該如何求解？
我原本是使用 eig(a) 後，再將特徵值除上 b
但發現結果是錯的，
我現在能想到的辦法是使用符號運算去求解，
但實在麻煩，請問版友們有沒有更方便的解法？
謝謝。
===============自問自答耍憨=============
剛剛無意間在大陸網站查到 eig(a,b) 這用法
廣義的特徵值解法，其解為正確，謝謝各位。
========================================
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.134.241.90
※ 編輯: magaiml 來自: 220.134.241.90 (01/02 22:30)