function?Aburptchangetime=mk（arraytimealphadrawing）?%突變檢驗(yàn)函數(shù)
%array為待檢驗(yàn)的時(shí)間序列
%time為與array對(duì)應(yīng)的時(shí)間，默認(rèn)按時(shí)間升序排列
%array與time均為行向量
%alpha為顯著水平
%drawing為判定條件，如果drawing=1，則繪圖
%Abruptchangetime為返回的突變時(shí)間

%正序計(jì)算
n=length（array）;%返回array元素個(gè)數(shù)
s=0;%初始化s
for?i=2:n
????for?j=1:i
????????if?array（i）>array（j）
????????????s=s+1;
????????end
????end
????Sk（i）=s;
????E=i*（i+1）/4;
????Var=i*（i-1）*（2*i+5）/72;
????UFk（i）=（Sk（i）-E）/Var^0.5;
end

%逆序計(jì)算
clear?Sk;
s=0;%再次初始化s
for?i=1:n
????arrayinverting（i）=array（n-i+1）;%反轉(zhuǎn)array為arrayinverting
end
for?i=2:n
????for?j=1:i
????????if?arrayinverting（i）>arrayinverting（j）
????????????s=s+1;
????????end
????end
????Sk（i）=s;
????E=i*（i+1）/4;
????Var=i*（i-1）*（2*i+5）/72;
????UB（i）=-（Sk（i）-E）/Var^0.5;
end
for?i=1:n
????UBk（i）=UB（n-i+1）;%反