版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

下一頁(yè)

heisekaishui

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 215
散金: 30
帖子: 52
在線: 59.2小時(shí)
蟲號(hào): 2445260
注冊(cè): 2013-05-02
性別: GG
專業(yè): 催化化學(xué)

[求助] 關(guān)于1STOPT擬合的小問(wèn)題已有2人參與

各位前輩，我最近用1stopt擬合了一組數(shù)據(jù)，但是發(fā)現(xiàn)擬合的數(shù)據(jù)使用不同階的龍格庫(kù)塔方程結(jié)果不同，但是精度都差不多，不知道是什么原因，希望大家不吝賜教。
代碼如下：
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í)，我在設(shè)定參數(shù)的范圍時(shí)，并沒有起到效果，而且1stopt的結(jié)果縱坐標(biāo)范圍是0-4而不是0-10不知道是什么原因

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有286人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有65人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

利用Excel數(shù)據(jù)前處理，結(jié)合1stopt數(shù)據(jù)自定義擬合初值，美美origin數(shù)據(jù)作圖分析！已經(jīng)有72人回復(fù)
關(guān)于用1stOpt進(jìn)行參數(shù)擬合問(wèn)題已經(jīng)有4人回復(fù)
求助！求助！小女子快崩潰了，關(guān)于1STOPT的擬合問(wèn)題，謝謝各位大神啦，求解答��！已經(jīng)有4人回復(fù)
個(gè)人發(fā)現(xiàn)軟件1stopt的隱函數(shù)擬合功能存在嚴(yán)重問(wèn)題已經(jīng)有8人回復(fù)
求助一個(gè)多元函數(shù)擬合的問(wèn)題。matlab或1stopt 已經(jīng)有2人回復(fù)
[求助] 1stopt 微分方程代碼求助已經(jīng)有5人回復(fù)
求助大神，關(guān)于1stopt帶積分的函數(shù)擬合已經(jīng)有7人回復(fù)
1stopt 非線性約束擬合已經(jīng)有6人回復(fù)
關(guān)于用debye公式擬合熱熔數(shù)據(jù) 求解debye溫度的問(wèn)題已經(jīng)有6人回復(fù)
關(guān)于1stopt擬合隱函數(shù)方程時(shí)如何設(shè)置目標(biāo)函數(shù)的問(wèn)題已經(jīng)有2人回復(fù)
關(guān)于1stOpt使用的問(wèn)題已經(jīng)有1人回復(fù)
1stopt軟件進(jìn)行非線性擬合的問(wèn)題已經(jīng)有1人回復(fù)
學(xué)1stopt要多久？已經(jīng)有12人回復(fù)
求助關(guān)于1stopt擬合~··~~ 已經(jīng)有4人回復(fù)

Don'tPanic!

1樓 2015-07-17 19:41:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ybkooo

至尊木蟲 (著名寫手)

應(yīng)助: 32 (小學(xué)生)
金幣: 11360
散金: 178
紅花: 15
帖子: 2310
在線: 543.5小時(shí)
蟲號(hào): 2710519
注冊(cè): 2013-10-09
性別: GG
專業(yè): 控制理論與方法

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
heisekaishui: 金幣+10, ★★★很有幫助, 非常感謝，我還想問(wèn)一下如果我要設(shè)置參數(shù)的取值范圍應(yīng)該怎么設(shè)置，比如我想讓所有的 2015-07-18 13:32:29
heisekaishui: 金幣+10, ★★★很有幫助, 非常感謝，我還想問(wèn)一下如果我要設(shè)置參數(shù)的取值范圍應(yīng)該怎么設(shè)置，比如我想讓所有的 2015-07-18 13:32:48

看看這樣能不能解決你的問(wèn)題
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;

贊一下

回復(fù)此樓

abcd

2樓2015-07-17 21:15:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ybkooo

至尊木蟲 (著名寫手)

應(yīng)助: 32 (小學(xué)生)
金幣: 11360
散金: 178
紅花: 15
帖子: 2310
在線: 543.5小時(shí)
蟲號(hào): 2710519
注冊(cè): 2013-10-09
性別: GG
專業(yè): 控制理論與方法

Parameters k1,k2,k3,k4;
修改方程，這個(gè)可以得到穩(wěn)定結(jié)果
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;

贊一下

回復(fù)此樓

abcd

3樓2015-07-17 21:41:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ybkooo

至尊木蟲 (著名寫手)

應(yīng)助: 32 (小學(xué)生)
金幣: 11360
散金: 178
紅花: 15
帖子: 2310
在線: 543.5小時(shí)
蟲號(hào): 2710519
注冊(cè): 2013-10-09
性別: GG
專業(yè): 控制理論與方法

