关于lyapunov指数计算的一点困惑

俊博

最近在算时间序列的lyapunov指数，先用C_C法算出了时间延迟tau和M,接下来想用小数量法计算最大lyapunov指数
目前遇到的问题是C-C法算出来的时间延迟tau是个小数，带到小数量法程序里算不了，这是怎么回事呢，请高人指教，谢谢

无水1324

这个问题还没有仔细想，有时间看看

gghhjj

你用的是混沌时间序列分析与预测工具箱？

俊博

原帖由 gghhjj 于 2007-7-3 05:52 发表

                               
登录/注册后可看大图

你用的是混沌时间序列分析与预测工具箱？

我分别用的是C-Cmethod 和小数量法的程序，在研学上找的，

gghhjj

能够贴一下你找得小数量法的程序

俊博

我已经换了方法计算tau,现在tau值已经不是小数了，但是用小数量法计算出来的lyapunov指数有好几十，请大家帮忙看看
可能是taumax出了问题，我不知道这是怎么定义的，就用了原数，这个程序是计算别的方程的，我把数据部分改成自己的了，

%输入---------------------------------------------------

%-----------------------0.008----------------------------------------
load hr_a1_s0.008.tim%
x=hr(:,9)


m=3;
N=length(x);
tau=4;
%求平均周期P---------------------------------------------------------
dt=herui(2,1)-herui(1,1);
fs=1/dt; % 采样频率
dflag='mean';
nfft=1024;
window=256;
noverlap=window/2;
[p,f]=psd(x,nfft,fs,window,noverlap,dflag);
%mpf:
mpf=trapz(f,f.*p)/trapz(f,p);
P=1/mpf;


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



taumax = 30;            % 最大离散步进时间

           

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

Y1 = Lyapunov_rosenstein_2(x,tau,m,taumax,P,fs);

figure
subplot(211); plot(Y1); grid; xlabel('i'); ylabel('y(i)'); title('Logistic')
subplot(212); plot(diff(Y1)); grid; xlabel('n'); ylabel('slope');

linear_zone = [9:30]';  % 线性区域
F1 = polyfit(linear_zone,Y1(linear_zone),1);

Lyapunov_2 = F1(1)                  % 以 2 为底的指数
Lyapunov_e = Lyapunov_2/1.4427      % 以 e 为底的指数

octopussheng

最好用别的小数量算法程序算一下，比较一下结果！

小菜

请发到邮箱里
zhaodeqin2009@126.com

cleverxue

我用互信息法算的延迟时间ｔａｕ也是小数阿，应该怎么处理阿？