版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

zhaoshazhu

新蟲 (小有名氣)

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

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

我的反應(yīng)是在固定床反應(yīng)器中，每一個(gè)條件下只知道停留時(shí)間和反應(yīng)進(jìn)出口的濃度，求九個(gè)參數(shù)�？偣沧隽�12組實(shí)驗(yàn)點(diǎn)進(jìn)行擬合。
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

回復(fù)此樓

» 猜你喜歡

中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有76人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有1人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

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

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

zhaoshazhu

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 125.7
散金: 10
帖子: 51
在線: 38小時(shí)
蟲號(hào): 1812564
注冊(cè): 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
相關(guān)系數(shù)(R): 0
相關(guān)系數(shù)之平方(R^2): 0
確定系數(shù)(DC): 0
F統(tǒng)計(jì)(F-Statistic): 0

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

非常感謝，可以幫我算另外一個(gè)模型的嗎？數(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

贊一下

回復(fù)此樓

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

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

查看全部 10 個(gè)回答

dingd

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

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
zhaoshazhu: 金幣+10, ★★★很有幫助 2015-01-04 16:59:48

不知道理解的對(duì)否：

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吧？感覺模型可簡(jiǎn)化下。

均方差(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

贊一下

回復(fù)此樓

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

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

zhaoshazhu

新蟲 (小有名氣)

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

引用回帖:

2樓: Originally posted by dingd at 2015-01-04 14:56:16
不知道理解的對(duì)否：

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+ ...

實(shí)在不好意思。第三個(gè)方程是 r3 = k9*k1*k5*x1*x5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;可以再幫我算一次嗎？參數(shù)全部需要大于零的。

贊一下

回復(fù)此樓

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

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

zhaoshazhu

新蟲 (小有名氣)

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

引用回帖:

非感謝您幫我。就是這樣理解的。我這幾個(gè)參數(shù)都需要計(jì)算出來所以模型應(yīng)該是不能簡(jiǎn)化。

贊一下

回復(fù)此樓

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

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

查看全部 10 個(gè)回答

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

[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +10	leaflight 2026-03-24	10/500	2026-03-29 20:58 by goldfish51
[考研] 291求調(diào)劑 +7	Y-cap 2026-03-29	8/400	2026-03-29 19:53 by klasasda
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8	mikasa的圍巾 2026-03-28	8/400	2026-03-29 12:48 by 無際的草原
[考研] 348求調(diào)劑 +5	小懶蟲不懶了 2026-03-28	5/250	2026-03-29 10:34 by 唐沐兒
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 086502化學(xué)工程342求調(diào)劑 +6	阿姨復(fù)古不過 2026-03-27	6/300	2026-03-28 07:06 by wangy0907
[考研] 266求調(diào)劑 +11	陽陽哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4	Ripcord上岸 2026-03-27	4/200	2026-03-27 17:00 by zhyzzh
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +5	SYA！ 2026-03-23	5/250	2026-03-27 09:29 by 不吃魚的貓
[考研] 317求調(diào)劑 +7	蛋黃咸肉粽 2026-03-26	7/350	2026-03-27 02:29 by fmesaito
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 打過很多競(jìng)賽，085406控制工程300分，求調(diào)劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑！�。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 一志愿國(guó)科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

zhaoshazhu

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

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

zhaoshazhu

dingd

【答案】應(yīng)助回帖

zhaoshazhu

zhaoshazhu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: