请教关联维数计算问题

buiesea · 发表于 2008-3-25 10:21

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

N=1000;
n=[0:N-1];
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;
以上是计算的时间序列和各参数，计算结果图如下:
我想问问怎么会出现这样的结果，根据图形怎么分析错误啊！谢谢了！

octopussheng · 发表于 2008-3-26 08:04

我怀疑是你的数据的问题，sin输入的点我从来没有算过，因此有此怀疑！呵呵！

buiesea · 发表于 2008-3-26 09:00

多谢主任！我到现在有些迷茫，关联维数计算的数据需要有什么要求吗？混沌时间序列应该是什么样的？比如说发动机在各种工作状态下的振动数据符合吗？

唳天鹄 · 发表于 2013-12-10 16:35

大仙，关联维数怎么求的，能教教我吗

gghhjj · 发表于 2014-2-21 05:08

唳天鹄发表于 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;

复制代码

账号		自动登录	找回密码
密码			我要加入

[分形与混沌] 请教关联维数计算问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。