| 2 | 1/1 | 返回列表 |
| 查看: 1234 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
zhangholmes至尊木蟲 (著名寫手)
逍遙門主
|
[求助]
matlab非線性最優(yōu)化求解,微分方程組的擬合,參數(shù)估計(jì)
|
|
|
如題,用微分方程組擬合4組不同溫度下的數(shù)據(jù),求理論值a和b,并使a和b的值分別與實(shí)驗(yàn)值y1和y2的差方和最小。并求此時(shí)的各個(gè)參數(shù)值(k1,k2,k3,k4,k5,k6),求各位大俠略施援手。 0<k1<30, 0<k2<10, 0<k3<5, 0<k4<20, 0<k5<20, 0<k6<5 t=[0.494660000000000,0.502760000000000,0.618930000000000,0.621470000000000,0.623470000000000,0.819730000000000,0.825080000000000,0.828350000000000]; T=463 K時(shí), y1= [0.00800000000000000,0.00990000000000000,0.00782000000000000,0.0108400000000000,0.00937000000000000,0.00882000000000000,0.00614000000000000,0.0108200000000000]; y2=[0.0625100000000000,0.0727400000000000,0.0540300000000000,0.0667600000000000,0.0620200000000000,0.0494400000000000,0.0452000000000000,0.0509400000000000]; %y1和y2分別為a和b的實(shí)驗(yàn)值 T=473 K時(shí), y1=[0.00829000000000000,0.0107700000000000,0.00938000000000000,0.0172000000000000,0.0128100000000000,0.0139200000000000,0.0108200000000000,0.0179000000000000]; y2=[0.0622600000000000,0.0719200000000000,0.0525700000000000,0.0604900000000000,0.0586200000000000,0.0444300000000000,0.0405900000000000,0.0439200000000000]; T=483 K時(shí), y1=[0.0169300000000000,0.0229900000000000,0.0188900000000000,0.0313900000000000,0.0265800000000000,0.0258100000000000,0.0194200000000000,0.0293300000000000]; y2=[0.0537100000000000,0.0598300000000000,0.0431500000000000,0.0465000000000000,0.0450300000000000,0.0326700000000000,0.0321200000000000,0.0327000000000000]; T=493 K時(shí), y1=[0.0328500000000000,0.0418500000000000,0.0330300000000000,0.0512500000000000,0.0424800000000000,0.0389500000000000,0.0300900000000000,0.0443800000000000]; y2=[0.0381700000000000,0.0415400000000000,0.0292800000000000,0.0273400000000000,0.0295700000000000,0.0197800000000000,0.0216900000000000,0.0180800000000000]; pA = 101*(1 - (2 b + a)); pW = 101*(a+ b); da/dt=k(4).*k(1).*k(2).*k(3).*pA.*pW./(k(1).*k(2).*k(4).*k(5).*pA^2+k(1).*k(2).*k(3).*pA.*pW+k(4)+k(1).*k(2).*k(4).*pA.*pW+k(1).*k(4).*pW; db/dt=k(6).*k(1).*k(2).*k(4).*k(5).*pA^2./(k(1).*k(2).*k(4).*k(5).*pA^2+k(1).*k(2).*k(3).*pA.*pW+k(4)+k(1).*k(2).*k(4).*pA.*pW+k(1).*k(4).*pW; [ 來自科研家族 木蟲催化 ] |

鐵桿木蟲 (職業(yè)作家)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[論文投稿]
申請回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|