function?CDstead
%?===================一維穩態問題===================
%?===================采用中心差分并與解析解比較===================
t0=1;t1=0;
r=1;L=1;k=0.1;u=2.5;
n1=100;
dx=L/n1;
x=linspace（0Ln1+1）;
for?i=1:n1
????xx（i）=（x（i）+x（i+1））/2;
end
F=r*u;
D=k/dx;
pe=F/D;
A=sparse（n1n1）;
b=zeros（n11）;
%?===================影響系數賦值===================
ae=D-F/2;
aw=D+F/2;
ap=ae+aw;
%?===================源項賦值===================
??sp1=-（2*D+F）;??sc1=2*D+F;???
??sp2=-（2*D-F）;sc2=2*D-F;
??b（1）=sc1*t0;
??b（n1）=sc2*t1;
%?===================對矩陣A主對角線賦值===================
for?i=2:n1-1
????A（ii）=ap;
end
%?===================兩邊界賦值===================
???A（11）=ae-sp1;
???A（n1n1）=aw-sp2;
%?===================對角線賦值===================
for?i=1:n1-1
???A（ii+1）=-ae;
???A（i+1i）=-aw;
end
T=A\b;
plot（T‘-*‘）;
hold?on
%===================對應該問題參數下的精確解===================
d=（1-exp（F*xx/k））/（exp（F*L/k）-1）+1;
plot（d‘r‘）
title（‘一維穩態對流擴散問題‘）;
legend（‘中心差分解‘‘解析解‘）
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????743??2011-04-27?22:29??CDstead.asv
?????文件????????1060??2011-04-28?19:59??CDstead.m
?????文件????????1031??2011-05-09?15:51??CFstead.m
?????文件????????2125??2011-04-26?21:10??Copy_of_onedicon.m
?????文件????????4045??2011-04-29?12:19??Copy_of_secdim.asv
?????文件????????4030??2011-04-29?12:23??Copy_of_secdim.m
?????文件??????????36??2011-04-26?18:28??matrixa.m
?????文件??????473647??2011-04-27?10:47??mccExcludedFiles.log
?????文件????????1978??2011-05-10?21:36??odim.m
?????文件????????2388??2011-04-27?22:39??onedicon.asv
?????文件????????2124??2011-05-09?15:53??onedicon.m
?????文件????????3990??2011-04-30?10:18??secdim.m
?????文件????????2932??2011-04-26?16:05??secdim2.asv
?????文件????????3952??2011-04-29?13:31??secdim2.m
?????文件????????3741??2011-05-10?22:16??sedim.m
?????文件????????1917??2011-05-10?22:22??yiwei.m