| 查看: 2442 | 回复: 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 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有219人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有30人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有56人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助用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 (初中生)
- 金币: 2109.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 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
