quuq432

新蟲 (初入文壇)

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

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

真心希望有人能幫我一下，這個(gè)問題已經(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來回算（初始值是0,0），然后帶入到LL里面�，F(xiàn)在反過來了不知道M、tsd和L的值我就不知道怎么寫了，用MATLAB的ode45和lsqnonlin寫了個(gè)程序，最后求出來不是很理想，而且跟初始值那些影響蠻大。用1stopt寫了一個(gè)，運(yùn)行不了（我覺得肯定是錯(cuò)的），不知道怎么寫，能幫我改改，然后算出來嗎？真的麻煩了，在線等啊。。。
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ù)此樓

» 猜你喜歡

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

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

求助�。�！1stopt進(jìn)行常微分方程組參數(shù)模擬已經(jīng)有5人回復(fù)
1stOpt軟件微分方程參數(shù)擬合，運(yùn)行后沒有任何反應(yīng)。已經(jīng)有5人回復(fù)
1stopt 微分方程數(shù)據(jù)擬合求助已經(jīng)有2人回復(fù)
1stopt微分方程組的參數(shù)估計(jì)的請教已經(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ù)來求方程中的參數(shù) 已經(jīng)有9人回復(fù)
求助1stopt或者matlab 微分方程擬合已經(jīng)有0人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有1人回復(fù)
求助，matlab/ 1stOpt 自定義微分方程擬合實(shí)驗(yàn)數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
幫忙看看能否用1stopt計(jì)算出微分方程組已經(jīng)有0人回復(fù)
1stOpt求解常微分方程邊值問題已經(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í)
蟲號: 291104
注冊: 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í)
蟲號: 4110317
注冊: 2015-09-29
專業(yè): 恒星與星際物質(zhì)

沒有人愿意幫幫忙嗎

贊一下

回復(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í)
蟲號: 291104
注冊: 2006-10-28

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

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

微分方程項(xiàng)v'和E'的初值應(yīng)該知道吧？還有數(shù)據(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í)
蟲號: 4110317
注冊: 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ù)的量級有點(diǎn)離譜吧。

初值都是0，量級是超新星爆發(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í)
蟲號: 291104
注冊: 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í)
蟲號: 4110317
注冊: 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ì)算

能幫我試試看嗎，我試過了，不知道是盜版的原因還是什么運(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í)
蟲號: 4110317
注冊: 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í)
蟲號: 4110317
注冊: 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í)
蟲號: 4110317
注冊: 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ā)表

[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 321求調(diào)劑 +8	何潤采123 2026-03-18	10/500	2026-03-19 16:46 by 何潤采123
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 274求調(diào)劑 +6	S.H1 2026-03-18	6/300	2026-03-19 09:34 by 花店相見
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu