|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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)
请问这个图说明什么?画的是什么图? |
|