求大牛帮忙怎么分岔图就是直线
tic;uu=100:1:1000;
t0=0;tf=10;
hold on
for j=1:length(uu)
=ode45('zf1',,,[],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=;
end
运算后只是一条直线，没得到分岔啊。求大牛们帮忙看看。不胜感激。
zf1.m
function rdot1=zf1(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;c4=15347*uu/497;d1=-4.2127e12/78.5;d2=2.1022e12/78.5;d3=-2.1325e11/78.5;e1=2.1022e12/78.5;e2=-6.4405e12/78.5;e3=1.132e12/78.5;f1=2.1325e11/162.269;f2=-1.132e12/162.269;
f3=(2.2857e11+237390*uu^2)/162.269;f4=-6420.5*uu/162.269;
rdot1=; plot(uu(j),y(1:1:10,6),'linewidth',5)
想让y坐标值更小，怎么设置？y(1:110)
页:
[1]