修正r s分析法matlab程序,关于R/S分析程序用法

看到有人用下面的源代码，但为何有时候求出的值是负数，请各位帮忙，我是新接触最好能稍微具体一点讲解，谢谢。

function [logRS,logERS,V]=RSana(x,n,method,q)

% 用 R/S 方法分析间序列

% logRS 是 log(R/S).

% logERS 是 log(R/S)期望.

% V 是统计量.

% x 是时间序列.

% n 是这个数列的子集.

% method 可以取下列值

%  'Hurst' 为了Hurst-Mandelbrot变量

%  'Lo' 是Lo变量.

%  'MW' 是Moody-Wu变量.

%  'Parzen' 是Parzen变量.

% q 可以是任意值

%  a 是非0整数.

%  'auto' 是 Lo的默认值.

if nargin<1 | isempty(x)==1

error('你应该给出一个时间序列.');

else

% x 必须是变量

if min(size(x))>1

error('时间序列无效.');

end

x=x(:);

% N 是时间序列的长度

N=length(x);

end

if nargin<2 | isempty(n)==1

n=1;

else

% n 必须是一个变化的标量或矢量

if min(size(n))>