[fanxing39]

哈哈,经过你的提醒,我明白了, 调整下" 相对误差容许上限" 就可以了
原程序中"ode45"命令前面先进行求解参数的修改,使得求解保证足够的精度
插入的命令如下:
options=odeset('RelTol',1e-7)%默认是0.001,现在改为10^(-7)
[time1,State1] = ode45(@EquationEx5,t,state0,options,Parameters);