应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 请问这个算是反应动力学方程吗？怎么求得k值？
81/1返回列表
查看: 3213  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

飞鸿印雪jay

银虫 (小有名气)

[求助] 请问这个算是反应动力学方程吗？怎么求得k值？

自己建的数学模型，反应速率方程算不算反应动力学方程？怎么用matlab求解参数k呢？假设浓度C可以试验测定。
请问这个算是反应动力学方程吗？怎么求得k值？
OE%JHWBEA]D{_`_H1%%8ZI8.jpg
回复此楼

» 收录本帖的淘帖专辑推荐

动力学拟合

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-07 15:01:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
微分方程拟合问题，用1stOpt很好解决，搜一下有很多例子。
赞一下(1人) 回复此楼
2楼2013-11-07 18:39:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞鸿印雪jay

银虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2013-11-07 18:39:26
微分方程拟合问题，用1stOpt很好解决，搜一下有很多例子。

请问怎么具体求解啊？现在参数k比方程的数目多，不知道怎么求解啊？想用matlab
赞一下 回复此楼
3楼2013-11-09 21:16:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
可采用逐步线性回归方法。
赞一下(1人) 回复此楼
4楼2013-11-09 21:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
飞鸿印雪jay: 金币+5, ★★★很有帮助 2014-01-08 23:58:45
以下是MATLAB拟合含2个方程的常微分方程组的实例：
http://muchong.com/bbs/viewthread.php?tid=6425538&authorid=1122189

你的是包含5个方程的，程序可类比。
赞一下(1人) 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2013-11-10 09:46:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞鸿印雪jay

银虫 (小有名气)
引用回帖:
5楼: Originally posted by 月只蓝 at 2013-11-10 09:46:19
以下是MATLAB拟合含2个方程的常微分方程组的实例：
http://muchong.com/bbs/viewthread.php?tid=6425538&authorid=1122189

你的是包含5个方程的，程序可类比。

function k1k2k3
format long
clear all
clc
tspan = [0  30  50  80 140 170 200 230 260 290 320 360 400 460 520 580 600];
x0 = [9;0;0;0;0];
k0 = [1  0  0  0  0  0  0  0  0  0];   
lb = [0  0  0  0  0  0  0  0  0  0];
ub = [1  1  1  1  1  1  1  1  1  1];

data=[30        7.939        1.458     7.939        1.458     7.939      
      50        7.687        1.535     7.687        1.535     7.687
      80        7.289        1.602     7.289        1.602     7.289
      140        6.658        1.717    6.658        1.717     6.658
      170        6.531        1.722    6.531        1.722     6.531
      200        6.218        1.671    6.218        1.671     6.218
      230        5.979        1.620    5.979        1.620     5.979
      260        5.591        1.550    5.591        1.550     5.591
      290        5.414        1.488    5.414        1.488     5.414
      320        4.968        1.433    4.968        1.433     4.968
      360        4.692        1.350    4.692        1.350     4.692  
      400        4.438        1.319    4.438        1.319     4.438
      460        4.144        1.294    4.144        1.294     4.144
      520        4.041        1.294    4.041        1.294     4.041
      580        4.052        1.287    4.052        1.287     4.052     
      600        4.052        1.287    4.052        1.287     4.052      
      ];
yexp = data(:,2:6);

[k,resnorm,residual,exitflag,output,lambda,jacobian]=...
lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp);      
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n')
fprintf('\tk1 = %.9f ± %.9f\n',k(1),ci(1,2)-k(1))
fprintf('\tk2 = %.9f ± %.9f\n',k(2),ci(2,2)-k(2))
fprintf('\tk3 = %.9f ± %.9f\n',k(3),ci(3,2)-k(3))
fprintf('\tk3 = %.9f ± %.9f\n',k(4),ci(3,2)-k(4))
fprintf('\tk3 = %.9f ± %.9f\n',k(5),ci(3,2)-k(5))
fprintf('\tk3 = %.9f ± %.9f\n',k(6),ci(3,2)-k(6))
fprintf('\tk3 = %.9f ± %.9f\n',k(7),ci(3,2)-k(7))
fprintf('\tk3 = %.9f ± %.9f\n',k(8),ci(3,2)-k(8))
fprintf('\tk3 = %.9f ± %.9f\n',k(9),ci(3,2)-k(9))
fprintf('\tk3 = %.9f ± %.9f\n',k(10),ci(3,2)-k(10))
fprintf('  The sum of the squares is: %.9e\n\n',resnorm)

function f = ObjFunc(k,tspan,x0,yexp)           % 目标函数
[t, Xsim] = ode45(@KineticsEqs,tspan,x0,[],k);
Xsim1=Xsim(:,1);
Xsim2=Xsim(:,2);
ysim(:,1) = Xsim1(2:end);
ysim(:,2) = Xsim2(2:end);

size(ysim(:,1));
size(ysim(:,2));
size(yexp(:,1));
size(yexp(:,2));

f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2))];

function dCdt = KineticsEqs(t,C,k)              % ODE模型方程
dCAdt = (k(1)+k(5)+k(6)+k(7))*C(1);
dCBdt = k(1)*C(1)-(k(2)+k(8)+k(9))*C(2);
dCCdt = k(5)*C(1)+k(2)*C(2)-(k(3)+k(10))*C(3);
dCDdt = k(6)*C(1)+k(9)*C(2)+k(3)*C(3)-k(4)*C(4);
dCEdt = k(7)*C(1)+k(8)*C(2)+k(10)*C(3)+k(4)*C(4);
dCdt = [dCAdt; dCBdt;dCCdt;dCDdt;dCEdt];
大神帮忙改一下啊
赞一下 回复此楼
6楼2013-11-10 21:30:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

278647873

新虫 (小有名气)
十个参数这么多，你得做多少次试验啊，你是没找到关键的地方吧，即使拟合出来也不会准确吧
赞一下 回复此楼
7楼2013-11-12 11:01:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zj34331197

木虫 (正式写手)

木虫
引用回帖:
3楼: Originally posted by 飞鸿印雪jay at 2013-11-09 21:16:29
请问怎么具体求解啊？现在参数k比方程的数目多，不知道怎么求解啊？想用matlab...

不知道问题解决没有，如果以反应进度表示动力学方程，可能要简单明了的多。
赞一下(1人) 回复此楼
路漫漫其修远兮
8楼2015-06-28 20:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 飞鸿印雪jay 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +17 木水思豆 2026-08-25 20/1000 2026-08-26 01:26 by chengyan1220
[基金申请] 今天务委会开完了，明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 今天基金会出结果吗？20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
关闭
请填处理意见
关闭