转子碰摩分岔，附程序，请教各位指错

xiaoshihanlan

请教各位前辈老师，刚接触分岔这块，做一个Jeffcott转子系统碰摩的分岔图，跟着书上方程编写的matlab程序，就是出不来正确的图，请教各位帮忙看看是哪里错了？
方程：function dx=pengmo_1(t,x,flag,gamma)
beta=3.0;
alfa=0.5;
epsilon=0.12;
f=0.12;
deta=0.16;
u=0.136;
w0=50*pi;
g=9.8;
w=w0*gamma;
dx=[x(2);
-2*epsilon*x(2)-x(1)-beta*(1-1/sqrt(x(1)^2+x(3)^2))*(x(1)-f*x(3))-alfa*(x(1)^2+x(3)^2)*x(1)+u*(gamma)^2*cos(w0*gamma*t);
x(4);
-2*epsilon*x(4)-x(3)-beta*(1-1/sqrt(x(1)^2+x(3)^2))*(f*x(1)+x(3))-alfa*(x(1)^2+x(3)^2)*x(3)+u*(gamma)^2*sin(w0*gamma*t)-g/((w^2)*deta)];
主程序：w0=50*pi;
gamma=1.10:0.01:1.90;
x0=[0.1;0;0.1;0];
for h=1:length(gamma)
     disp(gamma(h));
      T=2*pi/(w0*gamma(h));
[t x]=ode45('pengmo_1',[0:T/100:200*T],x0,[],gamma(h));
plot(gamma(h),x(19500:200:end,1),'k.');hold on
end

xiaoshihanlan

hsfy919

回复 2 # xiaoshihanlan 的帖子

1、既然是碰磨问题，那么判断是否接触的程序呢？
2、主程序中应该是用闪频法画的庞加莱截面，那么周期取点也不对。
3、你是在那本书上看到的？

xiaoshihanlan

回复 3 # hsfy919 的帖子

书上只给了微分方程，其它是自己写的，还要先判断是否发生碰摩？不懂。

伤痕累累

回复 4 # xiaoshihanlan 的帖子

你这是哪本书上的？碰摩的程序根本就不对

xiaoshihanlan

回复 5 # 伤痕累累 的帖子

在一本转子故障实验分析书上看到的，刚开始接触就跟着人家的方程做了,毕设设计做转子碰摩，现在就是程序搞不懂啊，请前辈帮帮忙

伤痕累累

回复 7 # xiaoshihanlan 的帖子

你是本科是吧？

xiaoshihanlan

回复 8 # 伤痕累累 的帖子

是的，在做毕业设计