exe_cengci.m

对层次分析法的程序实现

clear;
a=[1 5 2 1/3 1/2 1/5];
b1=[5 3 1/2];
b2=[1/7 1/3 5 1/7 1/5 1/7 1/5 2 1/5 3];%[1/7 1/3 5 1/7 1/5 1/7 1/5 3 1/5 5]
b3=[0.2 0.5 2];%[0.2 3 0.2]
b4=[7 3 0.2];
c1=[4 6 4];
c2=[3 5 1 3 3 3];%[3 5 5 3 5 3]
c3=[0.2];
c4=[5];
c5=[1/3 1/3 0.2 1/3 1/3 1.5];%[1/3 1/3 0.2 1/3 1/3 3]
c6=[0.5 3 5 3 5 1 3 5 1 1 0.5 1 0.25 0.25 0.25 5 8 2 2 2 7];%[1/7 1/3 7 1/3 7 0.5 0.2 7 0.2 1/3 1/7 3 1/7 1/7 1/7 5 8 5 5 5 8]
c7=[1/3];
c8=[0.5 0.5 0.5 3 3 5 0.5 0.5 1 0.2 2 2 3 0.5 3];%[0.2 0.2 1/3 3 3 5 0.2 0.2 0.25 1/3 0.2 2 3 0.2 3]
c9=[0.5	3	3];%[1/3 3 3]
c10=[3	3	1	0.5	0.2	0.2	1/3	1/3	1/3	5];%[3 3 1/3 1/3 0.2 0.2 1/3 1/3 0.2 5]
c11=[0.5	0.2	1/3	0.2	0.2	5	0.5	0.5	2	0.5	0.5	0.5	5	1	3];%[0.2 1/7 1/7 1/7 1/7 5 1/9 1/9 3 0.2 1/9 1/9 5 1 3]
c12=[5];
c13=[0.5	1/7	1/7	0.2	0.2	5	0.5	0.5	5	0.5];%[0.2 1/7 1/7 0.2 0.2 5 0.2 0.2 5 0.2]
c14=[2	6	4];%[3 5 5]
D1 = [ 0.2,0.8,0,0,0 ]; ; D2 = [ 0.3,0.7,0,0,0 ];  ; D3 = [ 0,0.3,0.7,0,0 ];
D4 = [ 0,0,0.6,0.4,0 ];  D5 = [ 0.3,0.7,0,0,0 ];   D6 = [ 0,0,0.3,0.7,0 ];
D7 = [ 0,0,0.4,0.6,0 ];  D8 = [ 0,0,0,0.6,0.4 ];   D9 = [ 0,0.6,0.4,0,0 ];
D10 = [ 0.2,0.8,0,0,0 ]; D11 = [ 0,0,0.2,0.4,0.4 ]; D12 = [ 0,0.8,0.2,0,0 ];
D13 = [ 0,0.3,0.4,0.3,0 ]; D14 = [ 0.4,0.6,0,0,0 ]; D15 = [ 0,0,0.5,0.5,0 ];
D16 = [ 0,0.7,0.3,0,0 ]; D17 = [ 0,0,0.2,0.4,0.4 ]; D18 = [ 0,0,0.7,0.3,0 ];
D19 = [ 0,0.7,0.3,0,0 ]; D20 = [ 0,0.8,0.2,0,0 ];  D21 = [ 0,0,0.8,0.2,0 ];
D22 = [ 0,0,0,0.3,0.7 ]; D23 = [ 0.3,0.7,0,0,0 ];  D24 = [ 0,0,0.3,0.5,0.2 ];
D25 = [ 0,0.7,0.3,0,0 ]; D26 = [ 0,0.2,0.5,0.3,0 ]; D27 = [ 0.3,0.7,0,0,0 ];
D28 = [ 0,0,0,0.3,0.7 ]; D29 = [ 0,0.7,0.3,0,0 ]; D30 = [ 0,0,0.2,0.5,0.3 ];
D31 = [ 0,0.2,0.6,0.2,0 ]; D32 = [ 0,0.6,0.4,0,0 ]; D33 = [ 0,0.2,0.6,0.2,0 ];
D34 = [ 0,0.2,0.5,0.3,0 ]; D35 = [ 0,0,0.3,0.5,0.2 ]; D36 = [ 0,0,0.5,0.3,0.2 ];
D37 = [ 0,0,0.4,0.6,0 ]; D38 = [ 0,0.2,0.6,0.2,0 ]; D39 = [ 0,0,0,0.3,0.7 ];
D40 = [ 0,0,0,0.8,0.2 ]; D41 = [ 0,0,0,0.6,0.4 ]; D42  = [ 0.2,0.8,0,0,0 ];
D43 = [ 0,0,0,0.8,0.2 ]; D44 = [ 0,0.2,0.5,0.3,0 ]; D45  = [ 0,0.3,0.7,0,0 ];
D46 = [ 0,0.5,0.5,0,0 ]; D47 = [ 0,0,0.2,0.8,0 ]; D48  = [ 0,0,0.3,0.7,0 ];
D49 = [ 0,0.7,0.3,0,0 ]; D50 = [ 0,0,0.3,0.7,0 ]; D51  = [ 0.2,0.5,0.3,0,0 ];
D52 = [ 0.2,0.8,0,0,0 ]; D53 = [ 0.2,0.8,0,0,0 ]; D54  = [ 0,0.4,0.6,0,0 ];
d=D1;
for i=2:54
    d=[d;eval(['D' num2str(i)])];
end
num=0;cc=[];eigv=[];crv=[];wvc=[];wvb=[];
for i=1:14
    [w eigenvalue cr matrix]=cengci(eval(['c' num2str(i)]));eigv=[eigv;eigenvalue];crv=[crv;cr];wvc=[wvc;w]
    eval(['cc' num2str(i) '=w''*d((num+1):(num+length(w)),:);'])
    cc=[cc;eval(['cc' num2str(i)])];
    matrix=[matrix w];
    save(['F:\ѧϰ2\matlab\cengci\c' num2str(i) '.txt'],'matrix','-ascii')
    num=num+length(w);
end
num=0;bb=[];
for i=1:4
    [w eigenvalue cr matrix]=cengci(eval(['b' num2str(i)]));eigv=[eigv;eigenvalue];crv=[crv;cr];wvb=[wvb;w]
    eval(['bb' num2str(i) '=w''*cc((num+1):(num+length(w)),:);'])
    bb=[bb;eval(['bb' num2str(i)])];
    matrix=[matrix w];
    save(['F:\ѧϰ2\matlab\cengci\b' num2str(i) '.txt'],'matrix','-ascii')
    num=num+length(w);
end
[w eigenvalue cr matrix]=cengci(a);eigv=[eigv;eigenvalue];crv=[crv;cr];
aa=w'*bb
matrix=[matrix w];
save('F:\ѧϰ2\matlab\cengci\a.txt','matrix','-ascii')