q3_1.asv

高速公路路面质量改进方案 采用matlab编程 达到了稳定点

function Q3_1
[N M G]=data2;

%插值法
[m n]=size(M);
for k=1:m
    if(sum(M(k,:))==0)
        continue;
    end
    line=1;
    x=1;
    y=1;
    for j=1:n
        if(M(k,j))
            x(line)=N(j);
            y(line)=M(k,j);
            line=line+1;
        end   
    end
    %display(k);
     MM(k,:)=interp1(x,y,N,'linear'); 
     %clear(x);
     %clear(y);
end

%删除NaN
MMlen = length(MM);
for i=1:MMlen
    if(MM(i,:)==0)
        continue;
    end
    for j=1:13
       if(isnan(MM(i,j)))
           MM(i,j)=100;
       end
    end
end

%plot(N,MM);%少了最后一行
%N=log(N);
plot(N,MM(1:10,:));

%计算每一个部分的百分比
for i=1:n-1
    MM(:,i)=MM(:,i)-MM(:,i+1);
end
%clc;
%[m n]=size(MM);
%for i=1:m
    %for j=1:n
        %if(MM(i,j)<0)
            %display(i);
            %display(j);
        %end
    %end
%end

%计算平均直径
D=MM*N'/100;

%计算方差
for i=1:MMlen
    mean=(N-D(i)*ones(1,n));
    %.*M(i,:)/100;
    A(i)=dot(mean,mean.*MM(i,:)/100);
    
    %A(i)=(M(i,:)-D(i)*ones(1,n))'*(M(i,:)-D(i)*ones(1,n));
end
A = A';


%根据Gmm挑出数据完全的组合
[m n]=size(G);
Glen=0;
for i =1:m
    if(G(i,1)~=0 & G(i,2)~=0)
        Glen=Glen+1;
        R(Glen,:)=[MM(i,:) D(i) A(i) G(i,:)];
    end
end


plot(N,R(:,1:13));
%plot3(R(:,17),R(:,14),R(:,15),'*');
%plot(N,R(:,1:13));  
%clear
%x = min(N):0.02:max(N);
%m=length(R);
%for i=1:m
    %draw(i,:)=interp1(N,R(i,1:13),x,'cubic'); 
%end
%y=interp1(N,R(1,1:13),x,'cubic');    
%plot(x,draw);
clc;