偏微分方程习题及答案.pdfVIP

偏微分方程习题及答案 【篇一：偏微分方程数值解法期末考试题答案】 题 答 案 及 评 分 标 准 学年学期：专 业：班 级：课 程：教学大纲：使用教材：教材作者： 出 版 社： 数学与应用数学 数学 偏微分方程数值解法《偏微分方程数值解法》 教学大纲（自编， 2006 ） 《偏微分方程数值解法》陆金甫、关治 清华大学出版社 一、判断题（每小题 1 分，共 10 分） 1、（ o ） 2 、（o ） 3 、（ x ） 4 、（ x ） 5 、（ o ） 6 、（ o ） 7 、（o ） 8 、（ x ） 9 、（x ） 10 、 （o ） 二、选择题（每小题 2 分，共 10 分） 11 、（ d ） 12 、（a ） 13 、 （c ） 14 、（ b ）15 、（c ） 三、填空题（每小题 2 分，共 20 分） ?2?2 16 、2?2? ?x1?x2?2 ?2 17 、a=[4 5 9;23 5 17;11 23 1] 18 、y=exp(-t/3)*sin(3*t) ?xn 19 、help 20 、zeros(m,n)21 、inva(a)*b 或者 a/b22 、 a=sym([cos(x-y),sin(x+y);exp(x-y),(x-1)^3]) ?(s)?1?(s)?c[? ?(s)]2?023 、a[?2(s)]2?2b?2 24 ?? ?? v(?)ed? 25 、 i?x u(xj,tn?1)?u(xj,tn) ? 四、计算题：（每小题 12 分，共 36 分） ?u?u ?0 （x?r,t?0 ）的有限差分方程（两层显示 26 、写成对流方程 ?a ?t?x 格式，用第 n 层计算第 n+1 层），并把有限差分方程改写为便于计 算的迭代格式 ???/h 为网格比。 解：在点 (xj,tn) 处，差分方程为 ?1un?unjj ? ?a n unj?1?uj h ?0 （j?0,?1,?2, ，n?0,1,2, ）（8 分） 便于计算的形式为 ?1nnn ???/h （4 分） un?u?a?(u?ujjj?1j) ， ?u?2u ?a2 的有限差分方程（中心差分格式，用第 n 层 27 、写出扩散方 程?t?x 计算第 n+1 层），并把有限差分方程改写为便于计算的迭代格 式， ???/h2 为网格比。 解 所给对流扩散方程的近似差分方程为 ?1nnun?ununjjj?1?2uj?uj?1 ?a?0 （j?0,?1,?2, ，n?0,1,2, ）（8 分） ?h2 便于迭代计算的格式为 ?1nnnn2 ， （4 分） ???/hun?u?a?(u?2u?u)jjj?1jj?1 ?1nn 28 、计算差分格式 un （其中???/h ，a?0 ）的增 长?unjj?a?(uj?1?uj) ， 因子，并根据 von neumann 条件给出差分格式稳定性条件。 nnijkhn?1nnn 解 令 uj?ve ，代入 uj?uj?a?(uj?1?uj) ，得到 vn?1eijkh?vneijkh?a?vn(1?e?ikh)eijkh 消去公因子有 vn?1?[1?a?(1?e?ikh)]vn （6 分） 增长因子为 g(?,k)?1?a?(1?e?ikh)?1?a?(1?coskh)?a?isinkh 所以有 kh 2 如果 a??1 ，则有 |g(?,k)|?1 ，根据 von neumann 条件，格式是稳 定的。（ 6 分） |g(?,k)|2?[1?a?(1?coskh)]2?[a?isinkh]2?1?4a?(1?a?)sin2