俊博 发表于 2007-7-2 11:06

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

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

无水1324 发表于 2007-7-2 15:50

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

gghhjj 发表于 2007-7-3 05:52

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

俊博 发表于 2007-7-3 09:29

原帖由 gghhjj 于 2007-7-3 05:52 发表 http://www.chinavib.com/forum/images/common/back.gif
你用的是混沌时间序列分析与预测工具箱？

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

gghhjj 发表于 2007-7-4 07:45

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

俊博 发表于 2007-7-4 15:52

我已经换了方法计算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;
=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 = ';% 线性区域
F1 = polyfit(linear_zone,Y1(linear_zone),1);

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

octopussheng 发表于 2007-7-5 08:03

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

小菜 发表于 2009-4-24 23:36

谁有混沌时间序列分析与预测工具箱

请发到邮箱里
zhaodeqin2009@126.com

cleverxue 发表于 2009-4-29 15:11

我用互信息法算的延迟时间ｔａｕ也是小数阿，应该怎么处理阿？
页: [1]
查看完整版本: 关于lyapunov指数计算的一点困惑