求计算关联维数的MATLAB程序
那位大虾有谢谢!急.
我的邮箱ndhutu_1981@163.com
[ 本帖最后由 eight 于 2008-3-2 17:36 编辑 ] wo ye shi zai deng a
我的邮箱是keen_2004@163.com
回复：（ndhutu）[求助]计算关联维数的MATLAB程序
是不是这个？function =G_P(data,N,tau,min_m,max_m,ss)
% the function is used to calculate correlation dimention with G-P algorithm
% data:the time series
% N: the length of the time series
% tau: the time delay
% min_m:the least embedded dimention m
% max_m:the largest embedded dimention m
% ss:the stepsize of r
%skyhawk
for m=min_m:max_m
Y=reconstitution(data,N,m,tau);%reconstitute state space
M=N-(m-1)*tau;%the number of points in state space
for i=1:M-1
for j=i+1:M
d(i,j)=max(abs(Y(:,i)-Y(:,j)));%calculate the distance of each two
end %points in state space
end
max_d=max(max(d));%the max distance of all points
d(1,1)=max_d;
min_d=min(min(d));%the min distance of all points
delt=(max_d-min_d)/ss;%the stepsize of r
for k=1:ss
r=min_d+k*delt;
C(k)=correlation_integral(Y,M,r);%calculate the correlation integral
ln_C(m,k)=log(C(k));%lnC(r)
ln_r(m,k)=log(r);%lnr
fprintf('%d/%d/%d/%d\n',k,ss,m,max_m);
end
plot(ln_r(m,:),ln_C(m,:));
hold on;
end
fid=fopen('lnr.txt','w');
fprintf(fid,'%6.2f %6.2f\n',ln_r);
fclose(fid);
fid = fopen('lnC.txt','w');
fprintf(fid,'%6.2f %6.2f\n',ln_C); fclose(fid); 是的
非常谢谢happy！ plot(ln_r(m,:),ln_C(m,:));
那位大侠知道，里面的笑脸代面什么意思 原帖由 realhappy 于 2006-8-5 11:20 发表
plot(ln_r(m,:),ln_C(m,:));
那位大侠知道，里面的笑脸代面什么意思
笑脸表示的是
:)
这个是显示问题，代码没有引用或者取消smiles造成的
[ 本帖最后由 eight 于 2007-3-30 17:23 编辑 ] plot(ln_r(m,,ln_C(m,);
用冒号＋右括号后还是提示有问题，不知道怎么回事。 这是干什么的一个程序呀 求混沌时间序列的关联维。 这是利用GP算法实现混沌时间序列关联维的求解。先看看能否成行，谢谢 这个程序我也有，可是我不知道怎么用？从哪里开始读取数据啊？ 原帖由 bella 于 2006-11-7 16:53 发表
这个程序我也有，可是我不知道怎么用？从哪里开始读取数据啊？
这是一个运算的function，数据(data)的导入需要你自己写程序读入或者以其他方式获取 原帖由 bella 于 2006-11-7 16:53 发表
这个程序我也有，可是我不知道怎么用？从哪里开始读取数据啊？
在网上查混沌时间序列数据，在国外的网站上有，下面就是一个 请问一下函数中的ss参数是什么意思，怎么确定？ 这里面的reconstitution 和correlation_integral函数没有, 无法执行的
页:
[1]
2