 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1237 | 回復(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; [ 來自科研家族 木蟲催化 ] |
回復(fù)此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有160人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有30人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有63人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 用Matlab做多變量非線性擬合求參數(shù)值，求各位大俠鼎力相助！
已經(jīng)有3人回復(fù)
- 微分方程組求參數(shù)問題，求高人指教，感謝
已經(jīng)有11人回復(fù)
- 微分方程組參數(shù)擬合的問題（多參數(shù)動(dòng)力系統(tǒng)參數(shù)估計(jì)）
已經(jīng)有11人回復(fù)
- 關(guān)于matlab微分、及非線性擬合問題
已經(jīng)有5人回復(fù)
- 二階非線性微分方程求解
已經(jīng)有14人回復(fù)
- 請(qǐng)問這個(gè)4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。
已經(jīng)有17人回復(fù)
- 請(qǐng)教用matlab求解一個(gè)非線性偏微分方程組的數(shù)值解
已經(jīng)有8人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較
已經(jīng)有4人回復(fù)
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復(fù)
- 求助，origin7.5 中的非線性擬合參數(shù)為什么一直在變，而且誤差很大�。�！
已經(jīng)有4人回復(fù)
- 求助：origin8.0 自定義非線性擬合，參數(shù)擬合值與設(shè)定初始值一樣
已經(jīng)有18人回復(fù)
- 【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 【求助】非線性參數(shù)的擬合方法或者擬合軟件？（分母、指數(shù)和對(duì)數(shù)）
已經(jīng)有3人回復(fù)
- 【分享】PDECOL-經(jīng)典數(shù)值求解非線性偏微分方程組源程序
已經(jīng)有142人回復(fù)
zhangholmes
至尊木蟲 (著名寫手)
逍遙門主
- 應(yīng)助: 17 (小學(xué)生)
- 金幣: 17040
- 散金: 50
- 紅花: 12
- 帖子: 1242
- 在線: 226.3小時(shí)
- 蟲號(hào): 967701
- 注冊(cè): 2010-03-11
- 性別: GG
- 專業(yè): 能源化工
|
謝謝回復(fù)，公式最后忘記加了一個(gè)括號(hào)。不要求統(tǒng)一參數(shù)，每組數(shù)據(jù)都有一組獨(dú)立的參數(shù)，共四組參數(shù)。 公式如下： 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）; |
3樓2013-05-12 15:59:36
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
2樓2013-05-12 14:34:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +18 | rare12345 2026-03-18 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過，有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|