声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2249|回复: 4

[稳定性与分岔] 非线性支承转子系统分岔图,求教,指点迷津

[复制链接]
发表于 2014-12-1 16:23 | 显示全部楼层 |阅读模式

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

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

x

参考别人论文里面的程序绘制的分岔图,是什么原因?


%非线性弹性支承转子系统
function dy=diffequat1(t,y)
global theta e kz k w W r2 M m g D beta;
e=0.0002;
theta=1.415263;
r2=1000.00;
w=0.25;
W=500;
kz=19304448;
M=50;
m=10;
D=0.008;
g=9.800;
k=1930444.8;
beta=1;
dy=[y(2);
e*kz/M/g*cos(t-theta)-r2*kz*w^2/M/g*sin(theta)-kz/M/W^2*cos(theta)-2*D*(kz/M)^0.5/W*(y(2)-y(4))-kz/M/W^2*(y(1)-y(3));
y(4);
-r2*kz*w^2/W^2/m/g*sin(theta)-kz/m/W^2*cos(theta)+2*D*M*(kz/M)^0.5/W/m*(y(2)-y(4))+kz/m/W^2*(y(1)-y(3))-2*D*M*(kz/M)^0.5/W/m*y(4)-k/m/W^2*y(3)-beta*k/m/W^2*y(3)^3;
y(6);
e*kz/M/g*sin(t-theta)-r2*kz*w^2/M/g*cos(theta)-kz/M/W^2*cos(theta)-2*D*(kz/M)^0.5/W*(y(6)-y(8))-kz/M/W^2*(y(5)-y(7));
y(8);
-r2*kz*w^2/W^2/m/g*cos(theta)-kz/m/W^2*sin(theta)+2*D*M*(kz/M)^0.5/W/m*(y(6)-y(8))+kz/m/W^2*(y(5)-y(7))-2*D*M*(kz/M)^0.5/W/m*y(8)-k/m/W^2*y(7)-beta*k/m/W^2*y(7)^3];


绘制分岔图程序:

clear all;

clc;

W=100:10:1300;

for h=1:length(W)

T=2pi;

[t,x]=ode45(‘diffequat1’,[1:T/100:1000*T],[11 1 1 1 1 1 1]),[],W(h));

plot(W(h),x(80000:100:end,1),’k.’);

hold on;

end

绘制的图出来感觉是不对的。


                               
登录/注册后可看大图


本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2014-12-9 09:58 | 显示全部楼层
看不到图?
难道是网速的原因
 楼主| 发表于 2014-12-9 11:16 | 显示全部楼层
帖子有问题,又不能删除
发表于 2014-12-11 21:51 | 显示全部楼层
发表于 2014-12-22 11:09 | 显示全部楼层
同问同问同问同问
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-20 01:19 , Processed in 0.080321 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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