[問題]解矩陣正定已回收

看板MATLAB作者 (好聽木琴)時間14年前 (2011/08/27 12:18), 編輯推噓3(3推 0噓 9→)
留言12則, 3人參與, 最新討論串1/1
假設我有個6x6 矩陣。 這邊方便書寫 我就寫2x2 說是[1,2;3,4]好了... 左上角那個位子 減個t 所以變成 [1-t,2;3,4] 想要求最大的t 使得這個矩陣 還是正定 (positive definite) 有指令 可以直接解這個t嗎? solve好像 沒有解這種的吼? 有別的嬤? 還有我也嘗試 硬幹...就是把矩陣 eigenvalue解出來... 他的e.value會是t變數的形式.. 然後我在解這些eigenvalue >0 去解t.. 可是今天他是6x6 好像有所謂的Galois theory問題... 我打eig(A) 他解不出來... 至少我可以打poly(A) = A的特徵多項式 然後再用solve指令... 可是也是解不出來...請問有什麼好方法 解出這個最大的t使得這矩陣還是正定呢? 感恩... -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 68.48.173.107
→
08/28 09:14, , 1F
這不就是lmi嗎?
08/28 09:14, 1F
→
08/28 11:47, , 2F
google了matlab lmi 沒看懂啊..
08/28 11:47, 2F
→
08/28 11:47, , 3F
可以請C大指點像我的問題 如何用lmi解
08/28 11:47, 3F
推
08/30 02:11, , 4F
太久沒碰程式了
08/30 02:11, 4F
推
08/30 02:15, , 5F
一元六次沒有解析解...有點困難
08/30 02:15, 5F
→
08/30 02:16, , 6F
所有根要大於零 印象中自動控制有這個課題
08/30 02:16, 6F
→
08/30 02:17, , 7F
以及相對應的方法 是Routh table嗎?
08/30 02:17, 7F
→
08/30 02:18, , 8F
不太熟 拍謝
08/30 02:18, 8F
→
08/30 08:08, , 9F
還是感謝樓上的回應指教..
08/30 08:08, 9F
推
08/30 17:39, , 10F
我印象中正定矩陣是不是要行列是直大於零
08/30 17:39, 10F
→
08/30 17:40, , 11F
由左上到右下 1*1 2*2 3*3~ 6*6行列式值大於零
08/30 17:40, 11F
→
08/30 17:40, , 12F
可以翻翻線性代數可能會比較清楚
08/30 17:40, 12F
‣ 更多 snaredrum 的文章
文章代碼(AID): #1EM70Vhf (MATLAB)