数值分析部分

1、数值分析讲义1第1章数值分析中的误差一、重点内容误差设精确值x*的近似值x,差e=x-x*称为近似值x的误差(绝对误差)。误差限近似值x的误差限e是误差e的一个上界，即|e|=|xx*|?&。相对误差er是误差e与精确值x*的比值，。常用计算。相对误差限是相对误差的最大限度，常用计算相对误差限。绝对误差的运算：£(x1±x2)=£(x1)+£(x2)s(x1x2)=|x1|e(x2)+|x2|e(x1)有效数字如果近似值x的误差限£是它某一个数位的半个单位，我们就说x准确到该位。从这一位起到前面第一个非0数字为止的所有数字称为x的有效数

2、字。关于有效数字：(1)设精确值x*的近似值x,x=±0.a1a2anx10ma1,a2,，an是09之中的自然数，且a1w0,|x-x*|?S=0.5X10m-l,1?l?n则x有l位有效数字.(2)设近似值x=±0.a1a2anx10m有n位有效数字，则其相对误差限(3)设近似值x=±0.a1a2anx10m的相对误差限不大于则它至少有n位有效数字。(4)要求精确到103,取该数的近似值应保留4位小数。一个近似值的相对误差是与准确数字有关系的，准确数字是从一个数的第一位有效数字一直数到它的绝对误差的第一位有效数字的前一位，例如具有绝对误差e=0.0926的数x

3、=20.7426只有三位准确数字2,0,7。一般粗略地说，具有一位准确数字，相对于其相对误差为10%的量级；有二位准确数字，相对于其相对误差为1%的量级；有三位准确数字，相对于其相对误差为0.1%的量级。二、实例例1设x*=p=3.1415926-近似值x=3.14=0.314X101,即m=1,它的误差是0.001526，有|xx*|=0.001526？0.5X1013即l=3,故x=3.14有3位有效数字。x=3.14准确到小数点后第2位。又近似值x=3.1416,它的误差是0.0000074,有|x-x*|=0.0000074？0.5X1015即m=1,l=5,x=3.1416有5位有效

4、数字。而近似值x=3.1415,它的误差是0.0000926有|x-x*|=0.0000926?0.5X1014即m=1,l=4,x=3.1415有4位有效数字。这就是说某数有s位数，若末位数字是四舍五入得到的，那么该数有s位有效数字；若末位数字不是四舍五入得到的，那么该数有s位或s1位有效数字。例2指出下列各数具有几位有效数字，及其绝对误差限和相对误差限：2.0004-0.0020090009000.00解因为x1=2.0004=0.20004X101,它的误差限0.00005=0.5X1015,即m=1,l=5,故x1=2.0004有5位有效数字。相对误差限。数值分析讲义2x2=0.002

5、00,误差限0.000005,因为m=-2,l=3,x2=-0.00200有3位有效数字。相对误差限er=0.000005/0.00200=0.25%。x3=9000,绝对误差限为0.5,因为m=4,l=4,x3=9000有4位有效数字，相对误差限er=0.5/9000=0.0056%。x4=9000.00,绝对误差限0.005,因为m=4,l=6,x4=9000.00有6位有效数字，相对误差限为er=0.005/9000.00=0.000056%。由x3与x4可以看到小数点之后的0,不是可有可无的，它是有实际意义的。例3ln2=0.69314718，精确到10-3的近似值是多少？解精确到10

6、-3=0.001,即绝对误差限是e=0.05%,故至少要保留小数点后三位才可以。ln2=0.693。三、练习题1 .设某数x*,它的保留三位有效数字的近似值的绝对误差是。2 .设某数x*,它的精确到10-4的近似值应取小数点后位。3 .()的3位有效数字是0.236X102。(A)235.54X10-1(B)235.418(C)2354.82X10-2(D)0.0023549X1034 .设a*=2.718181828，取a=2.718,则有()，称a有四位有效数字。(A)|a-a*|?0.5X104(B)|aa*|?0.5X1014(C)|aa*|?10-4(D)|aa*|?0.00035

