| 5 | 1/1 | 返回列表 |
| 查看: 783 | 回復: 2 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
iloveyou130金蟲 (小有名氣)
|
[求助]
求高版本的1stopt跑一下以下代碼,謝謝! 已有1人參與
|
|
|
第一次發(fā)帖寫程序,也不知道有沒有寫對,懸賞金幣也不知是否合適。在以下程序中,變量y8沒有在微分方程中出現(xiàn)過,由于y12~y14不可測量,所以沒有相關數(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 |
金蟲 (小有名氣)
版主 (知名作家)
|
均方差(RMSE): 0.0199436505934012 殘差平方和(SSR): 0.021876205944542 相關系數(shù)(R): 0.406206827247245 相關系數(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 目標 y1 計算 y1 目標 y2 計算 y2 目標 y3 計算 y3 目標 y4 計算 y4 目標 y5 計算 y5 目標 y6 計算 y6 目標 y7 計算 y7 目標 y8 計算 y8 目標 y9 計算 y9 目標 y10 計算 y10 目標 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 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 紡織、生物、化學、材料相關專業(yè)招生了 +3 | 耶耶業(yè) 2026-03-09 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 化工學碩306求調(diào)劑 +7 | 42838695 2026-03-12 | 7/350 |
|
|
[考研] 085701環(huán)境工程專碩求調(diào)劑 +7 | 幼稚園大班打手 2026-03-08 | 7/350 |
|
|
[考研] 一志愿江南大學085701環(huán)境工程專碩總分287求調(diào)劑 +5 | 18266118446 2026-03-09 | 5/250 |
|
|
[考研] 293求調(diào)劑,一志愿陜師大生物學 +3 | ??????.?.??? 2026-03-09 | 3/150 |
|
|
[考研] 312求調(diào)劑 +5 | 陌宸希 2026-03-10 | 5/250 |
|
|
[考研] 420求調(diào)劑 +3 | 莫向外求11 2026-03-10 | 3/150 |
|
|
[考研] 0817學碩華東區(qū)求調(diào)劑 +3 | 30660438 2026-03-08 | 3/150 |
|
|
[考研] 328,0703考生求調(diào)劑,一志愿為東北師范大學 +3 | 觀素律 2026-03-09 | 4/200 |
|
|
[考博] 2026博士申請 +6 | 起泡酒 2026-03-08 | 6/300 |
|
|
[考研] 085600材料與化工 326 求調(diào)劑 +4 | 熱愛生活ing 2026-03-09 | 4/200 |
|
|
[考研] 085600材料與化工,一志愿廣州985,求調(diào)劑 +15 | qqyyaill 2026-03-05 | 15/750 |
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
[考研] 301求調(diào)劑 +11 | 朝天椒的雙馬尾 2026-03-05 | 11/550 |
|
|
[考研] 安徽農(nóng)業(yè)大學材料與化學學院0856材料招收調(diào)劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
[考研] 288求調(diào)劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
[考研]
|
程晴之 2026-03-06 | 6/300 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財 2026-03-05 | 5/250 |
|