| 5 | 1/1 | 返回列表 |
| 查看: 2995 | 回復: 16 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
1stopt微分方程組參數(shù)擬合 已有2人參與
|
|||
|
真心希望有人能幫我一下,這個問題已經(jīng)困擾我很多天了,這10個金幣是我所有的金幣了。 具體的內容是這樣 然后再得 c是光速為3*10^10,k=0.5,M、tsd和L是未知參數(shù),需要擬合求解,現(xiàn)在只知道能得到t和LL的實驗值。具體的計算應該是v和Eint都是不知道的,若知道M、tsd和L的值,我用的是ode45,參數(shù)v和Eint根據(jù)t來回算(初始值是0,0),然后帶入到LL里面。現(xiàn)在反過來了不知道M、tsd和L的值我就不知道怎么寫了,用MATLAB的ode45和lsqnonlin寫了個程序,最后求出來不是很理想,而且跟初始值那些影響蠻大。用1stopt寫了一個,運行不了(我覺得肯定是錯的),不知道怎么寫,能幫我改改,然后算出來嗎?真的麻煩了,在線等啊。。。 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 |
|
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è)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專碩調劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 290求調劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 260求調劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 289求調劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
[考研] 求調劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 08工學調劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 286分人工智能專業(yè)請求調劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 化學求調劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 304求調劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 301求調劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 275求調劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 302求調劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 0854控制工程 359求調劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 304求調劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 070300化學學碩求調劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
[考研] 327求調劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 一志愿哈工大材料324分求調劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|