各位学长帮忙

吉祥如意

各位学长，请帮帮忙， 我想得到 当a 从0.01变化到1时，此系统的分岔图程序，我在论坛已搜索过分岔图的相关程序，可我修改后要不就是运行不了，要不就是运行结果不对，:@( 请各位学长帮忙编一下此程序，非常感谢！请各位学长帮帮忙！

无水1324

将你修改的程序贴上来，以便别人帮你修改

吉祥如意

我的程序是：
function xdot=ddddd(t,x,flag, gamma)
omega = 1;
xdot=[x(2); -gamma*x(2)-x(1)-x(1)*x(1)+x(1)^3+100*cos(omega*t)];

clear;clc;close all;
gamma=0.1:0.001:1.0;
omega=1:0.01:2;
for h=1:length(gamma)
    T=2*pi/omega(h);
    [t,x]=ode45('ddddd',[0:T/100:220*T],[0 0],[], gamma(h));
    plot(gamma(h),x(10000:100:end,2),'k.');hold on
end

无水1324

function xdot=ddddd(t,x,flag, gamma)
omega = 1;
xdot=[x(2); -gamma*x(2)-x(1)-x(1)*x(1)+x(1)^3+100*cos(omega*t)];

clear;clc;close all;
gamma=0.1:0.001:1.0;
omega=1;
for h=1:length(gamma)
    T=2*pi/omega;
    [t,x]=ode45('ddddd',[0:T/100:220*T],[0 0],[], gamma(h));
    plot(gamma(h),x(10000:100:end,2),'k.');hold on
end

吉祥如意

谢谢学长！