(1) 累计距平均
data=xlsread(‘test.xlsx’);
data(:,1)=[];
mean_data=mean(data);
diff_streamflow=data(:,1)-mean_data(1);
diff_baseflow=data(:,2)-mean_data(2);
diff_surfq=data(:,3)-mean_data(3);
Sum_streamflow_anomaly=zeros(55,1);
Sum_baseflow_anomaly=zeros(55,1);
Sum_surfq_anomaly=zeros(55,1);
for i=1:55;
Sum_streamflow_anomaly(i)=sum(diff_streamflow(1:i));
Sum_baseflow_anomaly(i)=sum(diff_baseflow(1:i));
Sum_surfq_anomaly(i)=sum(diff_surfq(1:i));
end
years=(1960:2014)’;
xlswrite(‘cumulative_anomaly.xlsx’,[years,Sum_streamflow_anomaly,Sum_baseflow_anomaly,Sum_surfq_anomaly])
通过绘图找到分界点即得到break point
(2) Pettitt检测
% This code is used to find the change point in a univariate continuous time series
% using Pettitt Test.
%
%
% The test here assumed is two-tailed test. The hypothesis are as follow:
% H (Null Hypothesis): There is no change point in the series
% H(Alternative Hypothesis): There is a change point in the se[……]
