quuq432

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 13
在線: 9.5小時(shí)
蟲號(hào): 4110317
注冊(cè): 2015-09-29
專業(yè): 恒星與星際物質(zhì)

[求助] 1stopt微分方程組參數(shù)擬合已有2人參與

真心希望有人能幫我一下，這個(gè)問(wèn)題已經(jīng)困擾我很多天了，這10個(gè)金幣是我所有的金幣了。
具體的內(nèi)容是這樣 $dv/dt=Eint/(M*v*t)$
$dEint/t=(L/(1+t/tsd)^2)-4*pi*c*Eint*t*v/(3*k*M)-Eint/t$
然后再得
$LL=Eint*t*(4*pi*c*v)/(3*k*M)$
c是光速為3*10^10，k=0.5，M、tsd和L是未知參數(shù)，需要擬合求解，現(xiàn)在只知道能得到t和LL的實(shí)驗(yàn)值。具體的計(jì)算應(yīng)該是v和Eint都是不知道的，若知道M、tsd和L的值，我用的是ode45，參數(shù)v和Eint根據(jù)t來(lái)回算（初始值是0,0），然后帶入到LL里面�，F(xiàn)在反過(guò)來(lái)了不知道M、tsd和L的值我就不知道怎么寫了，用MATLAB的ode45和lsqnonlin寫了個(gè)程序，最后求出來(lái)不是很理想，而且跟初始值那些影響蠻大。用1stopt寫了一個(gè)，運(yùn)行不了（我覺(jué)得肯定是錯(cuò)的），不知道怎么寫，能幫我改改，然后算出來(lái)嗎？真的麻煩了，在線等啊。。。
Parameters M,tsd,L;
Constant c=3*10^10,k=0.5;
Variable t,v,E,LL;
Function v'=E/(M*v*t);
         E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t;
         LL=E*t*(4*pi*c*v)/(3*k*M);
Data;
//t          LL
502272       5.65432*10^42
762325       1.30331*10^43
1114354       9.00958*10^42
1364122       5.3601*10^42
1706100       3.56514*10^42
1964980       2.5611*10^42
2140916       2.27497*10^42
2313059       1.94825*10^42
2489034       1.50971*10^42
2742752       1.22059*10^42
3001044       1.04568*10^42

回復(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)有259人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(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ì)您同樣有幫助:

求助�。。�1stopt進(jìn)行常微分方程組參數(shù)模擬已經(jīng)有5人回復(fù)
1stOpt軟件微分方程參數(shù)擬合，運(yùn)行后沒(méi)有任何反應(yīng)。已經(jīng)有5人回復(fù)
1stopt 微分方程數(shù)據(jù)擬合求助已經(jīng)有2人回復(fù)
1stopt微分方程組的參數(shù)估計(jì)的請(qǐng)教已經(jīng)有3人回復(fù)
用1stOpt求解微分方程已經(jīng)有2人回復(fù)
[求助] 1stopt 微分方程代碼求助已經(jīng)有5人回復(fù)
求幫忙運(yùn)行一下下面的1stopt的程序，微分方程-參數(shù)估計(jì) 已經(jīng)有6人回復(fù)
用1stopt軟件做微分方程組的參數(shù)估計(jì) 已經(jīng)有13人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來(lái)求方程中的參數(shù) 已經(jīng)有9人回復(fù)
求助1stopt或者matlab 微分方程擬合已經(jīng)有0人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來(lái)求方程中的參數(shù) 已經(jīng)有1人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來(lái)求方程中的參數(shù) 已經(jīng)有4人回復(fù)
幫忙看看能否用1stopt計(jì)算出微分方程組已經(jīng)有0人回復(fù)
1stOpt求解常微分方程邊值問(wèn)題已經(jīng)有5人回復(fù)

1樓 2016-03-03 17:11:35

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

回帖支持 ( 顯示支持度最高的前 50 名 )

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ùn)行只是效果不好。

贊一下(1人)

回復(fù)此樓

10樓2016-03-04 17:26:28

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

普通回帖

quuq432

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 13
在線: 9.5小時(shí)
蟲號(hào): 4110317
注冊(cè): 2015-09-29
專業(yè): 恒星與星際物質(zhì)

沒(méi)有人愿意幫幫忙嗎

贊一下

回復(fù)此樓

2樓2016-03-03 19:41:21

