用Matlab进行MK检验Word格式.doc

1、1jix（m） s=s+1; elseif x（n）=x（m） s=s+0; else s=s-1; end end end %计算vars和e%vars=len*（len-1）*（2*len+5）/18;%计算zc %if s zc=（s-1）/sqrt（vars）;else zc=（s+1）/sqrt（vars）;end %计算zaza=var（x）;sign=0;zc1=abs（zc）;if zc1 = za sign=1;else sign=0;end %计算倾斜度 ndash = len * （ len - 1 ） / 2;slope1= zeros（ ndash, 1 ）;m=1;

2、for k = 1:len-1,for j = k+1:len,slope1（m） = （ x（j） - x（k） ） / （ j - k ） ; % m = m + 1;end;slope= median（ slope1 ）;2、M-K突变检验clc;y=data（1:50,6）;%输入数据n=length（y）;Sk=zeros（size（y）;UFk=zeros（size（y）;s1=0;for i=2:n % for j=1:i if y（i）y（j） s1=s1+1; s1=s1+0; % （） end; end;Sk（i）=s1; % E=i*（i-1）/4;Var=i*（i-1）

3、*（2*i+5）/72;UFk（i）=（Sk（i）-E）/sqrt（Var）; % y2=zeros（size（y）;Sk2=zeros（size（y）;UBk=zeros（size（y）; s2=0;for i=1:n y2（i）=y（n-i+1）; for j=1: if y2（i）y2（j） s2=s2+1; s2=s2+0; end;Sk2（i）=s2; E=i*（i-1）/4; % Sk2（i）的均值 Var=i*（i-1）*（2*i+5）/72; % Sk2（i）的方差 % UBk（i）=0-（Sk2（i）-E）/sqrt（Var）; UBk2=zeros（size（y）; UBk2（i）=UBk（n-i+1）;year=1961:2010;year=year;xlswrite（D:test.xls,year,Sheet1,A1）;,UFk,B1,UBk2,C1