碰摩转子分岔图matlab程序源

碰摩转子求分岔图的matlab程序源，该程序可以获得很好的分岔图，很好显示了系统的P-1、p-2……拟周期、混沌运动。整理下，希望对大家有所帮助

clear all;
clc;
global w;
global u;
%global E;
%global beta;
u=0.00007;
f0=250;
range=[2.5:0.05:4];
k=1;
YY=[];
z0=[0 0.001 1.27e-3 0.001];
for w=range
   period=2*pi/w;
   step=period/512;
    w
    %j=j+1;
     % discard the first *** periodic data;
     %除去前面2000个周期的数据，并将最后的结果作为下一次积分的初值
     tspan=[0:step:2000*period];
     %options=odeset（‘RelTol‘10^-3‘AbsTol‘10^-5）;
    [tY]=ode45（‘rubbing‘tspanz0）;
    z0=Y（end:）
    %j=1;
   % for i=2000:2200
        tspan=[2000*period:step:2100*period];
        %options=odeset（‘RelTol‘10^-3‘AbsTol‘10^-5）;
        [tY]=ode45（‘rubbing‘tspanz0）;
        YY（k:）=Y（1:512:end1）;   % get the omega data from every period end
        %j=j+1;               
        z0=Y（end:）;
        k=k+1;
    %end
end
bifdata=YY（:end-51:end）;
plot（rangebifdata‘k.‘‘LineWidth‘1）;
xlabel（‘频率比w/w0‘）;ylabel（‘x‘）;title（‘随频率比变化的分岔图‘）;
%xlabel（‘刚度比kc/k0‘）;ylabel（‘x‘）;title（‘随刚度比变化的分岔图‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1023  2012-10-28 17:04  碰摩转子分岔图matlab程序\Bifurcation.m
     文件         485  2012-10-28 16:34  碰摩转子分岔图matlab程序\rubbing.m
     目录           0  2012-11-07 14:12  碰摩转子分岔图matlab程序\