声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1905|回复: 11

[经典算法] 请教一个幂律非线性方程和代数方程混合微分方程组初值问题的解法

[复制链接]
发表于 2008-1-18 15:52 | 显示全部楼层 |阅读模式

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

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

x
非线性微分方程组中包括幂律方程和代数方程,很大,超过600个变量
幂律方程中的f 是形如a*(y1^3.2)*(y2^(-0.1))之类的线性和。
如附图所示,
我查到有用mixed Adams-BDF 方法作的,
但是为了保证数值稳定性,只采用二阶方法,这样会不会损失精度?
我想问有没有其他更为合适的方法呢,谢谢!

方程组

方程组
回复
分享到:

使用道具 举报

 楼主| 发表于 2008-1-18 15:54 | 显示全部楼层

错了

错了,方程中因变量是x,我在文中说成y了,对不起!
发表于 2008-1-29 15:46 | 显示全部楼层
看看文章《一种求解微分代数混合方程组的统一算法》希望对你有所启发
 楼主| 发表于 2008-1-29 19:16 | 显示全部楼层

mixed Adams-BDF methods

mixed Adams-BDF methods
 楼主| 发表于 2008-1-29 19:18 | 显示全部楼层
PAGE 2
page2.jpg
 楼主| 发表于 2008-1-29 19:23 | 显示全部楼层
怎么只有图。。。ft!

看过大虾提到的文章
另外以前看过一个解法
mixed Adams-BDF method
但是对于具体怎么解还是不大懂
比如怎么选择Adams 和 BDF
是在一个步骤中同时使用这两种方法计算吗?
是不是分别计算不同的变量
比如微分方程的变量就用Adams方法计算
而代数方程中的变量就用BDF算法计算
还望大虾赐教
另外也想用那篇中文的方法试试
谢谢!
 楼主| 发表于 2008-1-29 19:24 | 显示全部楼层
page3.jpg
 楼主| 发表于 2008-1-29 19:24 | 显示全部楼层
PAGE 4
page4.jpg
 楼主| 发表于 2008-1-29 19:26 | 显示全部楼层
STILL
page5.jpg
 楼主| 发表于 2008-1-29 19:27 | 显示全部楼层
PAGE 6
page6.jpg
 楼主| 发表于 2008-1-29 19:29 | 显示全部楼层
page7.jpg
发表于 2008-2-19 16:15 | 显示全部楼层
Adams 适用与求解非刚性方程
BDF    用于求解刚性方程
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-24 00:55 , Processed in 0.099219 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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