请教高人：这个分叉图的程序有问题吗

baiyun_0810475

我想画一个方程的分插图，但运行起来特困难，不知道是不是程序有问题，我对编程不是很懂，麻烦朋友指点一下，非常感谢！
主程序是
%clear;
for gamma=20:0.1:40
options=odeset('RelTol',1e-7);
tt=2*pi/1.2;
%[t,x]=ode45(@dafen,[0:tt/100:80*tt],[0,0],options,[],ff);
[t,x]=ode45(@c3eqdae6,[0:tt/100:500*tt],[1;1;1],options,[],gamma);
%pojialai
i=40000:100:50000;
plot(gamma,log(x(i,1)),'k.');
hold on;
end
% axis([0 3.2 -8 0]);


调用程序是
function dx=c3eqdae6(t,x,flag,gamma)
dx=[gamma*x(1)-0.5*x(1)^2-x(1)*x(2)+1e-4*x(2)^2-x(1)*x(3)-1e-3;
    x(1)*x(2)-1e-4*x(2)^2-10*x(2)+1e-3;
    -x(1)*x(3)+1e-3+16.5*x(3)-0.5*x(3)^2];

[ 本帖最后由 咕噜噜 于 2009-10-16 09:22 编辑 ]

无水1324

你的系统是自治系统的，而你采用的是非自治系统的频闪法，所以得到不合理的图示很正常的。
建议：1、从定义出发选取截面，2、参考论坛上的最大值法做图

baiyun_0810475

好的，谢谢您的指导，我再好好修改一下！