matlab 自动选择不标准桥梁跨径代码（源代码完全版）

%x存放data
%寻找函数，找到所有整数的组合
%53为一跨长度可任意取值，一般跨径为40~60m
%此函数仅考虑纵向分割
data = [];
for k=3:13
    for b=0:0.5:5
        a=(53-2*b)/k;
        if rem(a,1) ~= 0
            a;
            continue
        else
            a;b;k;
            x=[a,b,k];
            data = [data;x];
        end
    end
end
data

%筛选函数，选出更适合的组合
index=[];
A=data(:,1);
for i=1:length(A)
    if A(i)>=3 & A(i)<=7
        index=[index;i];
    else
        continue
    end
end
index

dataX=[];
for i=1:length(index)
    aa=data(index(i),:);
    dataX=[dataX;aa];
end
dataX

%平均性检测，选择比较平均的几组数据
Z=ones(1,size(dataX,1));
AA=5;BB=2.5;KK=8;
for i=1:size(dataX,1)
    Z(i)=abs(dataX(i,1)-AA)/AA+abs(dataX(i,2)-BB)/BB+abs(dataX(i,3)-KK)/KK;
end
Z

[F,ix]=sort(Z,'descend');
ixx=ix(1,fix(size(dataX,1)/2):size(dataX,1));
%将ixx倒叙
ixxx=fliplr(ixx)

dataXX=[];
for i=1:length(ixxx)
    p=dataX(ixxx(i),:);
    dataXX=[dataXX;p];
end
dataXX

%对a打分
for i=1:size(dataXX,1)
    A2(i)=abs(dataXX(i,1)-5)/5;
end
A2
ta=tabulate(A2(:))
NZ=size(ta,1)
AAA=[];
for i =1:size(ta,1)
    for j=1:length(A2)
        if A2(j) == ta(i,1)
            AAA(j)=1-(i-1)*(1/NZ);
        end
    end
end
AAA

%对b打分
for i=1:size(dataXX,1)
    B2(i)=abs(dataXX(i,2)-2.5)/2.5;
end
B2
ta2=tabulate(B2(:))
for i =1:size(ta2,1)
    for j=1:length(B2)
        if B2(j) == ta2(i,1)
            BBB(j)=1-(i-1)*(1/NZ);
        end
    end
end
BBB
%对c打分
for i=1:size(dataXX,1)
    K2(i)=abs(dataXX(i,3)-8)/8;
end
K2
ta3=tabulate(K2(:))
for i =1:size(ta3,1)
    for j=1:length(K2)
        if K2(j) == ta3(i,1)
            KKK(j)=1-(i-1)*(1/NZ);
        end
    end
end
KKK

%层次分析法
%为设置判断矩阵
%一致性指标查表，3阶Ri为0.58
%原则上CR小于0.1时，可以考虑
P=[1,5,3;
   1/5,1,1/3;
   1/3,3,1;];
n1=length(P);
ri=[0,0,0.58,0.9,1.12];
[x,y]=eig(P);
lamda=max(diag(y))
cr0=(lamda-n1)/(n1-1)/ri(n1)
ZZJG=[];
if cr0<0.1
    disp('一致性检测通过')
    %归一化处理，得出权重系数
    for i =1:n1
        for j=1:n1
            PP(i,j)=P(i,j)/sum(sum(P));
        end
    end
    PP
    H=sum(PP,2);
    for i =1:length(H)
        W(i)=H(i)/sum(H);
    end 
    %权重值
    W
    x_new=[AAA;BBB;KKK];
    PJZ=W*x_new
    [Z,xidx]=sort(PJZ,'descend');
    xidx
    for i=1:length(xidx)
        ZZJG=[ZZJG;dataXX(xidx(i),:)];
    end
    ZZJG2=ZZJG(2:size(ZZJG,1),:);
    disp('跨进评价结果（单位m）:')
    disp('第一推荐值:')
    disp('切割长度 两侧长度 跨径')
    disp(ZZJG(1,:))
    disp('其他推荐值:')
    disp((ZZJG2))
else
    disp('一致性检测失败')
end

这是增加了AHP层次分析法后的代码，考虑了三种A,B,K的等级分数。

建立了3*3的评价矩阵P

通过权重W相乘得到了最后的结果。

下一部是建立一个gui界面，将其表现的更加直观