用matlab画分岔图

用matlab画分岔图,第1张

给你画了一个,但是跟你要的不完全一样。我想思路应该是对销嫌粗的。我猜可能是你给的初值跟图不一样吧。把下面的程序保存成一个文件,运行即可。

程序的基本思路是:从初值出发,先计算y对时间t的导数,得到y的变化速度;者掘求出x对t的导数,得到x随时间的亏镇变化速度。然后时间前进一小步h=0.002,同时x和y也都发生一个小的变化,即在原来数值基础上加上时间步长h和变化速度的乘积。如果x到达1或者-1,强制让y发生变化。这个过程不断迭代,得到结果。

clear

clc

gamma=2

r=0.9

N=400000

x=zeros(1,N)

y=zeros(1,N)

t=zeros(1,N)

h=0.002

for i=2:N

t(i)=t(i-1)+h

y(i)=y(i-1)+gamma*sin(t(i-1))*h

x(i)=x(i-1)+y(i-1)*h

if abs(x(i))>=1

y(i)=-r*y(i)

x(i)=x(i-1)+y(i-1)*h

end

end

figure

plot(x,y,'.','markersize',1)

axis([-1 1 -3 3])

xlabel('Displacement')

ylabel('Velocity')

clcclear all

global r

r=0.11

[T,X]=ode15s(@eq,[0,10],[0 0])

plot(T,X(:,1),'r--')

hold on

plot(T,X(:,2)) function dx=eq(t,x)

global r

dx=[x(2),-0.5*x(2)+x(1)^3-x(1)^5-r*cos(t)]

end

:本程序介绍了一种衫燃胡分岔图的绘制或拦方法与对段敬象,...bifurcation diagramdrawing method with an object,...·用MATLAB编写的svm源程序,


欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/yw/12355047.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-24
下一篇 2023-05-24

发表评论

登录后才能评论

评论列表(0条)

保存