7、.设某数x*,对其进行四舍五入的近似值是()，则它有3位有效数字，绝对误差限是0.5X10一4o(A)0.315(B)0.03150(C)0.0315(D)0.003156 .以下近似值中，保留四位有效数字，相对误差限为0.25X103。(A)0.01234(B)T2.34(C)220(D)0.22007 .将下列各数舍入成三位有效数字，并确定近似值的绝对误差和相对误差。(1) 2.1514(2)-392.85(3)0.0039228 .已知各近似值的相对误差，试确定其绝对误差：(1) 13267er=0.1%(2)0.896er=10%9 .已知各近似值及其绝对误差，试确定各数的有效位数。(

8、1) 0.3941e=0.25X10-2(2)293.481e=0.1(3)0.00381e=0.1X10-410.已知各近似值及其相对误差，试确定各数的有效位数。(1)1.8921er=0.1X102(2)22.351er=0.15(3)48361er=1%四、练习题答案1 .该数有效数字第四位的一半。2 .五3.(A)4.(B)5.(C)6.(D)7. (1)2.15,e=-0.14X10-2,er=0.65X103;(2)393,e=-0.15,er=0.38X103;(3)0.00392,e=0.2X105,er=0.51x10-38. (1)e=0.13X102;(2)0.9X101

9、9. (1)2;(2)3;(3)210. (1)3;(2)1;(3)2第15章线性方程组的数值解法数值分析讲义3一、重点内容1 .高斯顺序消去法解线性方程组AX=b,对增广矩阵顺序作初等行变换，使矩阵A化为上三角形矩阵，再回代，从而得到线性方程组的解。要求作初等行变换消元过程中，。注意：本章讨论线性方程组的解的方法，不讨论解的存在性。2 .高斯列主元消去法在高斯顺序消去法中，每次消元之前，要确定主元，(k=1,2,3,，n-1)把第r行作为主方程，做第k次消元。把系数矩阵化为上三角形矩阵，从而得到线性方程组的解。3 .雅可比迭代法(简单迭代法)解线性方程组AX=b的雅可比迭代法公式为(k=0,

10、1,2,)4 .高斯一一赛德尔迭代法解线性方程组AX=b的高斯一一赛德尔迭代法公式为(i=1,2,，n；k=0,1,2,)5 .解的收敛性定理【定理1】高斯消去法消元过程能进行到底的充分必要条件是系数矩阵A的各阶顺序主子式不为0;AX=b能用高斯消去法求解的充分必要条件是A的各阶顺序主子式不为0。【定理4】(迭代法基本定理)设线性方程组X=BX+f对于任意初始向量X(0)及任意f,对应此方程组的迭代公式X(k+1)=B(k)X+f收敛的充分必要条件是，其中入i(i=1,2,，n)为迭代矩阵B的特征根。当入i为复数时，|入i|表示入i的模。【定理6】(迭代法收敛的充分条件)设线性方程组AX=b,

11、(1)若A是严格对角占优矩阵，则雅可比迭代法和高斯一一赛德尔迭代法收敛；(2)若A为对称正定矩阵，则高斯一一赛德尔迭代法收敛。注：设矩阵A=aijn,若则称矩阵A是严格对角占优矩阵。二、实例例1用顺序消去法解线性方程组解顺序消元于是有同解方程组回代得解x3=1,x2=1,x1=1,原线性方程组的解为X=(1,1,1)T。例2取初始向量X(0)=(0,0,0)T,用雅可比迭代法求解线性方程组解建立迭代格式(k=1,2,3,)数值分析讲义4第1次迭代，k=0X(0)=0,得到X(1)=(1,3,5)T第2次迭代，k=1X(2)=(5,3,-3)T第3次迭代，k=2X(3)=(1,1,1)T第4次迭

