声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2142|回复: 2

[稳定性与分岔] 求大牛帮忙怎么分岔图就是直线

[复制链接]
发表于 2013-5-14 16:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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
运算后只是一条直线,没得到分岔啊。求大牛们帮忙看看。不胜感激。
回复
分享到:

使用道具 举报

 楼主| 发表于 2013-5-14 16:20 | 显示全部楼层
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=[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)+c4*y(9);y(1);y(2);y(3);d1*y(10)+d2*y(11)+d3*y(12);e1*y(10)+e2*y(11)+e3*y(12);f1*y(10)+f2*y(11)+f3*y(12)+f4*y(3);y(7);y(8);y(9)];
 楼主| 发表于 2013-5-15 10:19 | 显示全部楼层
plot(uu(j),y(1:1:10,6),'linewidth',5)
想让y坐标值更小,怎么设置?y(1:110)
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 22:26 , Processed in 0.048000 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表