bg算法 matlab,求助大神，有关BG算法

程序如下

function [ FLAGS ] = BGSA( x,y,P0,L0 )

% [ FLAGS ] = BGSA( x,y,P0,L0 )

%  x: 序列的x坐标(仅用于绘图，如果不使用则设置为[])

%  y: 降雨量序列

%  P0: 显著性水平门限值，低于此值的不再分割

%  L0: 最小分割尺度，子段长度小于此值的不再分割

%  FLAGS: 返回的逻辑向量与y大小相同，值为1表示相应的位置定位分割点

% 通过Bernaola-Galvan分割算法(BGSA)找出并绘制突变点及其最大密度间隔。

%   如何解释这个数字？

%   蓝色细曲线是原始序列; 而粗一个表示每个段中原始序列的间隔平均值t; 红色垂直线表示分割点所在的位置;阴影指示最大突变密度间隔。\

n = length(y);

if isempty(x)

x = 1:n;

end

FLAGS = zeros(1,n);

FLAGS(n) = 1;

% 通过标记1在FLAGS中找到分割点的位置。

num_sp = sum(FLAGS);

while 1

I = find(FLAGS==1);

numsp_bf = num_sp;

for k = 1:num_sp

if k==1

subx = x(1:I(k));

suby = y(1:I(k));

else

subx = x(I(k-1)+1:I(k));

su