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

PAGE PAGE # 如果|e| 如果|e|的一个上界已知, 界，简称误差限或误差界。 e x a 记er x x 1.2误差知识与算法知识 1.2.2绝对误差、相对误差与有效数字 设a是准确值x的一个近似值，记 e x a，称e为近似值a的绝对误差，简称误差。 记为 ，即|e| ，则称 为近似值a的绝对误差限或绝对误差 e 称er为近似值a的相对误差。由于 x未知，实际上总把 作为a的 a er er的上界，即r |a| 相对误差，并且也记为er a x a ，相对误差一般用百分比表示。 a e(u)大，则有误差估计(u)f ⑹ e(u) 大，则有误差估计 (u) f ⑹(a)k! k e(a) f(k)(a) k! o k (a) 称为近似值a的相对误差限或相对误差界。 定义 设数a是数x的近似值。如果 a的绝对误差限时它的某一位的半个单位，并且从该位 到它的第一位非零数字共有 n位，则称用a近似x时具有n位有效数字。 1.2.3函数求值的误差估计 设u f(x)存在足够高阶的导数，a是x的近似值，则u f(a)是u f(x)的近似值。 若f (a) 0且| f(a)|/| f (a)| 若f (a) 0且| f(a)|/| f (a)|不很大，则有误差估计 。 (u) f (a) (a) f(k 1 (a) / f(k) (a)不很若 f (a) f (a) ... f f(k 1 (a) / f(k) (a)不很 e(u)对于n元函数，有误差估计(u) TOC \o 1-5 \h \z \o Current Document f(a1, a2,..., an) ) e(u) 对于n元函数，有误差估计 (u) Xi i ;若一阶偏导全为零或很 f ( a1 , a2 ,..., an ) / \ (ai) Xi 小，则要使用高阶项。 1.2.4算法及其计算复杂性 (1 )要有数值稳定性，即能控制舍入误差的传播。 (2 )两数相加要防止较小的数加不到较大的数中所引起的严重后果。 要尽量避免两个相近的近似值相减，以免严重损失有效数字。 除法运算中，要尽量避免除数的绝对值远远小于被除数的绝对值。 1.3向量范数与矩阵范数 1.3.1向量范数 定义 定义在Rn上的实值函数 ？称为向量范数，如果对于 Rn中的任意向量x和y满足: (1 )正定性：|x 0 ,当且仅当x 0时，x 0 ； （2）齐次性：对任一数 k R，有 kx k|x ； （3）成立三角不等式：|x y |x 1 |刖。 疋理1.1对Rn中的任一向量x （X1,X2，…，Xn），记 iix maxx则?i iix maxx 则?i , ?2 和？都是向量范数。 定理1.2设？和？是Rn上的任意两种向量范数，则存在与向量 M（OmM），使下列关系式成立 mix ix m\x , x Rn 1.3.2矩阵范数 定义 定义在Rn n上的实值函数||?称为矩阵范数，如果对于 Rn 足： （1） | A 0,当且仅当A 0时，A 0 ； （2） 对任一数 k R，有 ||kA |k | A ； （3） A B A B ； （4） AB AB。 定义 对于给定的向量范数 ？和矩阵范数 ？，如果对于任一个x 足 |Ax| ||A|||x| ，则称所给的矩阵范数与向量范数是相容的。 定理1.3设在Rn种给定了一种向量范数，对任一矩阵 A Rn n， 此定义的？是一种矩阵范数，并且它与所给定的向量范数相容。 定理 1.4 设 A [aij] Rnn，则 x无关的常数 m和 n中的任意矩阵A和B满 Rn和任一个A Rn n满 令A = ma^ Ax，则由 ||x|i x i 1 nmax (A A)max…j ia n max (A A) max …j i aij 其中max(ATA)表示矩阵ATA的最大特征值(ATA是正定或半正定矩阵， 它的全部特征值 非负)。 还有一种常见的矩阵范数IIAIf J ai2，且与向量范数||? 2相容，但是不从属于任何 向量范数。单位矩阵 i, j 1 I的任何一种算子范数11二max Ix] 1。 定理1.5设矩阵A Rn n的某种范数 A 1,则1 A为非奇异矩阵，并且当该范数为算子 成立。1 A范数时，还有 I A 1|| 成立。 1 A 2.1 Gauss消去法 2.1.1顺序Gauss消去法 定理2.1顺序Gauss消去法的前 定理2.1顺序Gauss消去法的前 n-1个主元素 akJ (k 1,2,..., n 1)均不为零的充分必要条 (1) anak； (1) an ak；) 0,(k 1,2,..., n 1)件是方程组的系数矩阵 A的前n-1个顺序主子式Dk 0,(k 1,2,..., n 1) ak