12、代，k=3X(4)=(1,1,1)T例3填空选择题：1 .用高斯列主元消去法解线性方程组作第1次消元后的第2,3个方程分别为。解选a21=2为主元，作行互换，第1个方程变为：2x1+2x2+3x3=3,消元得到是应填写的内容。2 .用选主元的方法解线性方程组AX=b,是为了()(A)提高计算速度(B)减少舍入误差(C)减少相对误差(D)方便计算答案：选择(B)3 .用高斯一一赛德尔迭代法解线性方程组的迭代格式中=(k=0,1,2,)答案：解答：高斯一一赛德尔迭代法就是充分利用已经得到的结果，求x2的值时应该用x1的新值。4 .当a()时，线性方程组的迭代解一定收敛。(A)>6(B)=6(

13、C)<6(D)>6或V6答案：(D)解答：当|a|>6时，线性方程组的系数矩阵是严格对角占优矩阵，由教材第10章定理6,迭代解一定收敛。三、练习题1 .用高斯列主元消去法解线性方程组2 .用高斯一一赛德尔迭代法求解线性方程组取初始值(4.67,7.62,9.05)T,求二次迭代值。3 .证明线性方程组的迭代解收敛。4 .用高斯顺序消去法解线性方程组，消元能进行到底的充分必要条件是5 .用列主元消去法解线性方程组，第1次消元，选择主元为()(A)3(B)4(C)-4(D)9数值分析讲义5四、练习题答案1 .X=(-4,1,2)T2 .(4.66619,7.61898,9.047

14、53)T3 .提示：系数矩阵是严格对角占优矩阵。4 .线性方程组的系数矩阵的各阶顺序主子式均不为0。5.(C)第2章函数插值与最小二乘拟合一、重点内容1 .函数插值已知函数f(x)的n个函数值yk=f(xk),k=0,1,2,，n。构造一个多项式P(x),使得P(xk)=yk。P(x)就是插值多项式，f(x)就是被插函数，xk就是插值节点。误差R(x)=f(x)P(x)。2 .拉格朗日多项式称n次多项式Pn(x)=y0l0+y1l1+ynln=为拉格朗日插值多项式，其中基函数(i=0,1,2,，n)当n=1时，线性插值P1(x)=yklk(x)+yk+1lk+1(x)其中基函数。当n=2时，得

15、到二次多项式，就是二次插值。拉格朗日插值多项式的余项为：，其中七C(a,b)注意：过n+1个互异点，所得的多项式应该是次数不超过n的多项式。3 .均差与牛顿插值多项式函数值与自变量的差商就是均差，一阶均差(或记作fx0,x1);二阶均差(或记作fx0,x1,x2)均差有两条常用性质：(1)均差用函数值的线性组合表示；(2)均差与插值节点顺序无关。用均差为系数构造多项式，就是牛顿插值多项式Nn(x)=f(x0)+fx0,x1(xx0)+fx0,x1,x2(xx0)(xx1)十+fx0,x1,x2,，xn(xx0)(xx1)(xx2)(xxn-1)牛顿插值多项式的余项为:Rn(x)=f(x)-Nn

16、(x)=fx,x0,x1,x2,，xn(xx0)(xx1)(xx2)(xxn1)(xxn)4 .分段线性插值已知n+1个互异节点x0,x1,，xn构造一个分段一次的多项式P(x),且满足：(1)P(x)在a,b上连续；(2)P(xk)=yk(k=0,1,2,，n);(3)P(x)在xk,xk+1上是线性函数。分段线性插值函数其中lk(x)(k=0,1,2,，n)是分段线性插值基函数。(i=1,2,，n-1)5 .三次样条插值函数(k=0,1,2,，n-1)(xk?x?xk+1)其中S2(xk)=mk(k=0,1,2,，n),hk=xk+1xk(k=0,1,2,，n1),m0,m1,，mn满足的

17、方程组是数值分析讲义6(*)其中：，(k=1,2,，n-1)(1)当已知S«x0)=yG,Sqxn)=yG时，(*)式中m0=1,ln=1,(2)当已知S2(x0)=y20=m0,S2(xn)=y2n=mn时，(*)式化为6 .最小二乘法用j(x)拟合数据(xk,yk)(k=1,2,，n),使得误差的平方和为最小，求j(x)的方法，称为最小二乘法。(1)直线拟合若，a0,a1满足法方程组(2)二次多项式拟合若，a0,a1,a2满足法方程组二、实例例1已知函数y=f(x)的观察数据为xk-2045yk51-31试构造拉格朗日多项式Pn(x),并计算P(1)。只给4对数据，求得的多项式不

18、超过3次解先构造基函数所求三次多项式为P3(x)=P3(1)=例2已知函数y=f(x)的数据如表中第1,2歹U。计算它的各阶均差。解依据均差计算公式，结果列表中。kxkf(xk)一阶均差二阶均差三阶均差四阶均差00.400.4107510.550.578151.116000.650.696751.168000.280000.800.888111.275730.358930.197330.901.201521.384100.433480.213000.03134计算公式为一阶均差(k=0,1,2,3)二阶均差(k=0,1,2)三阶均差(k=0,1)四阶均差例3设x0,x1,x2,，xn是n+1个

