asas.m

来自「用于mann-kendall趋势分析」· M 代码 · 共 43 行

function asas
a=[-0.62
0.09
-0.16
-1.9
1.06
-1.13
2.69
2.47
0.9
0.65
1.16
-0.88
-0.98
1.28
1.14
-0.84
];
b=[-0.62
0.09
-0.16
-1.9
1.06
-1.13
2.69
2.47
0.9
0.65
1.16
-0.88
-0.98
1.28
1.14
-0.84
];
c=a-b
for i=1:1:15
   for j=2:1:16 
    c=(b(j)-a(i))/(j-i)
   end  
end