| 6 | 1/1 | 返回列表 |
| 查看: 1691 | 回復(fù): 5 | |||
ghw_nit鐵桿木蟲 (正式寫手)
|
[求助]
一個偏微分方程組的求解問題 已有1人參與
|
|
在一篇文獻(xiàn)中看到的某系統(tǒng)的方程,及相應(yīng)的一些條件如圖所示,文章中說用pdepe就可以求解,但是我就是解不出來呀,請大家看看有沒有什么問題。 360截圖-3516750.jpg |
鐵桿木蟲 (正式寫手)
|
這是我編寫的matlab程序,大家都可以調(diào)一下,看看怎么樣,總是說步長不夠小,我已經(jīng)調(diào)到很小了,比文獻(xiàn)中說的要小好多,都是不能求解的,不知道是怎么回事。 function [c,f,s]=pdefun(x,t,u,dudx) Pe_C=1; Pe_T=1; B_C=1; B_T=2; gamma=20; beta_T=1; u_c=0; y=u(2)/(1+u(2)/gamma); R_11=B_T*B_C*(1-u(1))*exp(y); R_12=-beta_T*(u_c-u(2)); R_13=-B_C*(1-u(1))*exp(y); c=[1 1]'; f=[1/Pe_C 1/Pe_T]'.*dudx; s=[R_13-dudx(1) R_11+R_12-dudx(2)]'; function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t) pl=[-ul(1),-ul(2)]'; ql=[1 1]'; pr=[0 0]'; qr=[1 1]'; function u0=icfun(x) u0=[0 0]'; m=0; t=linspace(0,1,500); x=linspace(0,1,100); sol=pdepe(m,@pdefun,@icfun,@bcfun,x,t); |
鐵桿木蟲 (正式寫手)
鐵桿木蟲 (正式寫手)
新蟲 (初入文壇)
鐵桿木蟲 (正式寫手)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +4 | 谷雨上岸 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點實驗室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +3 | 我愛學(xué)電池 2026-03-23 | 3/150 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|