二维传热数值计算matlab程序

二维流体流动数值计算matlab程序，供学习参考！

%Temperature distribution in a rod
%See example 9
%Discretization method: Finite difference
% Solution method: SOR
clear all;
a=[];b=[];c=[];d=[];x=[];T=[];analytical=[];
nn = input（‘Number of increment = ‘）
n = nn+1;% number of  grid points
L = 0.6;
dx = L/nn; % size of increment
Qprim = 50000; % heat source
lambda = 20; % thermal conductivity
for k = 1:n % Set coefficients
   a（k）=2.0;
   b（k）=1.0;
   c（k）=1.0;
   d（k）=Qprim*dx^2/lambda;
   T（k） = 30.0; % start value
   x（k） = （k-1）*dx;
end; %
a（1） = 3;
b（1） = 4;
c（1） = 0;
% Solver SOR if omega = 1 Gauss-Siedel
omega=2/（1+sin（pi/nn））;
counter = 0;
maxit = 200;
sumres = 1.;
maxres = 1.0e-5;
while （（sumres > maxres）&（counter < maxit） ）
   for k = 2:n-1 % SOR
      T（k）=omega*（b（k）*T（k+1）+c（k）*T（k-1）+d（k））/a（k

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1382  2003-08-05 10:25  Example1a.m
     文件         966  2003-08-05 10:33  Example1b.m
     文件        1225  2003-08-05 10:34  Example1c.m
     文件        1085  2003-08-05 10:36  Example1FVM.m
     文件        1496  2003-08-05 12:40  Example2a.m
     文件        1531  2003-08-06 14:12  Example2aFVM.m
     文件        1891  2003-08-08 12:34  Example2b.m
     文件        1795  2003-08-07 14:21  Example2bFVM.m
     文件        2600  2003-08-08 13:40  Example4.m
     文件        3036  2003-08-07 14:21  Example_3_SOR.m
     文件        3483  2003-08-07 14:21  Example_3_TDMA.m
     文件         218  2007-01-22 19:53  www.pudn.com.txt