| 10 | 1/1 | 返回列表 |
| 查看: 1458 | 回復(fù): 9 | |||
dragonpu木蟲 (小有名氣)
|
[求助]
動力學(xué)參數(shù)擬合
|
|
有哪位大神能夠幫我擬合一下動力學(xué)參數(shù)。atlab程序我編好了,但是結(jié)果非常依賴于初值。小弟不勝感激! 動力學(xué)方程組: dy1/dw=-k(1).*exp(-k(2)./(R.*T)).*P.*y(1)-k(3).*exp(-k(4)./(R.*T)).*P.*y(1)-k(7).*exp(-k(8)./(R.*T)).*P.*y(1)-k(9).*exp(-k(10)./(R.*T)).*P.*y(1); dy2/dw=k(1).*exp(-k(2)./(R.*T)).*P.*y(1)+k(5).*exp(-k(6)./(R.*T)).*P.*y(3); dy3/dw=k(3).*exp(-k(4)./(R.*T)).*P.*y(1)-k(5).*exp(-k(6)./(R.*T)).*P.*y(3)-k(11).*exp(-k(12)./(R.*T)).*P.*y(3); dy4/dw=k(7).*exp(-k(8)./(R.*T)).*P.*y(1); dy5/dw=k(9).*exp(-k(10)./(R.*T)).*P.*y(1)+k(11).*exp(-k(12)./(R.*T)).*P.*y(3); 其中P=100; R=8.314e-3; k(1)~k(12)為要擬合的參數(shù) y的初值為[0.5 0 0 0 0] 實驗數(shù)據(jù)為 y1=[0.480097655 0.474124441 0.462780469 0.451193215 0.440723531 0.429448481 0.47974112 0.472322212 0.459982252 0.446630481 0.435406437 0.422382655 0.477436093 0.470092739 0.455851064 0.442346819 0.430467491 0.418103342]'; y2=[0.002842264 0.003831889 0.005528585 0.007890192 0.008760934 0.010302604 0.002846129 0.004434475 0.006468111 0.009272929 0.010353186 0.012069968 0.004068987 0.005286294 0.007419716 0.010652046 0.012228119 0.013803879]'; y3=[0.00347535 0.004900899 0.007194984 0.009158732 0.011510271 0.013775715 0.002932906 0.004427726 0.006607841 0.008459383 0.010480658 0.013100752 0.002395948 0.003717976 0.006043769 0.007540305 0.009972262 0.011736397]'; y4=[0.000368155 0.000553876 0.001417951 0.002258876 0.003408988 0.004321705 0.000816502 0.00160397 0.002537401 0.003319669 0.004791179 0.006212197 0.002013246 0.002162904 0.00365511 0.005083102 0.006746324 0.007918681 ]'; y5=[0.00146887 0.001936815 0.003335547 0.003735631 0.005773633 0.007057201 0.002564792 0.002891745 0.004969817 0.006502123 0.008844523 0.009879428 0.003020119 0.004579988 0.006845988 0.00900583 0.009735804 0.013336633 ]'; yout=[y1;y2;y3;y4;y5]; t=[375 375 375 375 375 375 400 400 400 400 400 400 425 425 425 425 425 425]; T=t+273.15; w=[2.531646 3.375527 5.063291 6.751055 8.438819 10.12658 2.531646 3.375527 5.063291 6.751055 8.438819 10.12658 2.531646 3.375527 5.063291 6.751055 8.438819 10.12658]; |
鐵桿木蟲 (職業(yè)作家)
|
1stOpt求解: 均方差(RMSE): 0.00515849113704793 殘差平方和(SSE): 0.00239490277299019 相關(guān)系數(shù)(R): 0.999823871742008 相關(guān)系數(shù)之平方(R^2): 0.999647774505179 決定系數(shù)(DC): 0.999850282157082 F統(tǒng)計(F-Statistic): 6.45801799379734 參數(shù) 最佳估算 -------------------- ------------- k1 4.89399891352048E-6 k2 -10.6601371015196 k3 3.67383743863496E-5 k4 -1.43251983409417 k5 2.79201648997203E-8 k6 -49.2290083090162 k7 3.52249551424626E-6 k8 -10.9413062151302 k9 2.95332341973282E-5 k10 -1.07700586267194 k11 1.05224743887147E-8 k12 2.06702669054423E-8 |
木蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
木蟲 (正式寫手)
|
均方差(RMSE): 0.00515849113704793 殘差平方和(SSE): 0.00239490277299019 相關(guān)系數(shù)(R): 0.999823871742008 相關(guān)系數(shù)之平方(R^2): 0.999647774505179 決定系數(shù)(DC): 0.999850282157082 F統(tǒng)計(F-Statistic): 6.45801799379734 請問這個是怎么擬合的?能否把擬合代碼發(fā)給我看看,wczheng@163.com |
木蟲 (正式寫手)
| 請問樓主擬合好了嗎?我現(xiàn)在也有一個擬合,有17個反應(yīng)速率常數(shù)需要擬合,想借鑒一下你的經(jīng)驗和matlab代碼,wczheng@163.com |
鐵桿木蟲 (職業(yè)作家)
木蟲 (正式寫手)
木蟲 (正式寫手)
木蟲 (正式寫手)
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +11 | 陌の森林 2026-03-18 | 11/550 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +6 | Ymlll 2026-03-18 | 9/450 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|