24小时热门版块排行榜    

查看: 777  |  回复: 1

gxjhoney

银虫 (小有名气)

[求助] 急需要帮助,如何计算我的状态方程中的参数 已有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
请问各位大牛,如何计算式子中的四个参数?需要详细过程!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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


急需要帮助,如何计算我的状态方程中的参数-1
c2.jpg

2楼2016-01-06 15:01:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gxjhoney 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +23 医学老男孩 2026-08-13 49/2450 2026-08-14 20:02 by zhaifei
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +6 启萌科技 2026-08-12 21/1050 2026-08-14 19:48 by 启萌科技
[基金申请] filecode +6 cratir 2026-08-14 10/500 2026-08-14 18:29 by 笑叹辞穷
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
请填处理意见