应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 溶解度用已有的模型方程模拟
71/1返回列表
查看: 1592  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-06-26 20:00:05
已阅   回复此楼   关注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朵）

912816088
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 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