原创力文档- 1、本文档共12页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
利用matlab程序解决热传导问题
哈佛大学能源与环境学院
课程作业报告
作业名称：传热学大作业——利用matlab程序解决热传导问题
院系：能源与环境学院
专业：建筑环境与设备工程
学号：5201314
姓名：盖茨比
2015年6月8日
题目及要求
原始题目及要求
各节点的离散化的代数方程
源程序
不同初值时的收敛快慢
上下边界的热流量（λ=1W/(m℃））
计算结果的等温线图
计算小结
题目：已知条件如下图所示：
二、各节点的离散化的代数方程
各温度节点的代数方程
ta=(300+b+e)/4 ; tb=(200+a+c+f)/4; tc=(200+b+d+g)/4; td=(2*c+200+h)/4
te=(100+a+f+i)/4; tf=(b+e+g+j)/4; tg=(c+f+h+k)/4 ; th=(2*g+d+l)/4
ti=(100+e+m+j)/4; tj=(f+i+k+n)/4; tk=(g+j+l+o)/4; tl=(2*k+h+q)/4
tm=(2*i+300+n)/24; tn=(2*j+m+p+200)/24; to=(2*k+p+n+200)/24; tp=(l+o+100)/12
三、源程序
【G-S迭代程序】
【方法一】
函数文件为：
function [y,n]=gauseidel(A,b,x0,eps)
D=diag(diag(A));
L=-tril(A,-1);
U=-triu(A,1);
G=(D-L)\U;
f=(D-L)\b;
y=G*x0+f;
n=1;
while norm(y-x0)=eps
x0=y;
y=G*x0+f;
n=n+1;
end
命令文件为：
A=[4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0;
-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0;
0,-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0;
0,0,-2,4,0,0,0,-1,0,0,0,0,0,0,0,0;
-1,0,0,0,4,-1,0,0,-1,0,0,0,0,0,0,0;
0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0,0;
0,0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0;
0,0,0,-1,0,0,-2,4,0,0,0,-1,0,0,0,0;
0,0,0,0,-1,0,-1,0,4,0,0,0,-1,0,0,0;
0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0,0;
0,0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0;
0,0,0,0,0,0,0,-1,0,0,-2,4,0,0,0,-1;
0,0,0,0,0,0,0,0,-2,0,0,0,24,-1,0,0;
0,0,0,0,0,0,0,0,0,-2,0,0,-1,24,-1,0;
0,0,0,0,0,0,0,0,0,0,-2,0,0,-1,24,-1;
0,0,0,0,0,0,0,0,0,0,0,-1,0,0,-1,12];
b=[300,200,200,200,100,0,0,0,100,0,0,0,300,200,200,100];
[x,n]=gauseidel(A,b,[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],1.0e-6)
xx=1:1:4;
yy=xx;
[X,Y]=meshgrid(xx,yy);
Z=reshape(x,4,4);
Z=Z
contour(X,Y,Z,30)
Z =
139.6088 150.3312 153.0517 153.5639
108.1040 108.6641 108.3119 108.1523
84.1429 67.9096 63.3793 62.4214
20.1557 15.4521 14.8744 14.7746
【方法2】 t=zeros(5,5);
t(1,1)=100;
t(1,2)=100;
t(1,3)=100;
t(1,4)=100;
t(1,5)=100;
t(2,1)=200;
t(3,1)=200;
t(4,1)=200;
t(5,1)=200;
for i=1:10
t(2,2)=(300+t(3,2
文档评论（0）