24小时热门版块排行榜    

查看: 1960  |  回复: 6

912816088

银虫 (小有名气)

[求助] 溶解度用已有的模型方程模拟

用此公式回归实验数据:log s = A + BT
                        其中:   A = a + bx + cx2 + dx3 + ex4
                                      B = f + gT + hT2
      目的:确定a,b,c,d,e,f,g,h的值,及相关性的好坏
实验数据如下:
    温度         x               s
298.15           0        0.1481
              0.0417       0.1256
              0.0891       0.1049
              0.1436       0.0855
              0.2069       0.0725
              0.2813        0.0565
              0.3698        0.0453
              0.4772        0.0318
              0.6100       0.0215
              0.7790       0.0094
                      1        0.0018

303.15          0        0.1534
             0.0417        0.1296
             0.0891        0.1103
             0.1436        0.091
             0.2069        0.0755
             0.2813        0.0604
             0.3698        0.0473
             0.4772        0.0342
             0.6100        0.0022
             0.7790           0.0108
                      1        0.0022


308.15           0        0.0158
             0.0417        0.1333
             0.0891        0.1135
             0.1436        0.0954
             0.2069        0.0787
             0.2813        0.0632
             0.3698        0.0493
             0.4772        0.0359
             0.6100        0.0245
             0.7790        0.0118
                      1        0.0016


  313.15         0        0.0162
             0.0417        0.0139
             0.0891        0.012
             0.1436        0.0994
             0.2069        0.082
             0.2813        0.066
             0.3698        0.0514
             0.4772        0.0379
             0.6100         0.0258
             0.7790        0.0128
                       1        0.0031
回复此楼

» 猜你喜欢

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

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

somomo91

专家顾问 (职业作家)

LZ 应该对你的方程做个起码的说明啊
把问题往这里一丢,理直气壮,
然后等别人猜你心思 ?
回答问题 ?
2楼2013-06-27 05:17:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

912816088

银虫 (小有名气)

引用回帖:
2楼: Originally posted by somomo91 at 2013-06-27 05:17:36
LZ 应该对你的方程做个起码的说明啊
把问题往这里一丢,理直气壮,
然后等别人猜你心思 ?
回答问题 ?

哦,我这是在测水-乙醇中碘化钾的溶解度
我就是希望通过这些实验数据,用所给的方程回归出溶解度 s 和温度T、乙醇摩尔分数 x 的一个关系式
3楼2013-06-27 09:47:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
912816088: 金币+6, ★★★很有帮助 2013-06-27 19:15:15
csgt0: 金币+2, 谢谢 2013-06-28 11:35:48
感觉楼主比较粗心,数据似乎有明显的错误:

1:T=303.15时的数据中,倒数第三的”0.0022“实际应该为”0.022“吧
2:T=308.15时的数据中,第一个数据”0.0158“实际应该为”0.158“吧
3:T=313.15时的数据中,前三个数据”0.0162,0.0139, 0.012“实际分别应该为”0.162,0.139, 0.12“吧

如果上述理解正确,可用1stOpt来求解:
CODE:
VarConstant T1=[298.15,303.15,308.15,313.15];
ConstStr A1 = a + b*x + c*x^2 + d*x^3 + e*x^4, B1 = f + g*T1 + h*T1^2;
Variable x,s;
Function s = 10^(A1 + B1*T1);
Data;
//x        s
0        0.1481
0.0417        0.1256
0.0891        0.1049
0.1436        0.0855
0.2069        0.0725
0.2813        0.0565
0.3698        0.0453
0.4772        0.0318
0.6100        0.0215
0.7790        0.0094
1        0.0018
data;
0        0.1534
0.0417        0.1296
0.0891        0.1103
0.1436        0.091
0.2069        0.0755
0.2813        0.0604
0.3698        0.0473
0.4772        0.0342
0.6100        0.022
0.7790        0.0108
1        0.0022
data;
0        0.158
0.0417        0.1333
0.0891        0.1135
0.1436        0.0954
0.2069        0.0787
0.2813        0.0632
0.3698        0.0493
0.4772        0.0359
0.6100        0.0245
0.7790        0.0118
1        0.0016
data;
0        0.162
0.0417        0.139
0.0891        0.12
0.1436        0.0994
0.2069        0.082
0.2813        0.066
0.3698        0.0514
0.4772        0.0379
0.6100        0.0258
0.7790        0.0128
1        0.0031

