??%我自己寫的mk突變點分析，已經調試正確
%s是數據序列；sy是起始年份或序號（做圖的時候會用到）
function?[ufkubk]=?mkabrpt（xsy）
sy=load（‘F:\MK\2.txt‘）%這個是年份
x?=load（‘F:\MK\1.txt‘）%這個是徑流數據
????N=length（x）;??%%F:\畢業論文參考文獻\數學方法\mk\閩江水文站M_K\MkUB\月max\竹岐.txt
????ufk=smk（x）;
????y=flipud（x）;
????ubk=smk（y）;
????ubk=-ubk;
????ubk=flipud（ubk）;???
????%繪圖
????ey=sy+N-1;
????htx=[sy:ey];
????
????a=ceil（abs（max（[ufk;ubk]））+3）;
????b=floor（min（[ufk;ubk]）-1）;
?????????
????if?a<5
???????a=5
????end
????
????if?b>-4
???????b=-4
????end
????
????figure（2）
????plot（htxufk‘r-‘‘linewidth‘1.5）;
????hold?on
????plot（htxubk‘b-‘‘linewidth‘1.5）;
????plot（htx1.96*ones（N1）‘:‘‘linewidth‘1）;
????legend（‘UF統計量‘‘UB統計量‘‘0.05顯著水平‘）;
????xlabel（‘年度‘‘FontName‘‘TimesNewRoman‘‘FontSize‘16）;
????ylabel（‘統計量‘‘Font

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1106??2014-08-06?10:18??mk程序\mkabrpt.m

?????文件???????1994??2015-08-12?22:43??mk程序\mk計算Z值.m

?????文件????????736??2014-08-06?10:18??mk程序\smk.m

?????目錄??????????0??2015-08-12?22:43??mk程序

-----------?---------??----------?-----??----

?????????????????3836????????????????????4