计算物理 第6章习题集答案.docVIP

第6章：偏微分方程数值解法 6.1对流方程 ====【6.1.1】=== 考虑边值问题： 如果取： ，，， 采用Crank-Nicolson方法求出 【解】采用Crank-Nicolson方法 由边界条件：，取， ， 由初始条件：，， ， ， ，， ， 6.2抛物形方程 ====【6.2.1 分别用下面方法求定解问题： （1）取，用显式格式计算； (2) 取用隐式格式计算两个时间步。 【解】 （1）, 令, ，，则 ， 计算程序计算程： p6_2_1_1.m, difu.m (2) 对于隐式格式 计算程序: p6_2_1_2.m 和tri.m ,idifu.m 和tri.m ====【6.2. 数值求解扩散方程： 边界条件：；初始条件： 建议显式格式取：，隐式格式：。 解析解为：， ====【6.2. 求下列热传导方程的数值解 【解】 （1）采用显式FTCS方法 ， 整理得： ，， 计算程序：EFTCS.m 和demo_EFTCS.m (2) 采用隱式FTCS方法 ， 整理得： , ， 若的值满足Dirichlet型边界条件，则得到矩阵形式方程 计算程序：IBTCS.m 和demo_IBTCS.m 对于给定牛顿边值问题（即一端给定一阶导数值，另一端给定函数值） ，， 得到： ， 由上面式得： ， 得到： ， 整理得： ， 这时矩阵形式方程为