版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (6443)
>考研 (1698)
>导师招生 (1354)
>文献求助 (495)
>虫友互识 (290)
>考博 (279)
>硕博家园 (155)
>博后之家 (116)
>基金申请 (109)
>招聘信息布告栏 (92)
>论文投稿 (90)
>休闲灌水 (86)
>第一性原理 (65)
>教师之家 (64)
>公派出国 (43)
>外文书籍求助 (32)

返回列表

gxjhoney

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 474.7
散金: 17
帖子: 107
在线: 29.5小时
虫号: 2906953
注册: 2014-01-01
专业: 化工热力学和基础数据

[求助] 急需要帮助，如何计算我的状态方程中的参数已有1人参与

8.314*T/(V-0.6468*m*(1+λ*(1-(T/m)^0.5))^2)/ n-(31.6053*m^2*(1+k*(1-( T/m)^0.5))^2/n)/(V^2+1.2963*V*m*(1+λ*(1-(T/m)^0.5))^2/n-0.4184*m^2*(1+λ*(1-(T/m)^0.5))^4/(n^2))-P=0
n, λ, m,k
体积V/cm3/mol 温度T/k 压力p/Mpa
298.1613713 318.15 2
297.5363097 318.15 5
296.3889003 318.15 10
295.174061 318.15 15
294.1480253 318.15 20
302.6598305 333.15 2
301.7178487 333.15 5
300.2872339 333.15 10
298.8161612 333.15 15
297.3855464 333.15 20
306.7065098 348.15 2
305.4476077 348.15 5
303.8821332 348.15 10
302.5998641 348.15 15
301.290623 348.15 20
310.2542082 363.15 2
309.3459413 363.15 5
307.7332659 363.15 10
306.6532863 363.15 15
305.472162 363.15 20
314.1255699 378.15 2
312.9880414 378.15 5
311.5776556 378.15 10
310.531391 378.15 15
308.9456875 378.15 20
请问各位大牛，如何计算式子中的四个参数？需要详细过程！

回复此楼

» 猜你喜欢

光学工程学硕调剂信息已经有32人回复
欢迎加入课题组已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有167人回复
散金币，求好运，祝面上顺利！已经有79人回复
华东师范大学芯片设计徐珑真实水平咋样已经有1人回复
调剂已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有0人回复
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职已经有0人回复
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职已经有0人回复
【新加坡】纳米电子器件项目组有“联合培养博士生”名额已经有2人回复
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。已经有0人回复

1楼 2016-01-06 14:41:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gxjhoney: 金币+10, ★★★★★最佳答案, 非常感谢！ 2016-01-08 20:22:09

是非线性拟合问题吧，用1stOpt求解吧，很简单高效：

CODE:

Parameter n,r,m,k;

Variable V,T,p;

Function p=8.314*T/(V-0.6468*m*(1+r*(1-(T/m)^0.5))^2)/n-(31.6053*m^2*(1+k*(1-(T/m)^0.5))^2/n)/(V^2+1.2963*V*m*(1+r*(1-(T/m)^0.5))^2/n-0.4184*m^2*(1+r*(1-(T/m)^0.5))^4/(n^2));

Data;

//体积V/cm3/mol        温度T/k        压力p/Mpa

298.1613713        318.15        2

297.5363097        318.15        5

296.3889003        318.15        10

295.174061        318.15        15

294.1480253        318.15        20

....

结果稳定唯一：
均方差(RMSE):0.767919758588003
残差平方和(SSE):14.7425188907464
相关系数(R): 0.993061794541294
相关系数之平方(R^2): 0.986171727777574
确定系数(DC): 0.986170244943015
卡方系数(Chi-Square): 0.966262223400765
F统计(F-Statistic): 506.155023798799

参数最佳估算
---------- -------------
n -0.118509025251561
r -0.0646358182721027
m 13.0818803799013
k 1.28729810110388

c1.jpg

c2.jpg

赞一下