disp（‘請輸入判斷矩陣A（n階）‘）;
A=input（‘A=‘）;
[nn]=size（A）;
x=ones（n100）;
y=ones（n100）;
m=zeros（1100）;
m（1）=max（x（:1））;
y（:1）=x（:1）;
x（:2）=A*y（:1）;
m（2）=max（x（:2））;
y（:2）=x（:2）/m（2）;
p=0.0001;i=2;k=abs（m（2）-m（1））;
while??k>p
??i=i+1;
??x（:i）=A*y（:i-1）;
??m（i）=max（x（:i））;
??y（:i）=x（:i）/m（i）;
??k=abs（m（i）-m（i-1））;
end
a=sum（y（:i））;
w=y（:i）/a;
t=m（i）;
disp（w）;disp（t）;
?????????%以下是一致性檢驗
CI=（t-n）/（n-1）;RI=[0?0?0.52?0.89?1.12?1.26?1.36?1.41?1.46?1.49?1.52?1.54?1.56?1.58?1.59];
CR=CI/RI（n）;
if?CR<0.10
????disp（‘此矩陣的一致性可以接受!‘）;
????disp（‘CI=‘）;disp（CI）;
????disp（‘CR=‘）;disp（CR）;
end

function?AHPInit1（xy）
%層次分析的初始化
%默認只有兩層?x為準則數，y為方案數
%CToT為準則對目標生成的比較陣
%EigOfCri為準則層的特征向量
%EigOfOpt為選項層的特征向量

EigOfCri=zeros（x1）;%準則層的特征向量
EigOfOpt=zeros（yx）;
dim=x;%維度
RI=[0?0?0.58?0.90?1.12?1.24?1.32?1.41?1.45?1.49?1.51];%RI標準
%生成成對比較陣
for?i=1:dim
????????CToT（i:）=input（‘請輸入數據:‘）;
end
CToT?%輸出
pause
tempmatrix=zeros（x+1）;
tempmatrix=AHP1（dimCToT）;
EigOfCri=tempmatrix（1:x）;
ci1=tempmatrix（1+x）;
EigOfCri
ci1
pause

matrix=cell（x）;%元胞數組
ci=zeros（1x）;