达芬方程的分叉图
2018-4-19 23:11
阅读:5429
以达芬振子为例
dafen = @(t,x,flag,force)[x(2);force*cos(1.2*t)-x(1)^3+x(1)-0.3*x(2)];
ff=0.23:0.002:0.40;%参数
axis([0.22,0.4,0.2,0.6]);
options=odeset('RelTol',1e-7);tt=2*pi/1.2;hold on
for j=1:length(ff)
[t,x]=ode45(dafen,[0:tt/100:80*tt],[0,0],options,[],ff(j));
i=2000:100:3000;
plot(ff(j),x(i,2),'.')
end
转载本文请联系原作者获取授权,同时请注明本文来自王福昌科学网博客。
链接地址:https://wap.sciencenet.cn/blog-292361-1109869.html?mobile=1
收藏
当前推荐数:0
推荐到博客首页
网友评论0 条评论