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))>1
error('n 必须是一个变化的标量或矢量.');
end
% n 必须是个整数
if n-round(n)~=0
error('n 必须是个整数.');
end
% n 必须是确定
if n<=0
error('n 必须是确定.');
end
end
评论9