| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1604 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
1.414木蟲(chóng) (正式寫(xiě)手)
入門(mén)
|
[求助]
matlab 邊界問(wèn)題
|
|
|
我最近在用matlab來(lái)進(jìn)行計(jì)算。但是現(xiàn)在遇到一個(gè)微分方程求解的問(wèn)題,不知道該怎么解決。請(qǐng)各位大俠給支個(gè)招,看該用什么方法來(lái)求解。方程的模型見(jiàn)下圖 方程模型.jpg |

木蟲(chóng) (正式寫(xiě)手)
入門(mén)

鐵桿木蟲(chóng) (職業(yè)作家)
木蟲(chóng) (正式寫(xiě)手)
|
你看,這個(gè)方程有2個(gè)特點(diǎn),一個(gè)是在x=0,有點(diǎn)奇異,第二個(gè)是邊界條件是個(gè)等式,可能會(huì)有多解存在; 為了處理奇異,可以簡(jiǎn)單分析下,按你的邊界條件,y(0)似乎有界。 這樣方程的第二項(xiàng),就有可能無(wú)窮大了,f(0)如果是非0有限值,那按常理y‘(0)應(yīng)該是無(wú)界的。 1. 為了便于實(shí)施求解,不如先假設(shè)y(0)=0,這樣y(1)按條件就有了。 經(jīng)典的兩點(diǎn)邊值問(wèn)題,不過(guò)要處理下方程; y(x)在x=0附近展開(kāi),y(0),y’(0)=f(0)都知道,再算個(gè)y‘’(0),這樣y/x,局部得用級(jí)數(shù)取代,方程被分割為兩部分,0<=x<0.01,0.01<=x<=1.0;再matlab求解。打靶法很容易實(shí)施。 2. 如果第一步,你求得解了,可以考慮第二步試試;類(lèi)似的做法 y(0)非零 級(jí)數(shù)展開(kāi),y’(0),y‘’(0)通過(guò)方程求得; 然后就是這樣了,祝好 |

木蟲(chóng) (正式寫(xiě)手)
入門(mén)

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