untitled3.m

用于mann-kendall趋势分析

clear;clc;
yangben=[
47.924375 25.168125 827.4105438 330.08875 1045.164375 261.374375
47.924375 25.168125 827.4105438 330.08875 1045.164375 261.374375]; %样本数据

fangzhen=[
36.27 14.59 836.15 420.41 1011.83 189.54 
47.924375 25.168125 827.4105438 330.08875 1045.164375 261.374375

]; %待判数据

[rows,cols]=size(fangzhen);
p=0.5;    %分辨系数
[m,n]=size(yangben);      
R=[];
for irow=1:rows
     yy=fangzhen(irow,:);
     data=[yy;yangben];
     data_gyh1=mean(yangben)
     for i=1:m+1
         for j=1:n
             data_gyh(i,j)=data(i,j)/data_gyh1(j);
         end
     end
     
     for i=2:m+1
         for j=1:n
             Dij(i-1,j)=abs(data_gyh(1,j)-data_gyh(i,j));
         end
     end
     Dijmax=max(max(Dij));
     Dijmin=min(min(Dij));
     
     for i=1:m
         for j=1:n
             Lij(i,j)=(Dijmin+p*Dijmax)/(Dij(i,j)+p*Dijmax);
         end
     end
     
     LijRowSum=sum(Lij');
     
     for i=1:m
         Rij(i)=LijRowSum(i)/n;
     end
     R=[R;Rij];
end
R