数值计算方法复习资料

1、?数值计算方法?复习资料代'K弟一早数值计算方法与误差分析弟早非线性方程的数值解法第三章线性方程组的数值解法第四章插值与曲线拟讨第五章数值积分与数值微分第六章常微分方程的数值解法自测题课程的性质与任务数值计算方法是一门应用性很强的根底课,在学习高等数学,线性代数和算法语言的根底上,通过本课程的学习及上机实习、使学生正确理解有关的根本概念和理论,掌握常用的根本数值方法,培养应用计算机从事科学与工程计算的水平,为以后的学习及应用打下良好根底.第一章数值计算方法与误差分析一考核知识点误差的来源类型；绝对误差和绝对误差限,相对误差和相对误差限,有效数字；绝对误差的传播.二复习要求1 .知道产生

2、误差的主要来源.2 .了解绝对误差和绝对误差限、相对误差和相对误差限和有效数字等概念以及它们之间的关系.3 .知道四那么运算中的误差传播公式.三例题例1设x*=JF3.1415926近似值x=3.14=0.314X101,即m=1,它的绝对误差是0.0015926,有1-引=0削15926£05x1产即n=3,故x=3.14有3位有效数字.x=3.14准确到小数点后第2位.又近似值x=3.1416,它的绝对误差是0.0000074有=0.0000074-<0,5XW1-5即m=1,n=5,x=3.1416有5位有效数字.而近似值x=3.1415,它的绝对误差是0.0000926

3、-；有a-X*|=0.0000926-<0,5xW1-4即m=1,n=4,x=3.1415有4位有效数字.这就是说某数有s位数,假设末位数字是四舍五入得到的,那么该数有s位有效数字；例2指出以下各数具有几位有效数字,及其绝对误差限和相对误差限：2.0004-0.0020210009000.00解由于x1=2.0004=0.20004X101,它的绝对误差限0.00005=0.5X1065,即m=1,n=5,故x=2.0004有5位有效数字.a=2,相对误差限1O'-2x9x2=-0.00200,绝对误差限0.000005,由于m=-2,n=3,x2=-0.00200有3位有效数字

4、.a1=2,相对误差限多=103=0.00252x9x3=9000,绝对误差限为0.5X10°,由于m=4,n=4,x3=9000有4位有效数字,a=9,相对误差限s=1/产=0.0000562x9X4=9000.00,绝对误差限0.005,由于m=4,n=6,X4=9000.00有6位有效数字,相对误差限为专=!10U=0.000000562x9由X3与X4可以看到小数点之后的0,不是可有可无的,它是有实际意义的.例3ln2=0.69314718丁精确到10一3的近似值是多少解精确到103=0.001,意旨两个近似值Xi,x2满足卜一jJWO.OOl,由于近似值都是四舍五入得到的,

5、要求满足一/t0,001,近似值的绝对误差限应是8=0.0005,故至少要保存小数点后三位才可以.故ln2ft0.693.第二章非线性方程的数值解法一考核知识点二分法；迭代法；牛顿法；弦截法.二复习要求1 .知道有根区间概念,和方程f(x)=0在区间(a,b)有根的充分条件.2 .掌握方程求根的二分法,知道其收敛性；掌握二分法二分次数公式,掌握迭代法,知道其收敛性.3 .熟练掌握牛顿法.掌握初始值的选择条件.4 .掌握弦截法.三例题例1证实方程1xsinx=0在区间0,1内有一个根,使用二分法求误差不超过0.5X104的根要迭代多少次证实令f(x)=1xsinx,f(0)=1>0,f(1

6、)=-sin1<0f(x)=1xsinx=0在0,1有根.又f(x)=1cosx>0(xW0.1),故f(x)=0在区间0,1内有唯一实根.、一4,.给定误差限0.5X10,有-lnOJ+41nlO165n7n>11=1/ln2ln2只要取n=14.例2用迭代法求方程x54x2=0的最小正根.计算过程保存4位小数.分析容易判断1,2是方程的有根区间.假设建立迭代格式工=,即租=.同卜丁>1(工,此时迭代发散4A'A建立迭代格式x=灯4在2,w(x)=地工+21|pV)|=此时迭代收敛.解建立迭代格式x=04计2,中二即4工+2产(初二<-(xeQ,2),取

7、初始值%=115M4#+25=加瓦+2=6总14310叼=卬4万+2=历海"5051=麻+2=80204阳15165为二加角+2二将15182七二涵海二酶西阳L5185取1.5185例3用弦截法求方程x3-x2-1=0,在x=1.5附近的根.计算中保存5位小数点.分析先确定有根区间.再代公式.解f(x)=x3-x2-1,f(1)=-1,f(2)=3,有根区间取1,2.取xi=1,迭代公式为/(演)/、(n=1,2,)3(H-)二2-/1fcil.251253-1253-1运=1,25-子口Jx.25-2)总137662-1.25、125'-2'+2?x4=1.3766

8、2-1.37662-137662J-11一376螃一L37663-l25?+125x(l.37662-1.25)1,48881莅=148881-L488B1-1.488B1-11.48881s-L4888P-L376623+137662x(L4888M37662)=1.463481463481-14634E2-1%=146348J0小.1x(146348-148881)洌1465531.463483-146348'-1,48/+L488B12取/翱1.46553,f(1.46553)%0.000145例4选择填空题1 .设函数f(x)在区间a,b上连续,假设满足,那么方程f(x)=0在

9、区间a,b一定有实根.答案：f(a)f(b)<0解答：由于f(x)在区间a,b上连续,在两端点函数值异号,由连续函数的介值定理,必存在c,使得f(c)=0,故f(x)=0一定有根.2 .用简单迭代法求方程f(x)=0的实根,把方程(x)=0表成x=%x),那么f(x)=0的根是()(A)y=x与y=9(x)的交点(B)y=x与y=*x)交点的横坐标(C)y=x与x轴的交点的横坐标(D)y=4x)与x轴交点的横坐标答案：(B)解答：把f(x)=0表成x=*x),满足x=%x)的x是方程的解,它正是y=x与y=%x)的交点的横坐标.3.为求方程x3-x2-1=0在区间1.3,1.6内的一个根

