matlab ode45画延时,ode45求解微分方程并画图

clc

clear all

%画系统的庞加莱图主程序；

T=2*pi;%周期；

tspan=0:T/100:200*T;%求解区间；

y0=[0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005];%初始值；

options=odeset('RelTol',10^-3,'AbsTol',10^-3);%选项优化；

[t,y]=ode45('fun',tspan,y0,options);

Y1=y(15000:end,1);

Y2=y(15000:end,2);

Y3=y(15000:end,3);

Y4=y(15000:end,4);

Y5=y(15000:end,5);

Y6=y(15000:end,6);

Y7=y(15000:end,7);

Y8=y(15000:end,8);

Y9=y(15000:end,9);

Y10=y(15000:end,10);

%绘制主动轮垂直方向的位移-速度庞加莱图；

figure(1)

plot(y(15000:end,2),y(15000:end,3),'.');

function d=fun(t,y)

m1=220;%轴承转子的质量；

m2=20;%主动轮的质量；

m3=440;%从动轮的质量；

R1=0