 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1235 | 回復: 2 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
zhangholmes至尊木蟲 (著名寫手)
逍遙門主
|
[求助]
matlab非線性最優(yōu)化求解，微分方程組的擬合，參數(shù)估計
|
|
|
如題，用微分方程組擬合4組不同溫度下的數(shù)據(jù)，求理論值a和b，并使a和b的值分別與實驗值y1和y2的差方和最小。并求此時的各個參數(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時， 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的實驗值 T=473 K時， 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時， 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時， 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; [ 來自科研家族 木蟲催化 ] |
回復此樓» 猜你喜歡
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復
-
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有6人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有269人回復
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復
-
求調(diào)劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復
-
光學工程學碩調(diào)劑信息
已經(jīng)有30人回復
-
歡迎加入課題組
已經(jīng)有0人回復
-
散金幣，求好運，祝面上順利！
已經(jīng)有63人回復
-
華東師范大學芯片設計徐瓏真實水平咋樣
已經(jīng)有1人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 用Matlab做多變量非線性擬合求參數(shù)值，求各位大俠鼎力相助！
已經(jīng)有3人回復
- 微分方程組求參數(shù)問題，求高人指教，感謝
已經(jīng)有11人回復
- 微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計）
已經(jīng)有11人回復
- 關(guān)于matlab微分、及非線性擬合問題
已經(jīng)有5人回復
- 二階非線性微分方程求解
已經(jīng)有14人回復
- 請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。
已經(jīng)有17人回復
- 請教用matlab求解一個非線性偏微分方程組的數(shù)值解
已經(jīng)有8人回復
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應曲線？
已經(jīng)有13人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復
- 求助，origin7.5 中的非線性擬合參數(shù)為什么一直在變，而且誤差很大�。�！
已經(jīng)有4人回復
- 求助：origin8.0 自定義非線性擬合，參數(shù)擬合值與設定初始值一樣
已經(jīng)有18人回復
- 【求助】用mathematica 5.0求解一個非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復
- 【求助】非線性參數(shù)的擬合方法或者擬合軟件？（分母、指數(shù)和對數(shù)）
已經(jīng)有3人回復
- 【分享】PDECOL-經(jīng)典數(shù)值求解非線性偏微分方程組源程序
已經(jīng)有142人回復
zhangholmes
至尊木蟲 (著名寫手)
逍遙門主
- 應助: 17 (小學生)
- 金幣: 17040
- 散金: 50
- 紅花: 12
- 帖子: 1242
- 在線: 226.3小時
- 蟲號: 967701
- 注冊: 2010-03-11
- 性別: GG
- 專業(yè): 能源化工
|
謝謝回復，公式最后忘記加了一個括號。不要求統(tǒng)一參數(shù)，每組數(shù)據(jù)都有一組獨立的參數(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ā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +6 | 墨墨漠 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 329求調(diào)劑 +5 | 想上學吖吖 2026-03-19 | 5/250 |
|
|
[考研] 081700化工學碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +14 | enenenhui 2026-03-13 | 15/750 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +3 | Ncdx123456 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 本科鄭州大學物理學院，一志愿華科070200學碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|