| 查看: 2257 | 回复: 14 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
拜求！matlab拟合动力学参数已有2人参与
|
||||
|
实验数据如下 t（h）[RH] [ROOH] 0 8.32016 0 1 7.37673 0.21115 2 6.74433 0.7379 3 5.92613 1.10562 4 5.68806 1.29586 5 5.01584 1.58327 6 6.09249 1.36404 8 5.34827 1.84065 10 4.9961 2.00986 动力学模型如图 其中[o2]浓度未测量但认为恒定，想把K4[o2]看成一个参数，并估计出其他的k1 k1‘ K4 K5 K8 f 值，matlab一窍不通，看了几天了还是不明所以，所以恳请各位大神给拟合一下，并给出代码，因为还要拟合其他的，所以想大神给出一个后面好跟着做。。。  捕获.JPG |
回复此楼» 猜你喜欢
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有119人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求助用matlab拟合动力学方程
已经有9人回复- 求助用matlab拟合动力学方程
已经有13人回复
- matlab 拟合方法
已经有4人回复
- 动力学拟合求助
已经有4人回复
- matlab 拟合求助
已经有8人回复
- 关于matlab拟合
已经有4人回复
- matlab拟合Sanchez-lacombe方程
已经有3人回复
- 怎样用一组参数同时拟合两个曲线--matlab
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 如何拟合得到动力学模型参数？？
已经有3人回复
- matlab数据拟合问题求助
已经有7人回复
- 动力学参数拟合
已经有26人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab拟合求值
已经有17人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- matlab非线性参数拟合问题
已经有7人回复
- 动力学数据拟合
已经有5人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
- 【求助】使用Matlab拟合反应动力学方程问题
已经有7人回复
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2107.8
- 红花: 23
- 帖子: 282
- 在线: 111.9小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
随便编了一下，拟合结果有点问题，估计是k的初值k0取得有问题或者数据太少或者实验数据有问题。。。 程序如下，其中k0是我随便设的初值，k1=k1,k2=k1',k3=k8,k4=k4*[O],k5=k5,k6=f function piadatfit2 clear all;clc tspan=[0 1 2 3 4 5 6 8 10]; cexp=[8.32016 0;7.37673 0.21115; 6.74433 0.7379;5.92613 1.10562;5.68806 1.29586; 5.01584 1.58327;6.09249 1.36404;5.34827 1.84065;4.9961 2.00986]; k0=[-0.1 1 -2.3 13 0.2 0.4]; c0=[8.32016 0]; LB=[0 0 0 0 0 0]; UB=[+inf +inf +inf +inf +inf +inf]; [k,resnorm,residual]=lsqnonlin(@objpia,k0,LB,UB,[],cexp,tspan,c0) [tplot cplot]=ode45(@piakin,tspan,c0,[],k); plot(tspan,cexp(:,1),'bx',tplot,cplot(:,1),'b-',tspan,cexp(:,2),'ko',tplot,cplot(:,2),'k-') function f=objpia(k,cexp,tspan,c0) [t c]=ode45(@piakin,tspan,c0,[],k); f1=c(:,1)-cexp(:,1); f2=c(:,2)-cexp(:,2); f=[f1;f2]; function dcdt=piakin(t,c,k) dc1dt=-(k(4)*k(5)*c(1)*sqrt(k(1)*c(2)/k(3))/(k(4)+k(5)*c(1))+2*k(6)*k(2)*c(2)); dc2dt=(k(4)*k(5)*c(1)*sqrt(k(1)*c(2)/k(3))/(k(4)+k(5)*c(1))+2*k(6)*k(2)*c(2)); dcdt=[dc1dt dc2dt]'; |
小小莓邮8楼2014-02-19 15:16:14
2楼2014-02-18 16:54:32
3楼2014-02-18 19:00:19
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 计算强帖: 1
- 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2014-02-18 22:54:21