请教关联维数计算问题
N=1000;n=;
Fs=500;
t=n/Fs;
Lag=100;
data=sin(2*pi*10*t);
tau=7;ss=50;
N=length(data)
min_m=2;max_m=10;
以上是计算的时间序列和各参数，计算结果图如下:
我想问问怎么会出现这样的结果，根据图形怎么分析错误啊！谢谢了！ 我怀疑是你的数据的问题，sin输入的点我从来没有算过，因此有此怀疑！呵呵！ 多谢主任！我到现在有些迷茫，关联维数计算的数据需要有什么要求吗？混沌时间序列应该是什么样的？比如说发动机在各种工作状态下的振动数据符合吗？ 大仙，关联维数怎么求的，能教教我吗{:{23}:} 唳天鹄 发表于 2013-12-10 16:35
大仙，关联维数怎么求的，能教教我吗
具体理论可以找找书籍或者文献，一般用G-P算法比较多
程序你可以找陆振波的混沌工具包一类的，也可以自己编程
下面是一个简单的关联维数计算程序：
% estimation of the correlation dimension
% Grassberger and Procaccia algorithm
r=0.50;
m=22;
t=5;
dr=0.10;
keyboard; % input r,m,t,dr
n=length(fname);
mm=n-(m-1)*t;
eeaa=zeros(mm,m);
for i=1:mm
for j=1:m
eeaa(i,j)=fname(i+(j-1)*t);
end;
end;
fenmu=mm*(mm-1);
mmm=250/r;
eebb=zeros(mmm,4);
h=1;
for ij=r:r:250
k=0;
for i=1:mm-1
for j=i+1:mm
hh=norm(eeaa(i,:)-eeaa(j,:),inf);
if hh<ij
k=k+2;
end;
end;
end;
crm=k/fenmu;
eebb(h,1)=k;
eebb(h,2)=crm;
eebb(h,3)=log2(crm);
eebb(h,4)=log2(ij);
h=h+1;
end;
页:
[1]