function?[?mk?]?=MK（?x?）
%UNtitleD?Summary?of?this?function?goes?here
%???Detailed?explanation?goes?here
X=x;
n=length（X）;
S=0;
????for?k=1:（n-1）
????????for?j=（k+1）:n
????????????S=S+sign（X（j）-X（k））;
????????end???
????end
?VarS=n*（n-1）*（2*n+5）/18;
????if?S>0
????????Z=（S+1）/sqrt（VarS）;
????elseif?S==0
????????Z=0;
????else
????????Z=（S-1）/sqrt（VarS）;
????end
?
Qi=zeros（n*（n-1）/21）;
counter=1;
for?k=1:（n-1）
????for?j=（k+1）:n
????????Qi（counter）=（X（j）-X（k））/（j-k）;
????????counter=counter+1;
????end
end
Alpha=1-normcdf（Z01）;
disp（‘趨勢的顯著水平為‘）;
disp（Alpha）;

Q=median（Qi）;
disp（‘整體趨勢的變化速率為‘）;
disp（Q）;
mk=[ZAlphaQ];
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????433??2014-05-03?12:45??Pettitt.asv
?????文件?????????503??2014-05-04?17:40??Pettitt.m
?????文件?????????503??2014-05-06?10:15??Pettitt2.asv
?????文件?????????537??2014-05-06?11:04??Pettitt2.m
?????文件?????????678??2014-05-06?11:43??MK.asv
?????文件?????????687??2014-05-06?12:55??MK.m