| 查看: 1030 | 回復: 14 | ||
heisekaishui銀蟲 (小有名氣)
|
[求助]
關于1STOPT擬合的小問題 已有2人參與
|
|
各位前輩,我最近用1stopt擬合了一組數(shù)據(jù),但是發(fā)現(xiàn)擬合的數(shù)據(jù)使用不同階的龍格庫塔方程結果不同,但是精度都差不多,不知道是什么原因,希望大家不吝賜教。 代碼如下: Parameters k1,k2,k3; Variable t,a,b,c; ODEFunction a'=-k1*a-k2*a; b'=k1*a; c'=k2*a; Data; t=0,2,4,6,8,10; a=1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893; b=0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415; c=0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787; 同時,我在設定參數(shù)的范圍時,并沒有起到效果,而且1stopt的結果縱坐標范圍是0-4而不是0-10不知道是什么原因 |

|
看看這樣能不能解決你的問題 Parameters k1,k2,k3; Variable t,a,b,c; SharedModel; ODEFunction a'=-k1*a-k2*a; b'=k1*a; c'=k2*a; RowData; 0,2,4,6,8,10; 1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893; 0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415; 0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787; |

|
Parameters k1,k2,k3,k4; 修改方程,這個可以得到穩(wěn)定結果 Variable t,a,b,c; ODEFunction a'=-k1*a-k2*a; b'=k3*a; c'=k4*a; RowData; 0,2,4,6,8,10; 1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893; 0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415; 0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787; |

|
坐標顯示0-4可能是他的作圖功能不好,或不是正版軟件,只看他的計算值就行了(計算值不是有,每時刻的計算值嗎) 2 1.098850041 .060153328 .002571727 1.08633992640023 6.97307002743706E-02 5.60446897778746E-03 For a: 4 1.022186084 .132977117 .006411895 1.01597773545576 .134944946131689 1.08459367441951E-02 DC: 0.99487234139559 6 .952491579 .197793007 .01129051 .950172900633618 .195935268051574 1.57479148657442E-02 R: 0.999959508427279 8 .885120223 .254614938 .021839935 .888630242171105 .25297524940825 2.03323920723293E-02 RMSE: 0.00692414160103415 10 .825879893 .306825415 .028869787 .831073699086233 .306320753592402 2.46199328847424E-02 SSE: 0.000239718684555859 |


鐵桿木蟲 (職業(yè)作家)
|
穩(wěn)定唯一: 均方差(RMSE):0.00523400858904149 殘差平方和(SSE):0.000410922688652401 相關系數(shù)(R): 0.99229618847388 相關系數(shù)之平方(R^2): 0.98465172565979 確定系數(shù)(DC): 0.949184328781669 F統(tǒng)計(F-Statistic): 585.146344863956 參數(shù) 最佳估算 -------------------- ------------- k1 0.0310189175554483 k2 0.0024794201096689 |
銀蟲 (小有名氣)

銀蟲 (小有名氣)

銀蟲 (小有名氣)

銀蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專碩調劑 +6 | Doleres 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 329求調劑 +6 | 想上學吖吖 2026-03-19 | 6/300 |
|
|
[考研] 307求調劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
[考研] 294求調劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 【考研調劑】化學專業(yè) 281分,一志愿四川大學,誠心求調劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研] 085600材料與化工調劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 274求調劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 08工科 320總分 求調劑 +5 | 梨花珞晚風 2026-03-17 | 5/250 |
|
|
[考研] 268求調劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
[考研] 290求調劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 297求調劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 330求調劑 +3 | ?醬給調劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 304求調劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 311求調劑 +3 | 冬十三 2026-03-13 | 3/150 |
|