#### Readme
- hw1_explict:程序入口
- 网格大小dx,dy,dt设定
- 初值设定 (平板网格编号成一列)
- 调用显式求温度矩阵函数
- 调用画图函数
- heatExplict:显式实现温度矩阵求解函数
----实现方案均来自于教材。简单说明一下matlab实现。
- Dirichlet boundary
- 圆形边界分8种情况讨论,两个临近网格点在圆内的情况有4种,只有一个临近网格点在圆内的情况有4种。
- 用标记矩阵tag记录点是否在圆内,在圆内计为0,在圆外计为1。
- 通过圆与网格点交点计算出交点坐标,按照书中公式用显式方法求解。
- Neumann boundary
- 圆形边界分8种情况讨论,基本与狄利克雷边界情况相同。但是这里更为复杂,在一上基础上再讨论法线延长线与网格交于垂直方向,或者交于水平方向。
- 通过法线斜率和圆与网格交点坐标可确定法线延长线与网格交点具体坐标。
- 法线延长线与网格交点温度可通过插值法确定。
- 用法线延长线与网格交点温度确定圆与网格交点温度。
- 其他即与狄利克雷边界处理类似,用显式方法求解。
- plotHeatFigure:画图函数
- hw1_result.pdf 运行结果,四个时刻点截图说明。