热传导方程的差分格式.pdf

如文档 对 你有用， 请 下 载 支持！ 如文档 对 你有用， 请 下 载 支持！ 一维抛物方程的初边值问题 分别用向前差分格式、向后差分格式、六点对称格式，求解下列问题： 2 在 t 0.05,0.1和0.2 时刻的数值解，并与解析解 u(x,t ) e t sin( x)进行比较 。 1 差分格式形式 2 设空间步长 h 1/ N , 时间步长 0 ，T M ，网比 r / h . （1）向前差分格式 该问题是第二类初边值问题（混合问题），我们要求出所需次数的偏微商的函数 2 u u u( x,t ) ，满足方程 a 2 , 0 x 1,和初始条件 u( x , 0) sin x , 0 x 1， t x 及边值条件 u(0, t ) u(1,t ) 0, t 0 。 已知 sin x 在相应区域光滑，并且在 x 0, l 与边值相容，使问题有唯一充分光滑的 解。 取空间步长 h 1/ N ，和时间步长 T / M ，其中 N , M 都是正整数。用两族平行 ( 0 , 1 , 直 线 x j x j h j 和 t N t k k (k 0,1, , M ) 将 矩 形 域 G {0 x 1,t 0} 分割成矩形网络，网络格节点为 (x ,t ) 。以 G 表示网格内点集合， j k h 即位于矩形 G 的网点集合； Gh 表示闭矩形 G 的网格集合； h Gh Gh 是网格界点的集 合。 向前差分格式，即 k 1 k k k k u u u 2u u j j j 1 j j 1 a 2 f i （1） h f f ( x ) ， i i 其中 ， j 1,2, , N 1, k 1,2, , M 1. 以 r a / h2 表示网比。（ 1）式可改写成如下 ： 此格式 为显 格式。 其矩 阵表达式如下： （2）向后差分格式 向后差分格式，即 如文档 对 你有用， 请 下 载 支持！