坐標(biāo)顯示0-4可能是他的作圖功能不好，或不是正版軟件，只看他的計(jì)算值就行了（計(jì)算值不是有，每時(shí)刻的計(jì)算值嗎）
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

贊一下

回復(fù)此樓

abcd

4樓2015-07-17 21:48:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ybkooo

至尊木蟲 (著名寫手)

應(yīng)助: 32 (小學(xué)生)
金幣: 11360
散金: 178
紅花: 15
帖子: 2310
在線: 543.5小時(shí)
蟲號(hào): 2710519
注冊(cè): 2013-10-09
性別: GG
專業(yè): 控制理論與方法

1111.png

贊一下

回復(fù)此樓

abcd

5樓2015-07-17 21:51:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

穩(wěn)定唯一：

均方差(RMSE):0.00523400858904149
殘差平方和(SSE):0.000410922688652401
相關(guān)系數(shù)(R): 0.99229618847388
相關(guān)系數(shù)之平方(R^2): 0.98465172565979
確定系數(shù)(DC): 0.949184328781669
F統(tǒng)計(jì)(F-Statistic): 585.146344863956

參數(shù) 最佳估算
-------------------- -------------
k1 0.0310189175554483
k2 0.0024794201096689

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

heisekaishui

6樓2015-07-17 23:23:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

heisekaishui

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 215
散金: 30
帖子: 52
在線: 59.2小時(shí)
蟲號(hào): 2445260
注冊(cè): 2013-05-02
性別: GG
專業(yè): 催化化學(xué)

引用回帖:

2樓: Originally posted by ybkooo at 2015-07-17 21:15:36
看看這樣能不能解決你的問(wèn)題
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.95249 ...

我還想問(wèn)一下如果我要設(shè)置參數(shù)的取值范圍應(yīng)該怎么設(shè)置，比如我想讓所有的k大于0，應(yīng)該怎樣設(shè)置。

贊一下

回復(fù)此樓

Don'tPanic!

7樓2015-07-18 13:33:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

heisekaishui

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 215
散金: 30
帖子: 52
在線: 59.2小時(shí)
蟲號(hào): 2445260
注冊(cè): 2013-05-02
性別: GG
專業(yè): 催化化學(xué)

送紅花一朵

引用回帖:

6樓: Originally posted by dingd at 2015-07-17 23:23:35
穩(wěn)定唯一：

均方差(RMSE):0.00523400858904149
殘差平方和(SSE):0.000410922688652401
相關(guān)系數(shù)(R): 0.99229618847388
相關(guān)系數(shù)之平方(R^2): 0.98465172565979
確定系數(shù)(DC): 0.949184328781669
F統(tǒng)計(jì)(F-St ...

非常感謝

回復(fù)此樓

Don'tPanic!

8樓2015-07-18 13:37:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

heisekaishui

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 215
散金: 30
帖子: 52
在線: 59.2小時(shí)
蟲號(hào): 2445260
注冊(cè): 2013-05-02
性別: GG
專業(yè): 催化化學(xué)

引用回帖:

3樓: Originally posted by ybkooo at 2015-07-17 21:41:22
Parameters k1,k2,k3,k4;
修改方程，這個(gè)可以得到穩(wěn)定結(jié)果
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. ...

我發(fā)現(xiàn)這樣也是不穩(wěn)定，采用龍格庫(kù)塔一階擬合的和版主的就一致，但是階數(shù)越高差別越大，可能跟版本有關(guān)吧。

贊一下

回復(fù)此樓

Don'tPanic!

9樓2015-07-18 13:39:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

heisekaishui

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 215
散金: 30
帖子: 52
在線: 59.2小時(shí)
蟲號(hào): 2445260
注冊(cè): 2013-05-02
性別: GG
專業(yè): 催化化學(xué)

引用回帖:

我采用龍格庫(kù)塔一階擬合的和版主的就一致，但是階數(shù)越高差別越大，四階結(jié)果就是k1=1.3，K2=0.1了，不知道是不是跟版本有關(guān)。

贊一下

回復(fù)此樓

Don'tPanic!

10樓2015-07-18 13:41:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 heisekaishui 的主題更新

返回列表

下一頁(yè)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 297求調(diào)劑 +9	戲精丹丹丹 2026-03-17	9/450	2026-03-21 01:49 by JourneyLucky
[考研] 299求調(diào)劑 +4	shxchem 2026-03-20	5/250	2026-03-21 00:24 by 星空星月
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 321求調(diào)劑 +9	何潤(rùn)采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 求調(diào)劑 +5	Mqqqqqq 2026-03-19	5/250	2026-03-20 20:46 by zhukairuo
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明

24小時(shí)熱門版塊排行榜

[求助] 關(guān)于1STOPT擬合的小問(wèn)題 已有2人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

[求助] 關(guān)于1STOPT擬合的小問(wèn)題已有2人參與