 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (6360)
- >導(dǎo)師招生 (1286)
- >考研 (1179)
- >基金申請(qǐng) (530)
- >文獻(xiàn)求助 (503)
- >蟲(chóng)友互識(shí) (321)
- >考博 (279)
- >休閑灌水 (192)
- >碩博家園 (160)
- >論文投稿 (151)
- >招聘信息布告欄 (104)
- >博后之家 (88)
- >公派出國(guó) (47)
- >論文道賀祈福 (35)
- >綠色求助(高懸賞) (31)
- >SciFinder/Reaxys (30)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 5036 | 回復(fù): 1 | |||
zhaohuiyuan金蟲(chóng) (著名寫(xiě)手)
|
[求助]
一維熱傳導(dǎo)偏微分方程的數(shù)值解法求助
|
|
各位蟲(chóng)友們，你們好！由于我對(duì)編程軟件運(yùn)用的不是很熟悉，希望得到各位的幫助，我編出來(lái)的程序怎么都運(yùn)行不出來(lái) 不知道是為什么，希望各位高手能指點(diǎn)一下，我要求解的方程是（見(jiàn)圖片），如果哪位能幫上我，讓我的程序能運(yùn)行出結(jié)果來(lái)的話，我將用所有的金幣來(lái)酬謝，若能給出有用意見(jiàn)的也行，我將酌情給金幣，謝謝！ 程序如下 function heat_conduction() %一維齊次熱傳導(dǎo)方程 options={'空間半徑ro','空間點(diǎn)數(shù)N' ,'時(shí)間點(diǎn)數(shù)M','密度p','比熱c','導(dǎo)熱系數(shù)s',‘'對(duì)流傳熱系數(shù)h',}; topic='seting'; lines=1; def={'0.032','100','1000','1714.8','1.48','0.535','61.27'}; b=inputdlg(options,topic,lines,def);%創(chuàng)建輸入對(duì)話 ro=eval(b{1});%eval執(zhí)行字符串 N=eval(b{2}); M=eval(b{3}); p=eval(b{4}); c=eval(b{5}); s=eval(b{6}); h=eval(b{7}); %*************************************************** Fo=s/(p*c); or=ro/N;%空間步長(zhǎng) r=0  r:ro;r=r'; ot=Fo*or^2/a;%時(shí)間步長(zhǎng) tm=M*ot;%熱傳導(dǎo)的總時(shí)間tm t=0 t:tm;t=t'; %計(jì)算初值和邊值 T=zeros(N+1,M+1); Ti=init_fun(r);%fun為目標(biāo)函數(shù)的表達(dá)式字符串或MATLAB自定義函數(shù)的函數(shù)柄 To=border_funo(t); Te=border_fune(t); T(:,1)=Ti;%t=0時(shí)的初始溫度 T(1,  =To;%r=0時(shí)的溫度值T(N+1, =Te;%r=r0時(shí)的溫度%用差分法求出溫度T與半徑r、時(shí)間t的關(guān)系 for k=1:M m=2; while m<=N&&Fo<=(ro/(2*r+or)) T(m,k+1)=(Fo+(or*Fo)/ro)*(T(m+1,k)+(1-(or*Fo)/r+2*Fo)*T(m,k))+Fo*T(m-1,k); m=m+1; end; end; %設(shè)置立體網(wǎng)格 for i=1:M+1 X(:,i)=r; end; for j=1:N+1 Y(j, =t;end mesh(X,Y,T); view([1 -1 1]); xlabel('r'); ylabel('t'); zlabel('T'); function y=init_fun(r)%初值條件t=0 y=14+t.*0 return function y=border_funo(t)%r=0的邊界條件 for k=1:M m=1; while Fo<=0.5 T(m+1,k+1)=(1-2Fo)*T(m+1,k)+2Fo*T(m,k); end; end; return function y=border_fune(t)%r=ro的邊界條件 Tf=40+t.*0; A=(h*ot)/(p*c*or); for k=1:M m=N; while (1-A-Fo)>=0 T(m+1,k+1)=A*Tf+(1-A-Fo)*T(m+1,k))+Fo*T(m,k); end; end; return 需求解的方程 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 好帖子 |
» 本帖已獲得的紅花（最新10朵）
zhaohuiyuan» 猜你喜歡
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有13人回復(fù)
-
085600材料與化工 求調(diào)劑
已經(jīng)有15人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑招人�。。�
已經(jīng)有4人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有17人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有7人回復(fù)
-
申博26年
已經(jīng)有3人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 跟大家請(qǐng)教一個(gè)一維熱傳導(dǎo)方程建立的問(wèn)題（高分求助）
已經(jīng)有3人回復(fù)
- matlab 求解三元方程組
已經(jīng)有9人回復(fù)
- 求助，計(jì)算方程最大值
已經(jīng)有5人回復(fù)
- 中心差分求解二維熱傳導(dǎo)的matlab程序
已經(jīng)有12人回復(fù)
- erfc的求解問(wèn)題。謝謝。
已經(jīng)有5人回復(fù)
- 求助-幫我用mathematics計(jì)算下下面的方程組
已經(jīng)有13人回復(fù)
- 問(wèn)個(gè)熱傳導(dǎo)問(wèn)題的計(jì)算
已經(jīng)有15人回復(fù)
- 求助matlab編程：一維熱傳導(dǎo)問(wèn)題
已經(jīng)有5人回復(fù)
- 數(shù)值求解N-S方程一般用什么編程�。�？？
已經(jīng)有9人回復(fù)
- 求解，熱傳遞方程
已經(jīng)有6人回復(fù)
- 利用有限元方法求解Maxwell方程組、流體力學(xué)方程以及熱傳導(dǎo)方程的軟件
已經(jīng)有10人回復(fù)
- 哪位高手能用matlab編程求解一維瞬態(tài)熱傳導(dǎo)方程，先謝謝啦
已經(jīng)有7人回復(fù)
- 【求助】緊束縛模型與薛定諤方程的區(qū)別與聯(lián)系
已經(jīng)有13人回復(fù)
- 【資源】數(shù)值解 光纖中 非線性薛定諤方程 軟件
已經(jīng)有27人回復(fù)
- 【求助】一道博士生面試題：薛定諤方程和熱傳導(dǎo)方程
已經(jīng)有10人回復(fù)
- 【求助】解含時(shí)薛定諤方程的問(wèn)題
已經(jīng)有5人回復(fù)
- 【求助】O3的薛定諤方程表達(dá)式，及每一項(xiàng)的物理意義
已經(jīng)有4人回復(fù)
zhaohuiyuan
金蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1370.5
- 散金: 3625
- 紅花: 21
- 帖子: 1136
- 在線: 412.6小時(shí)
- 蟲(chóng)號(hào): 2005943
- 注冊(cè): 2012-09-16
- 性別: MM
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
2樓2014-03-06 18:06:57
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhaohuiyuan 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +7 | Ymlll 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +8 | zlingli 2026-03-13 | 8/400 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|