19、互异的插值节点，lk(x)(k=0,1,2,，n)是拉格朗日插值基函数，证明：数值分析讲义7(1)；(2)(m=0,1,2,，n)证明(1)Pn(x)=y0l0+ylll+ynln=当f(x)三1时，1=由于，故有(2)对于f(x)=xm,m=0,1,2,，n,对固定xm(0?m?n),作拉格朗日插值多项式，有当n>m1时，f(n+1)(x)=0,Rn(x)=0,所以注意：对于次数不超过n的多项式，利用上结果，有可见，Qn(x)的拉格朗日插值多项式就是它自身，即次数不超过n的多项式在n+1个互异节点处的拉格朗日插值多项式就是它自身。例4已知函数e-x的下列数据，用分段线性插值法求x=0.

20、2的近似值。x0.100.150.250.30e-x0.9048370.8607080.7788010.740818解用分段线性插值，先求基函数。所求分段线性插值函数为所以，e-0.2=P(0.2)=-0.81907X0.2+0.983569=0.819755例5已知数据如表的第2,3歹U,试用直线拟合这组数据。解计算列入表中。kxkykxkyk11414224.5493369184481632558.52542.5S153155105.5n=5。a0,a1满足的法方程组是解得a0=2.45,a1=1.25。所求拟合直线方程为y=2.45+1.25x例6选择填空题1.设y=f(x),只要x0,

21、x1,x2是互不相同的3个值，那么满足P(xk)=yk(k=0,1,2)的f(x)的插值多项式P(x)是(就唯一性回答问题)答案：唯一的数值分析讲义8解答：因为过3个互异节点，插值多项式是不超过2次的。设P(x)=a2x2+a1x+a0,其中a2,al,a0是待定数。P(xk)=yk,即这是关于a2,al,a0的线性方程组，它的解唯一，因为系数行列式所以，不超过2次的多项式是唯一的。2 .通过四个互异节点的插值多项式P(x),只要满足(),则P(x)是不超过一次多项式。(A)初始值y0=0(B)一阶均差为0(C)二阶均差为0(D)三阶均差为0答案：(C)解答：因为二阶均差为0,那么牛顿插值多项

22、式为N(x)=f(x0)+fx0,x1(x-x0)它是不超过一次的多项式。3 .拉格朗日插值多项式的余项是()，牛顿插值多项式的余项是()(A)(B) fx,x0,x1,x2,xn(xx1)(xx2)(xxn1)(xxn)(C)(D) fx,x0,x1,x2,，xn(xx0)(xx1)(xx2)(xxn1)(xxn)答案：(A),(D)o见教材有关公式。4.数据拟合的直线方程为y=a0+a1x,如果记那么系数a0,a1满足的方程组是()(A)(B)(C)(D)答案：(B)解答：因为法方程组为由第1个方程得到，将其代入第2个方程得到整理得故(B)正确。三、练习题1 .已知函数y=f(x),过点(

23、2,5),(5,9),那么f(x)的线性插值多项式的基函数为。2 .过6个插值节点的拉格朗日插值多项式的基函数l4(x)=。3 .已知多项式P(x),过点(0,0),(2,8),(4,64),(11,1331),(15,3375),它的3阶均差为常数1, 一阶，二阶均差均不为0,那么P(x)>()(A)二次多项式(B)不超过二次的多项式(C)三次多项式(D)四次多项式4 .已知y=f(x)的均差，。那么fx4,x2,x0=()(A)5(B)9(C)14(D)85 .求数据拟合的直线方程y=a0+a1x的系数a0,a1是使最小。6 .求过这三个点(0,1),(1,2),(2,3)的拉格朗日

24、插值多项式。7 .构造例2的函数f(x)的牛顿插值多项式，并求f(0.596)的近似值。8 .设l0(x)是以n+1个互异点x0,x1,x2,，xn为节点的格朗日插值基函数数值分析讲义9试证明:9,已知插值条件如表所示，试求三次样条插值函数。x123y2412y01-110.已知数据对(7,3.1),(8,4.9),(9,5.3),(10,5.8),(11,6.1),(12,6.4),(13,5.9)。试用次多项式拟合这组数据。四、练习题答案1.2.3.C4,B5.6.x+17,给定五对点，牛顿多项式是不超过4次的多项式。N4(x)=0.41075+1.11600(x0.40)+0.28000

