| 查看: 864 | 回复: 2 | ||
xxz903金虫 (小有名气)
|
[求助]
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 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有75人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】关于三甲基苯用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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
★
dbb627: 金币+1, 欢迎交流 2012-07-12 09:56:20
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
xxz903
金虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 1045.6
- 散金: 190
- 红花: 4
- 帖子: 253
- 在线: 66.7小时
- 虫号: 1360336
- 注册: 2011-08-03
- 性别: GG
- 专业: 化工系统工程
|
谢谢了啊！上一个的值有点问题 不应该是负的，我修改了一个目标函数 能再帮我拟合一个吗？ 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 |
3楼2012-07-11 09:31:53