已閱回復(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

微分方程項(xiàng)v'和E'的初值應(yīng)該知道吧？還有數(shù)據(jù)的量級(jí)有點(diǎn)離譜吧。

贊一下

回復(fù)此樓

3樓2016-03-03 21:12:00

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

quuq432

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 13
在線: 9.5小時(shí)
蟲號(hào): 4110317
注冊(cè): 2015-09-29
專業(yè): 恒星與星際物質(zhì)

引用回帖:

3樓: Originally posted by dingd at 2016-03-03 21:12:00
微分方程項(xiàng)v'和E'的初值應(yīng)該知道吧？還有數(shù)據(jù)的量級(jí)有點(diǎn)離譜吧。

初值都是0，量級(jí)是超新星爆發(fā)的理論計(jì)算值，這個(gè)大小算是一般的，能幫我看看是怎么回事嗎？我加了初值都還是運(yùn)行不了。

贊一下

回復(fù)此樓

4樓2016-03-03 21:50:59

已閱回復(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)助回帖

建議將一般公式：
LL=E*t*(4*pi*c*v)/(3*k*M);

寫成微分方程形式：LL'=。。。。

這樣才能進(jìn)行微分方程計(jì)算

贊一下

回復(fù)此樓

5樓2016-03-04 11:37:12

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

quuq432

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 13
在線: 9.5小時(shí)
蟲號(hào): 4110317
注冊(cè): 2015-09-29
專業(yè): 恒星與星際物質(zhì)

引用回帖:

5樓: Originally posted by dingd at 2016-03-04 11:37:12
建議將一般公式：
LL=E*t*(4*pi*c*v)/(3*k*M);

寫成微分方程形式：LL'=。。。。

這樣才能進(jìn)行微分方程計(jì)算

能幫我試試看嗎，我試過(guò)了，不知道是盜版的原因還是什么運(yùn)行不了

贊一下

回復(fù)此樓

6樓2016-03-04 13:04:23

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

quuq432

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 13
在線: 9.5小時(shí)
蟲號(hào): 4110317
注冊(cè): 2015-09-29
專業(yè): 恒星與星際物質(zhì)

引用回帖:

Parameters M,tsd,L;
Constant c=3*10^10,k=0.5;
InitialODEValue t=0.1,E=0.1,v=0.1,LL=0;
Variable t,LL；
ODEFunction v'=E/(M*v*t);
         E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t;
         LL'=E*(4*pi*c*v)/(3*k*M)
Data;
//t          LL
0.1          0.1
502272       1.12575E+37
762325          1.70965E+37
1114354          8.08503E+36
1364122          3.92934E+36
1706100          2.08964E+36
1964980       1.30337E+36
2140916          1.06261E+36
2313059          8.42281E+35
2489034          6.06545E+35
2742752          4.45025E+35
3001044          3.4844E+35

改成這樣了，還是不能運(yùn)行

贊一下

回復(fù)此樓

7樓2016-03-04 16:15:27

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

quuq432

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 13
在線: 9.5小時(shí)
蟲號(hào): 4110317
注冊(cè): 2015-09-29
專業(yè): 恒星與星際物質(zhì)

數(shù)據(jù)還是上面那個(gè)，我搞錯(cuò)了

贊一下

回復(fù)此樓

8樓2016-03-04 16:17:08

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

quuq432

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 11
帖子: 13
在線: 9.5小時(shí)
蟲號(hào): 4110317
注冊(cè): 2015-09-29
專業(yè): 恒星與星際物質(zhì)

Parameters M,tsd,L;
Constant c=3*10^10,k=0.5;
InitialODEValue E=0.1,v=0.1,LL=0.1;
Variable t,LL；
ODEFunction v'=E/(M*v*t);
         E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t;
         LL'=E*(4*pi*c*v)/(3*k*M)+E'*t*(4*pi*c*v)/(3*k*M)+E*t*(4*pi*c*v')/(3*k*M);
Data;
0.1          0.1
502272       5.65432*10^42
762325       1.30331*10^43
1114354       9.00958*10^42
1364122       5.3601*10^42
1706100       3.56514*10^42
1964980       2.5611*10^42
2140916       2.27497*10^42
2313059       1.94825*10^42
2489034       1.50971*10^42
2742752       1.22059*10^42
3001044       1.04568*10^42

應(yīng)該是這樣，運(yùn)行不了。。。。

贊一下

回復(fù)此樓

9樓2016-03-04 16:48:12

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

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

[考研] 085600材料與化工 +6	安全上岸！ 2026-03-16	6/300	2026-03-20 15:07 by lbsjt
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +4	Charlieyq 2026-03-19	4/200	2026-03-20 14:26 by 無(wú)懈可擊111
[考研] 295復(fù)試調(diào)劑 +7	簡(jiǎn)木ChuFront 2026-03-19	7/350	2026-03-20 14:11 by 星空星月
[考研] 271材料工程求調(diào)劑 +7	.6lL 2026-03-18	7/350	2026-03-20 09:10 by xingguangj
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫郵件了，沒(méi)回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3	一瓶苯甲酸 2026-03-14	3/150	2026-03-19 15:17 by 盡舜堯1
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky