微分方程数值解-课件

微分方程数值解

12024/4/3微分方程数值解12024/4/2微分方程数值解今天的数值计算方法，无论从形式到内容，还是从工具到效果，已远非半世纪前VonNeumann、Lax等先驱们所处的环境和条件了，计算机技术和应用软件的发展，让计算数学展开了双翼。许多迅速发展的其他学科和社会进步给计算数学的发展开拓出更为广阔的新天地。2微分方程数值解今天的数值计算方法，无论从形式到内容，还是从工精品资料3精品资料3你怎么称呼老师？如果老师最后没有总结一节课的重点的难点，你是否会认为老师的教学方法需要改进？你所经历的课堂，是讲座式还是讨论式？教师的教鞭“不怕太阳晒，也不怕那风雨狂，只怕先生骂我笨，没有学问无颜见爹娘……”“太阳当空照，花儿对我笑，小鸟说早早早……”44随着计算机软件硬件的不断更新和计算方法的迅速发展，科学计算与实验以及理论研究成为现代科学研究的三大主要手段。科学计算还能解决实验及理论无法解决的问题，并由此发现一些新的物理现象，加深人们对物理机理的理解和认识，促进科学的发展。作为三种科学研究手段之一的科学计算是一门工具性、方法性、边缘性的新学科，发展迅速，它的物质基础是计算机(包括其软硬件系统)，其理论基础主要是计算数学。5随着计算机软件硬件的不断更新和计算方法的迅速发展，科学计算与微分方程数值解26微分方程数值解2677微分方程数值解28微分方程数值解28微分方程数值解计算数学的发展与科学工程计算是紧密相联的，计算数学的发展历史也就是与其他学科结合，利用计算机不断形成新的理论及数值方法并不断形成新的学科的历史，例如：“计算物理”。9微分方程数值解计算数学的发展与科学工程计算是紧密相联的，计算微分方程数值解10微分方程数值解10微分方程数值解计算物理的物质基础是计算机；计算物理的关键技术是“计算方法”和“程序设计”;计算物理发展的原始动力是美国核武器研制的刺激。11微分方程数值解11微分方程数值解12微分方程数值解12微分方程数值解13微分方程数值解13微分方程数值解14微分方程数值解14微分方程数值解15微分方程数值解15微分方程数值解16微分方程数值解16微分方程数值解17微分方程数值解17微分方程数值解18微分方程数值解18微分方程数值解19微分方程数值解19微分方程数值解20微分方程数值解20微分方程数值解21微分方程数值解21微分方程数值解22微分方程数值解22微分方程数值解23微分方程数值解23微分方程数值解24微分方程数值解24微分方程数值解25微分方程数值解25微分方程数值解“战略计算”一词首次出现在1995年美国为了确保核库存的性能、安全性、可靠性和更新需要而实施的“加速战略计算创新（ASCI）计划”。26微分方程数值解“战略计算”一词首次出现在1995年美国为了确微分方程数值解27微分方程数值解27微分方程数值解28微分方程数值解28微分方程数值解29微分方程数值解29微分方程数值解30微分方程数值解30微分方程数值解31微分方程数值解31微分方程数值解32微分方程数值解32

提问：数值计算方法是做什么用的？研究对象：数值问题——有限个输入数据（问题的自变量、原始数据）与有限个输出数据（待求解数据）之间函数关系的一个明确无歧义的描述。如一阶微分方程初值问题求函数解析表达式数学问题求函数在某些点的近似函数值数值问题33提问：数值计算方法是做什么用的？研究对象：数值问程序设计上机计算设计高效、可靠的数值方法数值问题求解近似结果输出重点讨论数值问题的来源：实际问题建立数学模型数值问题34程序上机设计高效、可数值求解近似结果输出重点讨论数值问题的来数值方法的设计原则收敛性：方法的可行性稳定性：初始数据等产生的误差对结果的影响便于编程实现：逻辑复杂度要小计算量要小：时间复杂度要小，运行时间要短存贮量要尽量小：空间复杂度要小可靠性分析计算复杂性误差估计：运算结果不能产生太大的偏差且能够控制误差35数值方法的设计原则收敛性：方法的可行性稳定性：初始数据等产生

误差

/*Error*/一、误差的来源与分类

/*Source&Classification*/

1、从实际问题中抽象出数学模型

——模型误差

/*ModelingError*/

2、通过观测得到模型中某些参数（或物理量）的值

——观测误差

/*MeasurementError*/

3、数学模型与数值算法之间的误差求近似解

——方法误差

(截断误差

/*TruncationError*/)

4、由于机器字长有限，原始数据和计算过程会产生新的误差

——舍入误差

/*RoundoffError*/36误差/*Error*/一、误差的来二、误差分析的基本概念

/*BasicConcepts*/设为真值（精确值），为的一个近似值称为近似值的绝对误差，简称误差。

注：

