| 查看: 737 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】matlab非线性拟合模型参数
|
|||
|
我是菜鸟，算个模型参数， lny=-0.3915*[m1/2/(1+1.2*m1/2)+1.6667*ln(1+1.2m1/2)]+2*c(1)m+ 0.5*c(2) [1-(1+2*m1/2-2*m)exp(-2*m1/2)] +c(3)*m2 我写的计算程序 function f=Pitzer(c,mdata) f=(-0.3915).*(mdata.^0.5./(1+1.2.*mdata.^0.5)+1.66667.*log(1+1.2.*mdata.^0.5))+2.*c(1).*mdata+0.5.*c(2).*(1-(1+2.*mdata.^0.5-2.*mdata).*exp(-2.*mdata.^0.5))+c(3).*mdata.^2; mdata=[0.0855 0.0879 0.1001 0.1051 0.1273 0.1398 0.1729]; ydata=[-0.2731 -0.2640 -0.2917 -0.2837 0.3202 -0.3369 -0.3411]; c0=[1 1 1]; for i=1:100000 c=lsqcurvefit('Pitzer',c0,mdata,ydata); c0=c; end c 答案应该是-0.0224 0.1749 0.0020 可是我的出来的是 107.1222 -147.8175 -365.6594 请高手指点下，可能公式表达出现了问题，非常感谢！！！ [ Last edited by dongdong3881 on 2009-12-17 at 14:55 ] |
回复此楼» 猜你喜欢
🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队
已经有0人回复
-
第一性原理计算方向2026级博士申请 PRB*1，四级484
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有64人回复
-
求助VISSIM破解版软件
已经有0人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有4人回复
2楼2009-03-14 19:08:51
3楼2009-03-31 19:42:25
4楼2009-04-01 12:55:38
tianlangxingaa
铁杆木虫 (著名写手)
- 计算强帖: 4
- 应助: 52 (初中生)
- 贵宾: 0.02
- 金币: 7464.2
- 红花: 19
- 帖子: 1283
- 在线: 735.1小时
- 虫号: 438857
- 注册: 2007-09-30
- 专业: 理论和计算化学
5楼2009-04-01 23:45:46
wujixian
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2003
- 帖子: 423
- 在线: 250.1小时
- 虫号: 126308
- 注册: 2005-12-07
- 性别: GG
- 专业: 化学反应工程
6楼2009-04-13 10:28:12