(完整)北京航空航天大学数值分析课程知识点总结,推荐文档

1、1.2 误差知识与算法知识1.2.2 绝对误差、相对误差与有效数字设 a 是准确值 x 的一个近似值，记 ex a ，称 e 为近似值 a 的绝对误差，简称误差。如果 |e |的一个上界已知，记为，即 | e |，则称 为近似值 a 的绝对误差限或绝对误差界，简称误差限或误差界。记 erex a ，称 er为近似值 a 的相对误差。由于 x 未知，实际上总把e 作为 a 的xxaex ae 的上界，即 r相对误差，并且也记为 er，相对误差一般用百分比表示。aar| a |称为近似值 a 的相对误差限或相对误差界。定义 设数 a 是数 x 的近似值。如果a 的绝对误差限时它的某一位的半个单位，

2、并且从该位到它的第一位非零数字共有n 位，则称用a 近似 x 时具有 n 位有效数字。函数求值的误差估计设 uf (x) 存在足够高阶的导数，a 是 x 的近似值， 则 uf (a) 是 uf (x) 的近似值。若 f'(a) 0 且 | f ''(a) | / | f '(a) |不很大，则有误差估计e(u)f '(a)e(a)。(u)f '(a)(a)若 f '(a) f ''(a) .f (k 1) (a)0,f ( k) (a)0 ，且比值f( k)(a)ke(u)k!e( a)大，则有误差估计。f( k)(a)k

3、(u)(a)k !nf (a1, a2 ,., an ) e(a )e(u)i1xii对于 n 元函数，有误差估计nf ( a1 ,a2,., an )(u)(ai )i1xif (k 1) (a) / f (k ) (a) 不很；若一阶偏导全为零或很小，则要使用高阶项。算法及其计算复杂性（ 1）要有数值稳定性，即能控制舍入误差的传播。（ 2）两数相加要防止较小的数加不到较大的数中所引起的严重后果。（ 3）要尽量避免两个相近的近似值相减，以免严重损失有效数字。（ 4）除法运算中，要尽量避免除数的绝对值远远小于被除数的绝对值。1.3 向量范数与矩阵范数向量范数定义 定义在 Rn 上的实值函数 ?

4、称为向量范数，如果对于Rn 中的任意向量x 和 y 满足：（1）正定性： x0 ，当且仅当x 0 时， x0 ；1（2）齐次性：对任一数kR ，有kxkx ；（3）成立三角不等式：xyxy。定理 1.1对 Rn 中的任一向量x( x1, x2 ,., xn )T，记nx 1xii 1nx 2xi2i 1xmax xi1 in则 ?1， ? 2和 ? 都是向量范数。定理 1.2设? 和?是 Rn 上的任意两种向量范数，则存在与向量x 无关的常数 m 和M(0<m<M) ，使下列关系式成立m xxMx , xRn1.3.2 矩阵范数定义 定义在 Rn n 上的实值函数? 称为矩阵范数，

5、如果对于Rn n 中的任意矩阵 A 和 B 满足：（1） A0 ，当且仅当 A0时， A0；（2）对任一数kR ，有 kAkA ；（3）ABAB；（4） ABA B 。定义 对于给定的向量范数? 和矩阵范数 ? ，如果对于任一个 xRn 和任一个 ARn n 满足 AxAx ，则称所给的矩阵范数与向量范数是相容的。定理 1.3设在 Rn 种给定了一种向量范数，对任一矩阵ARn n ，令 A = maxAx ，则由x 1此定义的? 是一种矩阵范数，并且它与所给定的向量范数相容。定理 1.4 设 A aij Rn n ，则2nA 1maxaij1 j n i 1A 2max ( AT A)nAma

6、xaij1 i nj1其中 max ( AT A) 表示矩阵 AT A 的最大特征值 （ AT A 是正定或半正定矩阵，它的全部特征值非负）。naij2还有一种常见的矩阵范数A Fi , j，且与向量范数 ? 2 相容，但是不从属于任何1向量范数。单位矩阵I 的任何一种算子范数I= max Ix1。x1定理 1.5 设矩阵ARn nA1IA的某种范数为非奇异矩阵， 并且当该范数为算子，则范数时，还有 IA11成立。1A2.1 Gauss 消去法2.1.1 顺序 Gauss 消去法定理 2.1 顺序 Gauss 消去法的前n-1 个主元素 akk(k ) ( k1,2,., n 1) 均不为零的

