cscc_method.m

CC_method:CC法

% 此程序用来测试CC_method
% 2008-12-01
% zhangli

clear all
clear all


%利用方程获得
% 产生 Lorenz 时间序列
% dx/dt = sigma*(y-x)
% dy/dt = r*x - y - x*z
% dz/dt = -b*z + x*y
sigma=16;          % Lorenz 方程参数
b=4;                              
r=45.92;                                      
y=[-1,0,1];        % 起始点 (1 x 3 的行向量)
h=0.01;            % 积分时间步长
k1=10000;          % 前面的迭代点数
k2=3000;           % 后面的迭代点数
Z=LorenzData(y,h,k1+k2,sigma,r,b);
X=Z(k1+1:end,1);

max_d=200;         % 最大延迟时间

% 调用C_CMethod_inf，求tau
tic
[Smean_inf,Sdeltmean_inf,Scor_inf,tau_inf,tw_inf]=C_CMethod_inf(X,max_d);
toc
tau_inf
tw_inf
% 相关作图
figure('name','CC法求时间延迟');
plot(1:max_d,Smean_inf,'-b');hold on;
plot(1:max_d,Sdeltmean_inf,'-*c');hold on;
plot(1:max_d,Scor_inf,'-m');hold on;
plot(1:max_d,zeros(1,max_d),'r');
title('C_CMethod_inf');xlabel('Lag');
legend('S(t)平均值','ΔS(t)平均值','Scor_inf');
% 将数据保持下来
fid=fopen('Smean_inf.txt','w');
fprintf(fid,'%f\n',Smean_inf);
fclose(fid);
fid=fopen('Sdeltmean_inf.txt','w');
fprintf(fid,'%f\n',Sdeltmean_inf);
fclose(fid);
fid=fopen('Scor_inf.txt','w');
fprintf(fid,'%f\n',Scor_inf);
fclose(fid);