[求助] 龙格-库塔（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 使用什么方法或者软件能求得含有未知参的解求大神指点

摘自一篇文献文献原文为：酸催化酯化法制备生物柴油动力学研究（化学反应工程与工）彭宝祥舒庆王光润王金福

回复此楼

» 猜你喜欢

天津工业大学郑柳春团队长期招化学化工、高分子化学或有机合成方向的博士生已经有0人回复
天津工业大学郑柳春团队长期招化学化工、高分子化学或有机合成方向的博士生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有103人回复
同济大学张海军教授2026年度博士招生已经有1人回复
同济大学张海军教授2026年度博士招生已经有0人回复
工作一年半了，突然分配到浮选药剂的合成，我想问问浮选药剂是不是夕阳产业了已经有18人回复
我用苯甲酸做醇酸树脂，苯甲酸在80度左右就开始升华了，各位虫友，有什么好办法吗已经有1人回复
定向冷冻聚酰亚胺气凝胶已经有0人回复
HMDI和PTMEG合成的热熔胶为什么固化不了一直是发粘的已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

静下心来！！！！

1楼2013-03-08 09:58:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

TEPI: 2
应助: 20282 (院士)
金币: 145683
红花: 1374
帖子: 93066
在线: 7693.7小时
虫号: 1482829
注册: 2011-11-08
性别: GG
专业: 功能陶瓷

【答案】应助回帖

这个方程可是得到解析解的。
将约束条件代入微分方程中后，根据各参数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为积分常数，由初始条件确定。