资源描述:
《科学计算复习要点.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数值分析复习要点引论1数值计算研究的对象与特点计算方法研究的对象是专门研究各种数学问题的计算机解法(数值解法),包括方法的构造和求解过程的理论分析及软件实现,包括方法的收敛性、稳定性以及误差分析等.计算方法即具有纯数学的抽象性与严密性的特点,又具有应用的广泛性与实验的技术性特点.2误差的概念2.1误差的来源模型误差:数学模型的解与实际问题的解之间出现的误差,称为模型误差.测量误差:在测量具体数据时产生的误差称为测量误差.截断误差:数学模型的准确解与数值方法的准确解之间的误差称为截断误差.舍入误差:由于计算机字长的限制而产生的误
2、差,称为舍入误差.2.2误差的度量(1).绝对误差与绝对误差限(2).相对误差与相对误差限(3).有效数字2.3误差的传播和、差的误差限不超过各误差限的和.积、商的相对误差限不超过各相对误差限的和.3数值计算的若干原则避免两相近数相减和绝对值太小的除数、简化计算步骤、使用数值稳定的算法方程求根1二分法用二分法求方程的实根的近似值,其主要思想是:将含有根的隔离区间二分,通过判断二分点与边界点函数值的符号,逐步对半缩小隔离区间,直到缩小到满足精度要求为止,然后取最后二分区间的中点为根的近似值.2迭代法一般地,为了求一元非线性方程的
3、根,可以先将其转换为如下的等价形式然后构造迭代公式.3收敛性和收敛速度(收敛性基本定理)的条件和结论收敛速度的快慢可用收敛阶来衡量.(收敛阶)设序列收敛到,并记误差.若存在常数和,使得:则称序列是阶收敛的,当时,称为线性收敛,当时,称为超线性收敛,当时,称为二次收敛或平方收敛.4牛顿迭代公式及其收敛性牛顿迭代公式牛顿法的收敛性设是方程的单根,并且在的邻域上连续,则牛顿迭代法(3.4.1)至少平方局部收敛.解线性方程组的直接法1高斯消去法消元过程为:对逐次计算:回代过程:逐步回代求得原方程组的解高斯消去法的乘除法总计算量为:2高
4、斯—约当消去法约当消去法的计算过程为:对于计算:乘除法的总次数为:.它比高斯消去法的计算量大,但不需要回代过程3向量和矩阵的范数、条件数向量范数:1范数2范数范数矩阵的范数设为维向量,为阶方阵,则算子范数:称为矩阵的行范数。称为矩阵的列范数。设为阶可逆矩阵,则称数为条件数:,分别称为的条件数,条件数,条件数解线性方程组的迭代法1雅克比迭代法的迭代公式:矩阵形式:,2高斯—赛德尔迭代法迭代公式为:,成矩阵形式 , 3迭代法的收敛性判断(迭代法收敛的基本定理)设有阶方程组,对于任意初始向量和右端项,迭代法收敛的充分必要条件是迭
5、代矩阵的谱半径(迭代法收敛的充分条件)若,则由迭代公式(5.1.3)所产生的向量序列收敛于方程组的精确解,且有误差估计式,(充分条件)若线性方程组的系数矩阵为严格对角占优或不可约弱对角占优矩阵,则雅克比和高斯—赛德尔迭代法收敛。函数插值1插值的基本概念包括线性插值、抛物插值和多项式插值的存在惟一性。2拉格朗日插值3插值余项与误差估计若在上的插值多项式为,则称为的插值余项(也称误差)。设在上的阶导数连续,记为且在互异节点的函数值为。若满足插值条件的插值多项式为,则对有:其中,4牛顿插值数值积分1代数精度的概念及其求法。若数值求积
6、公式对被积函数都能精确成立,而对被积函数不能精确成立,则称求积公式具有次代数精度。2牛顿-柯特斯公式梯形求积公式抛物线求积公式或辛普生求积公式梯形公式的截断误差,3复合梯形求积公式将区间等分,记分点为并在每个小区间上应用梯形公式得:复合梯形公式的截断误差,4复合辛普生求积公式在每个小区间上,用辛普生公式得:其中为的中点,即5高斯求积公式若有一组节点,使插值型求积公式(8.5.1)具有次代数精度,则称此组节点为高斯点,并称相应的求积公式为高斯型求积公式。常微分方程初值问题的数值解法1欧拉公式包括显式、隐式、两步、改进的欧拉公式和
7、梯形公式。欧拉公式隐式欧拉公式为梯形公式改进的欧拉公式两步欧拉公式2单步法的局部截断误差和方法的阶设是微分方程的精确解,则称为单步法的局部截断误差。如果求微分方程数值方法的局部截断误差是,其中为整数,则称该方法是阶的,或该方法具有阶精度。越大,方法的精度越高。含的项,称为该方法的局部截断误差主项。欧拉公式是一阶方法,其截断误差主项为。隐式欧拉公式也是一阶方法,它的主项是梯形公式是二阶方法,其局部截断误差为可以证明,改进的欧拉方法也为二阶方法。3四阶龙格—库塔方法4阶经典方法形式为:4单步法的收敛性和稳定性若求微分方程的一种数值
8、方法对于任意固定的,当(同时)时,有则称该方法是收敛的。Euler方法是收敛的梯形公式是收敛的改进的方法也是收敛的Y67若用某一数值方法计算时,所得到的实际计算结果为,且由扰动引起以后各节点的扰动为,如果总有则称该方法是稳定的。单步法的稳定区间方法稳定区间