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

qq880416 · 发表于 2011-11-19 12:53

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

为什么没有分叉啊下面是我的程序是否什么地方有错期待解答
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 · 发表于 2011-11-19 20:49

建议检查一下系数。

qq880416 · 发表于 2011-11-20 12:10

回复 2 # chenshumeionly 的帖子

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

lgdlina · 发表于 2011-11-22 10:58

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

kezairenjian · 发表于 2011-11-22 11:16

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

live101 · 发表于 2011-11-22 20:57

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

chenshumeionly · 发表于 2011-11-24 11:47

回复 3 # qq880416 的帖子

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

账号		自动登录	找回密码
密码			我要加入

[稳定性与分岔] 大家帮忙看下我的分岔图有什么问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。