线性方程组的敏度分析与消去法的舍入误差分析

数值计算中的误差估计与分析在数值计算中，误差是无法避免的。

无论是数值积分、求根、线性方程组求解还是常微分方程求解，我们都需要对误差进行估计与分析，以保证结果的可靠性。

1.舍入误差：计算机中数字的存储精度是有限的，常用的浮点数表示法只能表示有限位数的小数。

当进行计算时，由于舍入操作会使结果产生一定的误差。

舍入误差是由于浮点数计算机表示能力造成的，它依赖于计算机所采用的机器数系统。

2.截断误差：在数值计算方法中，我们通常会使用有限项的级数展开式或多项式插值来近似解析解。

但由于展开或插值时的截断限制，会导致结果与真实结果之间的误差。

3.近似误差：数值计算方法本身就是在对问题进行近似求解，所以解的精确性受到近似精度的限制。

比如，对于数值积分来说，选择积分点的个数、插值多项式的次数都会影响结果的准确性。

4.舍入误差传播：在多步计算的过程中，每一步的舍入误差都会传播到下一步计算中，进而影响最终结果。

舍入误差的传播是一个累积效应，有时即使每一步舍入误差非常小，但在多步计算的累加下，也会导致结果产生很大的误差。

二、误差估计方法1.精度估计：对于一些数值方法，可以通过理论分析推导出误差的范围。

例如，对于数值积分，可以通过误差估计公式进行分析。

这种方法需要对问题进行数学建模，并具备一定的数学推导能力。

2.实验估计：对于一些复杂问题，很难通过理论分析得到精确的误差范围。

此时可以通过实验的方式来估计误差。

实验方法可以是计算机模拟实验，也可以是通过比较数值方法与解析解的差异来估计误差。

3.改进方法：除了估计误差大小，我们还可以通过改进数值方法来减小误差。

比如，可以采用更高阶的数值积分公式、使用更精确的数值微分方法等。

这些改进方法在一定程度上可以提高数值计算的准确性，并减小误差。

三、误差分析策略1.迭代策略：很多数值方法都是通过迭代来逐步逼近真实解的。

在迭代过程中，我们可以通过观察迭代序列的变化情况来判断结果是否趋近真实解，以及误差的变化是否在可接受范围内。

线性方程组的四种数值解法(电子科技大学物理电子学院，四川 成都 610054)摘要：本文介绍了四种求解线性方程组的数值解法: 雅克比迭代法、高斯赛德尔迭代法、高斯消去法和改进的平方根法的基本原理和算法流程，通过求解具体方程，对四种求解方法进行了对比。

对于雅克比迭代法和高斯赛德尔迭代法，研究了两种算法对求解同一方程组的迭代效率差异，结果表明高斯赛德尔迭代法达到同样精度所需迭代次数较少。

对于高斯消去法，通过选择列主元的方法提高算法的准确度，计算结果表明高斯消去法计算精确，且运算复杂度也不是很高。

对于改进的平方根法，其运算复杂度低，但对于给定的方程组有着严苛的要求。

关键词：雅克比迭代法；高斯赛德尔迭代法；高斯消去法；改进的平方根法；线性方程组引言线性方程组的求解在日常生活和科研中有着极其重要的应用，但在实际运算中，当矩阵的维数较高时，用初等方法求解的计算复杂度随维数的增长非常快，因此，用数值方法求解线性方程组的重要性便显现出来。

经典的求解线性方程组的方法一般分为两类：直接法和迭代法。

前者例如高斯消去法,改进的平方根法等，后者的例子包括雅克比迭代法，高斯赛德尔迭代法等。

这些方法的计算复杂度在可以接受的范围内，因此被广泛采用。

一般来说，直接法对于阶数比较低的方程组比较有效；而后者对于比较大的方程组更有效。

在实际计算中，几十万甚至几百万个未知数的方程组并不少见。

在这些情况下，迭代法有无可比拟的优势。

另外，使用迭代法可以根据不同的精度要求选择终止时间，因此比较灵活。

在问题特别大的时候，计算机内存可能无法容纳被操作的矩阵，这给直接法带来很大的挑战。

而对于迭代法，则可以将矩阵的某一部分读入内存进行操作，然后再操作另外部分。

本文使用上述四种算法求解对应的方程组，验证各种算法的精确度和计算速度。

