隐式差分方程组差分法matlab,热传导方程几种差分格式的MATLAB数值解法比较

h2)，O(τ2+h2)和O(τ2+h2)［1］．1.2

初、边值条件的处理

1

可算出u在第一层各个节点处的近似值uj．重

可以逐层计算出所有的uj．复使用此式，

隐式格式：将(7)式与离散化的初边值条件联立，得差分方程组：

k+1k+1k+1kk

－rUj－1+(1+2r)Uj－rUj+1=Uj+τfj 2，…，N－1，j=0，1，…，M－1) (k=1，

(11)

0

U=(k=0，1，…，N)jj k

U0=g1k，Uk(j=0，1，…，M)L=g2k将上述方程组改写成矩阵形式

k

对定解条件进行离散化．由初始条件及第一

类边界条件，可直接得到：

U00)=j，(j=0，1，…，