数值计算复习要点

1、数值分析复习要点引论1数值计算研究的对象与特点计算方法研究的对象是专门研究各种数学问题的计算机解法(数值解法)，包括方法的构造和求解过程的理论分析及软件实现，包括方法的收敛性、稳定性以及误差分析等计算方法即具有纯数学的抽象性与严密性的特点，又具有应用的广泛性与实验的技术性特点.2误差的概念2.1 误差的来源模型误差：数学模型的解与实际问题的解之间出现的误差，称为模型误差.测量误差：在测量具体数据时产生的误差称为测量误差截断误差：数学模型的准确解与数值方法白准确解之间的误差称为截断误差.舍入误差：由于计算机字长的限制而产生的误差，称为舍入误差.2.2 误差的度量(1) .绝对误差与绝对误差限(2

2、) .相对误差与相对误差限(3) .有效数字2.3误差的传播和、差的误差限不超过各误差限的和.积、商的相对误差限不超过各相对误差限的和3数值计算的若干原则避免两相近数相减和绝对值太小的除数、简化计算步骤、使用数值稳定的算法方程求根1二分法用二分法求方程f(x)=0的实根X*的近似值，其主要思想是：将含有根X*的隔离区间二分通过判断二分点与边界点函数值的符号，逐步对半缩小隔离区间，直到缩小到满足精度要求为止，然后取最后二分区间的中点为根X*的近似值.2迭代法一般地，为了求一元非线性方程f(X)=0的根，可以先将其转换为如下的等价形式X斗X法后构造迭代公式.Xk由=%Xk)k=0,1,23收敛性和

3、收敛速度(收敛性基本定理)的条件和结论收敛速度的快慢可用收敛阶来衡量.(收敛阶)设序列tk匕收敛到X*,并记误差ek=|xkx*|.若存在常数p21和c#0,使得:lim*=c1k1jek则称序列&篇是p阶收敛的，当p=1时，称为线性收敛，当p>1时，称为超线性收敛,当p=2时，称为二次收敛或平方收敛.4牛顿迭代公式及其收敛性牛顿迭代公式Xk.1=Xk-fxklk002,2-f(Xk)牛顿法的收敛性、1*、.设X是方程f(x)=0的单根，并且f'(X)在x的邻域上连续，则牛顿迭代法(3.4.1)至少平方局部收敛.解线性方程组的直接法1高斯消去法消元过程为：又k=1,2，,

4、n1逐次计算：1=aikk9/akkn,(i=k+1n)细(k)=a产likakj",(i,j=k+1;-,n)Pi(k)=b("likbk"(i=k+1;：n)回代过程：逐步回代求得原方程组的解xnm前产八八nXk=(bk)-"akjXj)/akk),(k=n1,n-2；-,1)j生1高斯消去法的乘除法总计算量为:1n3n2n2/nn3n2n32522332高斯一约当消去法约当消去法的计算过程为:对于k=1,2，n计算：ak)二ak：Q/akk(j=k1,n1)14k)=a(3-a产Makjk)(i=1,2,n且i#k;j=k+1,k+2,n+1)乘除

5、法的总次数为1n3n222它比高斯消去法的计算量大，但不需要回代过程3向量和矩阵的范数、条件数向量范数：1范数n|x|1=xx2范数i1|凶2=(2|Xi)M=maxXi01g包i矩阵的范数设X为n维向量,A为n阶方阵,则算子范数称为矩阵A的行范数。|A|1=paxgaj称为矩阵A的列范数。1空空一设A为n阶可逆矩阵，则称数Cond(A)=|A1|A为条件数:Cond.A)二网恒,10a,Condi(A)=|A|jAL分别称为A的2条件数，1条件数，2条件数解线性方程组的迭代法1雅克比迭代法的迭代公式：n、(k-H)1(k).冲=Zaijxj+biaaiikjmj#JCond2(A)|A2AJ

6、2矩阵形式：x(k询=Bjx的+fjBj=ID,A,fj=D/b2高斯一赛德尔迭代法迭代公式为：“i_1n、Zaijxjk*-Na。x(k)十«,(iH,2,，n)j二j二中x(k1)=1aii成矩阵形式x(k1)=BG，x(k)fG3Bgc=D-L“U,fG，=(D-L)b3迭代法的收敛性判断(迭代法收敛的基本定理)设有n阶方程组x=Bx+f,对于任意初始向量x(0)和右端项f,迭代法收敛的充分必要条件是迭代矩阵的谱半径二(B)门.(迭代法收敛的充分条件)若|B|<1,则由迭代公式(5.1.3)所产生的向量序列«(k)收.、)I，一.»,一一，一*敛于万程

