q3_2.asv

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

function Q3_2
[N M G]=data2;
VMA=[
    12.48
12.08
12.42
13.03
21.84
13.05
14.68
12.47
12.12
13.71
12.31
14.18
12.24
13.08
14.29
12.24
13.1
14.52
12.31
13.4
12.34
13.1
13.67
12.18
13.28
13.48
12.77
13.27
14.07
12.51
12.51
12.4
13.07
15
12.77
13.63
15.08
13.19
14.68
14.03
16.5
13.19
13.21
14.51
13.31
12.5
12.14
13.45
13.11
13.03
13.9
13.17
14.01
16.6
11.8
12.1
13.63
12.11
13.13
13.04
12.2
13.36
13.34
13.59
14.76
13.47
12.54
12.28
12.6
16.7
12.8
13.2
12.8
12.71
12.25
12.58
12.93
14.03
13.2
13.87
16.5
16.58
16.6
13.7
13.1
13.04
13.04
16.7
16.6
14.27
14.61
13.3
14
13.47
13.47
13.21
22.6
14.4
13.45
12.7
17.07
14.1
16.97
16.93
13.43
16.5
17
16.64
16.8
12.07
13.55
16.84
16.6
13.14
14.33
16.6
17.15
14.49
16.5
16.6
14.14
13.71
15.32
16.79
14.06
12.7
17.1
16.9
16.5
16.8
17
13.58
16.6
16.51
16.63
16.55
14.33
13
14.98
16.79
14.71
13.89
12.2
16.5
0
13.3
13.3
13.3
13.1
13.3
17.3
17.3
17.1
17
12.9
16.9
13.5
14.9
14.6
12.4
14
14.6
13.39
12.2
12
21.5
12.13
12.37
12.85
12.71
13.13
14.3
14
14
13.5
13.9
14.4
0
13.07
13.5
15
15.41
14.8
0
0
13.7
0
18.14
13.64
14.6
16.94
15.33
13.02
13.15
13.23
13.08
13.25
12.01
12.39
12.35
12.31
12.95
12.51
12.4
16.94
14.91
14.61
14.31
14.24
12.29
13.31
12.63
12.48
13.74
14.21
12.37
14.57
13.63
13.05
12.34
14.74
13.09
12.19
13.08
15.22
16.6
13.08
12.23
12.15
13.03
16.65
13.31
15.07
14.77
14.87
16.84
12.1
12.1
14.49
16.62
13.27
14.28
14.64
13.06
14.61
13.15
12.86
12.66
21.93
17.22
13.13
13.42
17.06
13.4
13.66
14.89
13.58
14.53
13.36
13.13
14.4
13.01
13.07
13.42
14.53
14.6
16.74
14.06
16.5
13.67
12.03
15.15
15.22
14.74
14.02
14.37
13.09
13.48
14.32
13.18
13.06
13.55
12.26
16.52
16.73
13.2
13.23
16.77
13.1
0
];
%插值法
[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:5,:));

%计算每一个部分的百分比
%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=length(VMA);
VMAlen=0;
for i =1:m
    if(VMA(i)~=0)
        VMAlen=VMAlen+1;
        R(VMAlen,:)=[MM(i,:) VMA(i)];
    end
end

for i=1:length(N)
    plot(R(:,i),R(:,legth(N)+1));
    pause;
  
end
%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;