24小时热门版块排行榜

返回列表

xxz903

金虫 (小有名气)

应助: 33 (小学生)
金币: 1045.6
散金: 190
红花: 4
帖子: 253
在线: 66.7小时
虫号: 1360336
注册: 2011-08-03
性别: GG
专业: 化工系统工程

[求助] 1stOpt求助！无法计算~说是fatal:line too long(more than 1023 characters)

求高手解答啊，代码如下：

Title "单电解质参数回归";
Parameters r1,q1,r2,q2,t12,t21;
Variable m,y1,y2;
Constant z=10;
conststr x1=m/(m+h);
conststr x2=1-x1;
conststr o1=x1*r1/(x1*r1+x2*r2);
conststr o2=x2*r2/(x1*r1+x2*r2);
conststr b1=x1*q1/(x1*q1+x2*q2);
conststr b2=x2*q2/(x1*q1+x2*q2);
conststr l1=0.5*z*(r1-q1)-(r1-1);
conststr l2=0.5*z*(r2-q2)-(r2-1);
SharedModel;
Function y1=ln(o1/x1)+0.5*z*q1*ln(o1/b1)+o2*(l1-r1*l2/r2)-q1*ln(b1+b2*t21)+b2*q1*(t21/(b1+b2*t21)-t12/(b2+b1*t12));
Function y2=ln(o2/x2)+0.5*z*q2*ln(o2/b2)+o1*(l2-r2*l1/r1)-q2*ln(b2+b1*t12)+b1*q2*(t12/(b2+b1*t12)-t21/(b1+b2*t21));
data;
0 0 0
0.1 -0.673344553 -0.003214304
0.2 -0.783071888 -0.005447496
0.3 -0.841647189 -0.008713882
0.4 -0.869884359 -0.010990339
0.5 -0.884307686 -0.014309599
0.6 -0.889162064 -0.018687569
0.7 -0.889162064 -0.02207116
0.8 -0.88673193 -0.025484701
0.9 -0.879476759 -0.029975034
1 -0.869884359 -0.038731848
1.2 -0.850971266 -0.042636996
1.4 -0.830113036 -0.051995295
1.6 -0.807436327 -0.061533872
1.8 -0.783071888 -0.071257081
2 -0.759286983 -0.081169452
2.5 -0.701179352 -0.107972916
3 -0.650087691 -0.136153773
3.5 -0.599656837 -0.165810242
4 -0.553385238 -0.197052419
4.5 -0.512493681 -0.227286842
5 -0.47320876 -0.260553522

回复此楼

» 猜你喜欢

新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有145人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有1人回复
投稿chemical physical letters不送审？已经有2人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】关于三甲基苯用NBS溴代的反应已经有15人回复
【求助】MS的界面模型建立问题，旋转已经有14人回复
求助N2吸附脱附实验问题已经有15人回复
求助哌嗪单boc保护反应已经有6人回复
求助关于solvate的一些疑惑已经有10人回复
紧急求助:SDS-PAGE上样一定要调蛋白含量吗？已经有9人回复
求助，换了高保真酶扩不出来目的条带已经有26人回复
50金币求助电荷计算的高手，这句话如何理解？已经有3人回复
求助磷酸铵，磷酸二氢铵，磷酸一氢铵分解能量顺序已经有17人回复
中试求助加热设备已经有8人回复
紧急求助！！！Endnote X4导入PDF无法识别了已经有27人回复
Gaussian 09 计算出来的chk文件无法转化成fchk文件，求高手指教！已经有9人回复
计算投稿求助，大侠赶紧进已经有11人回复
求助第二相体积分数的计算已经有3人回复
求助，我安装auto CAD2007时，提示：错误 1402。无法打开主键：怎么解决呀！急~~~~ 已经有8人回复

1楼 2012-07-10 20:38:59

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

dingd

铁杆木虫 (职业作家)

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

★
dbb627: 金币+1, 欢迎交流 2012-07-12 09:56:20

