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

%x存放data
%寻找函数，找到所有整数的组合
%53为一跨长度可任意取值，一般跨径为40~60m
%此函数仅考虑纵向分割
data = [];
for k=3:13for b=0:0.5:5a=(53-2*b)/k;if rem(a,1) ~= 0a;continueelsea;b;k;x=[a,b,k];data = [data;x];endend
end
data%筛选函数，选出更适合的组合
index=[];
A=data(:,1);
for i=1:length(A)if A(i)>=3 & A(i)<=7index=[index;i];elsecontinueend
end
indexdataX=[];
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);endend
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);endend
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);endend
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.1disp('一致性检测通过')%归一化处理，得出权重系数for i =1:n1for j=1:n1PP(i,j)=P(i,j)/sum(sum(P));endendPPH=sum(PP,2);for i =1:length(H)W(i)=H(i)/sum(H);end %权重值Wx_new=[AAA;BBB;KKK];PJZ=W*x_new[Z,xidx]=sort(PJZ,'descend');xidxfor i=1:length(xidx)ZZJG=[ZZJG;dataXX(xidx(i),:)];endZZJG2=ZZJG(2:size(ZZJG,1),:);disp('跨进评价结果（单位m）:')disp('第一推荐值:')disp('切割长度 两侧长度 跨径')disp(ZZJG(1,:))disp('其他推荐值:')disp((ZZJG2))
elsedisp('一致性检测失败')
end

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