 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 5037 | 回復(fù): 1 | |||
zhaohuiyuan金蟲 (著名寫手)
|
[求助]
一維熱傳導(dǎo)偏微分方程的數(shù)值解法求助
|
|
各位蟲友們，你們好！由于我對編程軟件運用的不是很熟悉，希望得到各位的幫助，我編出來的程序怎么都運行不出來 不知道是為什么，希望各位高手能指點一下，我要求解的方程是（見圖片），如果哪位能幫上我，讓我的程序能運行出結(jié)果來的話，我將用所有的金幣來酬謝，若能給出有用意見的也行，我將酌情給金幣，謝謝！ 程序如下 function heat_conduction() %一維齊次熱傳導(dǎo)方程 options={'空間半徑ro','空間點數(shù)N' ,'時間點數(shù)M','密度p','比熱c','導(dǎo)熱系數(shù)s',‘'對流傳熱系數(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)建輸入對話 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;%空間步長 r=0  r:ro;r=r'; ot=Fo*or^2/a;%時間步長 tm=M*ot;%熱傳導(dǎo)的總時間tm t=0 t:tm;t=t'; %計算初值和邊值 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時的初始溫度 T(1,  =To;%r=0時的溫度值T(N+1, =Te;%r=r0時的溫度%用差分法求出溫度T與半徑r、時間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» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 跟大家請教一個一維熱傳導(dǎo)方程建立的問題（高分求助）
已經(jīng)有3人回復(fù)
- matlab 求解三元方程組
已經(jīng)有9人回復(fù)
- 求助，計算方程最大值
已經(jīng)有5人回復(fù)
- 中心差分求解二維熱傳導(dǎo)的matlab程序
已經(jīng)有12人回復(fù)
- erfc的求解問題。謝謝。
已經(jīng)有5人回復(fù)
- 求助-幫我用mathematics計算下下面的方程組
已經(jīng)有13人回復(fù)
- 問個熱傳導(dǎo)問題的計算
已經(jīng)有15人回復(fù)
- 求助matlab編程：一維熱傳導(dǎo)問題
已經(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ù)
- 【求助】解含時薛定諤方程的問題
已經(jīng)有5人回復(fù)
- 【求助】O3的薛定諤方程表達(dá)式，及每一項的物理意義
已經(jīng)有4人回復(fù)
zhaohuiyuan
金蟲 (著名寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1370.5
- 散金: 3625
- 紅花: 21
- 帖子: 1136
- 在線: 412.6小時
- 蟲號: 2005943
- 注冊: 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ā)表 | |
|---|---|---|---|---|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +14 | enenenhui 2026-03-13 | 15/750 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +6 | sllhht 2026-03-19 | 7/350 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|