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

登录注册

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

24小时热门版块排行榜

返回列表

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

[求助] 求高版本的1stopt代跑一趟。谢谢已有1人参与

我的反应是在固定床反应器中，每一个条件下只知道停留时间和反应进出口的浓度，求九个参数。总共做了12组实验点进行拟合。
Title "Kinetics";
//Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9,;
//Variable t,x1,x2,x3,x4,x5,x6;
//ODEFunction
         r1 = k7*k1*k2*x1*x2 /(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;
         r2 = k8*k1*k4*x1*x4/((1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;
         r3 = k9*k1**k5x1*k5/((1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;
         dx1/dt = -2*r1-2*r2-r3;
         dx2/dt = -r1;
         dx3/dt = r1+r2;
         dx4/dt = r1-r2;
         dx5/dt = r2-r3;
         dx6/dt = r3;
Data
  t x1in    x2in x3in x4in    x5in x6in    x1out    x2out       x3out    x4out       x5out       x6out
2.05 0.92932 0.00258 0.06810 0 0 0 0.9282 0.0000088 0.0695 0.000031 0.001164 0.001137
1.37 0.93029 0.00155 0.06817 0 0 0 0.9283 0.0000142 0.0697 0.000078 0.001348 0.000573
1.62 0.91758 0.00173 0.08068 0 0 0 0.9150 0.0000168 0.0827 0.000104 0.001727 0.000475
1.98 0.89928 0.00187 0.09884 0 0 0 0.8952 0.0000488 0.1017 0.000289 0.002363 0.000444
1.21 0.91772 0.00158 0.08070 0 0 0 0.9142 0.0000739 0.0831 0.000414 0.001897 0.000280
1.49 0.89895 0.00224 0.09881 0 0 0 0.8952 0.0001940 0.1017 0.000641 0.002019 0.000292
2.70 0.81544 0.00531 0.17925 0 0 0 0.8116 0.0007259 0.1837 0.001304 0.002320 0.000368
1.02 0.92936 0.00254 0.06810 0 0 0 0.9274 0.0000399 0.0703 0.000206 0.001776 0.000361
0.82 0.92918 0.00273 0.06809 0 0 0 0.9272 0.0001846 0.0703 0.000465 0.001608 0.000189
0.97 0.91722 0.00213 0.08065 0 0 0 0.9144 0.0003093 0.0830 0.000790 0.001376 0.000149
1.19 0.89756 0.00379 0.09865 0 0 0 0.8947 0.0007324 0.1017 0.000972 0.001730 0.000210
2.16 0.81436 0.00663 0.17902 0 0 0 0.8103 0.0017511 0.1838 0.001744 0.002119 0.000349

回复此楼

» 猜你喜欢

瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有130人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
物理学I论文润色/翻译怎么收费? 已经有50人回复
PVA溶解已经有2人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复
什么时候开奖？已经有6人回复
散金币祈福已经有94人回复

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

1stopt 语句修改已经有23人回复
求助 1stopt低版本方程的拟合（含若干未知参数）已经有10人回复
关于用1stpro拟合数值的问题。已经有4人回复
1stopt自变量超多时是否有简单语法？AND 自变量多数据量大的话1stopt都会算的特慢么已经有7人回复
1stopt数据拟合求助已经有8人回复
求大神跑下1stOpt，一个sharemodel，谢谢！已经有7人回复
跪求高手帮忙看看这组数据如何进行拟合已经有7人回复
求高版本1stopt跑一个程序已经有7人回复
1stopt解方程组已经有12人回复
用1stOpt进行多模型多变量非线性拟合得不到结果已经有4人回复
求用高版本的1stOpt拟合方程已经有8人回复
1stopt 拟合计算已经有6人回复
求高版本1stopt拟合，已经有11人回复
1stOpt程序的一个小问题已经有5人回复
跪求1stOpt2.5以上版本，或者请高人帮忙拟合包含复数的函数的参数已经有8人回复

1楼 2015-01-03 20:31:42

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

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhaoshazhu: 金币+10, ★★★很有帮助 2015-01-04 16:59:48

不知道理解的对否：

CODE:

Title "Kinetics";

Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;

ParameterDomain = [0,];

Variable t,x1,x2,x3,x4,x5,x6;

ConstStr r1 = k7*k1*k2*x1*x2/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2,

         r2 = k8*k1*k4*x1*x4/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2,

         r3 = k9*k1*k5*x1*k5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;

ODEFunction

           x1' = -2*r1-2*r2-r3;

           x2' = -r1;

           x3' = r1+r2;

           x4' = r1-r2;

           x5' = r2-r3;

           x6' = r3;

Data;

//t        x1in        x2in        x3in        x4in        x5in        x6in

0        0.92932        0.00258        0.06810        0        0        0

2.05        0.9282        0.0000088        0.0695        0.000031        0.001164        0.001137

Data;

0        0.93029        0.00155        0.06817        0        0        0

1.37        0.9283        0.0000142        0.0697        0.000078        0.001348        0.000573

Data;

0        0.91758        0.00173        0.08068        0        0        0

1.62        0.9150        0.0000168        0.0827        0.000104        0.001727        0.000475

Data;

0        0.89928        0.00187        0.09884        0        0        0

1.98        0.8952        0.0000488        0.1017        0.000289        0.002363        0.000444

Data;

0        0.91772        0.00158        0.08070        0        0        0

1.21        0.9142        0.0000739        0.0831        0.000414        0.001897        0.000280

Data;

0        0.89895        0.00224        0.09881        0        0        0

1.49        0.8952        0.0001940        0.1017        0.000641        0.002019        0.000292

Data;

0        0.81544        0.00531        0.17925        0        0        0

2.70        0.8116        0.0007259        0.1837        0.001304        0.002320        0.000368

Data;

0        0.92936        0.00254        0.06810        0        0        0

1.02        0.9274        0.0000399        0.0703        0.000206        0.001776        0.000361

Data;

0        0.92918        0.00273        0.06809        0        0        0

0.82        0.9272        0.0001846        0.0703        0.000465        0.001608        0.000189

Data;

0        0.91722        0.00213        0.08065        0        0        0

0.97        0.9144        0.0003093        0.0830        0.000790        0.001376        0.000149

Data;

0        0.89756        0.00379        0.09865        0        0        0

1.19        0.8947        0.0007324        0.1017        0.000972        0.001730        0.000210

Data;

0        0.81436        0.00663        0.17902        0        0        0

2.16        0.8103        0.0017511        0.1838        0.001744        0.002119        0.000349

参数都要求大于0吧？感觉模型可简化下。

均方差(RMSE):0.00107423565746257
残差平方和(SSE):8.3086721839011E-5

参数最佳估算
-------------------- -------------
k1 1.90714390662338
k2 2448.35355734201
k3 2.70694132686204
k4 603360.71453709
k5 0.0435072616240413
k6 129731241697.305
k7 23.7865519895587
k8 0.143945603891662
k9 0.0565240940660268

赞一下

回复此楼

2楼2015-01-04 14:56:16

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

2楼: Originally posted by dingd at 2015-01-04 14:56:16
不知道理解的对否：

Title "Kinetics";
Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;
ParameterDomain = ;
Variable t,x1,x2,x3,x4,x5,x6;
ConstStr r1 = k7*k1*k2*x1*x2/(1+k1*x1+k2*x2+k3*x3+k4*x4+ ...

实在不好意思。第三个方程是 r3 = k9*k1*k5*x1*x5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;可以再帮我算一次吗？参数全部需要大于零的。

赞一下

回复此楼

3楼2015-01-04 16:44:04

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

非感谢您帮我。就是这样理解的。我这几个参数都需要计算出来所以模型应该是不能简化。

赞一下

回复此楼

4楼2015-01-04 16:59:24

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

可以再帮我算一次吗？实在不好意思我上次给的第三个动力学手误写错了，非常感谢。

赞一下

回复此楼

5楼2015-01-05 10:57:38

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

dingd

铁杆木虫 (职业作家)

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

引用回帖:

5楼: Originally posted by zhaoshazhu at 2015-01-05 10:57:38
可以再帮我算一次吗？实在不好意思我上次给的第三个动力学手误写错了，非常感谢。

...

求解代码里r3那一句已经改过了。

赞一下

回复此楼

6楼2015-01-05 14:37:12

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

6楼: Originally posted by dingd at 2015-01-05 14:37:12
求解代码里r3那一句已经改过了。...

我第一次是写的r3 = k9*k1*k5*x1*k5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;，您的代码里也是这样的，分子有两个k5，需要把其中的一个改为x5.

赞一下

回复此楼

7楼2015-01-05 16:03:57

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

dingd

铁杆木虫 (职业作家)

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

引用回帖:

7楼: Originally posted by zhaoshazhu at 2015-01-05 16:03:57
我第一次是写的r3 = k9*k1*k5*x1*k5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;，您的代码里也是这样的，分子有两个k5，需要把其中的一个改为x5....

均方差(RMSE):0.00106616879873234
残差平方和(SSE):8.18435453321056E-5
相关系数(R): 0
相关系数之平方(R^2): 0
确定系数(DC): 0
F统计(F-Statistic): 0

参数最佳估算
-------------------- -------------
k1 2.65934369026499
k2 58012.3883745
k3 12.3953925501171
k4 910.506489281551
k5 1139098.80593599
k6 310.722381455093
k7 0.630067032131075
k8 77.4677237219543
k9 3.42700408019413E-16

赞一下

回复此楼

8楼2015-01-05 16:52:09

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

8楼: Originally posted by dingd at 2015-01-05 16:52:09
均方差(RMSE):0.00106616879873234
残差平方和(SSE):8.18435453321056E-5
相关系数(R): 0
相关系数之平方(R^2): 0
确定系数(DC): 0
F统计(F-Statistic): 0

参数最佳估算
-------------------- ...

非常感谢，可以帮我算另外一个模型的吗？数据一样的，只是没有x6组分，模型不一样
Title "Kinetics";
Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;
ParameterDomain = [0,];
Variable t,x1,x2,x3,x4,x5;
ConstStr
r1 = k(6).*((k(1).*x(1).*k(2).*x(2)).^0.5)./...
(1+2.*(k(2).*x(2)).^0.5+2.*(k(4).*x(4)).^0.5+(k(1).*x(1)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;
r2 = k(7).*((k(1).*x(1).*k(4).*x(4)).^0.5)./...
(1+2.*(k(2).*x(2)).^0.5+2.*(k(4).*x(4)).^0.5+(k(1).*x(1)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;

r3 = k(8).*k(5).*k(1).*x(5).*x(1)./...
(1+2.*(k(1).*x(1)).^0.5+2.*(k(4).*x(4)).^0.5+(k(2).*x(2)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;
ODEFunction
         x1' = -2*r1-2*r2-r3;
         x2' = -r1;
         x3' = r1+r2;
         x4' = r1-r2;
         x5' = r2-r3;
Data;
//t       x1in       x2in       x3in       x4in       x5in
0       0.92932       0.00258       0.06810       0       0
2.05       0.9282       0.0000088       0.0695       0.000031       0.001164
Data;
0       0.93029       0.00155       0.06817       0       0
1.37       0.9283       0.0000142       0.0697       0.000078       0.001348
Data;
0       0.91758       0.00173       0.08068       0       0
1.62       0.9150       0.0000168       0.0827       0.000104       0.001727
Data;
0       0.89928       0.00187       0.09884       0       0
1.98       0.8952       0.0000488       0.1017       0.000289       0.002363
Data;
0       0.91772       0.00158       0.08070       0       0
1.21       0.9142       0.0000739       0.0831       0.000414       0.001897
Data;
0       0.89895       0.00224       0.09881       0       0
1.49       0.8952       0.0001940       0.1017       0.000641       0.002019
Data;
0       0.81544       0.00531       0.17925       0       0
2.70       0.8116       0.0007259       0.1837       0.001304       0.002320
Data;
0       0.92936       0.00254       0.06810       0       0
1.02       0.9274       0.0000399       0.0703       0.000206       0.001776
Data;
0       0.92918       0.00273       0.06809       0       0
0.82       0.9272       0.0001846       0.0703       0.000465       0.001608
Data;
0       0.91722       0.00213       0.08065       0       0
0.97       0.9144       0.0003093       0.0830       0.000790       0.001376
Data;
0       0.89756       0.00379       0.09865       0       0
1.19       0.8947       0.0007324       0.1017       0.000972       0.001730       0.000210
Data;
0       0.81436       0.00663       0.17902       0       0
2.16       0.8103       0.0017511       0.1838       0.001744       0.002119

赞一下

回复此楼

9楼2015-01-05 19:16:57

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

很想得到您的帮助，可以帮我再算另一个模型的参数估算吗？真的很感谢

赞一下

回复此楼

10楼2015-01-07 21:26:53

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

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

返回列表

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

[基金申请] 估计是周四 +3	archvillain 2026-08-18	3/150	2026-08-21 01:48 by jnhyjjm
[基金申请] 基金啊基金 +3	longfie172 2026-08-20	3/150	2026-08-21 00:04 by wuchongjun
[基金申请] 科研孤儿太难了 +14	我4大白菜 2026-08-20	15/750	2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧！ +8	大萍1987 2026-08-20	10/500	2026-08-20 21:48 by 大萍1987
[基金申请] 言之凿凿，明天官网会通知放榜？ +15	医学老男孩 2026-08-20	17/850	2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12	wuchongjun 2026-08-20	17/850	2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护，明天很可能放榜！ +10	zju2000 2026-08-16	11/550	2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18	基诺咪客 2026-08-17	23/1150	2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21	yaoyewhu2008 2026-08-20	22/1100	2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3	snowwithsea 2026-08-19	12/600	2026-08-20 15:56 by lfy8008
[基金申请] 今天基金会出结果吗？20260819 +14	kkkl_v 2026-08-19	15/750	2026-08-20 08:21 by sunzitan
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5	瞬息宇宙 2026-08-15	5/250	2026-08-19 20:23 by super2002521
[基金申请] filecode=后面第一个是大写字母 +10	wangze12014 2026-08-14	12/600	2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6	wangzilk 2026-08-18	8/400	2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖？ +6	CrisMessi 2026-08-18	6/300	2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜？ +5	Shxjjxjkx 2026-08-18	5/250	2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4	wwm09 2026-08-17	6/300	2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14	archvillain 2026-08-15	26/1300	2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护祝所有人高中 +8	gjjjzhong 2026-08-18	9/450	2026-08-18 13:01 by 家与远方
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7	Tide man 2026-08-14	8/400	2026-08-15 01:11 by arzu_hma

24小时热门版块排行榜

[求助] 求高版本的1stopt代跑一趟。谢谢 已有1人参与

» 猜你喜欢

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

【答案】应助回帖

[求助] 求高版本的1stopt代跑一趟。谢谢已有1人参与