[問題] 二維熱傳問題已回收
這是我目前在解的二維熱傳問題,
三邊固定是T1,一邊是T2,以下是我編寫的程式碼,
可以請大家幫我看看哪裡出錯嗎?
因為跑出來的結果很奇怪~~ 謝謝
clc
clear all
T1=70;
T2=200;
W=2;
L=2;
for x=0:0.2:2
for y=0:0.2:2
error1=100.0;
n=0;
while error1>0.01
n=n+1;
w1(n)=2/n*sin(n*pi*x/L).*(sinh(n*pi*y/L)./sinh(n*pi*W/L));
w(n)=(2/pi)*(sum(w1(1:n)));
T(n)=(T2-T1)*w(n)+T1;
if n>1
error1=abs((T(n-1)-T(n))/T(n-1));
end
end
T_final(floor(x/0.2)+1,floor(y/0.2)+1)=T(n);
w1(:)=0;
end
end
--
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.39.133
→
12/07 23:48, , 1F
12/07 23:48, 1F
→
12/08 00:04, , 2F
12/08 00:04, 2F
→
12/08 00:05, , 3F
12/08 00:05, 3F
→
12/08 00:07, , 4F
12/08 00:07, 4F
→
12/08 00:09, , 5F
12/08 00:09, 5F
討論串 (同標題文章)