隐式微分方程一阶广义-α求解
本帖最后由 yanyongju 于 2012-1-4 14:06 编辑在进行刚性微分方程求解时，MATLAB函数库中的ode15i 经常会遇到这样的问题
Warning: Failure at t=？？？， Unable to meet integration tolerances without reducing the step size below
the smallest value allowed
广义-α算法是一种适用于刚性微分方程求解数值方法，但 一般多用于二阶显式微分方程，
这里编写了适用于一阶隐式微分方程的M文件，供大家分享，算例去Van der Pol方程。
Lambda的取值范围为[-1,0]，表示算法阻尼。
页:
[1]