 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1236 | 回復(fù): 2 | |||
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年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有235人回復(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ù)
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
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
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhangholmes 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9 | Ymlll 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +5 | 想上學(xué)吖吖 2026-03-19 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +18 | rare12345 2026-03-18 | 18/900 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5 | Doleres 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 12/600 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+6
|
.6lL 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|