小弟最近刚接触EMD,想用Rilling的程序试验一下Huang的98年的文章中所给的duffing方程例子。求解此微分方程时遇到了问题。运行解方程程序得到
Unable to meet integration tolerances without reducing the step size below the smallest value allowed (3.552714e-015) at time t.
In ode23 at 369
In xx at 1
因此我怀疑得到的数值解是发散的,不是周期性的,但这与文章结果明显不符,请高手指点我错在何处?另外用Rilling程序运行此例需要更改参数吗? 希望哪位做过类似试验的大侠指点我该怎麽做,不胜感激!(方程及程序如图)