[原创]Matlab气候突变检测程序合集（更新了）

linwei

示范点发射点发射点犯得上犯得上犯得上犯得上发士大夫阿斯顿发射点

linwei

十大发射点犯得上范德萨范德萨范德萨发发士大夫s'f

青风海魅

???  In an assignment  A(:) = B, the number of elements in A and B
must be the same.
跑lepage突变程序的时候出现这个，怎么解？

青风海魅

n1 = 25;
n2 = n1;
N = n1+n2;
for i = 1:length(y)-N+1
y1(:) = y(i:i+n1-1);
y2(:) = y(i+n1:i+N-1);
[~,~,Wstats] = ranksum(y1,y2);
W = Wstats.ranksum;
EW = n1*(N+1)/2;
VW = n1*n2*(N+1)/12;
[~,~,Astats] = ansaribradley(y1,y2);
A = Astats.W;
EA = n1*(N+2)/4;
VA = n1*n2*(N^2-4)/48/(N-1);
WA(i) = (W-EW)^2/VW+(A-EA)^2/VA;
end
???  In an assignment  A(:) = B, the number of elements in A and B
must be the same.
跑lepage突变程序的时候出现这个，怎么解？
谢谢楼主

青风海魅

恳请回复。谢谢

sfhstcn2

青风海魅 发表于 2017-3-16 15:58
恳请回复。谢谢

这是由于做循环的时候用于赋值的数组的size超出了被赋值的数组，把y1和y2后面的(:)去掉吧。

Quincy

感谢分享！

nabula

这么厉害！！！

hydrocharlene

非常感谢

青风海魅

sfhstcn2 发表于 2017-3-17 22:14
这是由于做循环的时候用于赋值的数组的size超出了被赋值的数组，把y1和y2后面的(:)去掉吧。

谢谢楼主回复，目前已能跑出程序，但是还是有一个问题：
chip1(1:length(WA)) = chi2inv(1-p/2,2);
chip2(1:length(WA)) = chi2inv(p/2,2);

%return
figure
if isempty(x)
    plot(WA,'color','b');
    hold on
    plot(chip1,'color','r');
    hold on
    plot(chip2,'color','r');
else
    xt = x(n1:length(y)-n2);
    plot(xt,WA,'color','b');
    hold on
    plot(xt,chip1,'color','r');
    hold on
    plot(xt,chip2,'color','r');
end
title('Le Page method')
end

??? Undefined function or variable 'p'.
若是不加第一个end后也能跑出来，但是这个问题提示怎么破？不用图的话，自己导出数据也能做出图