%C-C法子函数

function [Smean,Sdeltmean,Scor,tau,tw]=C_CMethod(data,max_d)               %原理参考：吕金虎-混沌时间序列分析及其应用——P67-70

N=length(data);                                                            %时间序列的长度

Smean=zeros(1,max_d);                                                      %初始化矩阵

Scmean=zeros(1,max_d);                                                     %C_CMethod 经使用后证实可以运行，不过其中有一个小错误，在初始化矩阵的位置要把Scmean改为Sdeltmean

Scor=zeros(1,max_d);

sigma=std(data);                                                           %计算序列的标准差

for t=1:max_d                                                              %计算Smean,Sdeltmean,Scor

S=zeros(4,4);

Sdelt=zeros(1,4);

for m=2:5

for j=1:4

r=sigma*j/2;

Xdt=disjoint(data,t);                                          %将时间序列data分解成t个不相交的时间序列

s=0;

for tau=1:t

N_t=floor(N/t);                                            %分成的子序列长度

Y=Xdt(:,tau);                                              %每个子序列

Cs1(tau)=0;                                                %计算C(1,N/t,r,t),相当于调用Cs1(tau)=correlation_integral1(Y,r)

for ii=1:N_t-1

for jj=ii+1:N_t

d1=abs(Y(ii)-Y(jj));                               %计算状态空间中每两点之间的距离,取无穷范数

if r>d1

Cs1(tau)=Cs1(tau)+1;

end

end

end

Cs1(tau)=2*Cs1(tau)/(N_t*(N_t-1));

Z=reconstitution(Y,m,1);                                   %相空间重构

M=N_t-(m-1);

Cs(tau)=correlation_integral(Z,M,r);                       %计算C(m,N/t,r,t)

s=s+(Cs(tau)-Cs1(tau)^m);                                  %对t个不相关的时间序列求和

end

S(m-1,j)=s/tau;

end

Sdelt(m-1)=max(S(m-1,  :  ))-min(S(m-1,  :  ));                            %差量计算

end

Smean(t)=mean(mean(S));                                                %计算平均值

Sdeltmean(t)=mean(Sdelt);                                              %计算平均值

Scor(t)=abs(Smean(t))+Sdeltmean(t);

end

for i=2:length(Sdeltmean)-1                                                %寻找时间延迟tau：即Sdeltmean第一个极小值点对应的t

if Sdeltmean(i)

tau=i;

break;

end

end

for i=1:length(Scor)                                                       %寻找时间窗口tw：即Scor最小值对应的t

if Scor(i)==min(Scor)

tw=i;

break;

end

end

%**************************************************************************

我的程序：

clc;clear;close all;

data=xlsread('Sunspot_17531_200112.xlsx');

d0=data(:,3);

figure;plot(d0,'r');grid on;xlabel('时间/月');

ylabel('太阳黑字数');title('1753年1月--2001年12月太阳黑子曲线图');

%--------------------------------------------------------------------------

[Smean,Sdeltmean,Scor,tau,tw]=C_CMethod(d0,400);                           %tau=31,tw=81,m=4这样的结果与书上m=8不相符合？？？？？

这个怎么办呢？

附件是所有涉及到的程序(其中xlsx为太阳黑子月度数据，数据的第三列是月度数据，第四列是方差)