误差可正可负，常常是无限位的

绝对误差限/*accuracy*/——绝对值的上界如：

绝对误差还不能完全表示近似值的好坏（绝对误差/*absoluteerror*/）37二、误差分析的基本概念/*BasicConcept近似值的误差与准确值的比值：称为近似值的相对误差，记作

注：

实际计算时，相对误差通常取因为（相对误差/*relativeerror*/）38近似值的误差与准确值的比一个算法如果输入数据有扰动（即误差），而计算过程中舍入误差不增长,则称此算法是数值稳定的，否则此算法就称为不稳定的。（数值稳定性/*NumericalStability*/）对数学问题本身如果输入数据有微小扰动，引起输出数据（即问题真解）的很大扰动，这就是病态问题。（病态问题/*ill-posedproblem*/）

它是数学问题本身性质所决定的，与算法无关，也就是说对病态问题，用任何算法（或方法）直接计算都将产生不稳定性。

39一个算法如果输入数据有扰动（即误差），而计算过程中舍此公式精确成立记为则初始误差????!!!Whathappened?!例计算

公式一：40此公式精确成立记为则初始误差????!!!Whatha考察第n步的误差我们有责任改变。造成这种情况的是不稳定的算法

/*unstablealgorithm*/迅速积累，误差呈递增趋势。初始的小扰动

公式二：注意此公式与公式一在理论上等价。方法：先估计一个IN

,再反推要求的In(n<<N)。可取41考察第n步的误差我们有责任改变。造成这种情况的是不稳定的算法取42取42考察反推一步的误差：以此类推，对n<N

有：误差逐步递减,这样的算法称为稳定的算法

/*stablealgorithm*/

在我们今后的讨论中，误差将不可回避，算法的稳定性将会是一个非常重要的话题。43考察反推一步的误差：以此类推，对n<N有：误差逐步递例：蝴蝶效应

——纽约的一只蝴蝶翅膀一拍，风和日丽的北京就刮起台风来了？！纽约北京这是一个病态问题44例：蝴蝶效应——纽约的一只蝴蝶翅膀一拍，风和日蝴蝶效应在社会学界用来说明：一个坏的微小的机制，如果不加以及时地引导、调节，会给社会带来非常大的危害，戏称为“龙卷风”或“风暴”；一个好的微小的机制，只要正确指引，经过一段时间的努力，将会产生轰动效应，称为“革命”。蝴蝶效应是气象学家洛伦兹1963年提出来的。其大意为：一只南美洲亚马孙河流域热带雨林中的蝴蝶，偶尔扇动几下翅膀，可能在两周后引起美国德克萨斯引起一场龙卷风。45蝴蝶效应在社会学界用来说明：一个坏的微小的机制，如果不加以几点注意事项

/*Remarks*/1、

避免相近二数相减例：a1=0.12345，a2=0.12346，各有5位有效数字。而a2

a1=0.00001，只剩下1位有效数字。

几种经验性避免方法：当|x|<<1时：46几点注意事项/*Remarks*/1、避免相近二数2、避免小分母:分母小会造成浮点溢出

/*overflow*/3、避免大数吃小数例：用单精度计算的根。精确解为

算法1：利用求根公式在计算机内，109存为0.11010，1存为0.1101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1的指数部分须变为1010，则：1=0.011010，取单精度时就成为：

