matlab实现m-k突变的,用matlab进行mk趋势分析与突变检验.pdf

% M-K 趋势检定

clear, close all,clc % clear:清变数 close all:清图面 clc:清画面

% define and assign the full file path using "file open" dialog

[filename filepath]=uigetfile('data1.xls');

full_filepath=[filepath filename];

[X,TXTX,RAWX]=xlsread(full_filepath,1); % 数据全部读入,数据缺失不影响结果

x=X(:,1); % x 时间

y=X(:,2); % y 数据

% 计算 S

n=size(y,1); % 数据个数

S=0;

for i=1:n-1

S = S + sum(sign(y(i+1:n) - y(i))); % S 计算式

end

% 计算 VarS

VarS=n*(n-1)*(2*n+5)/18;

%计算 Z

if S>0

Z=(S-1)/sqrt(VarS);

else

Z=(S+1)/sqrt(VarS);

end

% 计算 Zabs

alpha1=0.05; % 信度 95% 的显著水平

alpha2=0.01; % 信度 99% 的显著水平

PZ1=norminv(1-alpha1/2,0,1);

PZ2=norminv(1-alpha2/2,0,1);

H=