clc
clear
%A=round（255*rand（5））;%創建一個在0~255取值的5階整數隨機矩陣
%?A=100*rand（5）;%創建一個在0~100取值的實數矩陣
A=[2-1-1;2-1-2;-112];%用一道例題檢驗
L0=eig（A）;%求矩陣A的特征值
L=single（L0）;%降低精度方便unique函數判斷相等關系
tzz=unique（L（:））;
J=[];
for?i=1:length（tzz）
????num=length（find（L==tzz（i）））;%第i個特征值的重復次數，即特征值的階數
????R=rank（eye（size（A））*tzz（i）-A）;
????[a~]=size（A）;
????x=a-R;%第i個特征值有x個線性無關的特征向量
????jie=num-x+1;%第i個特征值對應的最高階Jordan塊的階數
????%將第i個特征值對應的最高階Jordan塊寫入結果矩陣
????for?j=1:jie
????????[mn]=size（J）;
????????J（m+1n+1）?=?tzz（i）;
????????if?j>1
????????????J（m+1n）?=?1;
????????end
????end
????%將第i個特征值對應的剩余Jordan塊依次寫入結果矩陣
????if?x>1
????????[mn]=size（J）;
????????for?k=1:（x-1）
????????????J（m+kn+k）?=?tzz（i）;
????????end
????end
end
A
%將結果與matlab內置函數結果進行對照
J
J2?=?jordan（A）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????940??2015-12-29?16:55??Jordan001.m
?????文件?????????575??2015-11-10?22:32??LUfenjie.m
?????文件?????????477??2015-11-10?21:08??QRfenjie.m