當(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ù)此樓

» 猜你喜歡

計(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)有4人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有197人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名已經(jīng)有0人回復(fù)

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

1stopt 語(yǔ)句修改已經(jīng)有23人回復(fù)
求助 1stopt低版本方程的擬合（含若干未知參數(shù)）已經(jīng)有10人回復(fù)
關(guān)于用1stpro擬合數(shù)值的問題。已經(jīng)有4人回復(fù)
1stopt自變量超多時(shí)是否有簡(jiǎn)單語(yǔ)法？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)境化工

引用回帖:

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

非感謝您幫我。就是這樣理解的。我這幾個(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è)回答

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)境化工

引用回帖:

實(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)境化工

引用回帖:

可以再幫我算一次嗎？實(shí)在不好意思我上次給的第三個(gè)動(dòng)力學(xué)手誤寫錯(cuò)了，非常感謝。

贊一下

回復(fù)此樓

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

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

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

[考研] 336求調(diào)劑 +5	kiyy 2026-04-01	5/250	2026-04-02 07:14 by chixmc
[考研] 085801 總分275 本科新能源求調(diào)劑 +10	bradoner 2026-04-01	11/550	2026-04-01 22:57 by 來看流星雨10
[考研] 建環(huán)，能源，土木老師路過看一看�。。� +4	嘿嘿uu 2026-04-01	4/200	2026-04-01 20:42 by 無懈可擊111
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +14	一定上岸aaa 2026-03-27	18/900	2026-04-01 20:10 by 積極調(diào)劑的小學(xué)?/a>
[考研] 一志愿西安交大材料學(xué)碩（英一數(shù)二）347，求調(diào)劑到高分子/材料相關(guān)專業(yè) +7	zju51 2026-03-31	9/450	2026-04-01 19:35 by CFQZAFU
[考研] 348求調(diào)劑 +9	zzzzyk123 2026-04-01	9/450	2026-04-01 14:24 by w蟲蟲123
[考研] 材料專業(yè)求調(diào)劑 +7	月月鳥木 2026-04-01	7/350	2026-04-01 11:36 by wangjy2002
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +4	元元青青 2026-03-31	6/300	2026-04-01 11:13 by syh9288
[考研] 一志愿華南師范361分，化學(xué)求調(diào)劑 +4	Nicole88888 2026-04-01	4/200	2026-04-01 10:08 by 唐沐兒
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +10	@taotao 2026-03-29	10/500	2026-04-01 09:30 by oooqiao
[考研] 求調(diào)劑：一志愿：南京大學(xué) 專業(yè)：0705 總分320 ，本科985，四六級(jí)已過 +3	lfy760306 2026-03-31	3/150	2026-04-01 01:57 by Creta
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +4	1段時(shí)間 2026-03-31	5/250	2026-03-31 19:56 by wxiongid
[考研] 375求調(diào)劑 +7	雨夏整夜 2026-03-29	7/350	2026-03-31 18:52 by xhai2011
[考研] 286求調(diào)劑 +6	Faune 2026-03-30	6/300	2026-03-31 14:37 by jp9609
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	7/350	2026-03-30 18:29 by 兔兔兔111223314
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 2026年華南師范大學(xué)歡迎化學(xué)，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入！ +3	llss0711 2026-03-28	6/300	2026-03-29 10:26 by llss0711
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 藥學(xué)105500求調(diào)劑 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry