[討論] 取出 Jordan blocks

[math99 (新世界)]

我有一個 Jordan form , 我想要取出各別的 Jordan blocks
請教 matlab 有什麼方便的做法嗎？

目前我做法只能用 for loop 一個一個取, 希望能夠不要用到 loop

例

Input:
A = [2 0 0 0;
     0 1 1 0;
     0 0 1 0;
     0 0 0 2]

Output:

J1 = [2], J2=[1,1;0,1], J3=[2];

感謝
※ 編輯: math99 (140.115.25.108), 08/24/2015 15:11:56
推 sunev: 如果你已經知道各個block的大小，可參考這篇
→ sunev: http://0rz.tw/3i4pY 下面前三行，[2 3 5 1]為各個block之
→ sunev: 大小
→ math99: 感謝 block size 我是用 diag(A,1) 判斷 應該是可以