关于基于最大Lyapunov指数的一步预测的几点疑问？

winterdij

我寻找到了Chaos Toolbox Ver.2.0中MainPre_by_Lya_1.m(基于最大Lyapunov指数的一步预测)源程序，但有几点不懂，在这里向高人请教？
源程序如下：
%skyhawk
clear all;

m=6;     %嵌入维数
N=80;    %预测后N个点

A=load('kj.txt');
P=26; % 北空的平均循环周期＝26

whl=A(:,4);
[whsl,lll]=size(whl);  

% lmd_1=lyapunov(m,m,whl,whsl);%求lyapunov指数
% lmd_mm=lmd_1(m);
for j=1:whsl            
    whlsj(j)=whl(j);
end   

fch=0;
for i=whsl-N+1:whsl         %预测后N个点
    [lmd_m,idx,min_d,idx1,min_d1]=lyapunov(m,whlsj,i-1,P);  
    [y(i),z(i)]=pre_by_lya(m,lmd_m,whlsj,i-1,idx,min_d);%预测第i+1个点  
   
    fch=fch+(y(i)-whl(i))*(y(i)-whl(i));
%     fch=fch+(z(i)-whl(i))*(z(i)-whl(i));
%     clear whlsj;

    iii=whsl-i     %显示进度
end

fch=sqrt(fch)/N

% for i=whsl-N+1:whsl
%     p(i-(whsl-N+1)+1)=y(i);
%     q(i-(whsl-N+1)+1)=z(i);
%     w(i-(whsl-N+1)+1)=whl(i);
% end

% kk=1:N;
% plot(kk,p,'r',kk,w)

yyy=[whl,y'];
save('kjyc.txt','yyy','-ASCII');

kk=1:whsl;
plot(kk,whl,'b',kk,y,'r')

问题如下：
1.A=load('kj.txt');
   这里的数据是原时间序列还是重构相空间以后的数据组？
   这个数列是否包括跟预测出的数据进行对比的原数据？
2.P=26; % 北空的平均循环周期＝26
   这里的P是否为进行FFT变换时预测的数列周期？
3.lmd_1=lyapunov(m,m,whl,whsl)
   这里的lyapunov函数从哪得到的，我的混沌工具包里没有这个函数，高人要是有的话，希望能给我一个，本人感激之至。
4.最后预测出的数据在哪个参数中表现出来？
   先谢谢了，我头一次做，自己不是很明白，让高人见笑了。

zhanman

怎么没人回帖.我也想知道为什么

柏莱

pre_by_lya和lyapunov(m,m,whl,whsl)函数都不知道,叫别人怎么看懂你的程序呀!还不如自己看看这个算法,自己编编看呢!
A=load('kj.txt'),感觉kj.txt里应该是原时间序列,而不是重构后的;
P就是时间序列的平均周期了;
plot(kk,whl,'b',kk,y,'r')表明预测值是y;

jhy14313

里面有很多源代码没有给出，需要自己编写

zao0123you

没有lyapunov函数，还得自己找，不知道lz解决没

水光潋滟

我个人认为A是重构后的数据组，P的值同意伯莱的看法，那两个没有的函数，我也很想知道怎么编写，唉

ym_qi

是呀，不知道怎么编写的呀