[讨论]R/S算法设计

[tasmhust]

最近在实验室做网络流量的自相似分析，其中有种R/S算法，需要用Matlab实现。以下是我的代码：

clear all;
clc;
X=1:1000;
L=length(X);
N=10;
while N<L&j<L/10
d=ceil(L/N);
for i=1:d
if(i*N<L)
Xni=X([(i-1)*N+1:i*N])
A=mean(Xni)
B=var(Xni)
Yi=cumsum(Xni)-(1:N)*A
rs(i)=(max(Yi)-min(Yi))/B^0.5
end
end
for j=1:L/10
RSn(j)=mean(rs);
end
N=N+10;
end
plot(log(10:10:L),log(RSn));

得到的图形是条水平线，我想问下，这个算法有没有问题，可以再优化吗？