HOPF分岔图

21172485

下午无聊,画了一个系统的HOPF分岔
大家看着玩吧

无水1324

这个市用什么东西画的？做件好事，把程序也发上来吧

21172485

以下是MATLAB程序,很简单.我就不写注释了
clc
clear

axis([-2 2 -2 2])
hold on
fai0=pi/12
fai2=pi*(2/3)
fai3=pi/3
fai4=pi
fai5=1.5*pi
fai6=pi*(4/3)
fai7=pi*(-1/6)
fai8=0.01*pi*(0.00002)

u=0.5
r=0.001
t=0.01:0.01:20
d=sqrt(r*r+(abs(u)-r*r)*exp(-2*u*t))
p=(sqrt(abs(u))*r)./d
x=p.*cos(t)
y=p.*sin(t)
plot(x,y,'k')
hold on


x=p.*cos(t+fai4)
y=p.*sin(t+fai4)
plot(x,y,'c')

u=0.5
r=2
d=sqrt(r*r+(abs(u)-r*r)*exp(-2*u*t))
p=(sqrt(abs(u))*r)./d
x=p.*cos(t)
y=p.*sin(t)
plot(x,y,'r')
hold on

%
x=p.*cos(t+fai5)
y=p.*sin(t+fai5)
plot(x,y,'m')

x=p.*cos(t+fai2)
y=p.*sin(t+fai2)
plot(x,y,'k')

x=p.*cos(t+fai6)
y=p.*sin(t+fai6)
plot(x,y,'g')

x=p.*cos(t+fai7)
y=p.*sin(t+fai7)
plot(x,y,'y')

x=p.*cos(t+fai4)
y=p.*sin(t+fai4)
plot(x,y)

x=p.*cos(t+fai3)
y=p.*sin(t+fai3)
plot(x,y,'m')

yzsldj

这个不像是分岔图啊，有点像相空间轨道图

无水1324

谢谢。我原以为是一个函数画出来的

octopussheng

呵呵，太强了，今天也算学习到了！

gxc8208

请问能把您的分岔图反映的振动系统方程给出来么？

kangarooli

这个在刘延柱编的非线性振动那本书里出现过，是一个二元系统的hopf分岔相轨迹

kangarooli

这个在刘延柱编的非线性振动那本书里出现过，是一个二元系统的hopf分岔相轨迹