 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 582 | 回復(fù): 1 | ||
藍諳至尊木蟲 (職業(yè)作家)
|
[求助]
自己寫了一個二維熱傳導(dǎo)方程的向前差分法，提示有兩個錯誤，但不知怎么改。幫幫忙啦~
|
|
function u=uu(x,y)%邊界函數(shù) u=(x.^2)*sin(y)-(y.^2)*sin(x); end function [x,y,u]=heat(x0,xn,y0,yn,t0,tn,h1,h2,a)%二維熱傳導(dǎo)方程的向前差分 n1=(xn-x0)/h1;%x y的點的個數(shù) n2=(tn-t0)/h2;%t 的點的個數(shù) for i=1:n1 for j=1:n1 u1(i,j,1)=0;%第一層初始化 end end x(1)=x0;y(1)=y0; for i=1:n1 x(i+1)=x(i)+h1; y(i+1)=y(i)+h1;%求出x y的坐標 end for i=1:n1 for j=1:n2 u1(1,i,j)=feval(uu,x(1),y(i));%x=x0的邊值 end end for i=1:n1 for j=1:n2 u1(n1,i,j)=feval(uu,x(n1),y(i));%x=xn的邊值 end end for i=1:n1 for j=1:n2 u1(i,1,j)=feval(uu,x(i),y(1));%y=y0的邊值 end end for i=1:n1 for j=1:n2 u1(i,n1,j)=feval(uu,x(i),y(xn));%y=yn的邊值 end end r=a*h2/(h1^2); for k=2  n2-1)for i=2 n1-1)for j=2 n1-1)u(i,j,k+1)=u(i,j,k)+r*(u(i+1,j,k)-4*u(i,j,k)+u(i-1,j,k)+u(i,j+1,k)+u(i,j-1,k)) end end end for i=1:n1 for j=1:n1 u(i,j)=u1(i,j,n2); end end |
回復(fù)此樓» 猜你喜歡
焦慮
已經(jīng)有13人回復(fù)
-
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有9人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有5人回復(fù)
-
285化工學(xué)碩求調(diào)劑（081700）
已經(jīng)有11人回復(fù)
-
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有9人回復(fù)
-
274求調(diào)劑
已經(jīng)有6人回復(fù)
-
求材料調(diào)劑
已經(jīng)有10人回復(fù)
-
本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑
已經(jīng)有4人回復(fù)
-
國自科面上基金字體
已經(jīng)有8人回復(fù)
藍諳
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 10510.4
- 散金: 1340
- 紅花: 20
- 沙發(fā): 7
- 帖子: 3250
- 在線: 280.3小時
- 蟲號: 3885761
- 注冊: 2015-05-23
- 性別: MM
- 專業(yè): 計算機應(yīng)用技術(shù)
2樓2015-11-28 22:51:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 藍諳 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料調(diào)劑 +10 | 隔壁陳先生 2026-03-12 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +5 | 陌の森林 2026-03-18 | 5/250 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 331求調(diào)劑（0703有機化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|