clear?all
data=load（‘E:\PhD?paper\Huai?River_Actual?evapotranspiration\MK_mutation?Test\ETa_AA?model.txt‘）;

%?Mann-Kandel?突變點(diǎn)分析MATLAB腳本文件，使用前預(yù)先導(dǎo)入目標(biāo)數(shù)據(jù)，變量名設(shè)定為data，再將內(nèi)容復(fù)制至命令框即可。
N=size（data1）;
n=2;
nb=1;
nB=2;
ri=zeros（N1）;
rbi=zeros（N1）;
SFk=zeros（N1）;
SBk=zeros（N1）;
UFk=zeros（N1）;
UBk=zeros（N1）;
while?n<=N
????i=1;
????r=zeros（1N）;
????while?i????????if?data（i）<=data（n）
????????????r（i）=1;
????????else
????????????r（i）=0;
????????end
?????????i=i+1;
????end
????ri（n）=sum（r）;
????SFk（n）=SFk（n-1）+ri（n）;
????UFk（n）=（SFk（n）-n*（n-1）/4）/sqrt（n*（n-1）*（2*n+5）/72）;
????n=n+1;
end
while?nb<=N-1
????i=N;
????rb=zeros（1N）;
????while?i>nb
????????if?data（i）<=data（nb）
????????????rb（i）=1;
????????else
????????????rb（i）=0;
????????e