10、,把方程改写成以下形式,并建立相应的迭代公式,迭代公式不收敛的是()(B)-二：.H'-1(D).:i:出+政+1答案：(A)解答：在(A)中产=占,力工,研划/二1°75瞄故迭代发散.1121在(B)中x=1+=1+,=-<-=0.901<1,故迭代收敛.xx#1.3在©中,破力二7,加攵)1=即：03,;二；他因0至15<1,故Jl1TAJ3(1T1.3J迭代收敛.在(D)中,类似证实,迭代收敛.第三章线性方程组的数值解法一、考核知识点高斯顺序消去法,列主元消去法；雅可比迭代法,高斯一一赛德尔迭代法,超松弛迭代法；消去法消元能进行到底的条件,迭

11、代解数列收敛的条件.二、复习要求1 .知道高斯消去法的根本思想,熟练掌握高斯顺序消去法和列主元消去法.2 .掌握线性方程组雅可比迭代法和高斯一一赛德尔迭代法.3 .知道解线性方程组的高斯消去法消元能进行到底的条件,知道迭代解数列收敛概念和上述两种迭代法的收敛性的充分条件.三、例题2xx+x2+4/二T例1用顺序消去法解线性方程组?3乙+2/+x3=4%+2勺+4&=-1解顺序消元4-T144-1214-1_00.5-5550152-0.5E-1月.,214-100.5-55.50017-17为+0.5/+2%=0.5于是有同解方程组：4-10马=1117%=-17I-回代得解：X3=-

12、1,X2=1,Xl=1.原线性方程组的解为X=1,1,一11.例2取初始向量XO=0,0,0T,用雅可比迭代法求解线性方程组4 +2/-2xs=14工1+K?+%=32xi+2x2+x3=5解建立迭代公式"产=-2碎+2承+1铲；-婿+穹+3k=1,2,3,铲=-2出+暇+5第1次迭代,k=0,XO=0,得到X=1,3,5t,第2次迭代,k=1,=(-2x3+2x5)+l=5q-)+3=_3,得到x=(533)tk?=2(14-3)+5=-3第3次迭代,k=2,=G2X(-3)+2x(-3)+l=l<五2=_(5+_3)+3=1,得到x=(1,1,1)T=-2(5-33+5=1

13、第4次迭代,k=3,2)=(-2xl+2xl)+l=l二一(1+1)+3=1,得到x(4)=(1,1,1)T4)=-2(1+1)+5=1例3填空选择题：1.用高斯列主元消去法解线性方程组+2x2+/=042%+2/+3xs=3作第1次消元后的第2,3个方程分别为1个方程变为：2x+2x2+3x3=3,消元得到解答1.选a21=2为主元,作行互换,第z,+-0-5及=-1.5423是应填写的内容.-2%+L5为=3.5xl+2x2-2x5=12.用高斯-赛德尔迭代法解线性方程组?工+j=3的迭代格式中2x1+2超+/=5=(k=0,1,2,)解答高斯一赛德尔迭代法就是充分利用已经得到的结果,求X

