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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

allengjn

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3772.5
紅花: 1
帖子: 487
在線: 28.3小時
蟲號: 3367570
注冊: 2014-08-15
專業(yè): 機(jī)械動力學(xué)

引用回帖:

10樓: Originally posted by 月只藍(lán) at 2014-11-06 21:39:29
代碼參見9樓給出的鏈接中第4樓 dingd 的回復(fù)。...

這個一階導(dǎo)可以用x'，如果是二階導(dǎo)，比如我標(biāo)題里的式子，那怎么寫？是x后邊跟兩個單引號，還是x后邊一個雙引號？
寫成x''=（Q-C*x'-K*x）/M 這種形式嗎？
還有，使用1stopt擬合，按那個鍵開始擬合？是”執(zhí)行“、”快速公式擬合搜索“、”全執(zhí)行“、”熱執(zhí)行“這些按鍵的哪個？需要先保存在運(yùn)行嗎？

贊一下

回復(fù)此樓

11樓2014-11-06 22:08:08

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

dingd

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

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

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
allengjn: 金幣+10, ★★★★★最佳答案 2014-11-06 22:27:33

1stOpt代碼：

CODE:

Constant M=0.3986,Q=100;

Variable t,x,x',x'';

ODEFunction x''= (Q-C*x'-K*x)/M;

Data;

//時間        位移x        速度dx        加速度ddx

0        0        0        26663.2

0.000500005        0.000137226        0.382478        281.033

0.00100001        0.000360148        0.505711        229.476

0.00150001        0.000640082        0.612797        212.787

0.00200002        0.000965359        0.664756        -9.53666

0.00250002        0.00130017        0.676495        -53.0477

0.00300003        0.0016186        0.587299        -131.129

....

均方差(RMSE):57.5493822281378
殘差平方和(SSE):1987158.83690418
相關(guān)系數(shù)(R): 0.291877273942345
相關(guān)系數(shù)之平方(R^2): 0.0851923430440146
確定系數(shù)(DC): -2.20653576106198
F統(tǒng)計(F-Statistic): 25.27741280379

參數(shù) 最佳估算
-------------------- -------------
c 49.6465860963742
k 70046293.1591487

其實(shí)如果不考慮加速度ddx的話，效果似乎更好：

CODE:

Constant M=0.3986,Q=100;

Variable t,x,x';

ODEFunction x''= (Q-C*x'-K*x)/M;

Data;

//時間        位移x        速度dx

0        0        0

0.000500005        0.000137226        0.382478

0.00100001        0.000360148        0.505711

0.00150001        0.000640082        0.612797

0.00200002        0.000965359        0.664756

0.00250002        0.00130017        0.676495

0.00300003        0.0016186        0.587299

0.00350003        0.00191261        0.610209

.....

均方差(RMSE):0.134482150106052
殘差平方和(SSE):7.23417947885864
相關(guān)系數(shù)(R): 0.905464687918125
相關(guān)系數(shù)之平方(R^2): 0.819866301066668
確定系數(shù)(DC): -0.231632439076251
F統(tǒng)計(F-Statistic): 90.6285375843916

參數(shù) 最佳估算
-------------------- -------------
c 24.0628939603168
k 11361.0933692284

c122.jpg

c123.jpg

贊一下

回復(fù)此樓

12樓2014-11-06 22:21:38

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

allengjn

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3772.5
紅花: 1
帖子: 487
在線: 28.3小時
蟲號: 3367570
注冊: 2014-08-15
專業(yè): 機(jī)械動力學(xué)

引用回帖:

12樓: Originally posted by dingd at 2014-11-06 22:21:38
1stOpt代碼：

Constant M=0.3986,Q=100;
Variable t,x,x',x'';
ODEFunction x''= (Q-C*x'-K*x)/M;
Data;
//時間位移x 速度dx 加速度ddx
0 0 0 26663.2
0.000500005 0.000137226 0.382478 281.033
0.0010 ...

哦，非常感謝，還有最后一個小問題，將代碼輸入到1stopt之后怎么運(yùn)行才是對它的函數(shù)擬合，我試了一下午總是不靈啊，您的1stopt是哪個版本的，是版本的問題嗎？我輸入的代碼，之后點(diǎn)”執(zhí)行“，就全都是0；點(diǎn)”快速公式擬合搜索“，就是擬合出來一個sin函數(shù)。這個怎么搞？
Parameter K,C;
Variable t, x;
ODEFunction x''= (100-C*x'-K-x)/0.39;
Data;
0 0
0.000500005 0.000137226
0.00100001 0.000360148
0.00150001 0.000640082
0.00200002 0.000965359
0.00250002 0.00130017
0.00300003 0.0016186
0.00350003 0.00191261
0.00400004 0.00220283
0.00450004 0.0024554
0.00500005 0.00266415
0.00550005 0.00288092
0.00600005 0.00314881
0.00650006 0.00344612

贊一下

回復(fù)此樓

13樓2014-11-06 22:27:15

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

allengjn

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3772.5
紅花: 1
帖子: 487
在線: 28.3小時
蟲號: 3367570
注冊: 2014-08-15
專業(yè): 機(jī)械動力學(xué)

送紅花一朵

懸賞分給了樓上的那位了，不好意思就只能送您一朵花了，非常感謝，希望以后還能得到您的幫助

贊一下

回復(fù)此樓

14樓2014-11-06 22:29:08

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

dingd

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

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

引用回帖:

13樓: Originally posted by allengjn at 2014-11-06 22:27:15
哦，非常感謝，還有最后一個小問題，將代碼輸入到1stopt之后怎么運(yùn)行才是對它的函數(shù)擬合，我試了一下午總是不靈啊，您的1stopt是哪個版本的，是版本的問題嗎？我輸入的代碼，之后點(diǎn)”執(zhí)行“，就全都是0；點(diǎn)”快速公 ...

如果用的是1.5破解版就別試了，沒法運(yùn)行的。

贊一下

回復(fù)此樓

15樓2014-11-06 23:02:33

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

allengjn

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3772.5
紅花: 1
帖子: 487
在線: 28.3小時
蟲號: 3367570
注冊: 2014-08-15
專業(yè): 機(jī)械動力學(xué)

內(nèi)容已刪除

回復(fù)此樓

16樓2014-11-06 23:04:18

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

allengjn

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3772.5
紅花: 1
帖子: 487
在線: 28.3小時
蟲號: 3367570
注冊: 2014-08-15
專業(yè): 機(jī)械動力學(xué)

內(nèi)容已刪除

回復(fù)此樓

17樓2014-11-06 23:34:01

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

dikeway

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 320.3
散金: 1
紅花: 1
帖子: 102
在線: 34.4小時
蟲號: 3414256
注冊: 2014-09-13
性別: GG
專業(yè): 經(jīng)濟(jì)地理學(xué)

引用回帖:

17樓: Originally posted by allengjn at 2014-11-06 23:34:01
您的意思是不是基本上網(wǎng)上能下到的，都不能擬合？您的是花錢買的嗎？...

都告訴你了，低版本的不能擬合微分方程，所以你輸入代碼結(jié)果全是0。大神是高版本的，據(jù)悉相當(dāng)之貴，1stopt6.0加密狗版1.5萬+

贊一下

回復(fù)此樓

18樓2014-11-07 09:04:12

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

allengjn

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 3772.5
紅花: 1
帖子: 487
在線: 28.3小時
蟲號: 3367570
注冊: 2014-08-15
專業(yè): 機(jī)械動力學(xué)

版主，我看了你回復(fù)的http://www.gaoyang168.com/bbs/viewthread.php?tid=6425538&authorid=1122189這個帖子，我現(xiàn)在有些混亂了，我的二階微分方程如何在KineticsEqs方程中寫出啊，因?yàn)槲业臄M合數(shù)據(jù)有三組，位移、速度、加速度，這三組是相關(guān)的，如何編寫這個程序，希望版主能夠指導(dǎo)一下

贊一下

回復(fù)此樓

19樓2014-11-07 10:44:24

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

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

返回列表

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

[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +3	材料學(xué)257求調(diào)劑 2026-03-23	4/200	2026-03-23 09:25 by 材料學(xué)257求調(diào)劑
[考研] 求調(diào)劑院校信息 +6	CX 330 2026-03-21	6/300	2026-03-22 15:25 by 無懈可擊111
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	17/850	2026-03-22 15:18 by 脫穎而出
[考研] 0856材料專碩353求調(diào)劑 +4	NIFFFfff 2026-03-20	4/200	2026-03-22 09:49 by 2026paper
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學(xué)員8dgXkO
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 22408 344分求調(diào)劑一志愿華電計算機(jī)技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678