cqupenghao的留言板 - 声振论坛

flybett3952 2013-7-21 10:13

老师你好！我最近在做分数阶系统的指数问题，在论坛上看见你有时间延迟tau和嵌入维数m的相关程序，能给我提供一份吗？跪谢！我的邮箱是fly0164@126.com.

bohua1208 2011-11-12 15:36

cqupenghao: 你说的那个问题，应该可以判断是混沌运动，关于你计算的Lyapunov指数值的大小，好像与你的程序设置的参数有一定关系，你能把你计算的程序借我参考一下吗？

以下是关于指数谱的程序

clear all;
clc;
tic,
LEmax=[];
f=2.85:0.001:3.25;
for j=1:length(f)
j
yinit = [-0.1,0.5];
orthmatrix = [1 0 ;
            0 1 ];
y = zeros(6,1);
% 初始化输入
IC(1:2) = yinit;
IC(3:6) = orthmatrix;
InitialTime = 0; % 时间初始值
TimeStep = 1e-2; % 时间步长
wholetimes = 2e4; % 总的循环次数
UpdateStepNum = 10; % 每次演化的步数
iteratetimes = wholetimes/UpdateStepNum; % 演化的次数
DiscardItr=200;
Iteration=UpdateStepNum*TimeStep;
DiscardTime=DiscardItr*Iteration+InitialTime;

T1=InitialTime;
T2=T1+Iteration;
TSpan=[T1:TimeStep:T2];
options=odeset('RelTol',1e-5,'AbsTol',1e-6);
mod = zeros(1,2);
d=2;
Sum=zeros(1,d);
n1=0;
k1=0;
xData=[];
yData=[];
F=f(j);
for i=1:iteratetimes
   n1=n1+1;
[t,X] = ode45('zhgz',TSpan,IC,options,F);
% 取积分得到的最后一个时刻的值
[rX,cX]=size(X);
L=cX-d*d;

for i=1:d
   m1=L+1+(i-1)*d;
   m2=m1+d-1;
   A(:,i)=(X(rX,m1:m2))';
end
y0 = TwoGS(A);

% 取两个向量的模
mod(1) = sqrt(y0(:,1)'*y0(:,1));
mod(2) = sqrt(y0(:,2)'*y0(:,2));
y0(:,1) = y0(:,1)/mod(1);
y0(:,2) = y0(:,2)/mod(2);
if T2>DiscardTime
      Q0=y0;
   else
      Q0=eye(d);
   end

  if (T2>DiscardTime )
   k1=k1+1;
   T=k1*Iteration;
   TT=n1*Iteration+InitialTime;
   Sum=Sum+log(abs(mod));
   lambda=Sum/T;

   Lambda=fliplr(sort(lambda));
   [rxD,cxD]=size(xData);
   [ryD,cyD]=size(yData);

      xData=[xData;TT];
      yData=[yData;lambda];
   end
   ic=X(rX,1:L);
   T1=T1+Iteration;
   T2=T2+Iteration;
   TSpan=[T1:TimeStep:T2];

IC=[ic(:);y0(:)];
end
LEmax=[LEmax;yData(end,1)];
end
plot(f,LEmax,'k-');
xlabel('F','fontsize',14);
ylabel('\lambda_m_a_x','fontsize',14);
box on
toc

bohua1208 2011-11-11 19:52

程序计算的三个Lyapunov分别为0.05，0 ，-0.55，分形维数为2.09，最大lyapunov指数才0.05很小啊，能说明是混沌运动吗，或者有可能是拟周期运动吗？

bohua1208 2011-4-16 16:11

有Lyapunov指数随参数变化的程序吗，给我发一下，672504435@qq.com 不胜感激

yanzi 2011-3-22 15:41

cqupenghao: 你好，请问你有求Kolmogorov熵的程序吗？最近在搞这方面的，急求啊

不好意思没有啊

553535912 2010-11-16 11:47

你好，我也在研究CC方法和最小数据量方法，可以加我Q吗？有些问题想请教一下，我的553535912，匹马西风

lilili606 2010-10-10 18:52

cqupenghao: 你QQ是好多啊？？昵称是什么啊？？？？？？？？、

394459138 独自等待

lilili606 2010-9-26 12:35

cqupenghao: 化工学院的

呵呵和我们学院都住四舍的~

lilili606 2010-9-23 12:38

我是资环力学系的你是？

cqupenghao的个人空间 http://home.vibunion.com/?157013 [收藏] [复制] [分享] [RSS]

留言板

cqupenghao