声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1410|回复: 0

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

[复制链接]
发表于 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],表示算法阻尼。






Generalized-alpha method.rar

2.27 KB, 下载次数: 1

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-26 03:42 , Processed in 0.113701 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表