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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 用matlab求动力学参数
71/1返回列表
查看: 2022  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者自制奶酪将赠送您 100 个金币

自制奶酪

铜虫 (小有名气)

[求助] 用matlab求动力学参数 已有3人参与

本人MATLAB小白，最近审稿人要求做个动力学模型，我不会用MATLAB拟合求速率常数，希望木虫的大神可以教我，谢谢。
(dC_A)/dt=-k1*CA-k7*CA
(dC_B)/dt=k1*CA-k2*CB- k5*CB
(dC_C)/dt=k2*CB-k3*CC- k6*CC
(dC_D)/dt=k3*CC-k4*CD
(dC_E)/dt= k4*CD
(dC_F)/dt=k5*CB+k6*CC+ k7*CA
此外还采用Nelder-Mead 方法通过最小化目标函数优化速率常数
F_obj=∑_m▒{∑_(n=1)^(n_p)▒[C_(m,n,calc)-C_(m,n,exp) ] ^2 }

C_(m,n,calc)是给定的时间n下组分的实验浓度；C_(m,n,exp)是给定的时间n下组分的计算浓度
第一组：
T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)
A= (18.7, 6.5, 0.7, 2.1, 1.7, 1.5, 1.2, 0.2)
B= (48.9, 31.3, 15.3, 6.9, 2.2, 1.8, 1.7, 0.5)
C= (19.3, 34.7, 41.2, 40.9, 37.2, 31.7, 26.1, 19.4)
D= (9.7, 19.6, 28.7, 34.5, 38.4, 38.2, 36.3, 33.4)
E= (0.2, 2.6, 5.7, 9.3, 17.1, 24.8, 30.5, 36.7)

第二组：
T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)
A= (23.7, 7.2, 1.7, 1.1, 0.7, 0.5, 0.2, 0.2)
B= (55.4, 44.3, 27.3, 11.9, 1.2, 1.1, 0.7, 0.5)
C= (14.3, 33.3, 43.2, 49.3, 51.5, 44.7,36.1, 30.4)
D= (3.7, 11.3, 19.6, 26.4, 33.4, 35.2, 34.1, 32.4)
E= (0.2, 1.2, 3.7, 7.5, 13.7, 18.3, 23.2, 27.3)

第三组：
T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)
A= (33.7, 15.2, 3.7, 1.1, 0.7, 0.5, 0.2, 0.2)
B= (49.7, 53.4, 44.3, 28.3, 8.7, 1.1, 0.7, 0.5)
C= (9.3, 20.8, 32.2, 43.7, 53.5, 55.7, 51.4, 46.9)
D= (2.7, 7.9, 13.7, 18.4, 25.4, 28.2, 29.4, 30.4)
E= (0, 0.2, 2.1, 4.3, 8.7, 12.3, 15.7, 19.3)
回复此楼

» 猜你喜欢
简单快乐
1楼 2021-01-05 15:31:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
https://blog.sina.com.cn/s/blog_c0cb8ce60102ysqu.html
赞一下 回复此楼
数值计算
2楼2021-01-05 21:04:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
楼主是想3组数据求出一组共同的参数？还是每一组数据分别求出一组参数？自己用1stOpt试试。
赞一下 回复此楼
3楼2021-01-05 21:21:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

自制奶酪

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by dingd at 2021-01-05 21:21:18
楼主是想3组数据求出一组共同的参数？还是每一组数据分别求出一组参数？自己用1stOpt试试。

每一组分别求出一组参数k
赞一下(1人) 回复此楼
简单快乐
4楼2021-01-05 21:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

自制奶酪

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 独孤神宇 at 2021-01-05 21:04:32
https://blog.sina.com.cn/s/blog_c0cb8ce60102ysqu.html

您好，链接打开加密了
赞一下 回复此楼
简单快乐
5楼2021-01-05 21:42:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

F_obj=∑_m▒{∑_(n=1)^(n_p)▒[C_(m,n,calc)-C_(m,n,exp) ] ^2 }这个关系式不明确，请重新书写。
赞一下 回复此楼
QQ：2120156492
6楼2021-01-11 19:23:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梅子炒米粉

新虫 (小有名气)
请问楼主问题解决了吗？方便交流交流吗？

发自小木虫Android客户端
赞一下 回复此楼
7楼2021-04-24 21:52:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 自制奶酪 的主题更新
71/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[考研] 化学工程321分求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-15 18:49 by a不易
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分，求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 298求调剂 +3 Vv呀！ 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
[考研] 337求调剂 +3 睡醒，。 2026-03-09 3/150 2026-03-09 10:02 by 求调剂zz
信息提示
关闭
请填处理意见
关闭