各位，谁用Matlab做过Pettitt突变分析，我刚刚气温的突变，程序和图附上

xiaolan16

但是，用累计距平验证，总觉得突变的时间点不对，我觉得这个突变时间应该在1997年左右，但是Pettitt做出来太靠后，而且不明显
程序如下：
function [Ut,chg_point]=pettitt_change(x,data)
%
% Pettitt test for detection of change point
A=textread('气温年值.txt')
x=A(:,1)%时间
data=A(:,2)%年平均气温
n=length(data);
s=0;
Ut=zeros(n-1,1);
for i=1:n-1;
  for j=i+1:n;
      s=s+sign(data(i)-data(j));
  end
  Ut(i)=s;
end
Kmax=max(abs(Ut));
chg_point=find(abs(Ut)==Kmax);%找出最大值对应的序列号
coef=2*exp(-6*Kmax^2/(n^3+n^2));
figure(1)
hold on
plot(x(1:n-1),Ut,'K-','linewidth',1.5);%画出整个曲线的变化趋势
plot([x(chg_point),x(chg_point)],[-Kmax,Kmax],'r--','linewidth',1.5);%画出该竖线即为突变点
Klim=Kmax-Kmax*0.1
plot([x(1),x(n)],[-Klim,-Klim],'b-');
plot([x(1),x(n)],[Klim,Klim],'b-');
end

andrewsoong

运行了一下，跟我用R写的的不一样啊。你的结果可能有问题。MK检验与你的差别也很大

topmad

andrewsoong 发表于 2017-1-9 23:15
运行了一下，跟我用R写的的不一样啊。你的结果可能有问题。MK检验与你的差别也很大

大神 请教我用R语言

andrewsoong

topmad 发表于 2017-1-10 08:25
大神 请教我用R语言

老五别闹(ー_ー)!!

xiaolan16

andrewsoong 发表于 2017-1-9 23:15
运行了一下，跟我用R写的的不一样啊。你的结果可能有问题。MK检验与你的差别也很大

我不太会编程，R语言也用过，能不能给我发一个你的R程序，我再试试，非常着急，拜托了,我邮箱1286637055@qq.com

luluhoy

http://download.csdn.net/detail/luluhoy/9644703

luluhoy

xiaolan16 发表于 2017-1-10 09:22
我不太会编程，R语言也用过，能不能给我发一个你的R程序，我再试试，非常着急，拜托了,我邮箱

能否给我也发一份，hoyjiaya@126.com，非常感谢~