一维热传导方程_有限差分

分层的一维热传导方程的有限差分法MATLAB程序，是简单的显式差分格式，尽量使用了矩阵运算，没有符号运算。 文件chafen1.m以路面为例，假设有多层建筑材料，热传导系数已经在代码中给出。用户使用时请按自己需要

%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