分岔图程序出图对称
function ydot=xiezhenqi(t,y,flag,vb)b=5e-8;
k1=5;
k3=15e12;
m=5e-12;
d=2e-6;
c0=1.875e-18;
w=0.5;
w0=sqrt(k1/m);
vac=0.21;
F=0.5*c0*(vb+vac*sin(w*w0*t)).^2/(d-y(1)).^2-0.5*c0*vb^2/(d+y(1)).^2;
ydot=[y(2)
(-b*y(2)-k1*y(1)-k3*y(1).^3+F)/m];
end
clear
vb=0:0.1:10;
c=0.5;
T=2*pi/c;
hold on
for j=1:length(vb)
=ode45('xiezhenqi',,,[],vb(j));
plot(vb(j),y(:,2),'linewidth',5)
end
不知道我的这个程序问题出在哪里？分岔图运行出的图关于x轴对称。。。
对于频闪法画分岔图还望各位大大请教！！
页:
[1]