7、充分必要条a11(1). a1(1)k件是方程组的系数矩阵A 的前 n-1 个顺序主子式 D k. 0,( k 1,2,., n 1)ak(1)1. akk(k )2.1.2 列主元素 Gauss 消去法定理 2.2 设方程组的系数矩阵A 非奇异，则用列主元素Gauss 消去法求解方程组时，各个列主元素 ai(kkk) (k1,2,., n1)均不为零。2.2 直接三角分解法2.2.1 Doolittle分解法 (单位下三角 +上三角 )与 Crout 分解法 (下三角 +单位上三角 )定理 2.3 矩阵A aij nn (n 2) 有唯一的 Doolittle 分解的充分必要条件是A 的前

8、n-1 个顺序主子式 Dk0,( k 1,2,., n 1) 。推论矩阵Aaij n n ( n2) 有唯一的 Crout 分解的充分必要条件是A 的前 n-1 个顺序主子式 Dk 0,( k1,2,., n1) 。2.2.2 选主元的Doolittle分解法3定理 2.4若矩阵 ARn n 非奇异，则存在置换矩阵Q，使 QA 可做 Doolittle分解。2.2.3 三角分解法解带状线性方程组定理 2.5(保带状结构的三角分解) 设 A aij n n 是上半带宽为 s、下半带宽为 r 的带状矩阵，且 A 的前 n-1 个顺序主子式均不为零，则A 有唯一的 Doolittle分解a11a1,

9、s1. .ar1,1.A.an s, n. .an,nr.a1u11u12.u1,s 1l 2,11. . . .uns,nlr 1,1. . . . .un1,nl n,n r. ln ,n11unn为节省空间，用C(m,n) 存储 A的带内元素，其中m=r+s+1 ，并且 aijci j s 1, j 。拟三对角线性方程组的求解方法a1c1d1d2a2c2A. . . .dn 1an 1cn 1cndnanp11 q1s1d2p21 q2s2. . . .qn 2sn 2.dn 1pn 11sn 1r1r2. rn 2rn 1rn12.3 矩阵的条件数与病态线性方程组矩阵的条件数与线性方程

10、组的性态定义 对非奇异矩阵 A ，称量 | A | A 1 | 为矩阵 A 的条件数，记作 cond( A) | A | A 1 |。矩阵 A 的条件数与所取的矩阵范数有关，常用的条件数是4cond( A)| A | | A 1 | ， cond( A)2| A |2| A 1 |2性质 1对任何非奇异矩阵 A， cond( A)1。性质 2设 A 是非奇异矩阵， k0 是常数，则有 cond( kA)cond( A) 。性质 3设 A 是非奇异的是对称矩阵，则有cond( A)21，其中 1和n 分别是矩阵 A 的n模为最大和模为最小的特征值。性质 4设 A 是正交矩阵，则有cond( A)

11、21 。2.3.2 关于病态线性方程组的求解问题（1）采用高精度的算术运算。（2）平衡方法（行平衡，取每行绝对值最大数的倒数组成对角阵，乘在原方程左右两边）。（3）残差校正。2.4 迭代法迭代法的一般形式及其收敛性x( k 1)Gx( k)d( k0,1,.)定义 设 n n 矩阵 G 的特征值是 1 ,2 ,., n ，称(G)max |i |为矩阵 G 的谱半径。1 i n定理 2.9对任意的向量 d，迭代法收敛的充分必要条件是(G )1 。定理 2.9如果矩阵 G 的某种范数 |G|<1，则（1）方程组的解x* 存在且唯一；（2）对于迭代公式，有lim x(k )x* ,x(0)R

12、 ，且下列两式成立k| x( k)x*| G |k| x(1)x(0)|1|G| x( k)x*|G| x( k)x( k1) |1|G|2.4.2 Jacobi 迭代法A DLUx( k 1)D1 (LU ) x( k)D 1b(k0,1,.)GJD 1(LU )定理 2.10 Jacobi 迭代法收敛的充分必要条件是(GJ )1。定理 2.11 如果 |1，则 Jacobi 迭代法收敛。GJ引理 2.1 若矩阵 ARn n 是主对角线按行 (或按列 )严格占优阵，则A 是非奇异矩阵。5定理 2.12 如果方程组的系数矩阵式主对角线按行( 或按列 )严格占优阵，则用Jacobi 迭代法求解必