7、组x=Bx+f的精确解x,且有误差估计式x(k)*-x:二B(k)(k4)-1-bxx3(1)x(0)-x(充分条件)若线性方程组Ax=b的系数矩阵为严格对角占优或不可约弱对角占优矩阵，则雅克比和高斯一赛德尔迭代法收敛。函数插值1插值的基本概念包括线性插值、抛物插值和多项式插值的存在惟一性。-xj-'xj2拉格朗日插值nnxLn(x)八一i0j=0xijT3插值余项与误差估计若f(x)在a,b上的插值多项式为Ln(x),则称Rn(x)=f(x)Ln(x)为Ln(x)的插值余项(也称误差)。设f(x)在a,b上的n+1阶导数连续，记为f(x)Cn*a,b且f(x)在互异节点a4x0<

8、;xiL”<b的函数值为y0,yi，yn。若满足插值条件Ln(为)=yi(i=0,1,2,，n)的插值多项式为Ln(x),则对Vxa,b有：Rn(x)=f(x)Ln(x)(n1)!f(n1)()(n-1)!,-n1(x)n其中a<J<b,6n书(x)=口(xxj)j=04牛顿插值Nn(x)=f(x°)，f(x0,x1)(x-x。)匚户+f(x0,x，xn)(x-x0)(x-x1)一(x-xnj)数值积分1代数精度的概念及其求法若数值求积公式对被积函数f(x)=1,x,xm都能精确成立，而对被积函数f(x)=xm+不能精确成立，则称求积公式具有m次代数精度。2牛顿-柯

9、特斯公式bnI(f)=f(x)dx=(b-a)%Ci(n)f(xi)ai=0C(n)Ci(T严ni!(n-i)!0H(t-j)dtj=0j才梯形求积公式I(f)：T二吐aIf(a)-f(b)2抛物线求积公式或辛普生求积公式b-ab-aI(f):、S=f(a)4f()f(b)62_''.梯形公式的截断误差R1(f).l-L)b(x-a)(x-b)dx=a'击ea)3,”(a，b)3复合梯形求积公式将a,b反间n等分，记分点为Xi=aTh,b-a(h=,in=0,1,，n)并在每个小区间xi,xi1b1上应用梯形公式得:f(x)dxnxnh11f(x)dxXhnf(a)0f

10、(x,f(b)2一y复合梯形公式的截断误差Rn(f)=-bah2f''()"三(a,b)124复合辛普生求积公式在每个小区间如，为1上，用辛普生公式得：h|n1n1|Sn=1f(a)+4£f(x1)+2£f(xi)+f(b)6-i3日yJ其中x/为为，为扪的中点，即x=为12125高斯求积公式若有一组节点xo,x1，,xn虻-1,1,使插值型求积公式(8.5.1)具有2n+1次代数精度，则称此组节点为高斯点，并称相应的求积公式为高斯型求积公式。常微分方程初值问题的数值解法1欧拉公式包括显式、隐式、两步、改进的欧拉公式和梯形公式。欧拉公式'0

11、1=Vnhf(xn,Yn)隐式欧拉公式yn1小门,hf函.141)h为梯形公式yn1=ynf(xn,yn)，f(xn1,yn1)改进的欧拉公式yn1=yn.hf(xn/n),f(xn1必,hHxn,yn)2两步欧拉公式yn1-ynj-2hf(xn,yn)2单步法的局部截断误差和方法的阶设y(x)是微分方程的精确解，则Tn1y(xn1)-y(xn)-h,(xn，xn1,y(xn),y(xn1),h)称为单步法的局部截断误差。如果求微分方程数值方法的局部截断误差是Tn卡=O(hp*),其中p*为整数，则称该方法是p阶的，或该方法具有p阶精度。p越大，方法的精度越高。含hp卡的项，称为该方法的局部截

12、断误差主项。h2欧拉公式是一阶方法，其截断误差主项为三y"(xn)。隐式欧拉公式也是一阶方法，它的主项是-y*(xn)3梯形公式是二阶方法，其局部截断误差为-Xy'(xn).12可以证明，改进的欧拉方法也为二阶方法。3四阶龙格一库塔方法hyn1=yn(ki2k22k3-k4)ki=f(xn,yn),h.h4阶经典RK方法形式为:"k2f(xn,ynk1)22,h.hk3=f(xn+-,yn十一k2)22k4=f(xn+h,yn+hk3)4单步法的收敛性和稳定性若求微分方程的一种数值方法对于任意固定的xn=x0+nh，当hT0(同时nT)时，有ynTy(xn),则称该方法是收敛的。Euler方法是收敛的梯形公式是收敛的改进的Euler方法也是收敛的若用某一数值方法计算yn时，所得到的实际计算结果为n,且由扰动6nNyn-n|引起以后各节点ym(mn)的扰动为第，如果总有16m国加|,则称该方法是稳定的。单步法的稳定区间方法E(h九)稳定区间隐式欧拉法一二h.0梯形公式122“h12-：h<0考试题型判断题解同一个线性方程组，高斯消去法比约当消去法的计算量大。（）。选择题在牛顿柯特斯求积公式中，柯特斯系数的代数和为（）A.0B