109+1=0.100000001010+0.000000001010=0.100000001010大数吃小数472、避免小分母:分母小会造成浮点溢出/*over算法2：先解出再利用注：求和时从小到大相加，可使和的误差减小。例：按从小到大、以及从大到小的顺序分别计算1+2+3+…+40+1094、

先化简再计算，减少步骤，避免误差积累。一般来说，计算机处理下列运算的速度为5、选用稳定的算法。48算法2：先解出注：求和时从小到大相加，可使和的误差减小。例：微分方程数值解49微分方程数值解491.微分方程及其数值解计算机解决实际问题的步骤建立数学模型选择数值方法编写程序上机计算501.微分方程及其数值解计算机解决实际问题的步骤建立数学模型现实世界中绝大多数事物的内外联系是及其复杂的，其状态随着时间、地点、条件的不同而不同，我们只能通过对问题进行简化和作某些假定，从中找出其状态和状态的变化规律之间的关系，也即一个或一些函数与它们的导数之间的关系，这种关系的数学表达就是微分方程。

51现实世界中绝大多数事物的内外联系是及其复杂的，其状态随着时间偏微分方程数值解主要是有限差分法和有限元法。

52偏微分方程数值解主要是有限差分法和有限元法。52偏微分方程发展史：(1)十八世纪初,Taylor:(2)十九世纪中期，三类偏微分方程：

(3)十九世纪末到二十世纪初，其它方程：

高阶方程：KDV方程：53偏微分方程发展史：(1)十八世纪初,Taylor:(2如果能找到一个（或一族）具有所要求阶连续导数的解析函数，将它代入微分方程（组）中，恰好使得方程（组）的所有条件都得到满足，我们就将它称为这个方程（组）的解析解（也称古典解）。“微分方程的真解”或“微分方程的解”就是指解析解。寻找解析解的过程称为求解微分方程。

微分方程的解在数学意义上的存在性可以在非常一般的条件下得到证明，这已有许多重要的结论。但从实际上讲，人们需要并不是解在数学中的存在性，而是关心某个定义范围内，对应某些特定的自变量的解的取值或是近似值-这样一组数值称为这个微分方程在该范围内的数值解，寻找数值解的过程称为数值求解微分方程。

54如果能找到一个（或一族）具有所要求阶连续导数的解析函数，将它为什么要研究数值求解方法呢？

1）在实际问题中我们所能获取的或感兴趣的，往往只是一个特定点上的数据。如空间的温度分布只能一个点一个点地测定，火箭升空传回的控制信息只能以某个确定的时间为间隔，一个个地发送和接受，如此等等。这些离散点上的函数值对于解决实际问题，已经足够了，寻找解析解的一般形式未必必要。

55为什么要研究数值求解方法呢？1）在实际问题中我们所能获取的2）在很多情况下，寻找解析解也并无可能。现实问题中归结的微分方程不满足解析解的存在条件的比比皆是，方程中出现的有些函数连续性都无法保证，它们并不存在前述意义的解析解。于是，求数值解便成了在这种情况下解决问题的重要手段了。

562）在很多情况下，寻找解析解也并无可能。现实问题中归结的微分3）即使微分方程的解析解存在，以并不意味可以将它表示为初等函数，如多项式、对数函数、指数函数三角函数及它们的不定积分的有限组合形式——显式解。

事实上，有显式解的微分方程只占解析解存在的微分方程中的非常小的一部分。

573）即使微分方程的解析解存在，以并不意味可以将它表示为初等函5858由于数字电子计算机只能存储有限个数据和作有限次运算，所以任何一种适用于计算机解题的方法，都必须把连续问题离散化，最终化成有限形式的代数方程组。59由于数字电子计算机只能存储有限个数据和作有限次运算，所以任何微分方程数值解1）区域剖分把整个定义域分成若干个小块，以便对每小块上的点或片求出近似值，这样按一定规律对定义域分切的过程称为区域剖分。2）微分方程的离散区域剖分完毕后，依据原来的微分方程去形成关于这些离散点或片的函数值的递推公式或方程。这是它们的未知量已不是一个连续函数，而成了若干个离散的未知值的某种组合了，这个步骤称为微分方程离散。60微分方程数值解1）区域剖分2）微分方程的离散603）初始和边界条件处理离散后系统是一个递推公式，那它需要若干个初值才能启动。若是一个方程组，那它所含的方程个数一般少于未知量的个数，要想求解还需要补充若干个方程。这些需要补充的初值和方程往往可以通过微分方程的初始条件和边界条件来得到，这就是初始和边界条件处理过程。613）初始和边界条件处理614）离散系统的性态研究我们主要研究：这个系统是否可解，即解的存在性、唯一性问题；它与精确解的差距有多大，这个差距当区域剖分的尺寸趋于零时，是否也会趋于零，趋于零的速度多快，即解的收敛性和收敛速度问题；当外界对数据有所干扰时，所得的解是否会严重背离离散系统的固有的解，即解的稳定性问题。624）离散系统的性态研究62上述问题说道底是一个误差分析问题，因为如果从实际问题到得出数值解的每一步都没有任何误差的话（当然，这是不可能的），那么数值解就应该是离散点上的精确值，也就不用煞费苦心去讨论上面的问题了。63上述问题说道底是一个误差分析问题，因为如果从实际问题到得出数数值求解微分方程过程示意微分方程区域剖分离散系统的性态研究递推计算或解线性代数方程组微分方程离散初始和边界条件处理解的存在性、唯一性解的收敛性和收敛速度解的稳定性得到数值解64数值求解微分方程过程示意微分方程区域剖分离散系统的递推计算微分方程数值解现实问题数学模型离散格式模型误差建模离散舍入误差观测模型截断误差数值解计算65微分方程数值解现实问题数学模型离散格式模型误差建模离散舍入误微分方程数值解数值线性代数又称矩阵计算，它是科学与工程计算的核心。可以毫不夸张地讲，大部分科学与工程问题最终都要归结为一个矩阵计算问题，其中具有挑战性的问题是大规模矩阵计算问题。66微分方程数值解数值线性代数又称矩阵计算，它是科学与工程计算的微分方程数值解（1）求解线性方程组的问题。（2）矩阵特征值问题。即给定一个矩阵，求它的部分或全部特征值以及对应的特征向量（3）线性最小二乘问题。67微分方程数值解（1）求解线性方程组的问题。（2）矩阵特征值问微分方程数值解若