qqq.m

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

clear
[A]=Lee
[m n]=size(A)
len=283;
i1=1;
i2=1;
j=1;
k=1;
i3=1;
m=1;


for j=1:9
    for k=1:6

        for i1=1:len
            if(A(i1,7)==j&A(i1,10)==k)
                Matrix(i2,:)=A(i1,:);
                i2=i2+1;
            end
            
            i1=i1+1;
        end
        k=k+1;
    end
    j=j+1;
end

j=Matrix(1,7);
k=Matrix(1,10);
m=Matrix(1,8);
i4=1;

for i3=1: 233
    if(Matrix(i3,7)==j&Matrix(i3,8)==m&Matrix(i3,10)==k)
        dat1(i4)=Matrix(i3,1);
        dat2(i4)=Matrix(i3,2);
        dat3(i4)=Matrix(i3,3);
        dat4(i4)=Matrix(i3,4);
        i4=i4+1;
    end
    if(Matrix(i3,7)==j&Matrix(i3,8)~=m&Matrix(i3,10)==k)
        Matrix(i3-1,12)=mean(dat1);
        Matrix(i3-1,13)=mean(dat2);
        Matrix(i3-1,14)=mean(dat3);
        Matrix(i3-1,15)=mean(dat4);
        Matrix(i3-1,17)=cov(dat1);
        Matrix(i3-1,18)=cov(dat2);
        Matrix(i3-1,19)=cov(dat3);
        Matrix(i3-1,20)=cov(dat4);
        i4=1;
        dat1=0;
        dat2=0;
        dat3=0;
        dat4=0;
        dat1(i4)=Matrix(i3,1);
        dat2(i4)=Matrix(i3,2);
        dat3(i4)=Matrix(i3,3);
        dat4(i4)=Matrix(i3,4);
        i4=i4+1;
        m=Matrix(i3,8);
        
    end
    
    if(Matrix(i3,7)~=j|Matrix(i3,10)~=k)
        Matrix(i3-1,12)=mean(dat1);
        Matrix(i3-1,13)=mean(dat2);
        Matrix(i3-1,14)=mean(dat3);
        Matrix(i3-1,15)=mean(dat4);
        Matrix(i3-1,17)=cov(dat1);
        Matrix(i3-1,18)=cov(dat2);
        Matrix(i3-1,19)=cov(dat3);
        Matrix(i3-1,20)=cov(dat4);
        i4=1;
        dat1=0;
        dat2=0;
        dat3=0;
        dat4=0;
        dat1(i4)=Matrix(i3,1);
        dat2(i4)=Matrix(i3,2);
        dat3(i4)=Matrix(i3,3);
        dat4(i4)=Matrix(i3,4);
        i4=i4+1;
        
        j=Matrix(i3,7);
        k=Matrix(i3,10);
        m=Matrix(i3,8);
    end
    i3=i3+1;
end



j=Matrix(1,7);
k=Matrix(1,10);
m=Matrix(1,8);
i4=1;
i6=1;
%i7=1;

for i5=1: 233
    if(Matrix(i5,7)==j&Matrix(i5,10)==k)
        %if(Matrix(i5,12)~=0&Matrix(i5,13)~=0&Matrix(i5,14)~=0&Matrix(i5,15)~=0)
            dat1(i4)=Matrix(i5,12);  
            dat2(i4)=Matrix(i5,13);            
            dat3(i4)=Matrix(i5,14);    
            dat4(i4)=Matrix(i5,15);
            
            dat5(i4)=Matrix(i5,1);
            dat6(i4)=Matrix(i5,2);
            dat7(i4)=Matrix(i5,3);
            dat8(i4)=Matrix(i5,4);
            i4=i4+1;
            %i7=i7+1;
        
    
    
    else
         Matrix2(i6,1)=Matrix(i5-1,7);
         Matrix2(i6,2)=Matrix(i5-1,10);
         
         Matrix2(i6,4)=mean(dat1);
         Matrix2(i6,5)=mean(dat2);
         Matrix2(i6,6)=mean(dat3);
         Matrix2(i6,7)=mean(dat4);
         
         Matrix2(i6,9)=cov(dat1);
         Matrix2(i6,10)=cov(dat2);
         Matrix2(i6,11)=cov(dat3);
         Matrix2(i6,12)=cov(dat4);
         
         
         
         Matrix2(i6,21)=mean(dat5);
         Matrix2(i6,22)=mean(dat6);
         Matrix2(i6,23)=mean(dat7);
         Matrix2(i6,24)=mean(dat8);
         
         Matrix2(i6,25)=cov(dat5);
         Matrix2(i6,26)=cov(dat6);
         Matrix2(i6,27)=cov(dat7);
         Matrix2(i6,28)=cov(dat8);
         
         
         
         dat1=0;
         dat2=0;
         dat3=0;
         dat4=0;
         
         dat5=0;
         dat6=0;
         dat7=0;
         dat8=0;
         
         i4=1;
         
         dat1(i4)=Matrix(i5,12);  
         dat2(i4)=Matrix(i5,13);            
         dat3(i4)=Matrix(i5,14);    
         dat4(i4)=Matrix(i5,15);
         i4=i4+1;
         
         j=Matrix(i5,7);
         k=Matrix(i5,10);
         i6=i6+1;
    end
    
    i5=i5+1;
end

Y=data;
[m n]=size(Y);

for j=1:n
    len=0;
    for i=1:m
        if(Y(i,j))
           len=len+1;
           B(len)=Y(i,j);
        end
    end
    store(1,j)=mean(B);
    store(2,j)=var(B);
end

[xx xxx]=size(Matrix2);
for i=1:xx
    
    Matrix2(i,29:32)=Matrix2(i,9:12)./Matrix2(i,25:28);

end