24小时热门版块排行榜    

查看: 359  |  回复: 2

912816088

银虫 (小有名气)

[交流] 1stopt模拟溶解度,模拟结果总不是我想要的结果 已有1人参与

我想用这个公式:logS=logS0 - x*(a - b*T),A1=logS0来模拟溶解度S和x之间的关系,归出其中a,b两个参数是我想要回的系数。我用1stopt 做了一下,但是参数有四个,而且相关度也不好,哪位大神帮我看一下我的代码写的对不?如果哪位做了一下回归的结果不错的话,就把你写的代码留下吧,我做个参考!

这是我的结果:
均方差(RMSE): 0.67605216266479
残差平方和(SSE): 22.852326332187
相关系数(R): 0.889274910454564
相关系数之平方(R^2): 0.790809866363973
决定系数(DC): 0.790808566438579
卡方系数(Chi-Square): 0.334048062168422
F统计(F-Statistic): 181.456328392217

参数        最佳估算
----------        -------------
a                  0.0934029665911945
b                  0.299767785554259
a1                 1.5626372290429
t1                0.303424735051151
以下是我写的代码和数据
Parameters a,b;
Variables x,s;
VarConstant T1=[293.15,303.15,313.15,318.15,323,15],A1=[1.549346636,1.555384688,1.562768543,1.567167577,1.571137936];
ConstStr K = a - b*T1;
Function s=10^(A1-x*K);
Data;
0        35.428
3.645        34.738
5.244        34.413
7.485        34.043
9.952        33.513
12.087        33.189
14.268        32.7
16.91        32.172
20.154        31.605
24.387        30.83
Data;
0        35.924
3.499        35.249
4.924        35.023
7.168        34.538
9.503        34.469
11.772        33.74
13.925        33.312
16.15        32.937
19.253        32.338
23.379        31.544
Data;
0        36.54
3.337        35.916
4.689        35.634
6.84        35.223
9.126        34.814
11.294        34.456
13.563        34.044
15.345        33.8
18.404        33.039
22.47        32.345
Data;
0        36.912
3.216        36.293
4.426        36.076
6.678        35.654
9.014        35.208
10.82        34.897
13.29        34.427
14.743        34.207
18.145        33.533
21.682        32.931
Data;
0        37.251
3.153        36.598
4.251        36.476
6.501        36.027
8.742        35.666
10.472        35.32
13.006        34.903
    14        34.603
17.857        33.972
21.003        33.427
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
csgt0: 金币+2, 谢谢 2013-07-11 15:32:42
版本有问题吧,a1,T1也当成了参数?

用正版的试试,效果很好啊:

均方差(RMSE):0.0662538051010117
残差平方和(SSE):0.219478334518142
相关系数(R): 0.999037637472066
相关系数之平方(R^2): 0.998076201085766
确定系数(DC): 0.99799088343278
F统计(F-Statistic): 7112.47419746166

参数                  最佳估算
--------------------        -------------
a        0.00437812259881813
b        6.63154488021291E-6
2楼2013-07-11 14:25:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

912816088

银虫 (小有名气)


fegg7502: 金币+1, 鼓励交流 2013-07-16 13:22:42
引用回帖:
2楼: Originally posted by dingd at 2013-07-11 14:25:40
版本有问题吧,a1,T1也当成了参数?

用正版的试试,效果很好啊:

均方差(RMSE):0.0662538051010117
残差平方和(SSE):0.219478334518142
相关系数(R): 0.999037637472066
相关系数之平方(R^2): 0.99807620 ...

我能下载下来的模拟的结果都不好,你用的是哪个版本,可否发给我呢,我的邮箱是912816088@qq.com
3楼2013-07-11 19:00:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 912816088 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见