Daniel趋势分析在matlab上的实现

版权声明：转载需写明出处 https://blog.csdn.net/qq_23860475/article/details/80669048

原理与计算步骤

数据集

    某县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	显著

结束语

该篇资料出资我发表的论文：

[1]赵嘉阳,郭福涛,梁慧玲,苏漳文,王文辉,林玉蕊.福建长汀红壤区1965—2013年气温和降水量的变化趋势[J].福建农林大学学报(自然科学版),2016,45(01):77-83.
觉得有用的话写论文的时候麻烦加入参考文献哦！

备注：承接模型、算法代码实现（支持python、matlab），有意请联系QQ947943645 ，非诚勿扰！