版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

912816088

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1365.4
红花: 2
帖子: 77
在线: 96.2小时
虫号: 1874146
注册: 2012-07-01
性别: MM
专业: 色谱分析

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

用此公式回归实验数据：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

回复此楼

» 猜你喜欢

2026年循环经济功能材料国际会议（ICFMCE 2026）已经有0人回复
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有282人回复
哈尔滨理工大学物理系招收物理学考研调剂已经有5人回复
0702一志愿吉大B区求调剂已经有5人回复
求调剂已经有0人回复
0702一志愿吉大B区求调剂有论文已经有0人回复
请问还有没有用Latex写文章的小伙伴们？已经有0人回复
光学工程学硕调剂信息已经有29人回复
欢迎加入课题组已经有0人回复
散金币，求好运，祝面上顺利！已经有53人回复

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

吸附动力学模型和吸附动力学方程的区别或者联系已经有9人回复
湍流模型的计算结果是如何反馈到动量方程、连续性方程中影响速度、压力的计算的？已经有6人回复
传统傅立叶热传模型与双温方程模型有什么区别？已经有4人回复
ANSYS CFX中应用零方程湍流模型已经有5人回复
fluent里对于轴边界是怎样处理的？对于uds解方程怎样体现模型是二维轴对称模型呢？已经有10人回复
求Matlab仿真模拟一级、二级倒立摆模型方程及matlab仿真模拟程序已经有6人回复
热压缩试验的流变应力模型和本构方程已经有7人回复
如何计算二元方程的模型参数已经有3人回复
怎样用calphad方法计算spinodal分解曲线？已经有5人回复
Gaussian计算和其他方法得到的热力学数据哪个更准确呢已经有11人回复
origin非线性拟合已经有4人回复
【求助】紧束缚模型与薛定谔方程的区别与联系已经有13人回复
【教程】结构方程模型原理及AMOS应用(2008) 已经有12人回复

1楼 2013-06-26 20:00:05

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

somomo91

专家顾问 (职业作家)

专家经验: +3
计算强帖: 1
应助: 144 (高中生)
金币: 19353
红花: 115
帖子: 3854
在线: 1053.2小时
虫号: 2399301
注册: 2013-04-03
专业: 动力学与控制
管辖: 计算模拟

LZ 应该对你的方程做个起码的说明啊
把问题往这里一丢，理直气壮，
然后等别人猜你心思？
回答问题？

赞一下

回复此楼

2楼2013-06-27 05:17:36

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

912816088

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1365.4
红花: 2
帖子: 77
在线: 96.2小时
虫号: 1874146
注册: 2012-07-01
性别: MM
专业: 色谱分析

引用回帖:

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

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1365.4
红花: 2
帖子: 77
在线: 96.2小时
虫号: 1874146
注册: 2012-07-01
性别: MM
专业: 色谱分析

送红花一朵

引用回帖:

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

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1365.4
红花: 2
帖子: 77
在线: 96.2小时
虫号: 1874146
注册: 2012-07-01
性别: MM
专业: 色谱分析

引用回帖:

急问，我下的1stopt怎么不能用呢，没有输入代码的地方呀

赞一下

回复此楼

6楼2013-07-02 21:59:21

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

wujingyu28

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17
帖子: 11
在线: 3.8小时
虫号: 4770523
注册: 2016-06-13
专业: 分离过程

引用回帖:

大神，能否指导一下，这个怎么编程！

赞一下

回复此楼

7楼2016-09-10 17:38:29

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

相关版块跳转我要订阅楼主 912816088 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 085410人工智能专硕317求调剂（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 274求调剂 +5	S.H1 2026-03-18	5/250	2026-03-18 21:27 by guosr9609
[考研] 0703化学调剂，求各位老师收留 +9	秋有木北 2026-03-14	9/450	2026-03-18 19:49 by macy2011
[考研] 化学工程321分求调剂 +15	大米饭！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 0817 化学工程 299分求调剂有科研经历有二区文章 +7	rare12345 2026-03-18	7/350	2026-03-18 14:31 by laoshidan
[考研] 303求调剂 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 268求调剂 +8	一定有学上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 材料与化工专硕调剂 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太监王公公
[考研] 一志愿，福州大学材料专硕339分求调剂 +3	木子momo青争 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 283求调剂 +3	听风就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 热情沙漠
[考研] 11408 一志愿西电，277分求调剂 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 药学383 求调剂 +3	药学chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 326求调剂 +4	诺贝尔化学奖觊� 2026-03-15	7/350	2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 318求调剂 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 复试调剂 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 招收0805（材料）调剂 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 281求调剂 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[考研] 材料专硕350 求调剂 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 274求调剂 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 081200-11408-276学硕求调剂 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求调剂zz