13、收敛。迭代法ADLUx( k1)( DL) 1Ux (k )(DL) 1 b(k 0,1,.)GG(D L) 1U定理 2.13 GS 迭代法收敛的充分必要条件是(GG )1。定理 2.14如果|1，则迭代法收敛。|JacobiGG定理 2.15如果方程组的系数矩阵式主对角线按行(或按列 )严格占优阵， 则用 GS 法求解必收敛。定理 2.16如果方程组的系数矩阵A 是正定矩阵，则用GS 法求解必收敛。2.4.4 逐次超松弛迭代法1DL(11)DUAx(k 1)(1DL) 1(11) DU x(k )(1D L) 1b(k0,1,.)GS( 1 DL) 1(11 ) DU 实际使用的形式x(

14、k1)D 1 Lx( k 1)(11)ID1U x(k )D 1b( k0,1,.)它的分量形式是x( k 1)i 1 aijx(k 1)1) x(k )n aijx(k )bi( k0,1,.)j(1jij 1 aiiij i 1 aiiaii定理 2.17 SOR 方法收敛的充分必要条件是(GS )1 。定理 2.18如果|1，则方法收敛。|SORGS定理 2.19 SOR 方法收敛的必要条件是02 。定理 2.20如果方程组的系数矩阵A 是主对角线按行 (或按列 )严格占优阵， 则用 01的SOR 方法求解必收敛。定理 2.21如果方程组的系数矩阵A 是正定矩阵，则用02的 SOR 方法

15、求解必收敛。* 实系数二次方程x2pxq0 的两个根之模均小于1 的充要条件是：| q |1,1pq0,1pq 0*A为正定矩阵A 的各阶顺序主子式全大于零。3.1 幂法和反幂法幂法 （用于计算矩阵的按模为最大的特征值和相应的特征向量）第一种幂法迭代格式：6u0Rn & u00k1ukT1uk 1yk1uk1 / k1ukAyk1kykT1uk(k1,2,.)第二种幂法迭代格式：u0| hr( kyk 1ukk(k=(h1(0) ,.,hn(0) )T & u01)( k 1) | max | hj |uk 1 / | hr( k 1) |(k )( k)Ayk 1(h1,.,

16、hnk 1( k)sgn(hr)hr1,2,.)0)Tk 作为1 的近似值， yk -1 作为 A 的属于1 的特征向量。反幂法第一种反幂法迭代格式：uRn & u000k 1ukT1uk 1yk 1uk1/ k1Aukyk1kyTk 1uk(k1,2,.)1 作为n 的近似值， yk -1 作为 A 的属于n 的特征向量。还可以用带原点平移的反幂法求k矩阵 A 的某个特征值s 。3.3 QR 方法矩阵的 QR 分解设 vRn 是单位向量，令HI2vvT ，则 H 是对称正交矩阵，称为Householder 矩阵。引理3.1 设有非零向量sRn 和单位向量eRn ，必存在Househo

17、lder 矩阵 H ，使得Hse ，其中是实数，并且|sT s 。 (可取 vse)( se)T ( se)7定理 3.2任何 nn 实矩阵 A 总可以分解为一个正交矩阵Q 与一个上三角矩阵 R 的乘积。设 ai 1(i2,3,., n) 不全为零，令s(a ,., a) T111n1c1sgn(a11)s1T s1 (取 sgn(0)1)u1s1c1e1H 1 I 2u1u1T / (u1Tu1)c1a12(2).a1(2)nA2H 1 A0 . .0an(2)2.ann(2)对第 j 列， a(ij1,.,n) 不全为零，令sj(j )( j )T，并继续计算。ij(0,.,0, ajj,

18、., anj )最终得到 AnH n 1 H n 2 .H 1 A 是一个上三角矩阵。则QH 1H 2 .H n 1, R An ，且AQR 。3.3.2矩阵的拟上三角化设 ai 1(i3,4,., n) 不全为零，令s1 (0, a21,., an1 )Tc1sgn(a21 )|s1|2 (取 sgn(0)1)u1s1c1e2H 1 I 2u1u1T / (u1Tu1)a11 a12(2).a1(2)nc1.A(2)H1 AH10 . .0an(2)2.ann(2)对第 j 列，(2,.,) 不全为零，令( j )( j )T，并继续计算。aijijnsj(0,.,0, aj 1, j ,.

19、, anj)最终得到 A( n 1)H n2.H 2 H 1 AH 1H 2 .H n 2 为拟上三角矩阵，令P H 1H 2.H n 2 ，则 A( n 1)PTAP 。8带双步位移的QR 方法基本 QR 方法的迭代公式是A1ARn nAkQk RkAk1RkQk(k1,2,.)4.1 非线性方程的迭代解法简单迭代法及其收敛性定理 4.1设函数(1)当 x a,b(2)当 x(a,b)( x)Ca, b ,在 (a,b)内可导 ,且满足两个条件 :时 ,( x) a, b ；时 , |'( x) | L 1 , 其中 L 为一常数。则有如下结论 :(1)方程 x=(x) 在区间 a,

20、b 上有唯一的根s ；(2)对任取 x0a, b ，简单迭代法x=( x ) 产生的序列 x a,b 且收敛于s；k 1kk(3) 成立误差估计式Lkx0 |s-xk | x11L| s-xk |L | xkxk 1 |1L定理 4.2 设 s=( s) ，'( x) 在包含 s 的某个开区间内连续。如果| '( s)|<1 ，则存在0 ，当 x0 s , s 时，由简单迭代法 xk1= (xk ) 产生的序列 xk s , s 且收敛于 s。简单迭代法的收敛速度定理 4.3 设函数( x)Ca, b ，'( x)C a, b ，且满足如下条件：（1）当 x a,

21、 b（2）当 x(a,b)时 ,( x)a, b ；时 ,'( x)0 ， | '( x) | L1 , 其中 L 为一常数。则对任取对任取x0 a,b ，简单迭代法 xk 1=(xk ) 产生的序列 x 收敛于方程 x=(x) 在k a,b 内的唯一的根s ，并且当 xs 时 xk 是线性收敛的。0定理 4.4 设 s=( s) ，( m) (x) 在包含 s 的某个开区间内连续 ( m 2 )。如果9(i) (s)0(i1,2,., m1)(m) ( s) 0则存在0 ，当 x0 s, s 但 x0s 时，由简单迭代法xk 1= ( xk ) 产生的序列 x k以 m 阶收

22、敛速度收敛于s 。迭代法yk( xk ), zk( yk )( yx) 2xk 1xkkk(k0,1,.)zk2 ykxk定理 4.5（局部）设 s=(s) ，( x) 在包含 s 的某个开区间内具有连续的二阶导数，并且'(s)1，则存在0，当 x0 s, s0s 时，由 Steffensen 迭代法产生的序但 x列 xk 至少以二阶收敛速度收敛于s 。4.1.5 Newton 法xk 1 xkf ( xk )(k 0,1,.)f'(xk )定理4.6（局部）设 s 是方程f (x)0 的根，在包含s 的某个开区间内f ''(x) 连续且f '(x)0，

23、则存在0，当 x0 s, s 时，由 Newton 法产生的序列 x 收敛于 s；k若 f ''(s)0 且 x0s，则序列 xk 是平方收敛的。定理 4.7（大范围）设函数 f ( x) 在区间 a,b 上存在二阶连续导数，且满足条件：（1） f ( a) f (b)0；（2） f ''(x) 在区间 a,b 上不变号；（3）当 x a, b 时， f '( x) 0；（4） x0 a,b ， f (x0 ) f ''(x0 )0 。则由 Newton 法产生的序列 xk 单调收敛于方程 f ( x)0 在a,b 内唯一的根 s ，并且

24、至少是平方收敛的。割线法xk 1xkf (xk )( xkxk 1) (k 0,1,.)f ( xk )f (xk 1)10定理 4.8设 f ( s)0 ， 在 s 地 某 邻域 内 f ''(x) 连续 且 f '(x)0，则存在0 ， 当x 1, x0I 时，由割线法产生的序列 xk 收敛于 s ，且收敛速度的阶至少为 1.618。4.1.8 单点割线法xk 1xkf (xk )( xkx0 ) (k 1,2,.)f ( xk )f ( x0 )定理 4.9设函数 f ( x) 在区间 a,b 上存在二阶连续导数，且满足条件：（1） f ( a) f (b)0；（

25、2） f ''(x) 在区间 a,b 上不变号；（3）当 x a, b 时， f'( x) 0 ；（4） x0 , x1 a, b 且 f ( x0 ) f ''(x0 )0 ， f ( x0 ) f ( x1 )0 。则有单点割线法产生的序列 xk 单调收敛于方程 f (x)0 在 a,b 内唯一的根 s ，并且收敛速度是一阶的。4.2 非线性方程组的迭代解法简单迭代法定理 4.13（压缩映像原理）设G:DRnRn 在闭区域 D0D 上满足两个条件：（1） G 把 D0 映入它自身，G(D0)D0 ；（2） G 在 D0 上是压缩映射，即存在常数L(0

26、,1) ，使对任意的x, y D0 ，有| G (x) G ( y) | L | x y |则有下列结论：（1）对任取的 x(0)D0 ，由迭代公式 x( k1)G ( x(k ) )(k0,1,.) 产生的序列 x( k) D0 ，且收敛于方程组 xG ( x) 在 D 内的唯一解 x*；0（2）成立误差估计式| x*x( k)|Lk| x(1)x(0)|1L| x*x( k)|L| x(k)x( k 1) |1L定理 4.14（局部）设G:DRnRn ， x*int( D ) 是方程组 x G ( x) 的解， G 在 x* 处11可微。若 G '(x* ) 的谱半径(G '

