| 5 | 1/1 | 返回列表 |
| 查看: 822 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
莫selina新蟲 (初入文壇)
|
[求助]
這個(gè)方程怎樣實(shí)現(xiàn)差分?(偶是新人,目前只有兩個(gè)金幣,不要嫌少。
|
||
|
這樣一個(gè)方程怎樣實(shí)現(xiàn)差分? 在n點(diǎn)t時(shí)刻,方程如圖: 未命名.jpg http://kuai.xunlei.com/d/GXKYUGYYBUXX?p=130497 |
木蟲 (著名寫手)
|
其實(shí)查分本身非常容易,就是因?yàn)樘菀祝杂锌赡芮蠼饽氵@個(gè)方程是不穩(wěn)定的,所以也不是很建議你直接用,不過(guò)你是不是應(yīng)該提供更多的信息? 你給出T為常數(shù),其實(shí)這個(gè)常數(shù)是n的函數(shù)吧?譬如第一輛車T應(yīng)該為零的,否則整個(gè)方程組沒(méi)法開(kāi)始。從第二輛車開(kāi)始,加速度可以用你給的式子計(jì)算,但是第一輛車的的加速度就不能了,所以在你的模擬中,第一輛車的速度其實(shí)應(yīng)該是你提前假定好的(當(dāng)然可以變化,但你是不是應(yīng)該提前給出個(gè)函數(shù)呢),然后也作為方程的初始條件對(duì)吧? 你把加速度寫成dv/dt是應(yīng)該的,也是求解微分方程組的第一步。 至于你的這個(gè)方程組,假設(shè)有n輛車,你需要求解的是一個(gè)有n-1個(gè)方程的方程組,本身很容易求解,如果你懂點(diǎn)計(jì)算機(jī)的話,你可以用matlab的ode45,或者在mathematica里直接敲出式子來(lái)給出初始條件就可以直接得出結(jié)果了;蛘呷绻阌胮ython,直接用scipy.integrate.odeint就可以了,很簡(jiǎn)單的,而且貌似前段時(shí)間有本一個(gè)家伙寫的python數(shù)值計(jì)算的中文教程,你拿來(lái)讀一讀(當(dāng)字典來(lái)查就行了),你這個(gè)問(wèn)題也就幾十行的代碼吧,加上漂亮的繪圖,應(yīng)該也不長(zhǎng)的。 |
木蟲 (著名寫手)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +8 | 新時(shí)代材料 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué) +6 | 妮妮ninicgb 2026-03-27 | 6/300 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
[考研] 339求調(diào)劑 +4 | 烤麥芽 2026-03-27 | 5/250 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-24 | 3/150 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|