 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1843 | 回復(fù): 4 | ||
[求助]
求助1stopt擬合動(dòng)力學(xué)參數(shù) 已有1人參與
|
|
最近做反應(yīng)動(dòng)力學(xué)模型，需要擬合微分方程組的幾個(gè)參數(shù)，之前用matlab，ode45求解，得到的數(shù)值解套用非線性最小二乘擬合，可是總也擬合不好。不知道是初值不合適還是模型有問題，但模型之前有文獻(xiàn)報(bào)道過，有問題的可能性較小。 想求助用1stopt擬合試試，看看能不能得到比較好的效果~ 我對(duì)1stopt 不太熟，試著寫了下程序，不知道有沒有問題，具體如下： Parameter a,b,c,d,e,f,g; Variable t,x,y,z,l,m,n; Function x’=-a*x*n*(0.000384369/(1+d*x+e*y+f*z+g*m))； y’=(a*x*n-b*y*n)*(0.000384369/(1+d*x+e*y+f*z+g*m))； z’=(b*y*n-c*z*n)*(0.000384369/(1+d*x+e*y+f*z+g*m))； l’=c*z*n*(0.000384369/(1+d*x+e*y+f*z+g*m))； m’=(a*x+b*y+c*z)*(0.000384369/(1+d*x+e*y+f*z+g*m))； n’= -m’； RowData； t=0 5 10 20 30 40 60; x=0.930682749 0.598029237 0.541233252 0.502127725 0.46364134 0.442517119 0.404660654; y=0 0.192347906 0.219628094 0.23819892 0.251976191 0.252181432 0.264979117; z=0 0.044174154 0.058131509 0.066640197 0.073187605 0.084552378 0.093852475; l=0 0.096131452 0.111689894 0.123715906 0.141877613 0.15143182 0.167190503; m=0 0.569090571 0.670960794 0.742627034 0.823984239 0.875581648 0.954255576; n=4.591368228 4.022277657 3.920407434 3.848741194 3.767383988 3.71578658 3.637112651; 就是已知x,y,z,l,m,n 隨時(shí)間變化的值，想擬合出模型中的參數(shù) a,b,c,d,e,f,g(這些參數(shù)的值都是大于等于0） 之前有下過低版本的1stopt，結(jié)果點(diǎn)運(yùn)行沒反應(yīng)，不知道是不是版本問題還是代碼寫的有問題。。 求各路牛人幫忙用高版本運(yùn)行下，非常非常感謝~~ |
回復(fù)此樓» 本帖@通知
» 猜你喜歡
2026年循環(huán)經(jīng)濟(jì)功能材料國際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有207人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(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)有26人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有34人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求大神幫忙擬合一個(gè)非線性方程，求出模型參數(shù)
已經(jīng)有15人回復(fù)
- 求1stopt擬合動(dòng)力學(xué)參數(shù)
已經(jīng)有7人回復(fù)
- 【求專家】根據(jù)一組數(shù)據(jù)，求其方程，尋最優(yōu)擬合方程。要求截距為0
已經(jīng)有3人回復(fù)
- 求助用matlab擬合動(dòng)力學(xué)方程
已經(jīng)有9人回復(fù)
- matlab 擬合方法
已經(jīng)有4人回復(fù)
- 求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù)
已經(jīng)有13人回復(fù)
- 通過實(shí)驗(yàn)數(shù)據(jù)擬合，求解公式中的參數(shù)
已經(jīng)有14人回復(fù)
- 微分方程組參數(shù)擬合的問題（多參數(shù)動(dòng)力系統(tǒng)參數(shù)估計(jì)）
已經(jīng)有11人回復(fù)
- 重金請(qǐng)教-如何用MAtlab對(duì)力學(xué)模型進(jìn)行參數(shù)識(shí)別
已經(jīng)有12人回復(fù)
- MATLAB微分方程參數(shù)擬合問題，求大神
已經(jīng)有7人回復(fù)
- 大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒有遇到過這種情況？
已經(jīng)有8人回復(fù)
- 用MATLAB怎么打開點(diǎn)云并把它擬合到一個(gè)平面上？
已經(jīng)有8人回復(fù)
- matlab數(shù)據(jù)擬合問題求助
已經(jīng)有7人回復(fù)
- 動(dòng)力學(xué)參數(shù)擬合
已經(jīng)有26人回復(fù)
- matlab 擬合反應(yīng)動(dòng)力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a>
已經(jīng)有14人回復(fù)
- matlab擬合擬均相動(dòng)力學(xué)參數(shù)
已經(jīng)有6人回復(fù)
- 如何用VTF方程擬合電導(dǎo)率與溫度曲線
已經(jīng)有8人回復(fù)
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復(fù)
- 【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合
已經(jīng)有17人回復(fù)
- 【求助】擬合動(dòng)力學(xué)方程求助
已經(jīng)有13人回復(fù)
- 【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場參數(shù)的程序？【已解決】
已經(jīng)有15人回復(fù)
- 【求助】遞推公式的參數(shù)擬合
已經(jīng)有5人回復(fù)
終極愛好者
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 132.9
- 帖子: 32
- 在線: 13小時(shí)
- 蟲號(hào): 2999962
- 注冊: 2014-02-27
- 性別: GG
2樓2014-04-01 09:45:56
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
如果都要求大于0： 均方差(RMSE):0.291509496670823 殘差平方和(SSE):3.05920031937395 相關(guān)系數(shù)(R): 0.41204372086458 相關(guān)系數(shù)之平方(R^2): 0.169780027903928 確定系數(shù)(DC): -32.2049580782807 F統(tǒng)計(jì)(F-Statistic): -0.0103799073026604 參數(shù) 最佳估算 -------------------- ------------- a 2379421492.52982 b 5477871570.00215 c 4827348354.29146 d 0.000157698484986776 e 28.2125368960529 f 9.53786176032263E-8 g 333524832.481022 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z 目標(biāo) l 計(jì)算 l 目標(biāo) m 計(jì)算 m 目標(biāo) n 計(jì)算 n 1 5 0.598029237 0.461228912051106 0.192347906 0.212216074803243 0.044174154 0.154698496352178 0.096131452 0.102539265793472 0.569090571 0.18430565779089 4.022277657 4.40706257020911 2 10 0.541233252 0.352387819112766 0.219628094 0.194209366200932 0.058131509 0.186948110939243 0.111689894 0.19713745274706 0.670960794 0.259898922792036 3.920407434 4.33146930520796 3 20 0.502127725 0.241106037906148 0.23819892 0.153262848683966 0.066640197 0.191130610932665 0.123715906 0.34518325147722 0.742627034 0.355979862816147 3.848741194 4.23538836518385 4 30 0.46364134 0.179231653280274 0.251976191 0.121526115166401 0.073187605 0.174250824742148 0.141877613 0.455674155811178 0.823984239 0.419251043228081 3.767383988 4.17211718477191 5 40 0.442517119 0.138636909389627 0.252181432 0.0975440616385669 0.084552378 0.153342396668236 0.15143182 0.541159381303569 0.875581648 0.465201935622196 3.71578658 4.1261662923778 6 60 0.404660654 0.0883552573120677 0.264979117 0.0646893622260779 0.093852475 0.114305811509068 0.167190503 0.663332317952786 0.954255576 0.527617939775541 3.637112651 4.06375028822446 |
3樓2014-04-03 14:01:19
|
真是太感謝了。不過這個(gè)擬合效果確實(shí)挺差的，不知道是不是模型不太合適。。 你能幫我再試試這組數(shù)據(jù)么？謝謝~ t=0 2 5 10 15 20 30 x=0.97929461 0.848102002 0.791547773 0.742120992 0.665706681 0.627442016 0.567283612 y=0 0.111151416 0.154799598 0.178808192 0.203933465 0.230873341 0.255998614 z=0 0.029197641 0.031321393 0.058558032 0.078165633 0.1017788 0.122183142 l=0 0.000965082 0.003370178 0.008463798 0.012727551 0.019200453 0.030677511 m=0 0.142077349 0.222319921 0.295346436 0.454731222 0.492032299 0.601852621 n=1.929382737 1.787305387 1.707062816 1.634036301 1.474651515 1.437350438 1.327530115 |
4樓2014-04-03 15:01:04
5樓2015-12-07 10:55:40
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 chen1215 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +3 | 簡單點(diǎn)0 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 7/350 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[基金申請(qǐng)]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|
|
|
[考博] 2026年博士申請(qǐng) +3 | QwQwQW10 2026-03-11 | 3/150 |
|