matlab 牛顿冷却定律,传热学MATLAB度分布大作业完整版.doc

传热学MATLAB度分布大作业完整版

传热学大作业(第四章)

姓名： 张宝琪 学号：

一、题目及要求

各节点的离散化的代数方程

源程序

不同初值时的收敛快慢

上下边界的热流量(λ=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&#