27、;( x* )1，则存在开球D0 x | xx* |,0D ，使对任意的 x(0)D0 ，由迭代法 x( k 1)G (x(k) )(k0,1,.) 产生的序列 x( k ) D0 且收敛于 x*。4.2.3 Newton法x(k 1)x( k )F '(x( k) ) 1 F ( x( k) )(k0,1,.)定理4.15 设 x*int( D ) 是方程组 F (x) 0 的解， F : DRnR n 在包含 x*的某个开区域 SD 内连续可微，且F '(x* ) 非奇异，则存在闭球D0 x | xx* |,0S ，使对任意的的x(0)D0 ，由 Newton 法产生的序列

28、 x(k ) D0 且超线性收敛于x* ；若更有F ( x) 在域 S 内二次连续可微，则序列 x(k) 至少是平方收敛的。5.1 代数插值一元函数插值（Lagrange 、 Newton ）定理 5.1设 x0 , x1 ,., xn 是互异的实数，对于给定的实数x ，实值函数f (t ) 在区间 I x 上具有n+1 阶导数，则插值公式f (x)pn ( x) 的余项为Rnf ( n 1) ()n 1(x)(n 1)!:其中I x 且依赖于 x ，n 1( x)( xx0 )(xx1).( xxn ) 。f(n 1)( ):* f x0,I x, x1,., xn , x(n 1)!定理