14、2的值时应该用Xi的新值.答案是：53产-2淄10±-/-3/=723.当卜卜)时,线性方程组,-再+7/+3%二&3的迭代解一定收敛.2五一4%/=9.2(A)>6(B)=6(C)<6(D)>6解答：当a|>6时,线性方程组的系数矩阵是严格对角占优矩阵,由教材第3章定理知,迭代解一定收敛.应选择(A).第四章插值与曲线拟合一考核知识点插值函数,插值多项式,被插值函数,节点；拉格朗日插值多项式：插值基函数；差商及其性质,牛顿插值多项式；分段线性插值、线性插值基函数,最小二乘法,直线拟合.二复习要求1 .了解插值函数,插值节点等概念.2 .熟练掌握拉格朗

15、日插值多项式的公式,知道拉格朗日插值多项式余项.3 .掌握牛顿插值多项式的公式,了解差商概念和性质,掌握差商表的计算,知道牛顿插值多项式的余项.4 .掌握分段线性插值的方法和线性插值基函数的构造.5 .了解曲线拟合最小二乘法的意义和推导过程,以及线性拟合和二次多项式拟合的方法,三例题例1函数y=f(x)的观察数据为Xk-2045yk51-31试构造f(x)的拉格朗日多项式Pn(x),并计算f(-1)o解先构造基函数=用0-4)(x-5)=x(x-4)(x-5)1=(x+2)(x-4)(x-5)1 (0-(-2)(0-4)(0-5)(x+2)(x4)(x-5)40j=(x+2)式工-5)2 -(

16、4+2)(4-0)(4-5)x(x+2Xx-5)24i=,+2)K工-2)("4)3(5+2)(5-0)(5-4).+2)工(工-4)35所求三次多项式为3P3(X)=2MJUO=一"r"X"5)+84(i+2)x(x-4)352)("4)("5)IcT("式x+2)(x-S)241T145x42515524P3(-1)=-h1=4214217例2函数y=f(x)的数据如表中第2,3歹U.计算它的各阶均差.解依据均差计算公式,结果列表中.kxkf(xk)一阶均差二阶均差三阶均差四阶均差00.400.4107510.550.5

17、78151.1160020.650.696751.168000.2800030.800.888111.275730.358930.1973340.901.202121.384100.433480.213000.03134计算公式为：一阶均差冗Jt一五41二阶均差/血,凝中.=三出土外生山心QL2)丸-无心工例3设工是n+1个互异的插值节点,峋值=32用是拉格朗»日插值基函数,证实：2?的巨1证实Pn(x)=y0l0(x)+yill(x)+ynln(x)=Jt-0f"团0M二任"+e/"+凡X.料,工)当砥时,=*)+»)?.附+再由于1,故有二

18、二,例4满足条件咖=加川(1)山(2"的插值多项式p(x)=解设所求的为p(x)=a0+aix+a2X2+a3X3由插值条件知p(O)=%=07/=R=Q广(1)-17?(2)二0.+2叮1+4鼻二+8a3=2解之得a2=3/2a3=-1/2所求的插值多项式为p(x)=-1/2x3+3/2x2例5选择填空题1 .通过四个互异节点的插值多项式P(x),只要满足(,那么P(x)是不超过一次的多项式.(A)初始值yo=0(B)一阶均差为0(C)二阶均差为0(D)三阶均差为0解答：由于二阶均差为0,那么牛顿插彳1多项式为N(x)=f(x0)+f(x0,xi)(xx0)它是不超过一次的多项式.

19、应选择(C)正确.2 .拉格朗日插值多项式的余项是(,牛顿插值多项式的余项是()/皿也)(A)二:''1-n-(B)f(x,x0,xi,x2,xn)(xxi)(xx2)仅一xni)(xxn)(盟+1j尸砌G)(C)八,工(D)f(x,x0,xi,x2,xn)(xXo)(XXi)(X&)仅一xni)(xxn)5+1)!解答：(A),(D).第五章数值积分与数值微分一考核知识点数值求积公式,求积节点,求积系数,代数精度；插值型求积公式,牛顿一一柯特斯求积公式,柯特斯系数及其性质,(复化)梯形求积公式,(复化)辛卜生求积公式；高斯型求积公式,高斯点,(二点、三点)高斯一一勒让