结果不错:
均方差(RMSE):0.000951040644926142
残差平方和(SSE):3.97970455652674E-5
相关系数(R): 0.999805921734548
相关系数之平方(R^2): 0.99961188113547
确定系数(DC): 0.999611031661381
F统计(F-Statistic): 876.693006971876

参数                  最佳估算
--------------------        -------------
a        -32.3571031158458
b        -1.69843503496331
c        1.23272933123442
d        -0.719608863005286
e        -0.619352305392383
f        0.298905867622026
g        -0.000953475904737373
h        1.02480424084104E-6

====== 结果输出 ======

文件: 数据文件-1
No        目标 s        计算 s
1        0.1481        0.146190725462772
2        0.1256        0.124791385981643
3        0.1049        0.105396846101324
4        0.0855        0.0879124462604165
5        0.0725        0.0722396031281754
6        0.0565        0.0581834610779382
7        0.0453        0.0454404834301919
8        0.0318        0.0334832668790505
9        0.0215        0.0218012281244924
10        0.0094        0.0105111285808733
11        0.0018        0.00229220232731093

文件: 数据文件-2
No        目标 s        计算 s
1        0.1534        0.152174961853353
2        0.1296        0.129899652260905
3        0.1103        0.109711207630726
4        0.091        0.0915110935646911
5        0.0755        0.0751966912780035
6        0.0604        0.0605651688368288
7        0.0473        0.0473005644557682
8        0.0342        0.03485388586667
9        0.022        0.0226936493248726
10        0.0108        0.0109413958085693
11        0.0022        0.00238603236022339

文件: 数据文件-3
No        目标 s        计算 s
1        0.158        0.158013137850682
2        0.1333        0.13488323840842
3        0.1135        0.113920266277607
4        0.0954        0.0950219067985633
5        0.0787        0.0780816042278808
6        0.0632        0.0628887450064626
7        0.0493        0.0491152455090936
8        0.0359        0.0361910514384711
9        0.0245        0.0235642887333974
10        0.0118        0.0113611612785842
11        0.0016        0.00247757223435676

文件: 数据文件-4
No        目标 s        计算 s
1        0.162        0.16396018970614
2        0.139        0.139959763209826
3        0.12        0.118207819452968
4        0.0994        0.098598193016419
5        0.082        0.0810203177780118
6        0.066        0.0652556534342385
7        0.0514        0.0509637684604707
8        0.0379        0.0375531537455052
9        0.0258        0.0244511646535941
10        0.0128        0.0117887549342827
11        0.0031        0.00257081923111779

» 本帖已获得的红花(最新10朵)

4楼2013-06-27 15:09:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

912816088

银虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by dingd at 2013-06-27 15:09:07
感觉楼主比较粗心,数据似乎有明显的错误:

1:T=303.15时的数据中,倒数第三的”0.0022“实际应该为”0.022“吧
2:T=308.15时的数据中,第一个数据”0.0158“实际应该为”0.158“吧
3:T=313.15时的数据中, ...

其实不是我的错,他们给的数据就是这样的,我也觉得有问题
回头在遇到这方面的问题请多多帮忙,留下个QQ号吧,方便求助
5楼2013-06-27 19:14:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

912816088

银虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2013-06-27 15:09:07
感觉楼主比较粗心,数据似乎有明显的错误:

1:T=303.15时的数据中,倒数第三的”0.0022“实际应该为”0.022“吧
2:T=308.15时的数据中,第一个数据”0.0158“实际应该为”0.158“吧
3:T=313.15时的数据中, ...

急问,我下的1stopt怎么不能用呢,没有输入代码的地方呀
6楼2013-07-02 21:59:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujingyu28

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by dingd at 2013-06-27 15:09:07
感觉楼主比较粗心,数据似乎有明显的错误:

1:T=303.15时的数据中,倒数第三的”0.0022“实际应该为”0.022“吧
2:T=308.15时的数据中,第一个数据”0.0158“实际应该为”0.158“吧
3:T=313.15时的数据中, ...

大神,能否指导一下,这个怎么编程!
7楼2016-09-10 17:38:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 912816088 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见