24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

CYYT

金虫 (小有名气)

应助: 5 (幼儿园)
金币: 4098.2
散金: 10
红花: 3
帖子: 161
在线: 195.7小时
虫号: 2244949
注册: 2013-01-16
专业: 能源化工

[求助] 求高版本 1stopt跑一下已有1人参与

Title "fun";
Variable f1,f2,f3,f4,T,k,r1;
Function r1=ln(k1)*exp(E/-8.314/T)*(f1*f3*f3-f4/k)/[1+f1*exp(a1-b1*(1/T-1/508.9))+f2*exp(a2-b2*(1/T-1/508.9))+f3*exp(a3-b3*(1/T-1/508.9))]^3;
data;
0.055401453 0.06568836 0.310940986 0.03388052 505.85 0.375502591 3.637487987
0.062438972 0.073957325 0.350203427 0.03983632 504.45 0.400287073 3.875227293
0.06951785 0.082275676 0.389613788 0.046072104 504.85 0.393030011 3.838005558
0.075342296 0.089100948 0.422011831 0.050927616 504.35 0.402123899 3.969231143
0.062428462 0.074049173 0.349970255 0.0432256 516.85 0.229853036 3.827809524
0.062421901 0.074093139 0.349849554 0.04253792 523.55 0.172092587 3.774803171
0.06116489 0.072336052 0.343314957 0.027342157 490.85 0.75820361 2.327431006
0.061168402 0.072258127 0.343479618 0.021084269 482.65 1.132812726 1.789155527
0.028574399 0.034715226 0.350085136 0.022047021 515.55 0.243329481 1.770600166
0.027986893 0.033932423 0.343158611 0.022257254 493.65 0.662970349 1.680596038
0.067324173 0.03064251 0.285081748 0.03654528 528.45 0.139877688 3.992327182
0.067340041 0.03060705 0.285279281 0.042628301 514.05 0.259951812 4.380735361

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有100人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

1楼 2015-11-01 12:50:07

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

CYYT

金虫 (小有名气)

应助: 5 (幼儿园)
金币: 4098.2
散金: 10
红花: 3
帖子: 161
在线: 195.7小时
虫号: 2244949
注册: 2013-01-16
专业: 能源化工

引用回帖:

3楼: Originally posted by lipenggg at 2015-11-02 15:13:29
是的，公式里只能有小括号

不知道为什么，我拟合的结果里出现了“nan”，求解答啊。。

捕获.JPG

赞一下

回复此楼

5楼2015-11-02 18:55:37

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

查看全部 11 个回答

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
CYYT: 金币+10, ★★★很有帮助, 谢谢你了。。 2015-11-02 18:51:59

把中括号“[]”改成小括号“()”低版本的也能运行吧。
公式仔细检查下，似乎有过拟合现象，参数值不唯一：
1：
均方差(RMSE):0.21698029314672
残差平方和(SSE):0.564965371368438
相关系数(R): 0.975158853698942
相关系数之平方(R^2): 0.950934789947434
确定系数(DC): 0.950779667454345
卡方系数(Chi-Square): 0.130254031694643
F统计(F-Statistic): 11.6096040370822

参数最佳估算
---------- -------------
k1 0.999999988673665
e -95862.8129665484
a1 -867833.168533956
b1 43.9156147666947
a2 -4982.47558760389
b2 0
a3 -0.00833194402515029
b3 1656.35069912403

2：
均方差(RMSE):0.216980293146719
残差平方和(SSE):0.564965371368435
相关系数(R): 0.975158853730298
相关系数之平方(R^2): 0.950934790008589
确定系数(DC): 0.950779667454345
卡方系数(Chi-Square): 0.130254031979546
F统计(F-Statistic): 11.6096040370823

参数最佳估算
---------- -------------
k1 0.999999988673666
e -95862.8132743037
a1 -17440154.537724
b1 4235845.00626166
a2 -1395.58639416365
b2 -1613072.39980911
a3 -0.00833195730676999
b3 1656.35070431295

赞一下

回复此楼

2楼2015-11-01 19:01:37

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