 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 777 | 回復(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ù)此樓» 猜你喜歡
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有224人回復(fù)
-
計算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 310 070300化學(xué)求調(diào)劑 +4 | 撲風(fēng)鈴的貓 2026-03-08 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 0703求調(diào)劑 +10 | c297914 2026-03-06 | 11/550 |
|
|
[考研] 327分求調(diào)劑086 +3 | 西紅柿?小帥 2026-03-09 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 鶴遨予卿 2026-03-09 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 9/450 |
|
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +9 | 1一11 2026-03-03 | 9/450 |
|
|
|
[考研] 0856求調(diào)劑 +3 | squirtle11 2026-03-07 | 3/150 |
|
|
|
[考研] 一志愿中科大能動297求調(diào)劑，本科川大 +4 | 邵11 2026-03-03 | 4/200 |
|
|
|
[考研] 083000，總分284，求調(diào)劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|
|
[基金申請] 進(jìn)入個人成果庫好難，一下午都沒進(jìn)去 +6 | mi_dilee 2026-03-05 | 6/300 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工，前三科206，總分283，求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研]
 材料工程調(diào)劑
+6
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
|
[考研] 材料與化工，291，求調(diào)劑 +7 | 咕嚕咕嚕123123 2026-03-05 | 8/400 |
|
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +8 | 拾玖壹 2026-03-04 | 8/400 |
|
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|