| 查看: 1249 | 回复: 2 | ||
zhangholmes至尊木虫 (著名写手)
逍遥门主
|
[求助]
matlab非线性最优化求解，微分方程组的拟合，参数估计
|
|
如题，用微分方程组拟合4组不同温度下的数据，求理论值a和b，并使a和b的值分别与实验值y1和y2的差方和最小。并求此时的各个参数值(k1,k2,k3,k4,k5,k6)，求各位大侠略施援手。 0<k1<30, 0<k2<10, 0<k3<5, 0<k4<20, 0<k5<20, 0<k6<5 t=[0.494660000000000,0.502760000000000,0.618930000000000,0.621470000000000,0.623470000000000,0.819730000000000,0.825080000000000,0.828350000000000]; T=463 K时， y1= [0.00800000000000000,0.00990000000000000,0.00782000000000000,0.0108400000000000,0.00937000000000000,0.00882000000000000,0.00614000000000000,0.0108200000000000]; y2=[0.0625100000000000,0.0727400000000000,0.0540300000000000,0.0667600000000000,0.0620200000000000,0.0494400000000000,0.0452000000000000,0.0509400000000000]; %y1和y2分别为a和b的实验值 T=473 K时， y1=[0.00829000000000000,0.0107700000000000,0.00938000000000000,0.0172000000000000,0.0128100000000000,0.0139200000000000,0.0108200000000000,0.0179000000000000]; y2=[0.0622600000000000,0.0719200000000000,0.0525700000000000,0.0604900000000000,0.0586200000000000,0.0444300000000000,0.0405900000000000,0.0439200000000000]; T=483 K时， y1=[0.0169300000000000,0.0229900000000000,0.0188900000000000,0.0313900000000000,0.0265800000000000,0.0258100000000000,0.0194200000000000,0.0293300000000000]; y2=[0.0537100000000000,0.0598300000000000,0.0431500000000000,0.0465000000000000,0.0450300000000000,0.0326700000000000,0.0321200000000000,0.0327000000000000]; T=493 K时， y1=[0.0328500000000000,0.0418500000000000,0.0330300000000000,0.0512500000000000,0.0424800000000000,0.0389500000000000,0.0300900000000000,0.0443800000000000]; y2=[0.0381700000000000,0.0415400000000000,0.0292800000000000,0.0273400000000000,0.0295700000000000,0.0197800000000000,0.0216900000000000,0.0180800000000000]; pA = 101*(1 - (2 b + a)); pW = 101*(a+ b); da/dt=k(4).*k(1).*k(2).*k(3).*pA.*pW./(k(1).*k(2).*k(4).*k(5).*pA^2+k(1).*k(2).*k(3).*pA.*pW+k(4)+k(1).*k(2).*k(4).*pA.*pW+k(1).*k(4).*pW; db/dt=k(6).*k(1).*k(2).*k(4).*k(5).*pA^2./(k(1).*k(2).*k(4).*k(5).*pA^2+k(1).*k(2).*k(3).*pA.*pW+k(4)+k(1).*k(2).*k(4).*pA.*pW+k(1).*k(4).*pW; [ 来自科研家族 木虫催化 ] |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有300人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- 二阶非线性微分方程求解
已经有14人回复
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 请教用matlab求解一个非线性偏微分方程组的数值解
已经有8人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab非线性参数拟合问题
已经有7人回复
- 求助，origin7.5 中的非线性拟合参数为什么一直在变，而且误差很大！！！
已经有4人回复
- 求助：origin8.0 自定义非线性拟合，参数拟合值与设定初始值一样
已经有18人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】非线性参数的拟合方法或者拟合软件？（分母、指数和对数）
已经有3人回复
- 【分享】PDECOL-经典数值求解非线性偏微分方程组源程序
已经有142人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-05-12 14:34:13
zhangholmes
至尊木虫 (著名写手)
逍遥门主
- 应助: 17 (小学生)
- 金币: 17040
- 散金: 50
- 红花: 12
- 帖子: 1242
- 在线: 226.3小时
- 虫号: 967701
- 注册: 2010-03-11
- 性别: GG
- 专业: 能源化工
|
谢谢回复，公式最后忘记加了一个括号。不要求统一参数，每组数据都有一组独立的参数，共四组参数。 公式如下： da/dt=k(4).*k(1).*k(2).*k(3).*pA.*pW./(k(1).*k(2).*k(4).*k(5).*pA^2+k(1).*k(2).*k(3).*pA.*pW+k(4)+k(1).*k(2).*k(4).*pA.*pW+k(1).*k(4).*pW）; db/dt=k(6).*k(1).*k(2).*k(4).*k(5).*pA^2./(k(1).*k(2).*k(4).*k(5).*pA^2+k(1).*k(2).*k(3).*pA.*pW+k(4)+k(1).*k(2).*k(4).*pA.*pW+k(1).*k(4).*pW）; |
3楼2013-05-12 15:59:36
