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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab求解动力学 方程中的参数
141/2返回列表12下一页
查看: 2333  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yupeifeng

金虫 (小有名气)

[交流] matlab求解动力学 方程中的参数 已有3人参与

实验数据位（t,c）=(0,6.94),(2,6.4),(4,6.3),(8,6.2),(48,6.1),动力学方程为：dc/dt=-k*(c0-c)^(2/3)*(c-c~)^p，其中，c0指初始浓度（6.9），c~指平衡浓度（6.1）。
如何求方程中的k与p.文献中给出的方法是先积分求出一个含有k、p的方程，然后拟合求出。如何写程序呢？？谢谢大家帮忙了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
憧憬未来
1楼 2012-06-24 09:53:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wptom

铁虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
这个貌似数据处理吧，不是编程的问题。
1.令dc/dt=a; 动力学方程两边取log，变形为线性方程。
2.将实验数据整理成dc/dt；c0-c;c-c~;
3.在excel 列出上述3列，做图，线性趋势线，就可以了。
赞一下 回复此楼
2楼2012-06-24 14:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
2楼: Originally posted by wptom at 2012-06-24 14:00:09
这个貌似数据处理吧，不是编程的问题。
1.令dc/dt=a; 动力学方程两边取log，变形为线性方程。
2.将实验数据整理成dc/dt；c0-c;c-c~;
3.在excel 列出上述3列，做图，线性趋势线，就可以了。

你的思路是对的，关键是微分方程比较复杂，应先积分后去对数吧？
赞一下 回复此楼
憧憬未来
3楼2012-06-24 15:18:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihaiyong826

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
求积分用maple来做，不用编程，很简单。
赞一下 回复此楼
4楼2012-06-24 21:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wptom

铁虫 (正式写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
贺仪: 金币+2, 谢谢交流 2012-06-26 17:32:09
引用回帖:
3楼: Originally posted by yupeifeng at 2012-06-24 15:18:45
你的思路是对的，关键是微分方程比较复杂，应先积分后去对数吧？...

不需要积分。
你的目的是求 K和 P。
两边取log后， log(a)=log(K)+ (2/3)long(c0-c)+P log(c-c~)
令 y=log(a)-(2/3) log(C0-C); x=log(C-C~);
就是一般线性方程 y = m* x + b
其中 m=P; b=log(K)
如果数据点是近似在一条直线上不就求出m、b了？
赞一下 回复此楼
5楼2012-06-25 09:57:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wptom

铁虫 (正式写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
dbb627: 金币+1, 欢迎交流 2012-06-30 17:46:54
注意数据的变换；a=[C(t)-C(t0)]/(t-t0) 或者= c(t+1)-c（t）

[ Last edited by wptom on 2012-6-25 at 10:02 ]
赞一下 回复此楼
6楼2012-06-25 10:01:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
6楼: Originally posted by wptom at 2012-06-25 10:01:25
注意数据的变换；a=/(t-t0) 或者= c(t+1)-c（t）

你好，那个a=[C(t)-C(t0)]/(t-t0) 或者= c(t+1)-c（t）,这个是因为啥？原方程两边取对数后方程左侧是dc/dt,这个怎么取对数？
赞一下 回复此楼
憧憬未来
7楼2012-06-25 19:36:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
6楼: Originally posted by wptom at 2012-06-25 10:01:25
注意数据的变换；a=/(t-t0) 或者= c(t+1)-c（t）

我明白你的意思了。但是dc/dt可以用两个相邻点的纵坐标差减去横坐标差来代替？这样岂不是误差很大啊？有没有积分可以算出来的，
赞一下 回复此楼
憧憬未来
8楼2012-06-25 19:40:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wptom

铁虫 (正式写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
dbb627: 金币+1, 欢迎交流 2012-06-30 17:47:04

dc/dt的物理意义是浓度随时间的变化速率。
你老惦记着积分干嘛？？
问一问你老师吧。特别是教分析化学的老师。
赞一下 回复此楼
9楼2012-06-26 10:50:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
9楼: Originally posted by wptom at 2012-06-26 10:50:29

dc/dt的物理意义是浓度随时间的变化速率。
你老惦记着积分干嘛？？
问一问你老师吧。特别是教分析化学的老师。

这样得出的结果和别人积分所得的差别很大，拟合出来的相关系数相差很大。呵呵，不过也是一种方法，呵呵，谢谢了。
赞一下 回复此楼
憧憬未来
10楼2012-06-26 11:21:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yupeifeng 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:59 by zNcEhTcH7ihq
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:50 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:23 by CgyNCDVNhGVg
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:15 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 14:14 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