拉依达找异常值

    occupy=xlsread('C:\Users\Hasee\Desktop\rs\占有率.xlsx');
    follow=xlsread('C:\Users\Hasee\Desktop\rs\流量.xlsx');
    A=occupy;
    B=follow;

    %std(A,flag,dim) A矩阵，flag0为除N-1的标准差，dim1按列求标准差
    M=zeros(1,23);
    M(1,1:end)=std(A(2:end,2:end));
    N=zeros(1,23);
    N(1,1:end)=std(B(2:end,2:end));
   
    k =1;
    P=zeros(50,3);
    Q=zeros(50,3);
    for j=2:24;
        for i=2:288;
            if (abs(A(i,j)-mean(A(2:288,j)))>3*M(1,j-1))
            %将异常数据先保存到P表中
            P(k,1)=A(i,j);
            P(k,2)=i;P(k,3)=j;
            A(i,j)=(A(i,j)+A(i,j-1))/2;
            %将异常数据保存到Q表中
            Q(k,1)=A(i,j);
            Q(k,2)=i;Q(k,3)=j;
            k=k+1;
            end
        end
    end
    
    k=1;
    Y=zeros(50,3);
    Z=zeros(50,3);
    for j=2:24;
        for i=2:288;
            if (abs(B(i,j)-mean(B(2:288,j)))>3*N(1,j-1))
            %将异常数据先保存到Y表中
            Y(k,1)=B(i,j);
            Y(k,2)=i;Y(k,3)=j;
            B(i,j)=(B(i,j)+B(i,j-1))/2;
            %将异常数据保存到Z表中
            Z(k,1)=B(i,j);
            Z(k,2)=i;Z(k,3)=j;
            k=k+1;
            end
        end
    end