隐式微分方程一阶广义-α求解

yanyongju · 发表于 2012-1-4 13:59

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

本帖最后由 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方程。

Generalized-alpha method.rar (2.27 KB, 下载次数: 1)
Lambda的取值范围为[-1,0]，表示算法阻尼。

账号		自动登录	找回密码
密码			我要加入

[共享资源] 隐式微分方程一阶广义-α求解

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。