1 算法介绍1.1 雅克比迭代法 1.1.1 算法理论设线性方程组(1)b Ax的系数矩阵A 可逆且主对角元素 均不为零,令并将A 分解成 (2)从而(1)可写成令其中. (3)以B 1为迭代矩阵的迭代法(公式)(4)称为雅克比(Jacobi)迭代法(公式),用向量的分量来表示,(4)为(5)其中为初始向量.1.1.2 算法描述 1给定迭代初始向量X 0以及误差要求delta 2根据雅克比迭代公式计算出下一组向量 3判断X 是否满足误差要求，即||X k+1 – X k || < delta4若误差满足要求，则停止迭代返回结果；若否，则返回第二步进行下一轮迭代1.2 高斯赛德尔迭代法nna ,...,a ,a 2211()nna ,...,a ,a diag D 2211=()D D A A +-=()b x A D Dx +-=11f x B x +=b D f ,A D I B 1111--=-=()()111f x B x k k +=+⎩⎨⎧[],...,,k ,n ,...,i x a ba xnij j )k (j j i iii)k (i21021111==∑-=≠=+()()()()()Tn x ,...x ,x x 002010=1.2.1 算法理论由雅克比迭代公式可知,在迭代的每一步计算过程中是用的全部分量来计算的所有分量,显然在计算第i 个分量时,已经计算出的最新分量没有被利用，从直观上看,最新计算出的分量可能比旧的分量要好些.因此，对这些最新计算出来的第次近似的分量加以利用,就得到所谓解方程组的高斯—塞德尔（Gauss-Seidel ）迭代法.把矩阵A 分解成(6)其中,分别为的主对角元除外的下三角和上三角部分,于是,方程组(1)便可以写成即其中(7)以为迭代矩阵构成的迭代法(公式)(8)称为高斯—塞德尔迭代法(公式),用变量表示的形式为(9)1.2.2 算法描述 1给定迭代初始向量X 0以及误差要求delta2根据高斯赛德尔迭代公式计算出下一组向量()k x ()1+k x ()1+k ix ()()1111+-+k i k x ,...,x 1+k()1+k x()1+k jx U L D A --=()nna ,...,a ,a diag D 2211=U ,L --A ()b Ux x L D +=-22f x B x +=()()b L D f ,U L D B 1212---=-=2B ()()221f x B x k k +=+⎩⎨⎧[],...,,k ,n ,,i x a x a b a xi j n i j )k (j ij )k (j ij i ii)k (i21021111111==∑∑--=-=+=++3判断X 是否满足误差要求，即||X k+1 – X k || < delta4若误差满足要求，则停止迭代返回结果；若否，则返回第二步进行下一轮迭代1.3 高斯消去法 1.3.1 算法理论下面三种变换称为初等行变换：1.对调两行；2.以数k ≠0乘某一行中的所有元素；3.把某一行所有元素的k 倍加到另一行对应的元素上去。

数值分析讲义第三章线性方程组的解法§3.0 引言§3.1 雅可比(Jacobi)迭代法§3.2 高斯-塞德尔(Gauss-Seidel)迭代法§3.3 超松驰迭代法§3.7 三角分解法§3.4 迭代法的收敛性§3.8 追赶法§3.5 高斯消去法§3.9 其它应用§3.6 高斯主元素消去法§3.10 误差分析§3 作业讲评3 §3.11 总结§3.0 引言重要性：解线性代数方程组的有效方法在计算数学和科学计算中具有特殊的地位和作用.如弹性力学、电路分析、热传导和振动、以及社会科学及定量分析商业经济中的各种问题.分类：线性方程组的解法可分为直接法和迭代法两种方法.(a) 直接法：对于给定的方程组，在没有舍入误差的假设下，能在预定的运算次数内求得精确解.最基本的直接法是Gauss消去法，重要的直接法全都受到Gauss消去法的启发.计算代价高.(b) 迭代法：基于一定的递推格式，产生逼近方程组精确解的近似序列.收敛性是其为迭代法的前提，此外，存在收敛速度与误差估计问题.简单实用,诱人.§3.1 雅可比Jacobi 迭代法 (AX =b )1基本思想:与解f (x )=0 的不动点迭代相类似,将AX =b 改写为X =BX +f 的形式,建立雅可比方法的迭代格式：X k +1=BX (k )+f ,其中，B 称为迭代矩阵.其计算精度可控,特别适用于求解系数为大型稀疏矩阵(sparse matrices)的方程组. 2问题:(a) 如何建立迭代格式？(b) 向量序列{X k }是否收敛以及收敛条件? 3 例题分析:考虑解方程组⎪⎩⎪⎨⎧=+--=-+-=--2.453.82102.7210321321321x x x x x x x x x (1)其准确解为X *={1, 1.2, 1.3}. 建立与式(1)相等价的形式：⎪⎩⎪⎨⎧++=++=++=84.02.01.083.02.01.072.02.01.0213312321x x x x x x x x x (2) 据此建立迭代公式：⎪⎩⎪⎨⎧++=++=++=+++84.02.01.083.02.01.072.02.01.0)(2)(1)1(3)(3)(1)1(23)(2)1(1k k k k k k kk k x x x x x x x x x (3) 取迭代初值0)0(3)0(2)0(1===xxx,迭代结果如下表.JocabiMethodP31.cpp迭代次数 x1 x2 x30 0 0 01 0.72 0.83 0.842 0.971 1.07 1.153 1.057 1.1571 1.24824 1.08535 1.18534 1.282825 1.095098 1.195099 1.2941386 1.098338 1.198337 1.2980397 1.099442 1.199442 1.2993358 1.099811 1.199811 1.2997779 1.099936 1.199936 1.29992410 1.099979 1.199979 1.29997511 1.099993 1.199993 1.29999112 1.099998 1.199998 1.29999713 1.099999 1.199999 1.29999914 1.1 1.2 1.315 1.1 1.2 1.34Jocobi迭代公式:设方程组AX=b, 通过分离变量的过程建立Jocobi迭代公式，即),,2,1()(1),,2,1(0,11n i x a b a x n i a b x a n ij j j ij i iii ii ni i j ij =∑-==≠∑=≠== 由此我们可以得到Jacobi 迭代公式：),,2,1()(11)1(n i x a b a xn ij j k i ij i iik i=∑-=≠=+[Jacobi 迭代公式的算法] 1: 初始化. n , (a ij ), (b j ), (x 1) , M . 2: 执行k =1直到M 为止. ① 执行i =1直到n 为止.ii nij j j ij i i a x a b u /)(1∑-←≠= ;② 执行i =1直到n 为止.i i u x ← ;③输出k , (x i ).另外，我们也可以建立Jacobi 迭代公式的矩阵形式. 设方程组AX =b ,其中，A =(a ij )n 为非奇异阵，X =(x 1,x 2,…,x n )T , b =(b 1,b 2,…,b n )T将系数阵A 分解为: A =U +D +L ,U 为上三角矩阵，D 为对角矩阵，L 为下三角矩阵.于是AX =b 可改写为 (U +D +L )X =b⇔ X =D -1b -D -1(U +L )X由此可得矩阵形式的Jocobi 迭代公式： X k +1=BX (k )+f □§3.2 高斯-塞德尔Gauss-Seidel 迭代法注意到利用Jocobi 迭代公式计算)1(+k ix 时，已经计算好)(1)(2)(1,,,k i k k x x x - 的值，而Jocobi 迭代公式并不利用这些最新的近似值计算，仍用)(1)(2)(1,,,k i k k x x x - .这启发我们可以对其加以改进,即在每个分量的计算中尽量利用最新的迭代值，得到),,2,1()(1111)1()1(n i x a x a b a xn i j k jij i j k j ij i iik i=∑-∑-=+=-=++上式称为Gauss-Seidel 迭代法. 其矩阵形式是X =-(D +L )-1UX +(D +L )-1b , X k +1=BX (k )+f .迭代次数 x1 x2 x3 0 0 0 0 1 0.72 0.902 1.1644 2 1.04308 1.167188 1.282054 3 1.09313 1.195724 1.2977714 1.099126 1.199467 1.2997195 1.09989 1.199933 1.2999656 1.099986 1.199992 1.2999967 1.099998 1.199999 1.2999998 1.1 1.2 1.3§3.3 超松驰迭代法SOR 方法1基本思想:逐次超松弛迭代法(Successive Over Relaxation Method,简写为SOR)可以看作带参数ω的高斯-塞德尔迭代法，是G-S 方法的一种修正或加速.是求解大型稀疏矩阵方程组的有效方法之一. 2 SOR 算法的构造:设方程组AX =b , 其中，A =(a ij )n 为非奇异阵，X =(x 1,x 2,…,x n )T , b =(b 1,b 2,…,b n )T . 假设已算出x (k )，),,2,1()(1111)1()1(n i x a x a b a xn i j k j ij i j k j ij i iik i=∑-∑-=+=-=++ (1)相当于用高斯-塞德尔方法计算一个分量的公式. 若对某个参数ω，作)1(+k ix与)(k i x 加权的平均,即)()1()()1()()1()(1k i k ik i k ik ik ix xx xxx-+=+-=+++ωωω (2)其中，ω称为松弛因子.用(1)式代入(2)式，就得到解方程组AX =b 的逐次超松弛迭代公式：⎪⎩⎪⎨⎧=∑-∑-=∆∆+==-=++),,2,1()()(11)1()()1(n i x a x a b a x x x x n ij k j ij i j k j ij i iii i k i k i ω (3) 显然，当取ω=1时，式(3)就是高斯-塞德尔迭代公式. 3 例题分析:利用SOR 方法解方程组⎪⎩⎪⎨⎧=+---=-+-=--3322242024321321321x x x x x x x x x (1) 其准确解为X *={1, 1, 2}. 建立与式(1)相等价的形式：⎪⎪⎩⎪⎪⎨⎧++=-+=+=132315.05.05.025.05.021*******x x x x x x x x x (2) 据此建立迭代公式：⎪⎪⎩⎪⎪⎨⎧++=-+=+=+++132315.05.05.025.05.0)(2)(1)1(3)(3)(1)1(23)(2)1(1k k k k k k kk k x x x x x x x x x (3)利用SOR 算法，取迭代初值1)0(3)0(2)0(1===x x x ,ω=1.5,迭代结果如下表.逐次超松弛迭代法次数 x1 x2 x3 1 0.625000 0.062500 1.750000 2 0.390625 0.882813 1.468750 3 1.017578 0.516602 1.8085944 0.556885 0.880981 1.7104495 1.023712 0.743423 1.8681036 0.746250 0.908419 1.8387377 0.997715 0.860264 1.9138948 0.864050 0.936742 1.9086059 0.986259 0.922225 1.94552310 0.928110 0.958649 1.94749311 0.985242 0.955944 1.96619812 0.961661 0.973818 1.96952113 0.988103 0.974699 1.97928914 0.979206 0.983746 1.98217215 0.991521 0.985318 1.98741616 0.988509 0.990038 1.98951317 0.994341 0.991414 1.99239718 0.993538 0.993946 1.99380619 0.996367 0.994950 1.99542420 0.996313 0.996342 1.99633121 0.997724 0.997018 1.99725422 0.997871 0.997798 1.99782223 0.998596 0.998234 1.998355GS迭代法须迭代85次得到准确值X*={1, 1, 2};而SOR方法只须55次即得准确值.由此可见，适当地选择松弛因子ω，SOR法具有明显的加速收敛效果. □§3.4 迭代法的收敛性1. 向量和矩阵范数 (a) 向量范数R n 空间的向量范数 || · || ,对任意n R y x ∈,, 满足下列条件：00||||;0||||)1(=⇔=≥x x x (正定性)||||||||||)2(x x⋅=αα (齐次性)||||||||||||)3(y x y x+≤+ （三角不等式）常见的向量范数有： (1) 列范数：(2) 谱范数：(欧几里德范数或向量的长度,模)(3) 行范数：(4) p 范数:上述范数的几何意义是:∞||||x =max(|x 2-x 1|,|y 2-y 1|) ; 1||||x =|x 2-x 1|+|y 2-y 1| ;2122122)()(||||y y x x x -+-=.向量序列}{)(k x依坐标收敛于向量x * 的充要条件是向量序列}{)(k x 依范数收敛于向量x *,即0||||lim *)(=-∞→x x k k .(b) 矩阵范数n m R ⨯空间的向量范数 || ·|| ,对任意 n m R B A ⨯∈,, 满足下列条件：|||||||| || AB || (4)||||||||||||)3(||||||||||)2(00||||;0||||)1(B A B A B A A A A A A ≤+≤+⋅==⇔=≥αα常见的矩阵范数有：∑==∞≤≤nj ij a A ni 1||max ||||1 (行和范数)∑==≤≤ni ij a A nj 11||max ||||1 (列和范数))(||||m ax 2A A A T λ= (谱范数)若A 对称，则有)()(2m ax m ax A A A T λλ=.矩阵A 的谱半径记为)(||||2A A ρ=,ρ(A ) =||m ax 1i ni λ≤≤,其中λi 为A 的特征根。

消去法等解线性方程组（一）中文摘要Gauss消去法是目前求解中小规模线性方程组（即阶数不要太高，例如不超过1000）常用的方法，它一般用于系数矩阵稠密（即矩阵的绝大多数元素都是非零的）而没有任何特殊结构的线性方程组。

本文对GAUSS、列主元及完全主元消去法这几种方法进行了比较。

求解n阶线性方程组Ax=b的基本思想是在逐步消元的过程中，把方程组的系数矩阵化为上三角矩阵，从而将原方程组约化为容易求解的等价三角方程组，然后进行回代求解。

关键词：Gauss消去法； C++编程ELIMINATION METHOD FOR SOLVING LINEAREQUATIONS (一)ENGLISH ABSTRACTGauss elimination method is the small size of solving linear equations (i.e. the order is not too high, for example, not more than 1000) commonly used method, it is generally used for the coefficient matrix dense (ie, most of the matrix elements are nonzero) without any special structure of linear equations. In this paper, GAUSS, the main element and completely Elimination PCA these methods were compared. Solving n linear equations Ax = b is the basic idea of the gradual elimination process, the equations of the coefficient matrix into upper triangular matrix, thus reduced to the original equations easy to solve trigonometric equations equivalence group then back substitution to solve.Keywords: Gauss elimination method; C + + Programming目录消去法等解线性方程组（一） (1)目录 (3)一、问题背景 (4)二、问题提出 (5)三、数学原理与算法分析 (5)3.1不选主元Gauss消去法 (5)3.1.1不选主元Gauss消去法数学原理 (5)3.1.2不选主元Gauss消去法算法步骤 (6)3.2列主元Gauss消去法 (7)3.2.1列主元Gauss消去法数学原理 (7)3.2.1列主元Gauss消去法算法描述 (8)3.3完全主元Gauss消去法 (8)3.3完全主元Gauss消去法数学原理 (8)3.3.1算法描述如下 (10)3.4计算结果及敏度分析 (11)四、其他方法 (13)五、参考文献 (14)六、附录 (15)一、问题背景在古代数学和近代数学数值计算和工程应用中，求解线性方程组都是重要的课题。

解线性方程组的迭代法Haha送给需要的学弟学妹摘要:因为理论的分析表明，求解病态的线性方程组是困难的，但是实际情况是否如此，需要我们来具体检验。

系数矩阵H 为Hilbert 矩阵，是著名的病态问题。

因而决定求解Hx b =此线性方程组来验证上述问题。

详细过程是通过用Gauss 消去法、J 迭代法、GS 迭代法和SOR 迭代法四种方法求解Hx b =线性方程组。

关键词：病态方程组、Gauss 消去法、J 迭代法、GS 迭代法、SOR 迭代法目录：一、问题背景介绍二、建立正确额数学模型 三、求解模型的数学原理1、Gauss 消去法求解原理2、Jacobi 迭代法求解原理3、G-S 迭代法求解原理4、SOR 迭代法求解原理5、Jacobi 和G-S 两种迭代法收敛的充要条件 四、计算过程（一）Hilbert 矩阵维数n=6时1、Gauss 消去法求解2、Jacobi 迭代法求解3、G-S 迭代法求解4、SOR 迭代法求解（二）Hilbert 矩阵维数n=20、50和100时1、G-S 迭代法求解图形2、SOR 迭代法求解图形 五、编写计算程序 六、解释计算结果1、Gauss 消去法误差分析2、G-S 迭代法误差分析3、SOR 迭代法误差分析G-S 迭代法与SOR 迭代法的误差比较 七、心得体会正文：一、问题背景介绍。

理论的分析表明，求解病态的线性方程组是困难的。

实际情况是否如此，会出现怎样的现象呢？二、建立正确的数学模型。

考虑方程组Hx b =的求解，其中系数矩阵H 为Hilbert 矩阵，,,1(), , ,1,2,,1i j n n i j H h h i j n i j ⨯===+-这是一个著名的病态问题。

通过首先给定解（为方便计算，笔者取x 的各个分量等于1），再计算出右端,b Hx =这样Hx b =的解就明确了，再用Gauss 消去法、J 迭代法、GS 迭代法和SOR 迭代法四种方法分别求解,Hx b =将求解结果与给定解比较，而后求出上述四种方法的误差，得出哪种方法比较好。

《数值代数》课程教学大纲Numerical Algebra课程代码：课程性质：专业方向理论课庞修适用专业：信息计算开课学期：7总学时数：48总学分数：3编写年月：200徉阴修订年月：200强阴执笔：徐圣兵一、课程的性质和目的本课程是信息与计算科学专业信息计算方向的一门专业方向选修课，它是科学与工程计算的核心。

本课程的主要内容就是，如何针对各类科学与工程问题所提出的矩阵计算问题的特点，设计出相应的快速可靠的算法。

要求学生根据计算机的特点研究计算时间最短、需要计算机内存最少的矩阵计算方法，并具备对矩阵计算问题的理论进行研究和探讨的能力，为参加大型科学工程计算实践打下必要的基础二、课程教学内容及学时分配第一章绪论（2学时）本章要求懂得数值代数研究的基本问题，了解研究数值方法的必要性，初步了解矩阵分解、敏度分析与误差分析、算法复杂性与收敛速度等数值代数的基本概念。

重点是要懂得数值代数所研究的基本问题。

本章知识点为：矩阵分解，敏度分析与误差分析，条件数，数值稳定性，算法复杂性与收敛速度。

_第二章线性方程组的直接解法（8学时，其中实验2学时）本章要求掌握解线性方程组的最基本的直接解法一Gauss消去法，懂得三角分解与选主元三角分解、平方根法、分块三角分解等算法。

重点掌握列主元三角分解算法。

本章知识点为：前代法，回代法，Gauss变换，三角分解，选主元三角分解，平方根法，改进平方根法，分块三角分解。

为进一步加强对线性方程组直接法求数值解的理解，安排实验一次:直接法求解线性方程组。

第三章线性方程组的敏度分析与消去法的舍入误差分析（8学时，其中实验2学时）本章要求掌握向量范数和矩阵范数的概念及其基本性质，懂得对线性方程组进行敏度分析和舍入误差讨论，并对列主元Gauss消去法进行详细的舍入误差分析，最后就是了解一种估计计算解的精度的实用方法以及改进其计算精度的迭代方法。

重点是对线性方程组进行敏度分析和舍入误差分析。

数值分析试题及答案一、选择题（每题3分，共30分）1. 下列关于数值分析的说法，错误的是（）。

A. 数值分析是研究数值方法的科学B. 数值分析是研究数值方法的数学理论C. 数值分析是研究数值方法的误差分析D. 数值分析是研究数值方法的数学理论、误差分析及数值方法的实现答案：B2. 在数值分析中，插值法主要用于（）。

A. 求解微分方程B. 求解积分方程C. 求解线性方程组D. 通过已知数据点构造一个多项式答案：D3. 线性方程组的解法中，高斯消元法属于（）。

A. 直接方法B. 迭代方法C. 矩阵分解方法D. 特征值方法答案：A4. 牛顿法（Newton's method）是一种（）。

A. 插值方法B. 拟合方法C. 迭代方法D. 优化方法答案：C5. 在数值分析中，下列哪种方法用于求解非线性方程的根？A. 高斯消元法B. 牛顿法C. 雅可比方法D. 斯托尔-温格尔方法答案：B6. 下列关于误差的说法，正确的是（）。

A. 绝对误差总是大于相对误差B. 相对误差总是小于绝对误差C. 误差是不可避免的D. 误差总是可以消除的答案：C7. 在数值分析中，下列哪个概念与数值稳定性无关？A. 条件数B. 截断误差C. 舍入误差D. 插值多项式的阶数答案：D8. 用泰勒级数展开函数f(x)=e^x，下列哪一项是正确的？A. f(x) = 1 + x + x^2/2! + x^3/3! + ...B. f(x) = 1 - x + x^2/2! - x^3/3! + ...C. f(x) = x + x^2/2 + x^3/6 + ...D. f(x) = x - x^2/2 + x^3/6 - ...答案：A9. 插值多项式的次数最多为（）。

A. n-1B. nC. n+1D. 2n答案：B10. 下列关于数值积分的说法，错误的是（）。

A. 梯形法则是一种数值积分方法B. 辛普森法则是一种数值积分方法C. 龙格法则是数值积分方法中的一种D. 数值积分方法总是精确的答案：D二、填空题（每题3分，共15分）1. 在数值分析中，条件数是衡量问题的______。

第2章线性方程组的解法--------学习小结一、本章学习体会本章主要学习的是线性方程组的解法。

而我们则主要学习了高斯消去法、直接三角分解法以及迭代法三种方法。

这三种方法的优缺点以及适用范围各有不同。

高斯消去法中，我们又学习了顺序高斯消去法以及列主元素高斯消去法。

顺序高斯消去法可以得到方程组的精确解，但要求系数矩阵的主对角线元素不为零，而且该方法的数值稳定性没有保证。

但列主元素高斯消去法因为方程顺序的调整，其有较好的数值稳定性。

直接三角分解法中，我们主要学习了Doolitte分解法与Crout分解法。

其思想主要是：令系数矩阵A=UL，其中L为下三角矩阵，U是上三角矩阵，为求AX=b 的解，则引进Ly=b,Ux=y两个方程，以求X得解向量。

这种方法计算量较小，但是条件苛刻，且不具有数值稳定性。

迭代法（逐次逼近法）是从一个初始向量出发，按照一定的计算格式，构造一个向量的无穷序列，其极限才是所求问题的精确解，只经过有限次运算得不到精确解。

该方法要求迭代收敛，而且只经过有限次迭代，减少了运算次数，但是该方法无法得到方程组的精确解。

二、本章知识梳理针对解线性方程组，求解线性方程组的方法可分为两大类：直接法和迭代法，直接法（精确法）：指在没有舍入误差的情况下经过有限次运算就能得到精确解。

迭代法（逐次逼近法）：从一个初始向量出发，按照一定的计算格式，构造一个向量的无穷序列，其极限才是所求问题的精确解，只经过有限次运算得不到精确解。

我们以前用的是克莱姆法则，对于计算机来说，这种方法运算量比较大，因此我们学习了几种减少运算次数的方法，有高斯消去法、直接三角分解法，同时针对病态方程组，也提出了几种不同的解法。

2.1 Gauss消去法Gauss消去法由消元和回代两个过程组成，消元过程是指针对方程组的增广矩阵，做有限次初等行变化，使它系数矩阵变为上三角矩阵。

2.1.1顺序Gauss消去法消元过程：对于K=1,2,3…,n-1执行(1)如果,则算法失效，停止计算；否则转（2）(2)对于计算回代过程：综上：顺序Gauss消去法的数值稳定性是没有保证的。

线性方程组的敏度分析与消去法的舍入误差分析我们已经在上一章详细讨论了如何用Ｇauss消去法求解线性方程组，但从数值计算的角度来看，这是不够的。

这是因为在实际计算时，我们不仅要面对如何求解的问题，而且也要面对数据不精确的问题和机器的有限精度问题。

只有讨论了后面这两个问题，我们才可能知道，所求得的解是否可靠。

而要讨论这两个问题，我们需要用范数来描述向量和矩阵的扰动的大小等概念。

所以在这一章中，我们首先介绍向量范数和矩阵范数的概念及其基本性质；然后对线性方程组进行敏度分析，讨论舍入误差问题，并对列主元Ｇauss消去法进行详细的舍入误差分析；最后介绍一种估计计算精度的实用方法以及改进其计算精度的迭代方法。

2.1 向量范数与矩阵范数2.1.1 向量范数向量范数的概念是复数模的概念的自然推广，其定义如下：定义2.1.1一个从到的非负函数叫做上的向量范数，如果满足：（1）正定性：对所有的有，而（2）且当且仅当；（3）齐次性：对所有的和有；（4）三角不（5）等式：对所有的有.定理 2.1.1 上的向量范数是连续函数。

证明对任意，因，由性质（3）得，即，及从而得到再注意到，由性质（2）和（3）得由此可知，作为上的实函数是连续的。

定理证毕最常用的向量范数是P范数（亦称H?lder范数）：其中是最重要的，即它们分别叫做1范数、2范数和∞范数。

这三个范数的正定性与齐次性是显然的，而且也容易证明和满足三角不等式。

对于2范数要证明三角不等式成立，需要用到Cauchy-Schwartz不等式（其证明参见习题2.1）这个不等式是H?lder不等式的特殊情形。

事实上，利用Cauchy-Schwartz不等，我们有由此即知2范数满足三角不等式。

尽管在上可以引进各种各样的范数，但在下面定理所述的意义下所有这些范数都是等价的。

定理2.1.2设和是上任意二个范数，则存在正常数和，使得对一切有证明分两步来证明这个结论。

首先来证明和等价，即证明存在常数，使对一切且有考虑泛函记，则是一个有界闭集。