大家帮忙看下我的分岔图有什么问题

qq880416

为什么没有分叉啊 下面是我的程序 是否什么地方有错 期待解答
function  dy=dddd(t,y)
global  ff
dy=[y(2);-0.1*y(2)-48.704*y(1)-25.35e4*y(1)^3+ff*cos(2*t)];

clc
clear all
global f ww ff
ww=2;
f=0.1:0.1:10;
www=ww;
               
for i=1:length(f)
ff=f(i);

period=1/www;

tspan=[0:period/100:100*period];
[t,y]=ode45('dddd',tspan,[1e-8 0]);

plot(f(i),y(1000:100:end,2),'k.','MarkerSize',2);
hold on;
end

chenshumeionly

建议检查一下系数。

qq880416

回复 2 # chenshumeionly 的帖子

我也觉得可能系数有问题 但这个系数是试验参数，那是否说明这个系统不会发生分叉呢，主要想知道我的程序是否有问题

lgdlina

period=2*pi/www,你再试试看

kezairenjian

这个程序有问题，没有设初值。

live101

应该是周期取错了，同意4L的

chenshumeionly

回复 3 # qq880416 的帖子

你再检查一下有没有算错，或者你试着调整一下数据。