马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
tic;
uu=100:1:1000;
t0=0;tf=10;
hold on
for j=1:length(uu)
[t,y]=ode45('zf1',[t0,tf],[0,0,0,0,0,0,0,0,0,0,0,0],[],uu(j));
plot(uu(j),y(1:1:10,6),'linewidth',5)
title('bifurcation')
xlabel('参数 uf')
ylabel('w位移')
end
toc
zf1.m
function rdot=zf(t,y,flag,uu)
a1=-3.2341*10^11;a2=1/78.5*(10.511e11);a3=-1/78.5*(1.066e11);b1=1/78.5*(1.0511e12);b2=1/78.5*(5.8546e12);b3=-1/78.5*(1.1319e12);c1=1/497*(1.0611e11);c2=1/497*(1.1319e12);
c3=(2.277e11+141850*(uu)^2)/497;
rdot=[a1*y(4)+a2*y(5)+a3*y(6);b1*y(4)-b2*y(5)-b3*y(6);c1*y(4)-c2*y(5)+c3*y(6);y(1);y(2);y(3)];
end
运算后只是一条直线,没得到分岔啊。求大牛们帮忙看看。不胜感激。
|