有懂裂纹转子分岔图的吗？请帮忙看一下出现这种情况是什么原因。
http://a1.qpic.cn/psb?/ca0a0372-67fc-412a-acf1-6ad530fef523/m23V*GgW6U*gmx6XtIJTs0gm2wbnMSnorPojlY3d3iQ!/b/dOfWF4.zIwAA&bo=ewL*AQAAAAABAKA!&rf=viewer_4&t=5这是书上的正确的图。
http://a2.qpic.cn/psb?/ca0a0372-67fc-412a-acf1-6ad530fef523/debk0ZWCzMgkyHBGE.fAWDCZqGFNb3E8O3aZdGkLbzw!/b/dEvLqI8bJAAA&bo=MAKkAQAAAAABALA!&rf=viewer_4&t=5
这是我编程出来得到的图。
两张图差太多了，不知道是什么原因，会仅仅只是参数初值的影响吗？
下面附一下我的程序。
------------------------------------------------------------
function dx=jeffliewen(t,x)
global w
m=32.1;
c=2100;
e=0.002;
k=2.5e7;
g=9.8;
s0=m*g/k;
omega=sqrt(k/m);
beta=0;
theta=0;
k1=1-beta/pi+sin(4*beta)/(4*pi);%法向刚度
k2=1-beta/pi+2*sin(2*beta)/(3*pi)-sin(4*beta)/(12*pi);%切向刚度
k3=0;%耦合刚度u
k11=k1*(cos(theta))^2+k2*(sin(theta))^2-k3*sin(2*theta);%x 方向刚度
k22=k1*(sin(theta))^2+k2*(cos(theta))^2+k3*sin(2*theta);%y方向刚度
k33=(k1-k2)*sin(theta)*cos(theta)+k3*cos(2*theta);%x,y方向耦合刚度（以上刚度来自中性轴理论）
dx=[x(2);
-x(2)*c/(m*omega)-x(1)*k11/k-x(3)*k33/k+1+cos(theta+beta)*e*w^2/(s0*omega^2);
x(4);
-x(4)*c/(m*omega)-x(3)*k22/k-x(1)*k33/k+sin(theta+beta)*e*w^2/(s0*omega^2)];
----------------------------------------------------------------------
clear;
clc;
hold on;
global w
f=10:3:1200
for i=1:length(f)
disp(f(i));
w=f(i);
T=2*pi;
x0=;
tspan=0:T/100:100*T;
=ode45('jeffliewen',tspan,x0);
x0=x(end,:);
plot(f(i),x(4000:100:end,1),'markersize',5);%画分岔图
xlabel ('转速w');
ylabel ('位移x1');
end http://a1.qpic.cn/psb?/V12ilAMs3XYRKq/UfxdsvOcs0aZp.*fjiyF2Z9qcx0l8OtJMD1NaUEpO60!/b/dOfWF4.zIwAA&bo=ewL*AQAAAAABAKA!&rf=viewer_4&t=5
正确的。 这图片怎么老发不上来。。。反正我的图根本就没有分岔，出来像一个牛角，由细到宽呈放射性扩大，这只跟参数初值有关吗》？ 我现在也在研究裂纹转子的非线性，请问你看的是什么书啊，发现关于裂纹转子方面的资料很少，可以分享一下吗
页:
[1]