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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 请问这个算是反应动力学方程吗？怎么求得k值？
81/1返回列表
查看: 3221  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金不中，共勉 +5 eulota 2026-08-26 5/250 2026-08-26 23:33 by young90
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +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
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