24小时热门版块排行榜    

CyRhmU.jpeg
查看: 961  |  回复: 3

xys1109

铜虫 (小有名气)

[求助] 龙格-库塔(runger-kuta)求解微分方程已有2人参与

微分方程为 d(A)/dt=k1[A][B]-k2[C][D] 初始条件:t=0 [A]=A0 [B]=B0 [C]=0 [D]=0
约束条件为[A]+[C]=A0  [B]+[C]=B0 [C]=[D]  采用四阶rungerkuta求解微分方程 该方程有两个位置参数k1 k2 使用什么方法或者软件能求得含有未知参的解 求大神 指点

摘自一篇文献 文献原文为: 酸催化酯化法制备生物柴油动力学研究(化学反应工程与工)彭宝祥 舒庆 王光润 王金福
回复此楼
静下心来!!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jlzeng

木虫 (正式写手)

直接找清华的老师问啊,呵呵
沉默坚守
2楼2014-01-02 10:23:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bhyxw

金虫 (正式写手)

【答案】应助回帖

使用龙格-库塔算法,用C语言编程求解
3楼2014-01-24 14:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

这个方程可是得到解析解的。
将约束条件代入微分方程中后,根据各参数A0、B0、k1、k2等数值的取值的不同,可以有以下两种解的形式:
Arctg{([A]+b)/b}/b+Const 或1/(2*b)*Ln{([A]+a-b)/([A]+a+b)}+Const。
a、b均为由A0、B0、k1、k2等已知数表达的函数;Const为积分常数,由初始条件确定。
4楼2014-01-25 23:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xys1109 的主题更新
信息提示
请填处理意见