微分方程数值解教学大纲
《微分方程数值解》教学大纲
一、课程的基本信息
课程名称：《微分方程数值解》
英文名称：Numerical solution for differential equaiton
课程性质：专业方向选修课
课程编号：1623313002
周学时：3学时
总学时：48学时（理论40+实验8）
学分：3学分
适用专业：
适用于信息与计算科学专业
预备知识：数值计算、常微分方程、数值逼近、数理方程
课程教材：
李立康主编，《微分方程数值解法》，复旦大学出版社出版、1999年
参考书目：
[1] 戴嘉尊主编，《微分方程数值解法》，东南大学出版社、2008年．
[2] 李荣华主编，《微分方程数值解法》(第四版)，高等教育出版社、2009年.
考核方式：考查
制定时间：2013年10月制定
二、课程的目的与任务
《微分方程数值解》是高等院校信息与计算科学专业的专业选修课之一。本课程主要内容为常微分方程和偏微分方程的数值求解问题，包括各种差分方法，有限元方法等的基本理论。通过微分方程数值解的教学，使学生了解和掌握微分方程数值解这一学科的基本概念、理论，培养学生的理论思维能力，为从事信息与计算科学学科的教学和研究打下一定的理论基础。
通过本课程的学习，学生应熟练掌握常微分方程和偏微分方程的常用数值求解方法和分析手段，从能力方面，应使学生初步认识如何从实际问题出发，建立微分方程数学模型，将连续问题离散化，由微分方程转化为差分方程，利用计算机实现数值方法求解一个微分方程的定解问题，并对结果给以几何解释。从教学方法上，着重体现思维方式，注重解决实际问题的方法以及利用计算机进行科学计算的能力
培养。
第一章微分方程数值解法（10学时）
一、本章基本要求
1．掌握线性多步方法，Runge-Kutta方法，Gear方法等计算常微分方程的计算格式；2．掌握相容性，稳定性，绝对稳定性概念和相互关系；
3．了解刚性问题和辛计算格式。
二、教学内容
1．微分方程模型和定性理论
2．计算格式：线性多步方法和高阶单步方法
3．稳定性和收敛性分析
4．刚性问题和其他
第二章椭圆方程差分方法（8学时）
一、本章基本要求
1．掌握椭圆型方程的五点、九点差分格式和有限体积法；
2．掌握极值原理，收敛性分析和误差估计。
二、教学内容
1．椭圆方程模型和定性理论
2．椭圆边值问题的差分方法
3．椭圆差分方程的形态研究
第三章发展方程差分方法（12学时）
一、本章基本要求
1．掌握抛物型方程和双曲型方程的差分方法；
2．掌握稳定性分析包括直接法，分离变量法，最大模方法，传播因子法； 3．掌握Courant-Friedrichs-Lewy 条件和V on Neumann 分析方法。 二、教学内容
1．发展方程模型和定性理论 2．抛物型方程的差分方法
3．抛物型方程差分方法的稳定性分析 4．双曲型方程的差分方法 第四章 有限元方法（10学时） 一、本章基本要求
1．掌握椭圆型方程的变分法； 2．掌握有限元计算格式的推导； 3．掌握基本的有限元误差分析。 二、教学内容 1．变分法概述
2．椭圆型方程的有限元方法 3．有限元方法分析
课内实验项目及基本要求（8学时） 实验一 Euler 法与改进Euler 法 一、目的和要求
1．通过计算比较两种方法的收敛速度和精度； 2．初步认识数值解法的重要性。 二、实验内容
考虑下列常微分方程的初值问题
5(01)u u t '=-<< ， (0)1u =
步长：h=0.1,0.05
采用Euler 法与改进Euler 法进行求解 实验二 线性多步法和Runge-Kutta 法 一、目的和要求
1．采用任何一种线性多步法（如admas 内插，外插法），求解微分方程；
2．采用4阶Runge-Kutta 法，求解微分方程； 3．用实验一的方法求解评论几种方法的优劣。 二、实验内容
考虑下列常微分方程的初值问题
57cos 2(03)y y y x x '=-++≤≤ ， y(0)1=
采用线性多步法、Euler 法、改进Euler 法、Runge-Kutta 求解 实验三 两点边值问题的线性有限元法 一、目的和要求
1. 用线性元求解两点边值问题；
2. 掌握有限元数值解的求解过程。 二、实验内容
用线性元求解两点边值问题：
2,01u u x x ''-+=<< ，(0)(1)0u u ==
执 笔 人：
审 核 者：
教学院长：
偏微分方程数值解期末试题及标准答案
偏微分方程数值解试题(０6B ） 参考答案与评分标准 信息与计算科学专业 一(10分）、设矩阵A 对称，定义)(),(),(2 1)(n R x x b x Ax x J ∈-=，)()(0x x J λλ?+=.若0)0('=?，则称称0x 是)(x J 的驻点(或稳定点).矩阵A 对称(不必正定)，求证0x 是)(x J 的驻点的充要条件是:0x 是方程组 b Ax =的解 解: 设n R x ∈0是)(x J 的驻点，对于任意的n R x ∈,令 ),(2),()()()(2 000x Ax x b Ax x J x x J λλλλ?+-+=+=, （3分) 0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0，则有0||||),(2000=-=--b Ax b Ax b Ax ，得到b Ax =0． (3分) 反之,若n R x ∈0满足b Ax =0,则对于任意的x ,)(),(2 1)0()1()(00x J x Ax x x J >+==+??,因此0x 是)(x J 的最小值点. (４分） 评分标准：)(λ?的展开式3分, 每问3分，推理逻辑性１分 二（10分）、 对于两点边值问题:?????==∈=+-=0 )(,0)(),()('b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ],[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的Ritz 形式和Galerkin 形式的变分方程。 解: 设}0)(),,(|{11=∈=a u b a H u u H E 为求解函数空间,检验函数空间．取),(1b a H v E ∈,乘方程两端,积分应用分部积分得到 (３分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(1 b a H v E ∈? 即变分问题的Galerkin 形式. （3分)
偏微分方程数值解法试题与答案
一．填空（1553=?分） 1．若步长趋于零时，差分方程的截断误差0→lm R ，则差分方程的解lm U 趋近于微分方 程的解lm u . 此结论_______（错或对）； 2．一阶Sobolev 空间{} )(,,),()(21 Ω∈''=ΩL f f f y x f H y x 关于内积=1),( g f _____________________是Hilbert 空间； 3．对非线性（变系数）差分格式，常用 _______系数法讨论差分格式的_______稳定性； 4．写出3 x y =在区间]2,1[上的两个一阶广义导数：_________________________________， ________________________________________； 5．隐式差分格式关于初值是无条件稳定的. 此结论_______（错或对）。 二．（13分）设有椭圆型方程边值问题 用1.0=h 作正方形网格剖分 。 （1）用五点菱形差分格式将微分方程在内点离散化； （2）用截断误差为)(2 h O 的差分法将第三边界条件离散化； （3）整理后的差分方程组为 三．（12）给定初值问题 x u t u ??=?? ， ()10,+=x x u 取时间步长1.0=τ，空间步长2.0=h 。试合理选用一阶偏心差分格式（最简显格式）, 并以此格式求出解函数),(t x u 在2.0,2.0=-=t x 处的近似值。 1．所选用的差分格式是： 2．计算所求近似值： 四．（12分）试讨论差分方程 ()h a h a r u u r u u k l k l k l k l ττ + - = -+=++++11,111 1 逼近微分方程 0=??+??x u a t u 的截断误差阶R 。 思路一：将r 带入到原式，展开后可得格式是在点（l+1/2,k+1/2）展开的。 思路二：差分格式的用到的四个点刚好是矩形区域的四个顶点，可由此构造中心点的差分格 式。
微分方程数值解试题库2011(试题参考)
---------------------------------------------------------------------------------------------------------------------- 《常分方程数值解法》试题一及答案 ---------------------------------------------------------------------------------------------------------------------- 1．用欧拉法解初值问题???1 =060≤≤0--='2)() .(y x xy y y ，取步长 h =0.2.计算 过程保留4位小数。 解：h =0.2, f (x )=－y －xy 2.首先建立欧拉迭代公式 ),,k )(y x (y .y hx hy y )y ,x (hf y y k k k k k k k k k k k 21042021=-=--=+=+ 当k =0，x 1=0.2时，已知x 0=0,y 0=1，有 y (0.2)≈y 1=0.2×1(4－0×1)＝0.800 0 当k ＝1，x 2=0.4时，已知x 1=0.2, y 1=0.8，有 y (0.4)≈y 2=0.2×0.8×(4－0.2×0.8)＝0.614 4 当k =2，x 3=0.6时，已知x 2=0.4,y 2=0.614 4，有 y (0.6)≈y 3=0.2×0.614 4×(4－0.4×0.4613)＝0.800 0 2．对于初值问题? ??1=0='2 )(y xy y 试用(1)欧拉法；(2)欧拉预报－校正公式； (3)四阶龙格－库塔法分别计算y (0.2),y (0.4)的近似值. 3．证明求解初值问题的梯形公式是 y k +1=y k +)],(),([2 11+++k k k k y x f y x f h , h =x k +1－x k (k =0,1,2,…,n －1)，
微分方程数值解--大纲
偏微分方程数值解 （Numerical Methods for Partial Differential Equations） 课程代码：10210801 学位课程/非学位课程：非学位课程 学时/学分：46/3 课程简介： 《偏微分方程数值解》是数学类专业必修的一门专业课。主要内容包括：变分形式和Galerkin有限元法、椭圆型方程的差分方法、抛物型方程的差分方法、双曲型方程的差分方法、离散方程的解法。通过本课程的学习，使学生掌握求解偏微分方程数值解的基本方法，能够根据具体的微分方程使用合适的计算方法。 一、教学目标 1、知识水平教学目标 偏微分方程数值解课程的教学，要使学生掌握椭圆型微分方程、抛物型微分方程、双曲型微分方程等典型方程的差分方法，了解与之相关的理论问题，理解变分原理、有限元方法以及离散方程的解法，理解各种计算方法的收敛条件和收敛速度。 2、能力培养目标 通过偏微分方程数值解课程教学,应注意培养学生以下能力： （1）连续问题离散化能力——掌握科学的思维方法，能够使用差分方法和有限元方法的各种格式对三类典型方程进行离散化处理。 （2）算法分析与设计能力——结合各类偏微分方程的特点，设计各种计算方法，对计算方法的收敛条件和收敛速度等进行分析，具体设计易于上机实现的算法。（3）离散方程组的快速求解能力——理解离散方程组的特点，使用数学软件编程，具体上机实现，进行数值模拟的动手能力。 3、素质培养目标 通过数学物理方程课程教学,应注重培养学生以下素质： （1）具体问题有限化——善于对现实世界中得到的偏微分方程进行有限差分、有限元分析的有限化思想素养。 （2）数值解法定性化——通过学习，引导学生树立偏微分方程数值求解的基本原则，培养学生对数值方法中的稳定性、收敛性和误差等进行定性分析的素质。（3）算法实现程序化——培养学生的创造性和具体实现程序化的思维，使学生学会用数学中算法的观点思考实际问题，用程序和计算机解决数学问题。 二、教学重点与难点 1、教学重点：椭圆型、抛物型、双曲型等微分方程的差分方法，有限元方法。 2、教学难点：各种计算方法的稳定性、收敛性和误差分析，变分形式。 三、教学方法与手段 以教师讲授为主，安排上机实验，辅以习题课、课堂讨论、小论文，注重理论联系实际。 四、教学内容与目标 教学内容教学目标课时分配 （46学时） 1. 边值问题的变分形式 6 二次函数的极值掌握 两点边值问题掌握
偏微分方程数值解复习题(2011硕士)
偏微分方程数值解期末复习（2011硕士） 一、考题类型 本次试卷共六道题目，题型及其所占比例分别为： 填空题20%；计算题80% 二、按章节复习内容 第一章 知识点：Euler法、向前差商、向后差商、中心差商、局部截断误差、整体截断误差、相容性、收敛性、阶、稳定性、显格式、隐格式、线性多步法、第一特征多项式、第二特征多项式、稳定多项式、绝对稳定等； 要求: 会辨认差分格式, 判断线性多步法的误差和阶; 第二章 知识点：矩形网格、(正则,非正则)内点、边界点、偏向前(向后,中心)差商、五点差分格式、增设虚点法、积分插值法、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和、稳定性等； 要求: 建立椭圆型方程边值问题的差分格式, 极值原理; 第四章 知识点：最简显格式、最简隐格式、CN格式、双层加权格式、Richardson 格式、网格比、传播因子法(分离变量法) 、传播因子、传播矩阵、谱半径、von Neumann条件、跳点格式、ADI格式、线性椭圆型差分格式、极值原理、比较定理、五点差分格式的相容收敛和稳定性等； 要求: 建立抛物型方程边值问题的差分格式, 计算局部截断误差; 第五章 知识点：左偏心格式、右偏心格式、中心格式、LF格式、LW格式、Wendroff 格式、跳蛙格式、特征线、CFL条件等； 要求: 建立双曲型方程边值问题的差分格式, 计算局部截断误差; 第七章 要求: 会用线性元(线性基)建立常微分方程边值问题的有限元格式
三 练习题 1、 已知显格式21131()22 n n n n u u h f f +++-=-，试证明格式是相容的，并求它的阶。 P39+P41 2、用Taylor 展开原理构造一元函数一阶导数和二阶导数的数值微分公式。 提示：向前、向后和中心差商与一阶导数间关系，二阶中心差商与二阶导数 之间的关系 课件 3、用数值微分方法或数值积分方法建立椭圆型方程 2222(,),(,),u u f x y x y x y ??--=?∈Ω?? :01,01x y Ω≤≤≤≤ 内点差分格式。 P75+课件 4、构造椭圆型方程边值问题的差分格式. P101 (4)题 5、构建一维热传导方程220,(0)u u Lu a a t x ??=-=>??的数值差分格式(显隐格式等)。 参考P132-135相关知识点 6、设有逼近热传导方程22(0)u u Lu a f a const t x ??≡-==>??的带权双层格式 ()()1111111122(1)2k k j j k k k k k k j j j j j j u u a u u u u u u h θθτ++++-+-+-??=-++--+?? 其中[0,1]θ∈，试求其截断误差。并证明当2 1212h a θτ=-时，截断误差的阶最 高阶为24()O h τ+。 P135+P165+课件 7、传播因子法证明抛物型方程22(0)u u Lu a f a const t x ??≡-==>??的最简显隐和六点CN 格式稳定性。 P156+课件 8、对一阶常系数双曲型方程的初边值问题 0,0,0,0,(,0)(),0,(0,)(),0, u u a t T x a t x u x x x u t t t T φψ???+=<≤<<∞>?????=≤<∞??=≤≤?
偏微分方程数值解试题06B答案
专业班级 姓名 学号 开课系室数学与计算科学学院 考试日期
偏微分方程数值解试卷 一（15分）、（1）简述用差分方法求解抛物型方程初边值问题的数值解的一般步骤.(2)写出近似一阶偏导数 n m x u |??的三种有限差分逼近及其误差阶,写出近似 n m x u |22 ??的差分逼近及其误差阶. 评分标准： （1） 7分，三个离散4分，其他步骤3分 （2） 8分，每个格式及误差2分。 二（15分）、（1）以抛物型方程的差分格式为例，解释差分格式的相容性,稳定性和收敛性概念,分析相容性,稳定性和收敛性与误差的关系，简述 Lax 等价性定理。(2) 简述差分格式稳定性分析的Fourier 级数法（或称为Neumann Von 方法，分离变量法）的一般步骤。 （1）8分，解释概念6分，等价关系2分 （2）7分，典型波2分，放大因子与条件3分，其他2分 三（20分）、对于边值问题 ?? ???=?=∈=??+???0 |) 1,0()1,0(),(,92 222G u G y x y u x u （1）建立该边值问题的五点差分格式（五点棱形格式又称正五点格式），推导截 断误差的阶。 （2）取3/1=h ，求边值问题的数值解（写出对应的方程组的矩阵形式并求解） （3）就取5/1=h 的情况写出对应方程组的系数矩阵（用分块矩阵表示）。 解：（1）7分，离过程与格式
第二页（共五页） 四（20分）、对于初边值问题??? ????≤≤==<<=≤<<
微分方程数值解试卷
中国矿业大学2008～2009学年第 1 学期 《微分方程数值解法》试卷（B ）卷 考试时间：100 分钟 考试方式：半开卷 学院 班级 姓名 序号 1、下面关于Euler 公式的结论哪些是正确的（打√）？哪些是错误的（打×）？ （1）二阶方法；（2）一阶方法；（3）显式公式；（4）隐式公式；（5）是数值稳定的。 2、如果微分方程为,(0)1u tu u '==，则用Taylor 级数法求()u h 时，它的前两项为： 。 3、二阶差商 11 2 2i i i u u u h +--+近似二阶导数()i u x ''局部截断误差为 。 4、算术平均11 2 i i u u +-+近似函数值()i u x 的局部截断误差为 。 5、在课本P98差分方程(3.10)中，第二个方程的局部误差是什么？ 。 6、函数空间0()C I ∞ 中函数满足什么性质？ 。 二、（10分）求解常系数齐次差分方程21120,1,2, 1,1 i i i u u u i u u ++-+==?? =-=?的解。 三、（25分）已知数值解公式21132(2)m m m m m u u u h f f +++-+=- （1）写出与它们对应的特征多项式。 （2）这个多步法相容吗？ （3）利用课本P47公式（2.66）求公式的局部截断误差的主项。 （4）讨论这个算法的零稳定性。 （5）求这个算法的绝对稳定区间。 四、（10分）试利用初值问题的数值解公式 11 11(,) (,)n n n n n n n n u u hf x u u u hf x u ++++=+?? =+? （1）构造一个PECE 预测校正系统；
偏微分方程数值解期末试题及答案(内容参考)
偏微分方程数值解试题(06B) 参考答案与评分标准 信息与计算科学专业 一（10分）、设矩阵A 对称，定义)(),(),(2 1 )(n R x x b x Ax x J ∈-= ，)()(0x x J λλ?+=.若0)0('=?,则称称0x 是)(x J 的驻点（或稳定点）.矩阵A 对称（不必正定），求证0x 是)(x J 的驻点的充要条件是：0x 是方程组 b Ax =的解 解: 设n R x ∈0是)(x J 的驻点,对于任意的n R x ∈,令 ),(2 ),()()()(2 000x Ax x b Ax x J x x J λλλλ?+ -+=+=, (3分) 0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有 0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若 n R x ∈0满足 b Ax =0,则对于任意的 x ,)(),(2 1 )0()1()(00x J x Ax x x J >+ ==+??,因此0x 是)(x J 的最小值点. (4分) 评分标准:)(λ?的展开式3分, 每问3分,推理逻辑性1分 二（10分）、 对于两点边值问题：????? ==∈=+-=0 )(,0)() ,()(' b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ] ,[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的Ritz 形式和 Galerkin 形式的变分方程。 解: 设}0)(),,(|{11 =∈=a u b a H u u H E 为求解函数空间,检验函数空间.取),(1 b a H v E ∈,乘方程两端,积分应用分部积分得到 (3分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(1 b a H v E ∈? 即变分问题的Galerkin 形式. (3分)
微分方程数值解
浅谈微分方程数值解法（双语）课堂教学模式 姓名：肖录明 学号：11301010232 摘要：微分方程数值解是高等院校信息与计算科学专业的一门重要专业基础课。这是一门本具有较强实际背景，专门研究科学计算的课程。这门课程理论性较强，公式多而且难记。我们还需要通过一门语言（比如MATLAB语言）来实现我们数值计算算法。由于解微分方程在科学计算中极为常见，故学好这门课程就非常有用且能为以后的学习打下基础。在我国双语教学正在慢慢的被倡导，且益处明显。本文主要探讨该课程的双语教学模式,并对在学习过程中出现的一些问题进行了思考。 关键词：微分方程数值解法双语教学科学计算 1引言 微分方程数值解法在数值分析中占有重要的地位，它以逼近论，数值代数等学科为基础，反过来又推动这些学科的发展。微分方程数值解法就主要研究如何通过离散算法将连续形式的微分方程转化为有限维问题，如代数方程组，进而来求解其近似解[1]。主要包括求解区域网格划分、离散方程的建立、方程性能分析、近似解收敛性分析等环节。微分方程数值解法在科学计算、工程技术等领域有极其广泛的应用，比如在计算物理、化学、流体力学航空航天等很多工程领域都有用到。目前已发展成为一门计算技术学科，其核心理论内容也成为高校计算数学和应用数学等专业的核心基础专业课程之一[2]。
2双语教学的必要性 双语教学主要指中英双语教学，是一种重要的教学模式，具有特殊效果和意义。 1．双语教学可丰富教学模式，转变教学理念，促进教育改革和开放。双语教学提倡用原版教材和国外的教学方式。其语言文字原汁原味，叙述合情合理，注重启发性，内容安排适合学生。这不仅使学生学到专业知识，且有助于提高英语水平，特别是专业英语阅读和写作能力。国外的教学模式以人为本，有助于转变以教师为中心、以学习知识体系为主的教育理念，促进教育改革。 2．双语教学有助于提高学生的人文素质。多学习和运用英语可以让我们发现和扬弃汉语中那些带有落后的人文价值观念和行为方式的词汇和句子，批判地接受一些思想观念和做法，使人的思维灵活有深度，个性得以发展，创新能力不断提高。大范围开展双语教学，有助于培养出具有世界主流人文素质且能很好地参与国际交流和合作的人才。 3．双语教学有助于学生以后在国内外学习、工作、考研和国际合作等带来很多方便。 微分方程数值解法既有数学上严密的逻辑性、独特的理论结构体系，又在各种工程计算中有着重要的应用，因此是联系纯数学理论和工程应用的桥梁和纽带。很多工业应用软件是利用数值方法开发成的，并且大都用英语写成。因此，有必要用双语的形式讲授这门课，让学生在学习专业知识的同时，还掌握专业英语词汇，有助于学生以后的学习和发展。从课程的体系和内容衔接上看，这门课一般安排在大学三年级。这时侯，学生对于数学分析、常微分方程、数学物理方程和计算方法等课程有了很好的基础，其中的很多概念如:导数、定积分、
偏微分方程数值解(试题)
偏微分方程数值解试题 1、考虑一维的抛物型方程： 2200, [0,], 0t T (,), (,)(,0)() x x u u x t x u x t u u x t u u x x ππνπ?==??=∈≤≤??=== （1）导出时间离散是一阶向前Euler 格式，空间离散是二阶精度的差分格式； （2）讨论（1）中导出的格式的稳定性； （3）若时间离散为二阶精度的蛙跳格式， 11 2n n n t t u u u t t +-=?-= ?? 空间离散是二阶精度的中心差分，问所导出的格式稳定吗？为什么？ 2、考虑Poission 方程 2(,)1, (,)0, in AB and AD (,)0, in BC and CD u x y x y u n u x y -?=∈Ω ?=?= 其中Ω是图1中的梯形。 使用差分方法来离散该方程。由于梯形的对称性，可以考虑梯形的一半，如图2， 图2 从物理空间到计算区域的几何变换 图1 梯形
为了求解本问题，采用如下方法：将Ω的一半投影到正方形区域?Ω ，然后在?Ω上使用差分方法来离散该方程。在计算区域?Ω 上用N N ?个网格点，空间步长为1/(1)N ξη?=?=-。 （1）引入一个映射T 将原区域Ω（带有坐标,x y ）变换到单位正方形?Ω（带有坐标,ξη）。 同时导出在新区域上的方程和边界条件。 （2）在变换区域，使用泰勒展开导出各导数项在区域内部和边界点上的差分格式。 3、对线性对流方程0 constant >0u u a a t x ??+=??，其一阶迎风有限体积法离散格式为 1?n j u +=?n j u a t x ?-?（?n j u 1?n j u --） （1）写出0a <时的一阶迎风有限体积法的离散格式； （2）写出a 为任意符号的常数的一阶迎风有限体积法的守恒形式。 （3）使用0 u u u t x ??+=??说明一阶迎风有限体积法不是熵保持的格式。 4、对一维Poission 方程 , (0,1) (0)(1)0 x xx u xe x u u ?-=∈? ==? 将[]01，分成(1)n +等分，写出用中心差分离散上述方程的差分格式，并问： （1）该差分格式与原微分方程相容吗？为什么？ （2）该差分格式稳定吗？为什么？ （3）该差分格式是否收敛到原微分方程的解？为什么？ （4）取(1)6n +=，写出该差分格式的矩阵表示。 5、叙述二重网格方法的执行过程，并对一维常微分方程边值问题 2 25, (0,1) (0)(1)0 xx u x x x u u πππ?-=∈? ==?（sin(5)+9sin(15)） 给出限制算子和延拓算子矩阵（以细网格h ：7n =，粗网格2h ：3n =为例）。 6、对一阶波动方程 01(,0)sin(), (0,1)2(0,)(1,)u u t x u x x x u t u t π???+=???? ? =∈?? =??? （1）写出用中心差分进行空间离散，用一阶向后Euler 进行时间离散的差分格式；
微分方程的分类及其数值解法
微分方程的分类及其数值解法 微分方程的分类： 含有未知函数的导数，如dy/dx=2x 、ds/dt=0.4都是微分方程。 一般的凡是表示未知函数、未知函数的导数与自变量之间的关系的方程，叫做微分方程。未知函数是一元函数的，叫常微分方程；未知函数是多元函数的叫做偏微分方程。微分方程有时也简称方程。 一、常微分方程的数值解法： 1、Euler 法： 00d (,), (1.1)d (), (1.2) y f x y x y x y ?=???=? 001 (),(,),0,1,,1n n n n y y x y y hf x y n N +=??=+=-? (1.4) 其中0,n b a x x nh h N -=+=. 用(1.4)求解(1.1)的方法称为Euler 方法。 后退Euler 公式???+==+++),,(),(111 00n n n n y x hf y y x y y 梯形方法公式 )].,(),([2 111+++++=n n n n n n y x f y x f h y y 改进的Euler 方法11(,),(,),1().2p n n n c n n p n p c y y hf x y y y hf x y y y y ++?=+??=+???=+??? 2、Runge-Kutta 方法： p 阶方法 ： 1()O h -=?总体截断误差局部截断误差 二阶Runge-Kutta 方法 ??? ????++==++=+),,(),,(,2212 1211hk y h x f k y x f k k h k h y y n n n n n n
偏微分方程数值解(试题)
1 / 7 偏微分方程数值解试题 1、考虑一维的抛物型方程： 2200, [0, ], 0t T (,), (,)(,0)() x x u u x t x u x t u u x t u u x x ππνπ?==??=∈≤≤??=== （1）导出时间离散是一阶向前Euler 格式，空间离散是二阶精度的差分格式； （2）讨论（1）中导出的格式的稳定性； （3）若时间离散为二阶精度的蛙跳格式， 11 2n n n t t u u u t t +-=?-=?? 空间离散是二阶精度的中心差分，问所导出的格式稳定吗？为什么？ 2、考虑Poission 方程 2(,)1, (,)0, in AB and AD (,)0, in BC and CD u x y x y u n u x y -?=∈Ω ?=?= 其中Ω是图1中的梯形。 使用差分方法来离散该方程。由于梯形的对称性，可以考虑梯形的一半，如图2， 图2 从物理空间到计算区域的几何变换 图1 梯形
2 / 7 为了求解本问题，采用如下方法：将Ω的一半投影到正方形区域?Ω ，然后在?Ω上使用差分方法来离散该方程。在计算区域?Ω 上用N N ?个网格点，空间步长为1/(1)N ξη?=?=-。 （1）引入一个映射T 将原区域Ω（带有坐标,x y ）变换到单位正方形?Ω（带有坐标,ξη）。 同时导出在新区域上的方程和边界条件。 （2）在变换区域，使用泰勒展开导出各导数项在区域内部和边界点上的差分格式。 3、对线性对流方程 0 constant >0u u a a t x ??+=??，其一阶迎风有限体积法离散格式为 1?n j u +=?n j u a t x ?-?（?n j u 1?n j u --） （1）写出0a <时的一阶迎风有限体积法的离散格式； （2）写出a 为任意符号的常数的一阶迎风有限体积法的守恒形式。 （3）使用0 u u u t x ??+=??说明一阶迎风有限体积法不是熵保持的格式。 4、对一维Poission 方程 , (0,1)(0)(1)0 x xx u xe x u u ?-=∈?==? 将[]01，分成(1)n +等分，写出用中心差分离散上述方程的差分格式，并问： （1）该差分格式与原微分方程相容吗？为什么？ （2）该差分格式稳定吗？为什么？ （3）该差分格式是否收敛到原微分方程的解？为什么？ （4）取(1)6n +=，写出该差分格式的矩阵表示。 5、叙述二重网格方法的执行过程，并对一维常微分方程边值问题 225, (0,1)(0)(1)0 xx u x x x u u πππ?-=∈?==?（sin(5)+9sin(15)） 给出限制算子和延拓算子矩阵（以细网格h ：7n =，粗网格2h ：3n =为例）。 6、对一阶波动方程 1(,0)sin(), (0,1)2(0,)(1,)u u t x u x x x u t u t π???+=?????=∈??=??? （1）写出用中心差分进行空间离散，用一阶向后Euler 进行时间离散的差分格式；
数值分析试题及答案.
一、单项选择题（每小题3分，共15分） 1. 和分别作为π的近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式 ()()2 1 121 1()(2)636f x dx f Af f ≈ ++? ，则A ＝（ ） A ． 16 B ．13 C ．12 D ．2 3 3. 通过点 ()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足（ ） A ． ()00l x ＝0， ()110 l x = B ． () 00l x ＝0， ()111 l x = C ． () 00l x ＝1， ()111 l x = D ． () 00l x ＝1， ()111 l x = 4. 设求方程 ()0 f x =的根的牛顿法收敛，则它具有（ ）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组12312312 20223332 x x x x x x x x ++=?? ++=??--=? 作第一次消元后得到的第3个方程（ ）. A ．232 x x -+= B ．232 1.5 3.5 x x -+= C ． 2323 x x -+= D ． 230.5 1.5 x x -=- 单项选择题答案
二、填空题（每小题3分，共15分） 1. 设T X )4,3,2(-=, 则=1||||X ，2||||X = . 2. 一阶均差 ()01,f x x = 3. 已知3n =时，科茨系数 ()()() 33301213,88C C C === ，那么() 3 3C = 4. 因为方程 ()420 x f x x =-+=在区间 []1,2上满 足 ，所以 ()0 f x =在区间内有根。 5. 取步长0.1h =，用欧拉法解初值问题 ()211y y y x y ?'=+?? ?=? 的计算公 式 . 填空题答案
数值分析试题及答案
数值分析试题及答案 一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为的近似数具有（）和（）位有效数字. A．4和3 B．3和2 C．3和4 D．4和4 2. 已知求积公式，则＝（） A． B．C．D． 3. 通过点的拉格朗日插值基函数满足（） A．＝0，B．＝0， C．＝1，D．＝1， 4. 设求方程的根的牛顿法收敛，则它具有（）敛速。 A．超线性B．平方C．线性D．三次 5. 用列主元消元法解线性方程组作第一次消元后得到的第3个方程（）. A．B． C．D． 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得分评卷 人 二、填空题（每小题3分，共15分） 1. 设, 则， . 2. 一阶均差 3. 已知时，科茨系数，那么 4. 因为方程在区间上满足，所以在区间内有根。 5. 取步长，用欧拉法解初值问题的计算公式.填空题答案
1. 9和 2. 3. 4. 5. 得分评卷 人 三、计算题（每题15分，共60分） 1. 已知函数的一组数据：求分段线性插值函数，并计算的近似值. 计算题1.答案 1. 解， ， 所以分段线性插值函数为 2. 已知线性方程组 （1）写出雅可比迭代公式、高斯－塞德尔迭代公式； （2）对于初始值，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算（保留小数点后五位数字）. 计算题2.答案 1.解原方程组同解变形为 雅可比迭代公式为 高斯－塞德尔迭代法公式 用雅可比迭代公式得 用高斯－塞德尔迭代公式得 3. 用牛顿法求方程在之间的近似根 （1）请指出为什么初值应取2？ （2）请用牛顿法求出近似根，精确到0.0001. 计算题3.答案
偏微分方程数值解试题参考答案
偏微分方程数值解 一（10分）、设矩阵A 对称正定，定义)(),(),(2 1 )(n R x x b x Ax x J ∈-= ，证明下列两个问题等价：(1)求n R x ∈0使)(min )(0x J x J n R x ∈=;(2)求下列方程组的解：b Ax = 解: 设n R x ∈0是)(x J 的最小值点,对于任意的n R x ∈,令 ),(2 ),()()()(2 000x Ax x b Ax x J x x J λλλλ?+ -+=+=, (3分) 因此0=λ是)(λ?的极小值点,0)0('=?,即对于任意的n R x ∈,0),(0=-x b Ax ,特别取b Ax x -=0,则有0||||),(2000=-=--b Ax b Ax b Ax ,得到b Ax =0. (3分) 反之,若 n R x ∈0满足 b Ax =0,则对于任意的x , )(),(2 1 )0()1()(00x J x Ax x x J >+ ==+??,因此0x 是)(x J 的最小值点. (4分) 评分标准:)(λ?的表示式3分, 每问3分,推理逻辑性1分 二（10分）、对于两点边值问题：????? ==∈=+-=0 )(,0)() ,()(b u a u b a x f qu dx du p dx d Lu 其中]),([,0]),,([,0)(min )(]),,([0min ] ,[1b a H f q b a C q p x p x p b a C p b a x ∈≥∈>=≥∈∈ 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的Ritz 形式和 Galerkin 形式的变分方程。 解: 设}0)()(),,(|{11 0==∈=b u a u b a H u u H 为求解函数空间,检验函数空间.取),(10b a H v ∈,乘方程两端,积分应用分部积分得到 (3分) )().(),(v f fvdx dx quv dx dv dx du p v u a b a b a ==+=??,),(10 b a H v ∈? 即变分问题的Galerkin 形式. (3分)
微分方程数值解习题(李立康)
常微分方程习题 《李立康》 习题 1.用Euler 方法求初值问题 ? ? ?=-='0)0(21u tu u 在1=t 时的近似解（取4 1= h ）。 2.初值问题 1 3 00 u u u()??'=? ?=? 有解32 23/u(t )t ?? = ? ?? 。但若用Euler 方法求解，对一切N T ,和H T h = ,都只能得到N t u t ,...,2,1,0==，试解释此现象产生的原因。 3.用Euler 方法计算 ?? ?=='1 )0(u u u 在1=t 处的值，取16 1 和41= h ，将计算结果与精确值e =）1（u 相比较。 4.设),(u t f 满足定理2.1的条件，对改进Euler 法（2.10）式证明： （1）其局部截断误差为)()(12 43 h O t u h -'''- ； （2）当1 ?? ?=='1 )0(u u u 计算公式 m m h h u ??? ? ??-+=22 取4 1 = h 计算)1(u 的近似值，并与习题3的结果比较。 6.就初值问题 ?? ?=+='0 )0(u b at u 分别导出用Euler 方法和改进Euler 法求近似解的表达式，并与真解 bt t a u += 22 相比较。 7.证明改进Euler 法的绝对稳定区域是整个左半平面0)Re( 包括基本概念，差分格式的构造、截断误差和稳定性，这些内容是贯穿整个教材的主线。解答问题关键在过程，能够显示出你已经掌握了书上的内容，知道了解题方法。这次考试题目的类型：20分的选择题，主要是基本概念的理解，后面有五个大题，包括差分格式的构造、截断误差和稳定性。 习题一 1． 略 2． y y x f -=),(，梯形公式：n n n n n n y h h y y y h y y )121(),(2111+-+=+- =+++，所以0122)1(01])121[()121()121(y h h y h h y h h y h h n h h n n n +--+--+-+=+-+==+-+= ，当0→h 时， x n e y -→。 同理可以证明预报-校正法收敛到微分方程的解. 3． 局部截断误差的推导同欧拉公式; 整体截断误差: ? ++++++-++≤1 ),())(,(11111n n x x n n n n n n n dx y x f x y x f R εε 11)(++-++≤n n n y x y Lh R ε，这里R R n ≤ 而111)(+++-=n n n y x y ε，所以 R Lh n n += -+εε1)1(，不妨设1 x ∈R n 2 ( Ax, x) , J ( x + x) = ? (1) = ? (0) + ( Ax, x) > J ( x ) ,因此 x 是 J ( x ) 的最小值点. (4 分) 2 二（10 分）、对于两点边值问题： ? dx dx a(u , v) = ?b ( p . + q u v)dx = ?b fvdx = f (v) , ? v ∈ H 1 (a , b ) dx dx a a 偏微分方程数值解 一（10 分）、设矩阵 A 对称正定，定义 J ( x ) = 1 ( Ax , x ) - (b , x ) ( x ∈ R n ) ，证明下 2 列两个问题等价：(1)求 x ∈ R n 使 J ( x ) = min J ( x ) ;(2)求下列方程组的解：Ax = b 解: 设 x ∈ R n 是 J ( x ) 的最小值点,对于任意的 x ∈ R n ,令 ?(λ) = J ( x + λx) = J ( x ) + λ( Ax - b , x) + λ2 (3 分) 因此 λ = 0 是 ?(λ) 的极小值点 , ? ' (0) = 0 ,即对于任意的 x ∈ R n , ( Ax - b , x) = 0 ,特 0 别取 x = Ax - b ,则有 ( Ax - b , Ax - b ) =|| Ax - b || 2 = 0 ,得到 Ax = b . (3 分) 0 0 反 之 , 若 x ∈ R n 满 足 Ax = b , 则 对 于 任 意 的 x , 1 0 0 0 评分标准: ?(λ) 的表示式 3 分, 每问 3 分,推理逻辑性 1 分 ? d du ?Lu = - ( p ) + qu = f x ∈ (a, b ) ?? u (a) = 0, u (b ) = 0 其中 p ∈ C 1 ([a , b ]), p ( x ) ≥ min p ( x ) = p x ∈[a,b ] min > 0, q ∈ C ([a , b ]), q ≥ 0, f ∈ H 0 ([a , b ]) 建立与上述两点边值问题等价的变分问题的两种形式：求泛函极小的 Ritz 形式和 Galerkin 形式的变分方程。 解 : 设 H 1 = {u | u ∈ H 1 (a , b ), u (a ) = u (b ) = 0} 为求解函数空间 , 检验函数空间 . 取 v ∈ H 1 (a, b ) ,乘方程两端,积分应用分部积分得到 (3 分) du dv 即变分问题的 Galerkin 形式. (3 分) 微分方程数值解―― 第二章 习题 1． 设)('x f 为)(x f 的一阶广义导数，试用类似办法定义)(x f 的k 阶广义导数) () (x f k （ ,2,1=k ）。 解：对一维情形，函数的广义导数是通过分部积分来定义的。 我们知，)(x f 的一阶广义导数位)(x g ，如果满足 dx x x f dx x x g b a b a )()()()('?? -=?? 类似的，)(x f 的k 阶广义导数为)()() (x f x g k =，如果有 dx x x f dx x x g b a k k b a )()()1()()()(?? -=?? 2． 试建立与边值问题 ?????====<<=+=) 2.1(0)()(,0)()() 1.1(,''44b u b u a u a u b x a f u dx u d Lu 等价的变分问题。 证明： 设}0)()(,0)()(),(|{' '2====∈=b v a v b v a v I H v v V 对方程)1.1(两边同乘以v ，再关于x 在),(b a 上积分)(V v ∈，得 ??=+b a b a fvdx vdx u dx u d )(44 其中 dx dx dv dx u d dx dx dv dx u d dx u d v dx u d d v vdx dx u d b a b a b a b a b a ???? -=-==33 33333344|)( dx dx v d dx u d dx dv dx u d dx u d d dx dv b a b a b a ??+-=-=22222222|)( dx dx v d dx u d b a ? = 2 222 (*) 记dx uv dx v d dx u d v u a b a ?+=)(),(2 222，?=b a fvdx v f ),(。于是我们得到以下等价变分问题的提法：微分方程数值解法答案
偏微分方程数值解试题参考答案
微分方程数值解――