MATLAB-ODE求解器：无法满足积分容差

Question

我在MATLAB中有一个ode求解器的问题。我使用了所有的ode求解器，如ode23s、ode23、ode15s、ode45等等。我的代码无法计算，因为错误警告：

警告：在t=8.190397E+01时失败。在时间t时，步长不减小到允许的最小值(2.273737E-13)以下，就无法满足积分公差。

我想计算它，请帮助我直接在我的代码。谢谢你。

第一个脚本：

% floq.m
global c_alpha c_beta c_gama om ms ks bs mii

 % Parameters            
   c_alpha=1;
   c_beta=1.1;
   c_gama=1.1;  
   ms=1;                             
   ks=1;                                                       
   D=0.01;
   OM=sqrt(ks/ms);                     
   bs=2*D*ms*OM;                       

 % Solver    
   sourad = 0.1:0.13:10;
   for pom = 1:length(sourad)
       eta= sqrt(1/sourad(pom));
       om=eta*OM;                          
       T=2*pi/(eta*OM);    

       for mii=-10*(eta^2):0.13:10*(eta^2)  
           tspan=0:0.01:T;
           [tt,x1]=ode23(@fun,tspan,[1; 0]);    
       end
   end       

% fun.m

 function v=fun(tt,x1);
 global c_alpha c_beta c_gama om ms ks bs mii

 mt=ms*cos(om*tt);
 bt=bs*cos(2*om*tt);
 kt=ks*cos(2*om*tt);

 % Matrix A
 A=[(-bs+mii*c_beta*bt)/(ms-mii*c_gama*mt) (-ks+mii*c_alpha*kt)/(ms-mii*c_gama*mt); 1 0];

 % Method
 v=A*x1;

谢谢你。

Answer 1

您正在对ODE函数的一个极点进行积分。到了极点，每一个解都结束了。极点左右的轨迹不能合并成一个更大的轨迹。

极点是1=mii*c_gama*cos(om*t)的第一个正解t。如果mii*c_gama>=1总有这样的解决方案。