%b:人工鋪設的各層厚度（4元素列向量）
%lb:半凍土層厚度
%T:要求解的時間總長（秒）
%k:t的步長
%h:x的步長
%l0fun路面溫度關于時間的函數

%u:要求解的溫度分布u（xt）

%function?u=ufun（blbTkh）

function?chafen1（blbTkh）
%a：材料參數a（5元素列向量，從瀝青層算起）
???a=[0.114795918367347
???0.713610586011342
???0.630434782608696
???0.032258064516129
???0.358910891089109]/1000;

%k：導熱系數
???ka=[0.270000000000000
???1.510000000000000
???1.160000000000000
???0.020000000000000
???0.580000000000000];


%r:步長比：r=k/（h^2）
r=k/（h^2）;

%確定各邊界位置
bound=[0;cumsum（[b;lb]1）];

%初始化各層的u（xt）
%u=zeros（floor（bound（end）/h）floor（T/k））;

%確定各層各維度的長度
%并構建每一層材料參數的矩陣
u1=[floor（bound（2）/h）floor（T/k）];
a1=a（1）*ones（u1（1）1）;
u2=[floor（（bound（3）-bound（2））/h）floor（T/k）];
a2=a（2）*ones（u2（1）1）;
u3=[floor（（bound（4）-bound（3））/h）floor（T/k）];
a3=a（3）*ones

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-09-07?20:42??A3\
?????文件?????3584760??2018-09-07?20:42??A3\ans_to_q1.eps
?????文件??????295856??2018-09-07?20:42??A3\ans_to_q1.fig
?????文件??????150100??2018-09-07?20:42??A3\ans_to_q1.mat
?????文件?????3485734??2018-09-07?20:42??A3\ans_to_q4.eps
?????文件?????1066008??2018-09-07?20:42??A3\ans_to_q4.fig
?????文件??????547473??2018-09-07?20:42??A3\ans_to_q4.mat
?????文件??????184102??2018-09-07?20:42??A3\ans_to_q4_2.fig
?????文件????????1654??2018-09-07?20:42??A3\chafen1.m
?????文件??????????70??2018-09-07?20:42??A3\l0fun.m
?????目錄???????????0??2018-09-07?20:42??A3\q2\
?????文件????????2219??2018-09-07?20:42??A3\run1.m
?????文件????????2123??2018-09-07?20:42??A3\run2.m
?????文件??????358912??2018-09-07?20:42??A3\u（xt）數據.xls
?????文件?????1294336??2018-09-07?20:42??A3\模擬結果.xls