将你的多余空格删掉,把不是6*6的矩阵用0填成6*6矩阵运行结果如下:
>> a=[1,1.023,1.124,1.180,1.242,1.238
0.977,1,1.099,1.153,1.213,1.210
0.889,0.910,1,1.050,1.104,1.101
0.847,0.867,0.953,1,1.052,1.049
0.805,0.824,0.905,0.950,1,0.997
0.808,0.827,0.908,0.953,1.003,1];
[x,y]=eig(a);eigenvalue=diag(y);lamda=eigenvalue(1);
ci1=(lamda-6)/5;cr1=ci1/1.24
w1=x(:,1)/sum(x(:,1))
b1=[1,0.248,0,0,0,0;
4.027,1,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0];
[x,y]=eig(b1);eigenvalue=diag(y);lamda=eigenvalue(1);
ci21=(lamda-2)/1;cr21=ci21/0.00
w21=x(:,1)/sum(x(:,1))
b2=[1,0.366,0,0,0,0;
2.733,1,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0];
[x,y]=eig(b2);eigenvalue=diag(y);lamda=eigenvalue(1);
ci22=(lamda-2)/1;cr22=ci22/0.00
w22=x(:,1)/sum(x(:,1))
b3=[1,0.647,18.125,0,0,0;
1.545,1,28,0,0,0;
0.055,0.036,1,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0];
[x,y]=eig(b3);eigenvalue=diag(y);lamda=eigenvalue(1);
ci23=(lamda-3)/2;cr23=ci23/0.58
w23=x(:,1)/sum(x(:,1))
b4=[1,0.802,12.462,0,0,0;
1.247,1,15.538,0,0,0;
0.08,0.064,1,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0];
[x,y]=eig(b4);eigenvalue=diag(y);lamda=eigenvalue(1);
ci24=(lamda-3)/2;cr24=ci24/0.58
w24=x(:,1)/sum(x(:,1))
b5=[1,1.144,49.75,0,0,0;
0.874,1,43.5,0,0,0;
0.02,0.023,1,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0];
[x,y]=eig(b5);eigenvalue=diag(y);lamda=eigenvalue(2);
ci25=(lamda-3)/2;cr25=ci25/0.58
w25=x(:,2)/sum(x(:,2))
b6=[1,1.193,102,0,0,0;
0.838,1,85.5,0,0,0;
0.01,0.012,1,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0;
0,0,0,0,0,0];
[x,y]=eig(b6);eigenvalue=diag(y);lamda=eigenvalue(1);
ci26=(lamda-3)/2;cr26=ci26/0.58
w26=x(:,1)/sum(x(:,1))
w_sum=[w21,w22,w23,w24,w25,w26]*w1
ci=[ci21,ci22,ci23,ci24,ci25,ci26];
cr=ci*w1/sum(0.58*w1)
cr1 =
-8.6843e-005
w1 =
0.1877
0.1834
0.1670
0.1591
0.1512
0.1517
cr21 =
-Inf
w21 =
0.1988
0.8012
0
0
0
0
cr22 =
Inf
w22 =
0.2679
0.7321
0
0
0
0
cr23 =
0.0013
w23 =
0.3846
0.5942
0.0213
0
0
0
cr24 =
-0.0024
w24 =
0.4298
0.5359
0.0344
0
0
0
cr25 =
-2.5846
w25 =
0.6463
0.3753
-0.0216
0
0
0
cr26 =
0.0131
w26 =
0.5411
0.4535
0.0054
0
0
0
w_sum =
0.3988
0.5947
0.0066
0
0
0
cr =
-0.3890
w_sum=[w21,w22,w23,w24,w25,w26]*w1
这条语句里面的[w21,w22,w23,w24,w25,w26]矩阵组合里的子矩阵必须要是一样的大小。如果不够的话要补零。下面也是一样的