| 10 | 1/1 | 返回列表 |
| 查看: 1339 | 回復(fù): 9 | ||
zhaoshazhu新蟲 (小有名氣)
|
[求助]
求高版本的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 |
鐵桿木蟲 (職業(yè)作家)
|
不知道理解的對(duì)否: 參數(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 |
新蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
新蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
|
均方差(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ù) 最佳估算 -------------------- ------------- 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 |
新蟲 (小有名氣)
|
非常感謝,可以幫我算另外一個(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 |
新蟲 (小有名氣)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 286分 材料求調(diào)劑 +5 | 麻辣魷魚 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +6 | 白QF 2026-03-21 | 7/350 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-27 | 7/350 |
|
|
[考研] 317求調(diào)劑 +5 | 十閑wx 2026-03-24 | 5/250 |
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|