海关论文matlab程序以及其他的程序.docx
《海关论文matlab程序以及其他的程序.docx》由会员分享，可在线阅读，更多相关《海关论文matlab程序以及其他的程序.docx（10页珍藏版）》请在冰豆网上搜索。
海关论文matlab程序以及其他的程序
影响
X=[1,1/4
4,1];
[x,y]=eig（X）;
Jmax1=max（max（y））
w1=x（:
1）/sum（x（:
1））
[N,M]=size（X）;
RI=[000.580.901.121.241.321.411.451.491.51];
CI1=（Jmax1-N）/（N-1）
CR1=CI1/RI（1,N）
ifCR1<0.1
disp（'通过一致性检验'）
else
disp（'未通过一致性检验'）
end
A=[1,1/5,1/5,1,2
5,1,1,4,7
5,1,1,4,7
1,1/4,1/4,1,1/2
1/2,1/7,1/7,1/2,1];
[x,y]=eig（A）;
Jmax21=max（max（y））
w21=x（:
1）/sum（x（:
1））
[N,M]=size（A）;
RI=[000.580.901.121.241.321.411.451.491.51];
CI21=（Jmax21-N）/（N-1）
CR21=CI21/RI（1,N）
ifCR21<0.1
disp（'通过一致性检验'）
else
disp（'未通过一致性检验'）
end
B=[1,1,1/6,1/3,1/5
1,1,1/5,1/2,1/5
6,5,1,3,2
3,2,1/3,1,1/2
5,5,1/2,2,1];
[x,y]=eig（B）;
Jmax22=max（max（y））
w22=x（:
1）/sum（x（:
1））
[N,M]=size（B）;
RI=[000.580.901.121.241.321.411.451.491.51];
CI22=（Jmax22-N）/（N-1）
CR22=CI22/RI（1,N）
ifCR22<0.1
disp（'通过一致性检验'）
else
disp（'未通过一致性检验'）
end
运行结果：
Jmax1=
2
w1=
0.2000
0.8000
CI1=
NaN
CR1=
0
通过一致性检验
Jmax21=
4.8519
w21=
0.0877
0.3907
0.3907
0.0801
0.0508
CI21=
-0.0370
CR21=
-0.0331
通过一致性检验
Jmax22=
5.0524
w22=
0.0614
0.0692
0.4302
0.1506
0.2886
CI22=
0.0131
CR22=
0.0117
通过一致性检验
概率
X=[1,1/6
6,1];
[x,y]=eig（X）;
Jmax1=max（max（y））
w1=x（:
1）/sum（x（:
1））
[N,M]=size（X）;
RI=[000.580.901.121.241.321.411.451.491.51];
CI1=（Jmax1-N）/（N-1）
CR1=CI1/RI（1,N）
ifCR1<0.1
disp（'通过一致性检验'）
else
disp（'未通过一致性检验'）
end
A=[1,1/4,1/2,2,1/3
4,1,3,6,2
2,1/3,1,3,1/2
1/2,1/6,1/3,1,1/5
3,1/2,2,5,1];
[x,y]=eig（A）;
Jmax21=max（max（y））
w21=x（:
1）/sum（x（:
1））
[N,M]=size（A）;
RI=[000.580.901.121.241.321.411.451.491.51];
CI21=（Jmax21-N）/（N-1）
CR21=CI21/RI（1,N）
ifCR21<0.1
disp（'通过一致性检验'）
else
disp（'未通过一致性检验'）
end
B=[1,2,1/4,1/2,1/3
1/2,1,1/6,1/3,1/5
4,6,1,2,2
2,3,1/2,1,1/3
3,5,1/2,3,1];
[x,y]=eig（B）;
Jmax22=max（max（y））
w22=x（:
1）/sum（x（:
1））
[N,M]=size（B）;
RI=[000.580.901.121.241.321.411.451.491.51];
CI22=（Jmax22-N）/（N-1）
CR22=CI22/RI（1,N）
ifCR22<0.1
disp（'通过一致性检验'）
else
disp（'未通过一致性检验'）
end
运行结果：
Jmax1=
2
w1=
0.1429
0.8571
CI1=
NaN
CR1=
0
通过一致性检验
Jmax21=
5.0463
w21=
0.0951
0.4247
0.1562
0.0554
0.2686
CI21=
0.0116
CR21=
0.0103
通过一致性检验
Jmax22=
5.1086
w22=
0.0941
0.0548
0.3941
0.1581
0.2988
CI22=
0.0271
CR22=
0.0242
通过一致性检验
[x,y]=eig（A）;
Jmax=max（max（y））
[N,M]=size（A）;
RI=1.24;
CI=（Jmax-N）/（N-1）
CR=CI/RI
ifCR<0.1
disp（'通过一致性检验'）
B=zeros（N,M）;
fori=1:
M
b（:
i）=a（:
i）/sum（a（:
i））;
end
b;
fori=1:
N
W（1,i）=sum（b（i,:
））;
end
W;
s=sum（W）;
c=W./s
else
disp（'未通过一致性检验'）
end
clc
clear %修改对比矩阵、一致性检验就可以
a=[1,1,1,4,1,1/2
1,1,2,4,1,1/2
1,1/2,1,5,3,1/2
1/4,1/4,1/5,1,1/3,1/3
1,1,1/3,3,1,1
2,2,2,3,3,1];
[x,y]=eig（a）;
Jmax=max（max（y））
w1=x（:
1）/sum（x（:
1））
[N,M]=size（a）;
RI=[000.580.901.121.241.321.411.451.491.51];
ci1=（Jmax-N）/（N-1）
cr=ci1/RI（1,N）
ifcr<0.1
disp（'通过一致性检验'）
else
disp（'未通过一致性检验'）
end
b1=[1,1/4,1/2;4,1,3;2,1/3,1];
[x,y]=eig（b1）;eigenvalue=diag（y）;lamda=eigenvalue
（1）;
ci21=（lamda-3）/2;cr21=ci21/0.58
w21=x（:
1）/sum（x（:
1））
b2=[1 1/4 1/5;4 1 1/2;5 2 1];
[x,y]=eig（b2）;eigenvalue=diag（y）;lamda=eigenvalue
（1）;
ci22=（lamda-3）/2;cr22=ci22/0.58
w22=x（:
1）/sum（x（:
1））
b3=[1 3 1/3;1/3 1 1/7;3 7 1];
[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 1/3 5;3 1 7;1/5 1/7 1];
[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 7;1 1 7;1/7 1/7 1];
[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 7 9;1/7 1 1;1/9 1 1];
[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）
ri=[0,0,0.52,0.89,1.12,1.26,1.36,1.41,1.46,1.49];%随机一致性指标RI对照表
a=[1,1/3,1/2;3,1,3;2,1/3,1];%输入矩阵a
d=eig（a）;%获得特征根
[m,n]=size（a）;%获取矩阵的阶数
maxd=max（d）%获得最大特征根
ci=（maxd-m）/（m-1）%计算一致性指标CI
cr=ci/ri（1,m）%计算一致性比率CR
升级会员 