Daniel趋势分析在matlab上的实现

原理与计算步骤

数据集

    某县1965-2013年年际气温数据，对其各时期进行Daniel趋势分析。

matlab代码

clc;clear;
data=[18.484 	18.664 	18.280 	18.498 	17.995 	18.117 	18.154 	18.388 	18.619 	18.243 	18.401 	17.742 	18.236 	18.377 	18.322 	18.326 	18.319 	18.251 	18.345 	17.531 	17.849 	18.075 	18.724 	17.890 	18.140 	18.583 	18.490 	18.034 	18.110 	18.608 	18.055 	17.962 	18.439 	19.400 	18.700 	18.465 	18.689 	18.999 	18.756 	18.296 	18.672 	19.049 	18.974 	18.796 	19.526 	19.257 	18.678 	19.008 	19.358 
];
Y={};B=[];c=[];
Y{1,1}=data;
Y{1,2}=data(:,1:5);
Y{1,3}=data(:,6:15);
Y{1,4}=data(:,16:25);
Y{1,5}=data(:,26:35);
Y{1,6}=data(:,36:49);
[c,N]=size(Y);
for j=1:N
    y1=Y{1,j};
    y2=Y{1,j};
    [r,n]=size(y1);
    for i=1:n
       for j1=i+1:n
           if y2(1,i) > y2(1,j1) 
            t=y2(1,i);y2(1,i) = y2(1,j1);y2(1,j1)=t;
           end
       end
    end
    k=0;
    for i=1:n
        for j1=1:n
            if y1(1,i)==y2(1,j1)
                k=k+1;
                c(1,k)=j1;
                break
            end
        end
    end
    s=0;
    for i=1:n
       s=s+(c(1,i)-i)^2;
    end
    B(j,1)=1 - (6 * s / (n ^ 3 - n));
end
disp('各时期秩相关系数为：')
disp(B)

运行结果

 

年份	Rs	趋势	Wp	显著性
1965-2013	0.550	上升	0.238	显著
1965-1969	-0.500	下降	0.900	不显著
1970-1979	0.127	上升	0.564	不显著
1980-1989	-0.297	下降	0.564	不显著
1990-1999	0.285	上升	0.564	不显著
2000-2013	0.591	上升	0.464	显著

如果对你有帮助，请点下赞，予人玫瑰手有余香！

时时仰望天空，理想就会离现实越来越近！