www.pudn.com > Chaos_Prediction.rar > Main_BoxDimension_TS.m

% 计算时间序列的盒维数 - 主函数
% 使用平台 - Matlab7.0
% 作者：陆振波，海军工程大学
% 欢迎同行来信交流与合作，更多文章与程序下载请访问我的个人主页
% 电子邮件：luzhenbo@yahoo.com.cn
% 个人主页：http://luzhenbo.88uu.com.cn

clc
clear all
close all

%--------------------------------------------------------------------------
% produce the FBM time series，Set parameter H to 0.6 and sample length

rand('state',sum(100*clock))
H = 0.6; lg = 1024;
s = wfbm(H,lg);                    % 时间序列

figure; plot(s,'.-')

%--------------------------------------------------------------------------
% 计算盒维数

partition = 2^7;                   % 每一维坐标上的分割数
[log2C,log2r] = BoxDimension_TS(s,partition);

figure
plot(log2r,log2C,'bo'); xlabel('log2(r)'); ylabel('log2(C(r))'); hold on;

%--------------------------------------------------------------------------
% 确定线性区域

Linear = 1:length(log2C);
par = polyfit(log2r(Linear),log2C(Linear),1);
Df = par(1)              % 盒维数

log2C_estimate = polyval(par,log2r(Linear),1);
plot(log2r(Linear),log2C_estimate,'r-'); hold off;