%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????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