关于基于最大Lyapunov指数的一步预测的几点疑问？
我寻找到了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);
=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个点
=lyapunov(m,whlsj,i-1,P);
=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=;
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.最后预测出的数据在哪个参数中表现出来？
先谢谢了，我头一次做，自己不是很明白，让高人见笑了。 怎么没人回帖.我也想知道为什么 pre_by_lya和lyapunov(m,m,whl,whsl)函数都不知道,叫别人怎么看懂你的程序呀!还不如自己看看这个算法,自己编编看呢!
A=load('kj.txt'),感觉kj.txt里应该是原时间序列,而不是重构后的;
P就是时间序列的平均周期了;
plot(kk,whl,'b',kk,y,'r')表明预测值是y; 里面有很多源代码没有给出，需要自己编写 没有lyapunov函数，还得自己找，不知道lz解决没 我个人认为A是重构后的数据组，P的值同意伯莱的看法，那两个没有的函数，我也很想知道怎么编写，唉 是呀，不知道怎么编写的呀
页:
[1]