版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

[求助] 求高版本的1stopt代跑一趟。謝謝已有1人參與

我的反應是在固定床反應器中，每一個條件下只知道停留時間和反應進出口的濃度，求九個參數(shù)�？偣沧隽�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

回復此樓

» 猜你喜歡

調劑已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有225人回復
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。已經(jīng)有0人回復
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！已經(jīng)有1人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人已經(jīng)有0人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

1stopt 語句修改已經(jīng)有23人回復
求助 1stopt低版本方程的擬合（含若干未知參數(shù)）已經(jīng)有10人回復
關于用1stpro擬合數(shù)值的問題。已經(jīng)有4人回復
1stopt自變量超多時是否有簡單語法？AND 自變量多數(shù)據(jù)量大的話1stopt都會算的特慢么已經(jīng)有7人回復
1stopt數(shù)據(jù)擬合求助已經(jīng)有8人回復
求大神跑下1stOpt，一個sharemodel，謝謝！已經(jīng)有7人回復
跪求高手幫忙看看這組數(shù)據(jù)如何進行擬合已經(jīng)有7人回復
求高版本1stopt跑一個程序已經(jīng)有7人回復
1stopt解方程組已經(jīng)有12人回復
用1stOpt進行多模型多變量非線性擬合得不到結果已經(jīng)有4人回復
求用高版本的1stOpt擬合方程已經(jīng)有8人回復
1stopt 擬合計算已經(jīng)有6人回復
求高版本1stopt擬合，已經(jīng)有11人回復
1stOpt程序的一個小問題已經(jīng)有5人回復
跪求1stOpt2.5以上版本，或者請高人幫忙擬合包含復數(shù)的函數(shù)的參數(shù) 已經(jīng)有8人回復

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +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

參數(shù)都要求大于0吧？感覺模型可簡化下。

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

參數(shù) 最佳估算
-------------------- -------------
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發(fā)消息送TA紅花 TA的回帖

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

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;可以再幫我算一次嗎？參數(shù)全部需要大于零的。

贊一下

回復此樓

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

非感謝您幫我。就是這樣理解的。我這幾個參數(shù)都需要計算出來所以模型應該是不能簡化。

贊一下

回復此樓

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

可以再幫我算一次嗎？實在不好意思我上次給的第三個動力學手誤寫錯了，非常感謝。

贊一下

回復此樓

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 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那一句已經(jīng)改過了。

贊一下

回復此樓

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

6樓: Originally posted by dingd at 2015-01-05 14:37:12
求解代碼里r3那一句已經(jīng)改過了。...

我第一次是寫的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發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計算強帖: 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
相關系數(shù)(R): 0
相關系數(shù)之平方(R^2): 0
確定系數(shù)(DC): 0
F統(tǒng)計(F-Statistic): 0

參數(shù) 最佳估算
-------------------- -------------
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發(fā)消息送TA紅花 TA的回帖

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

8樓: Originally posted by dingd at 2015-01-05 16:52:09
均方差(RMSE):0.00106616879873234
殘差平方和(SSE):8.18435453321056E-5
相關系數(shù)(R): 0
相關系數(shù)之平方(R^2): 0
確定系數(shù)(DC): 0
F統(tǒng)計(F-Statistic): 0

參數(shù) 最佳估算
-------------------- ...

非常感謝，可以幫我算另外一個模型的嗎？數(shù)據(jù)一樣的，只是沒有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發(fā)消息送TA紅花 TA的回帖

zhaoshazhu

新蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時
蟲號: 1812564
注冊: 2012-05-12
專業(yè): 環(huán)境化工

引用回帖:

很想得到您的幫助，可以幫我再算另一個模型的參數(shù)估算嗎？真的很感謝

贊一下

回復此樓

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 zhaoshazhu 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 295求調劑 +4	wei-5 2026-03-26	4/200	2026-03-28 23:20 by 小木蟲tim
[考研] 一志愿211 初試270分求調劑 +7	谷雨上岸 2026-03-23	8/400	2026-03-28 21:47 by lbsjt
[考研] 0703化學 +10	妮妮ninicgb 2026-03-27	10/500	2026-03-28 14:07 by 噠噠噠呱呱呱
[考研] 求化學調劑 +4	wulanna 2026-03-28	4/200	2026-03-28 13:37 by 唐沐兒
[考研] 085600 286分材料求調劑 +7	麻辣魷魚 2026-03-27	8/400	2026-03-28 12:17 by zllcz
[考研] 291求調劑 +15	hhhhxn.. 2026-03-23	21/1050	2026-03-28 11:26 by self2008
[考研] 286求調劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 086000調劑 +3	7901117076 2026-03-26	3/150	2026-03-27 21:34 by Jianing_Mi
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調劑 +10	幸運的醬醬 2026-03-22	12/600	2026-03-27 16:08 by muchong357
[考研] 一志愿北京化工大學 070300 學碩 336分求調劑 +9	vv迷 2026-03-22	9/450	2026-03-27 15:59 by 不吃魚的貓
[考研] 一志愿華東理工大學081700，初試分數(shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 329求調劑 +7	鈕恩雪 2026-03-25	7/350	2026-03-27 04:28 by wxiongid
[考研] 321求調劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 325求調劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 297求調劑 +6	田洪有 2026-03-26	6/300	2026-03-26 15:55 by 不吃魚的貓
[考研] 打過很多競賽，085406控制工程300分，求調劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息
[考研] 求b區(qū)院校調劑 +4	周56 2026-03-24	5/250	2026-03-25 17:12 by yishunmin
[考研] 284求調劑 +15	Zhao anqi 2026-03-22	15/750	2026-03-25 12:51 by wht0531
[考研] 一志愿吉大化學322求調劑 +4	17501029541 2026-03-23	6/300	2026-03-24 10:21 by 戴圍脖的小蚊子
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3	lishucheng1 2026-03-22	5/250	2026-03-24 08:22 by gltch

24小時熱門版塊排行榜

[求助] 求高版本的1stopt代跑一趟。謝謝 已有1人參與

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

【答案】應助回帖

[求助] 求高版本的1stopt代跑一趟。謝謝已有1人參與

» 本主題相關價值貼推薦，對您同樣有幫助: