% C-C方法
clc
begin_time = clock;
[tt,y]=ode45(@LorenzDifEqn1,[0,75],[0.1,0.1,0.1]);
xx = y(:,3)';
delt_t = 0.01;
t = 0:delt_t:75;
x = spline(tt, xx, t);
x = x(1000:1:4000);
tmax = 200;
mean_delt_S = zeros(1,tmax);
mean_S = zeros(1,tmax);
S_cor = zeros(1,tmax);
sigma = std(x);
for t = 1:1:tmax
    [mean_S_t,mean_delt_S_t,S_cor_t] = Cal_S_ave(x,t,sigma);
    mean_delt_S(t) = mean_delt_S_t;
    mean_S(t) = mean_S_t;
    S_cor(t) = S_cor_t;
    disp(t);
end
figure
hold
plot(mean_delt_S);
figure
hold
plot(mean_S);
figure
hold
plot(S_cor);
end_time = clock;
run_time = end_time-begin_time;
text = sprintf('run time = %d:%d:%d', run_time(4),run_time(5),ceil(run_time(6)));
disp(text);