混沌理论 发表于 2010-4-23 09:40

请帮我看一下这个程序

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

Zh1=mu*Xh1*Yh1/b;

PingHenDian1=,
PingHenDian2=,
PingHenDian3=[-Xh1 -Yh1 Zh1],
Xh=Xh1;Yh=Yh1;Zh=Zh1;
A=[-sgmtao+eps*Zheps*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=[-sgmtao0
       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)
请问这个图说明什么？画的是什么图？

无水1324 发表于 2010-4-24 17:29

好难猜得............
页: [1]
查看完整版本: 请帮我看一下这个程序