ode45本身是变步长的呀,怎么设置最小步长?我改用其它ode函数计算都会出现这个警告Warning: Failure at t=5.373174e-001. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.908935e-015) at time t.,用ode4函数算了也出现这个问题,缩小tspan也不行,现贴上我的程序请高手指点。
clear
clc
miu=0.0178;
L=201.6e-3;
R=210e-3;
LD=L/(2*R);
C=0.699e-3;
pofai=C/R;
mg=104860;
m=mg/9.8;
ome=450;
pa=2*miu*ome/(pofai.^2);
Kxx = 6.453356740612708;
Kxy = -4.133413328965316;
Kyx =1.146014052781378; Kyy = 7.846056548809618;
Dxx =1.145752328144106;Dxy =1.200350375824454;
Dyx =1.200350375824454;Dyy = 2.277131377326176;
KXxx = 4.310137875942596;KXxy =1.621686159880267;KXyy=-32.99946362265302;
DXxx =