声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1299|回复: 1

[分形与混沌] 请帮我看一下这个程序

[复制链接]
发表于 2010-4-23 09:40 | 显示全部楼层 |阅读模式

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

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

x
Xh1=Xh_PingFang^(1/2)
Yh1=b*r*Xh1/(q*b-s*mu*Xh1^2);

Zh1=mu*Xh1*Yh1/b;
  
PingHenDian1=[0 0 0],
PingHenDian2=[Xh1 Yh1 Zh1],
PingHenDian3=[-Xh1 -Yh1 Zh1],
Xh=Xh1;Yh=Yh1;Zh=Zh1;
A=[  -sgm  tao+eps*Zh  eps*Yh
   r+s*Zh          -q    s*Xh
    mu*Yh       mu*Xh      -b],
v=eig(A),
d=det(A)
KL=(sgm+b+q)*(sgm*q+sgm*b+q*b-(mu*eps*Yh^2+(r+s*Zh)*(tao+eps*Zh)+mu*s*Xh^2))
det_A=-sgm*q*b+(tao+eps*Zh)*s*Xh*mu*Yh+eps*Yh*(r+s*Zh)*mu*Xh+q*eps*mu*Yh^2+b*(tao+eps*Zh)*(r+s*Zh)+sgm*mu*s*Xh^2
KLM=KL+det_A

A0=[-sgm  tao  0
       r  -q   0
       0   0  -b],
v0=eig(A0),
d0=det(A0)


gg=[];
ss=[];
for i=1:250
    ss(i)=i;
    s_temp=-200+(i-1)+0.44444444
    gg(i)=KL_M(s_temp);
end
figure(8);
plot(ss,gg,'-');                  %graph of KL_M
% xlabel('x'); ylabel('y'); zlabel('z');
%title('The Lorenz Chaos');

for sss=-3000 : 0.25 : 177
    if KL_M(sss)>=0
        sss
        KL_M(sss)
        break
    end
end
sss
KL_M(sss)                     %max sss=77

for sss=-3000 : -0.25 : -4000
    sss
    if KL_M(sss)>=0
        sss
        KL_M(sss)
        break
    end
end
sss
KL_M(sss)
请问这个图说明什么?画的是什么图?

图像.doc

37.5 KB, 下载次数: 10

回复
分享到:

使用道具 举报

发表于 2010-4-24 17:29 | 显示全部楼层
好难猜得............
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-19 13:51 , Processed in 0.057431 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表