function?onedicon

%?和一維非穩態擴散算例一樣的初始和邊界條件
%?采用乘方格式?時間步長為0.001s
%?初始溫度場為200，速度為2m/s，長度為2cm，t=0s時刻東側溫度突然降至0C。
%?時間差分采用全隱式格式
clc
disp（‘一維無內熱源非穩態全隱?乘方格式離散，初始場為200C，右側邊界溫度變為0，左側邊界條件可選擇‘）;
disp（‘------------------‘）
a=input（‘選擇左邊界條件：1—第一類邊界條件，2-第二類邊界條件?:???‘）;
q=0;?t1=0;
if（a==0）
????q=input（‘熱流為?：?‘）;
else
????t1=input（‘選擇左邊界第一類條件溫度，取值在0-200內：?‘）;
end
t2=input（‘選擇右邊界第一類條件溫度，取值在0-200內：??‘）;
u=input（‘please?input?the?velocity?:???‘）;
n=input（‘長度2cm上劃分網格數?:???‘）;
dt=input（‘時間步長?:???‘）;
%幾何條件，物性參數
L=0.02;r=1.0;k=0.03;dx=L/n;
F=r*u;D=k/dx;pe=F/D;
TT=ones（n+11）;
%沒有內熱源
sp=0;
%初始溫度場
tp0=200*ones（n1）;
x=linspace（00.02n+1）;
for?i=1:n
?????xx（i）=（x（i）+x（i+1））/2;
end
xxx（n+1）=0.02;
xxx（1:n）=xx;
%乘方格式賦值影響系數
ae=D*max（0（1-0.1*abs（pe^5）））+max（-F0）;
aw=D*max（0（1-0.1*a

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1758??2013-06-15?16:02??onedimensional.asv

?????文件???????1758??2013-06-15?16:05??onedimensional.m

?????文件???????3125??2013-06-15?16:02??twodimensional.asv

?????文件???????3125??2013-06-15?16:05??twodimensional.m

-----------?---------??----------?-----??----

?????????????????9766????????????????????4