【答案】应助回帖

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 新建文本文档.txt

2012-07-11 13:50:40, 1.63 K

» 猜你喜欢

第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有110人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有14人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复

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

求助哪位傅里叶红外高人帮忙看一下我的红外图谱已经有18人回复
求助，核磁解析高手，帮忙分析一下。已经有11人回复
【求助】求流变达人帮忙分析一下流变图谢谢已经有12人回复
【求助】求分析一个XRD物相组成，请高手帮忙，谢谢已经有8人回复
【求助】请高手帮忙分析一下这张固化DSC曲线图，非常感谢已经有18人回复
【求助】欧洲药典6.0版或者更新的版本在哪可以免费下载，晓得朋友帮帮忙哈！！已经有3人回复
【求助】求助高手帮忙已经有8人回复
求uvprobe 2.1, 2.2, 2.3版本安装文件(楼主取消求助) 已经有4人回复
【求助】CoMFA 等值图怎么从sybyl X 1.1 083 版本中导出？已经有15人回复
【求助】请高手帮忙设计一下理财计划已经有12人回复

1楼2012-07-11 13:50:44

回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

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

★
dbb627: 金币+1, 感谢应助 2012-07-14 18:47:27

去掉第一组有0的数据：

均方差(RMSE): 0.00476945824112568
残差平方和(SSE): 0.000955404740381351
相关系数(R): 0.99996260092014
相关系数之平方(R^2): 0.999925203238972
决定系数(DC): 0.999925202966366
F统计(F-Statistic): 11441.3414586701

参数最佳估算
-------------------- -------------
h 6.88822200813756
r1 0.306423464581
q1 0.235046609670791
r2 1.2339292576153
q2 1.00223437207462
t12 6.91775076188525
t21 1.88540175538118E-5

赞一下(1人)

6楼2012-07-12 14:47:40

普通回帖

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

感谢参与，应助指数 +1

第一行数据不对吧？ m=0 ---> x1=m/(m+h) --> ln(o1/x1)，会导致被0除的错误。

2楼2012-07-11 21:17:41

xxz903

金虫 (小有名气)

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

引用回帖:

2楼: Originally posted by dingd at 2012-07-11 21:17:41
第一行数据不对吧？ m=0 ---> x1=m/(m+h) --> ln(o1/x1)，会导致被0除的错误。

o1=x1*r1/(x1*r1+x2*r2);、
所以当x1=0时 O1/X1=r1/r2啊

3楼2012-07-11 21:46:40

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
dbb627: 金币+1, 感谢应助 2012-07-12 09:58:23
xxz903: 金币+6, ★★★★★最佳答案 2012-07-12 21:35:28

你的第一个公式：
Function y1=ln(o1/x1)....

里面有“o1/x1”吧？

4楼2012-07-12 09:49:40

xxz903

金虫 (小有名气)

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

引用回帖:

4楼: Originally posted by dingd at 2012-07-12 09:49:40
你的第一个公式：
Function y1=ln(o1/x1)....

里面有“o1/x1”吧？

嗯有的拟合不出结果吗？

5楼2012-07-12 11:52:21

xxz903

金虫 (小有名气)

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

引用回帖:

6楼: Originally posted by dingd at 2012-07-12 14:47:40
去掉第一组有0的数据：

均方差(RMSE): 0.00476945824112568
残差平方和(SSE): 0.000955404740381351
相关系数(R): 0.99996260092014
相关系数之平方(R^2): 0.999925203238972
决定系数(DC): 0.9999252029663 ...

谢了这下结果给力了~ 正在考虑去买个高版本的软件呢~

7楼2012-07-12 15:27:35

xxz903

金虫 (小有名气)

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

引用回帖:

你好我已经购买了5.0的了为啥跑不出这个结果啊？~

8楼2012-09-21 10:19:54

dingd

铁杆木虫 (职业作家)

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

有些难度，运行代码如下，注意第一行为0的数据去除：

CODE:

ParameterDomain = [0,];

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

基本上都能得到不错的结果。

9楼2012-09-21 11:21:53

xxz903

金虫 (小有名气)

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

引用回帖:

9楼: Originally posted by dingd at 2012-09-21 11:21:53
有些难度，运行代码如下，注意第一行为0的数据去除：

ParameterDomain = ;
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*r ...

我还是算不出这个结果你看看我算的就是刚才的编码

10楼2012-09-21 11:40:58