5.0版运行没任何问题啊：

均方差(RMSE): 0.0124262605314638
残差平方和(SSE): 0.0067941258350159
相关系数(R): 0.999734045868563
相关系数之平方(R^2): 0.999468162468726
决定系数(DC): 0.999468099291204
F统计(F-Statistic): 1113.25913600748

参数最佳估算
-------------------- -------------
h 37.8325607068761
r1 -1905.12207854584
q1 -0.256488526823719
r2 -279.726412704903
q2 -0.00834967640705846
t12 4.47682374015692
t21 -0.045782658788392

赞一下

回复此楼

2楼2012-07-11 09:18:22

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

xxz903

金虫 (小有名气)

应助: 33 (小学生)
金币: 1045.6
散金: 190
红花: 4
帖子: 253
在线: 66.7小时
虫号: 1360336
注册: 2011-08-03
性别: GG
专业: 化工系统工程

引用回帖:

2楼: Originally posted by dingd at 2012-07-11 09:18:22
5.0版运行没任何问题啊：

均方差(RMSE): 0.0124262605314638
残差平方和(SSE): 0.0067941258350159
相关系数(R): 0.999734045868563
相关系数之平方(R^2): 0.999468162468726
决定系数(DC): 0.9994680992912 ...

谢谢了啊！上一个的值有点问题不应该是负的，我修改了一个目标函数能再帮我拟合一个吗？
Title "单电解质参数回归";
Parameters r1,q1,r2,q2,t12,t21;
Variable m,y1,y2;
Constant z=10;
conststr x1=m/(m+h);
conststr x2=1-x1;
conststr o1=x1*r1/(x1*r1+x2*r2);
conststr o2=x2*r2/(x1*r1+x2*r2);
conststr b1=x1*q1/(x1*q1+x2*q2);
conststr b2=x2*q2/(x1*q1+x2*q2);
conststr l1=0.5*z*(r1-q1)-(r1-1);
conststr l2=0.5*z*(r2-q2)-(r2-1);
SharedModel;
Function y1=ln(o1/x1)+0.5*z*q1*ln(o1/b1)+o2*(l1-r1*l2/r2)-q1*ln(b1+b2*t21)+b2*q1*(t21/(b1+b2*t21)-t12/(b2+b1*t12))-ln(r1/r2)-05*z*q1*ln(q1*r2/(q2*r1))-(l1-r1*l2/r2)+q1*ln(t21)-q1*(1-t12);
Function y2=ln(o2/x2)+0.5*z*q2*ln(o2/b2)+o1*(l2-r2*l1/r1)-q2*ln(b2+b1*t12)+b1*q2*(t12/(b2+b1*t12)-t21/(b1+b2*t21))-ln(r2/r1)-05*z*q2*ln(q2*r1/(q1*r2))-(l2-r2*l1/r1)+q2*ln(t12)-q2*(1-t21);
data;
0 0 0
0.1 -0.673344553 -0.003214304
0.2 -0.783071888 -0.005447496
0.3 -0.841647189 -0.008713882
0.4 -0.869884359 -0.010990339
0.5 -0.884307686 -0.014309599
0.6 -0.889162064 -0.018687569
0.7 -0.889162064 -0.02207116
0.8 -0.88673193 -0.025484701
0.9 -0.879476759 -0.029975034
1 -0.869884359 -0.038731848
1.2 -0.850971266 -0.042636996
1.4 -0.830113036 -0.051995295
1.6 -0.807436327 -0.061533872
1.8 -0.783071888 -0.071257081
2 -0.759286983 -0.081169452
2.5 -0.701179352 -0.107972916
3 -0.650087691 -0.136153773
3.5 -0.599656837 -0.165810242
4 -0.553385238 -0.197052419
4.5 -0.512493681 -0.227286842
5 -0.47320876 -0.260553522

赞一下(1人)

回复此楼

3楼2012-07-11 09:31:53

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

相关版块跳转我要订阅楼主 xxz903 的主题更新

返回列表