版权声明:转载需写明出处 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 ,非诚勿扰!