29、5.2设 x0 , x1 ,., xn 是互异的实数，对于给定的实数x ，实值函数f (t ) 在区间 I x 上具有m+n+2 阶导数， H m n 1H m n 1( xi )yi ,(i0,1,., n)( x) 则是满足条件(x )y',( k 0,1,., m)H 'm n 1ikik的 Hermite 插值多项式，则用H m n 1 (x) 近似代替 f ( x)的余项为f( m n 2)( )mR( x)(m n2)! n 1( x) k 0( x xi k ):其中I x 且依赖于 x 。5.3 样条插值12样条函数定义步长为 1、内节点等距的k 次 B 样条k

30、 (x)1 k1( 1)jk 1( xk 1k,<x<k! j0j2j )+性质 1k ( x)k ( x) 。性质 2 当 |x|k1( x)k1k ( x)0 。时，k0 ；当 |x|2时，2定义 步长为 h、内节点等距的k 次 B 样条x xi (k1)/21k 1jk1k<x<k ()j 0 ( 1),hk ! hkj( x xi k j )+性质 1 对称轴为 x=xi(k1)/2。性质 2k ( xxi ( k 1)/2 )0, x( xi k , xi 1)h0, x ( xi k , xi 1)空间 Dk,常用的两组基底及表示L1,x,., xk ,( x x1)k ,.( xxn 1 )k s( x)kaj xj1 n1k, a x bk ! jc j ( x x j )j01L k( xxj(k 1)/2),axb( j0,1,., nk1)hnk 1( xxj ( k 1)/2 ), a x bs( x)c jkj0h样条为基底的三次样条插值函数第一种边界条件s''(x0 )y0'' , s''(xn )yn'' ：c02c1c2h2 y0''c