20、德求积公式；(二点、三点)插值型求导公式.复习要求1, 了解数值积分和代数精度等根本概念.2, 了解牛顿柯特斯求积公式和柯特斯系数的性质.熟练掌握并推导(复化)梯形求积公式和(复化)辛卜生求积公式.3, 知道高斯求积公式和高斯点概念.会用高斯一勒让德求积公式求定积分的近似值.4,知道插值型求导公式概念,掌握两点求导公式和三点求导公式.例1试确定求积公式JJ&哥/(-j=)解当f(x)取1,x,x,计算求积公式何时精确成立.取f(x)=1,有：左边=(2)取f(x)=x,有：左边=|xdx=0,右边=0(3)类似导出,取f(x)=x2,x3,有左边=右边(5)取f(x)=x4,有：左边=

21、2/5,右边=2/9当k金求积公式精确成立,而x4公式不成立,可见该求积公式具有3次代数精度.例2试用梯形公式、科茨公式和辛卜生公式计算定积分,4改(计算结果取5位有效数字)JOJ(1)用梯形公式计算/(0.5)+/0)=0.25x0.70711+1=0.4267873212327(2)用柯特斯公式系数为9090909090版人7,而1+32、阪J+12x园+加画+7x如2_x4.94975+25.29822+1039223+2993326+7=0.43096(3)如果要求精确到10一5,用复化辛卜生公式,截断误差为/网二去叱Ma=max|/4)(z)6醛Y制=max15-7/2A<1,

22、5Rn二二百；只需把0.5,14等分,分点为0.5,0.625,0.75,0.875,1G125(卢双超-/(0.5)+2/(075)+4(/(0碰»+/(0.875)+川)0,70711+2K0,86602分4x(0.79057+0.9354(+I-0.43096例3用三点高斯-勒让德求积公式计算积分sin(t+1)解做变量替换工有查表得节点为40.774596669和0;系数分另为0.5555555556和0.8888888889-1sin(i+1)一出sin-(-0.774S96669+1)x05555555556x_J-0.774596669+1sinl(0+l)sin-(0

23、.774596669+1)+0.888888889X2+0555555556X1=0.940831240+10774596669例4函数值f(1.0)=0.250000,f(1.1)=0.226757,f(1.2)=0.206612,用三点公式计算j在x=1.0,1.1,1.2处的导数值.解三点导数公式为K4仁北I+4九-匕2h/k白.i一4八+如*Qk=1,2,3,n-1本例取X0=1.0,x1=1.1,X2=1.2,y0=0.250000,y1=0.226757,y2=0.206612,h=0.1.于是有/q_0)冏*0,25000000+4x0.226757-0206612)=-0247

24、嵬-Ql)(-0,25000000+0.206612)-0.2169412x01(0.25000000-4x0.226757+3x0.206612)=-0185%例5选择填空题1.如果用复化梯形公式计算定积分e-!tdx,要求截断误差不超过0.5X104,试问n-()(A)41(B)42(C)43(D)40解答；复化的梯形公式的截断误差为二=104感P*1212制<0,5xl0-4,n=40.8,取n241.应选择(A).2.n=3时,柯特斯系数等,以那么呢0a解答：由柯特斯系数的归一性质,第六章常微分方程的数值解法考核知识点尤拉公式,梯形公式,改良尤拉法,局部截断误差；龙格一一库塔法,

25、局部截断误差.复习要求1 .掌握尤拉法和改良的尤拉法梯形公式、预报-校正公式,知道其局部截断误差.2 .知道龙格一库塔法的根本思想.知道二阶、三阶龙格一库塔法.掌握四阶龙格一-库塔法,知道龙格一库塔法的局部截断误差.三例题yr=-y-10工x£0.6例1用尤拉法解初值问题",取步长h=0.2o计算过程保70=1留6位小数.解h=0.2,f(x)=-yxy2.首先建立尤拉迭代格式加二M+版("/)二兀-机-融&=02九(4-1加已=0,1,2)当k=0,xi=0.2时,xo=0,yo=1,有y(0.2)；tyi=0.2X1(40X1)=0.8当k=1泾=0.4时,xi=0.2,yi=0.8,有y(0.4)即=0.2X0.8X(4-0.2X0.8)=0.6144当k=2,x3=0.6时,x2=0.4,y2=0.6144,有y(0.6)醐3=0.2X0.6144X(40.4X0.4613)=0.8一,y'+y+v疝>0例2用尤拉预报校正公式求解初值问题V,取步长h=0.2,刈=1