原微分方程可以求出原函数,不需要做常微分拟合:
dx/dt=a*(1-x)^2-b*(x/2)^2
dx/(a*(1-x)^2-b*(x/2)^2)=dt,两边积分得:
t = -0.2e1 * ((a * b) ^ (-0.1e1 / 0.2e1)) * atanh(((2 * (4 * a - b) * x - 8 * a) * (a * b) ^ (-0.1e1 / 0.2e1)) / 0.4e1)
得到上述代数式之后,直接作非线性拟合即可,1.5版本1stOpt即可实现。
积分结果.png