25、(x0.40)(x0.55)+0.19733(x0.40)(x0.55)(x0.65)+0.03134(x0.40)(x0.55)(x0.65)(x0.80)将x=0,596代入牛顿多项式N4(x)中，得到：f(0.596)N(0.596)=0.631928,提示：求l0(x)的牛顿插值多项式。9.10. y=-0,145x2+3.324x-12.794第4章数值积分与微分一、重点内容1, m次代数精度求积公式对于任意不超过m次的代数多项式都准确成立，而对某一个m+1次代数多项式不成立。2,牛顿一一科茨求积公式：截断误差科茨系数：(k=0,1,2,，n),有两条性质。(2)牛顿一一科茨求积公式

26、的求积系数：Ak=(k=0,1,2,，n)(3)常见牛顿一一科茨求积公式梯形公式截断误差：R1f=复化梯形公式截断误差：，M2=抛物线公式复化抛物线公式截断误差：，科茨公式3 .高斯一一勒让德求积公式数值分析讲义10节点为的零点(高斯点)其余项：4 .微分公式(1)等距节点两点求导公式：(k=0,1,2,，n-1)(2)等距节点三点求导公式：(k=1,2,，n-1)二、实例例1试确定求积公式的代数精度。依定义，对xk(k=0,1,2,3,)，找公式精确成立的k数彳1解当f(x)取1,x,x2,计算求积公式何时精确成立。(1)取f(x)=1,有左边=,右边二(2)取f(x)=x，有左边=,右边二

27、(3)取f(x)=x2,有左边=,右边二(4)取f(x)=x3,有左边=,右边二(5)取f(x)=x4,有左边=,右边二当k?3时求积公式精确成立，而x4公式不成立，可见该求积公式具有3次代数精度。例2试用梯形公式、抛物线公式和科茨公式计算定积分(计算结果取5位有效数字)(1)用梯形公式计算(2)用抛物线公式用科茨公式系数为如果要求精确到10-5,用复化抛物线公式，截断误差为，,N?2只需把0.5,14等分,分点为0.5,0.625,0.75,0.875,1数值分析讲义11例3用三点高斯一一勒让德求积公式计算积分高斯型求积公式只能计算1,1上的定积分解做变量替换，查表得节点±0.77

28、4596669和0;系数分别为0.5555555556和0.8888888889注：该积分准确到小数点后七位是0.9460831,可见高斯型求积公式的精度是很高的。教材的第12章12.2节，用多种方法计算过该积分，它们的精度请读者自行比较。例4用三点公式计算在x=1.0,1.1,1.2处的导数值。已知函数值f(1.0)=0.250000,f(1.1)=0.226757,f(1.2)=0.206612解三点导数公式为k=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。于是有计算例5选择填空题1

29、.牛顿一一科茨求积公式与高斯型求积公式的关键不同点是。解答：牛顿一一科茨求积公式的节点和求积系数确定后，再估计其精度；高斯型求积公式是由精度确定其节点和求积系数。2.如果用复化梯形公式计算定积分,要求截断误差的绝对值不超过0.5X104,试问n?()(A)41(B)42(C)43(D)40答案：(A)解答；复化的梯形公式的截断误差中，故,n=40.8,取n?41。故选择(A)。3.已知n=3时，科茨系数，那么=答案：1/8解答：由科茨系数的归一性质,三、练习题1 .试确定求积公式的待定参数，使求积公式=A0f(0)+A1f(1)+A2f(2)的代数精度尽可能的高。2 .用复化抛物线公式计算定积

30、分。取n=4,保留4位有效数字。3 .试用四点(n=3)高斯一一勒让德求积公式计算积分4 .已知条件见例4。用两点求导公式计算f0.0),f(1.1)。5 .若用复化抛物线公式计算积分，要求截断误差的绝对值不超过0.5X104,试问n?()(A)1(B)2(C)4(D)36 .当n=6时，=()代数精度的。7 .用三点高斯一一勒让德求积公式计算积分四、练习题答案数值分析讲义121. A0=A2=1/3,A1=4/32. 0.11093.3.1416244.0.23243;0.201455.(B)6.(D)7.5次第13章方程求根一、重点内容1 .二分法：设方程f(x)=0在区间a,b内有根，用二分有根区间的方法，得到有根区间序列:x*=xn=(aO=a,b0=b),n=0,1,2,有误差估计式：？x*-xn?,n=0,1,2,二分区间次数：2 .简单迭代法：若方程f(x)=0表成x=j(x),于是有迭代格式：xn=j(xn1)(n=1,2,)x*xn若存在0vlv1Jj(