 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 774 | 回復(fù): 2 | ||
iloveyou130金蟲 (小有名氣)
|
[求助]
求高版本的1stopt跑一下以下代碼，謝謝！ 已有1人參與
|
|
第一次發(fā)帖寫程序，也不知道有沒有寫對，懸賞金幣也不知是否合適。在以下程序中，變量y8沒有在微分方程中出現(xiàn)過，由于y12~y14不可測量，所以沒有相關(guān)數(shù)據(jù)，不知道這樣可以擬合k的值嗎？ --------------------------------------------------------------------------- Parameter k(17); k(17)>=0; Variable t,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11; ODEFunction y1'=-2*k1*y2*y1-2*k3*y2*y1-8*k16*y2*y1-k17*y7*y1+k17/2.2*y6*y9; y2'=-k1*y2*y1-k3*y2*y1-k5*y2*y6-k7*y2-k16*y2*y1; y3'=k2*y5*y6+k4*y4*y6+k4*y4*y6+k5*y2*y6-k6*y3; y4'=k3*y2*y6-k4*y4*y6; y5'=k1*y2-k2*y5*y6; y6'=-2*k1*y2*y1-k2*y5*y6-k4*y4*y6-4*k5*y2*y6+k9*y14+3*k15-y11+10*k16*y2*y1+k17*y7*y1-k17/2.2*y6*y9; y7'=2*k1*y2*y1+k10*y14+8*k16*y2*y1-k17*y7*y1+k17/2.2*y6*y9; y8'=k4*y4*y6+k11*y14; y9'=k3*y2*y1+k12*y14+k17*y7*y1-k17/2.2*y6*y9; y10'=k5*y2*y6+k13*y14; y11'=k2*y5*y6+k3*y2*y1+k5*y2*y6+k14*y14-2*k15*y11; y12'=k15*y11; y13'=k8*y14; y14'=k6*y3+k7*y2-k8*y14; data; 10 5.41 0.19 9.56E-05 0.000966 0.000549 0.00165 0.00119 0.00134 0.00162 0 0.00169 20 5.45 0.17 2.62E-05 0.0023 0.00112 0.00302 0.00221 0.00174 0.00338 0.000283 0.00452 30 5.46 0.15 0.000786 0.00427 2.90E-03 0.00403 0.00317 0.00209 0.000633 0.00071 0.00527 45 5.49 0.12 0.00109 0.00461 0.0026 0.00998 0.00599 0.00644 0.00282 0.000824 0.0052 60 5.51 0.11 0.00208 0.00633 0.00434 0.0211 0.0053 0.0184 0.00841 0.00221 0.0078 80 5.44 0.1 0.00361 0.00783 0.00651 0.0136 0.00622 0.0103 0.00308 0.00153 0.00844 |
回復(fù)此樓» 猜你喜歡
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有168人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有3人回復(fù)
-
計算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有4人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
獨孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動力學(xué)
- 管轄: 計算模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
iloveyou130: 金幣+20, ★★★很有幫助 2020-06-18 11:17:03
感謝參與，應(yīng)助指數(shù) +1
iloveyou130: 金幣+20, ★★★很有幫助 2020-06-18 11:17:03
|
均方差(RMSE): 0.0199436505934012 殘差平方和(SSR): 0.021876205944542 相關(guān)系數(shù)(R): 0.406206827247245 相關(guān)系數(shù)之平方(R^2): 0.165003986502274 修正R平方(Adj. R^2): 1.15044908256948 確定系數(shù)(DC): -0.68967324594123 F統(tǒng)計(F-Statistic): -0.66404334401412 參數(shù) 最佳估算 -------------------- ------------- k1 4.60883336079867E-8 k2 1.54247889887596 k3 5.1150136876246E-16 k4 0.000387115599728942 k5 3.30143387048788 k6 9.90173114206637E-12 k7 1.41860538976745 k8 19771.887540242 k9 0.341602464209552 k10 1125.69386161916 k11 2369.80489888995 k12 693.150581014912 k13 33.3865189378566 k14 1095.01407663084 k15 0.0176173700760016 k16 4.28454221172775E-13 k17 5.96905321164476E-13 y12 初值 0 y13 初值 0 y14 初值 0.000805185581681898 ====== 結(jié)果輸出 ====== 文件: 數(shù)據(jù)文件-1 No t 目標(biāo) y1 計算 y1 目標(biāo) y2 計算 y2 目標(biāo) y3 計算 y3 目標(biāo) y4 計算 y4 目標(biāo) y5 計算 y5 目標(biāo) y6 計算 y6 目標(biāo) y7 計算 y7 目標(biāo) y8 計算 y8 目標(biāo) y9 計算 y9 目標(biāo) y10 計算 y10 目標(biāo) y11 計算 y11 1 20 5.45 5.40999999188834 0.17 0.16673223415736 2.62E-5 0.000288890181847264 0.0023 0.000965999875847417 0.00112 0.00054871967692699 0.00302 0.0054791674073478 0.00221 0.00254889247969801 0.00174 0.00420071647181485 0.00338 0.00245673858528659 0.000283 0.000233311486798891 0.00452 0.0031972088649323 2 30 5.46 5.40999998463729 0.15 0.14576533273367 0.000786 0.000629754879635674 0.00427 0.00096599962755348 0.0029 0.000548158657804571 0.00403 0.00867615151602474 0.00317 0.00372337399212577 0.00209 0.00667321386445097 0.000633 0.00317992589617602 0.00071 0.000608447276701977 0.00527 0.00466761215373619 3 45 5.49 5.40999997830098 0.12 0.127319431959793 0.00109 0.00105189712260829 0.00461 0.000965999276436763 0.0026 0.000547365929143034 0.00998 0.0114145598811571 0.00599 0.00474963372253263 0.00644 0.00883367755620465 0.00282 0.00381184560421971 0.000824 0.00106023275575494 0.0052 0.00607038597326002 4 60 5.51 5.40999997276868 0.11 0.111119483171172 0.00208 0.00151543628490384 0.00633 0.000965998835359388 0.00434 0.000546371494207135 0.0211 0.0138612901082428 0.0053 0.0056457324530104 0.0184 0.0107201287662758 0.00841 0.00436361864140109 0.00221 0.00154935297132375 0.0078 0.00738348099987025 5 80 5.44 5.40999996794202 0.1 0.0969107562795792 0.00361 0.00199525344001146 0.00783 0.000965998312451843 0.00651 0.000545194749301669 0.0136 0.016125635375396 0.00622 0.00642747209680412 0.0103 0.012365833385529 0.00308 0.00484497503960053 0.00153 0.00205117706028891 0.00844 0.00859746204094071 |
2樓2020-06-16 21:59:52
iloveyou130
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1189.4
- 散金: 20
- 帖子: 80
- 在線: 18.8小時
- 蟲號: 2835371
- 注冊: 2013-11-29
- 專業(yè): 多相流熱物理學(xué)
3樓2020-06-18 11:09:00
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 iloveyou130 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 | +5 | 橙子cyx 2026-03-06 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風(fēng)了呢 2026-03-05 | 9/450 |
|
|
[考研] 一志愿中石油，化學(xué)工程，301求調(diào)劑 +4 | 豚饅枳枳 2026-03-04 | 4/200 |
|
|
|
[考研] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +3 | ipys00zm58 2026-03-06 | 3/150 |
|
|
|
[考研] 070300調(diào)劑一志愿南京大學(xué)化學(xué)專業(yè)，初試273分，有論文有專利有競賽，求調(diào)劑B區(qū)211 +3 | 82206202 2026-02-28 | 13/650 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工，前三科206，總分283，求調(diào)劑 +6 | 26考研求調(diào)劑 2026-03-06 | 6/300 |
|
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線，總分 281 求調(diào)劑 +6 | 香香 Lu 2026-03-03 | 11/550 |
|
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
 5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +5 | ccc! 2026-03-03 | 6/300 |
|
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
|
[考研] 0855機(jī)械工程324求調(diào)劑 +3 | 鄭SK 2026-03-02 | 7/350 |
|
|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+6
|
15779376950 2026-03-01 | 8/400 |
|
|
|
[考研] 能動297求調(diào)劑，本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|