数值分析复习资料

1、数值分析复习资料一、重点公式第一章 非线性方程和方程组的数值解法1）二分法的基本原理，误差：2）迭代法收敛阶：，若则要求3）单点迭代收敛定理：定理一：若当时，且，则迭代格式收敛于唯一的根；定理二：设满足：时，则对任意初值迭代收敛，且：定理三：设在的邻域内具有连续的一阶导数，且，则迭代格式具有局部收敛性；定理四：假设在根的邻域内充分可导，则迭代格式是P阶收敛的（Taylor展开证明）4）Newton迭代法：，平方收敛5）Newton迭代法收敛定理：设在有根区间上有二阶导数，且满足：；：；：初值使得；则Newton迭代法收敛于根。6）多点迭代法：收敛阶：7）Newton迭代法求重根（收敛仍为线性收

2、敛），对Newton法进行修改：已知根的重数r，（平方收敛）：未知根的重数：，为的重根，则为的单根。8）迭代加速收敛方法：当不动点迭代函数在的某个邻域内具有二阶导数，平方收敛9）确定根的重数：当Newton迭代法收敛较慢时，表明方程有重根10）拟Newton法其中11）秩1拟Newton法：Broyden秩1方法第二章 线性代数方程组数值解法1）向量范数：非负性：，且的充要条件是；：齐次性：三角不等式：1范数：2范数：范数：p范数：2）矩阵范数：非负性：，且的充要条件是；：齐次性：三角不等式：乘法不等式：F范数：1范数：，列和最大范数：，行和最大2范数：，其中，为的特征值，3）Gauss消元法

3、（上三角阵）：；Gauss-Jordan消元法（对角阵）：；列选主元消元法：在消元之前进行行变换，将该列最大元素换置对角线主元位置；（可用于求逆矩阵）全选主元消元法：全矩阵搜索矩阵最大元素进行行变换和列变换至其处于对角线主元位置；4）三角分解法：Doolittle分解法：A=LU，L单位下三角阵，U上三角阵：Crout分解法：A=LU，L下三角阵，U单位上三角阵：Cholesky分解法：A对称正定，L为单位下三角阵：改进的Cholesky分解法：A对称正定，L为单位下三角阵，D为对角阵：追赶法：Crout分解法解三对角方程5）矩阵的条件数，谱条件数：6）如果，则为非奇异阵，且7）迭代法基本原理

4、：迭代法：(，迭代格式收敛)：至少存在一种矩阵的从属范数，使8）Jacobi迭代：9）Gauss-Seidel迭代：10）超松弛迭代法11）二次函数的一维搜索：12）最速下降法：选择方向进行一维搜索：，其中13）共轭梯度法：第一步：最速下降法，第二步：过选择的共轭方向，其中，过以为方向的共轭直线为，进行二次函数的一维搜索14）一般的共轭梯度法：第三章 插值法与数值逼近1）Lagrange插值：，余项：2）Newton插值：差商表余项3）反插值4）Hermite插值（待定系数法）其中余项：5）分段线性插值：插值基函数：余项：分段余项6）有理逼近：反差商表有理逼近函数式：7）正交多项式的计算：定理

5、：在上带权函数的正交多项式序列，若最高项系数唯一，它便是唯一的，且由以下的递推公式确定其中定理3.88）连续函数的最佳平方逼近：在上，法方程为，其中，均方误差：最大误差：9）离散函数的最佳平方逼近（曲线的最小二乘拟合）：法方程其中第四章 数值积分1）代数精度的概念及应用：对r次多项式的精确成立，以及代入法求解系数。2）Lagrange插值代入Lagrange插值基函数，其中误差：定理：数值积分公式具至少有n次代数精度其是差值型的3）等距节点的Newton-Cotes公式将拉格朗日差值积分公式中的差值节点即可，其中；，令（Cotes系数）则：N-C公式的数值稳定性：当同号时是稳定的，否则不稳定，

6、（其中）N-C公式至少具有n次代数精度，若n为偶数，则其代数精度可提高到n+1次；余项：当n为偶数时，当n为奇数时，4）复化的N-C公式复化的梯形公式：将积分区间n等分，然后在每个区间上应用梯形公式复化的Simpson公式：将积分区间n等分，然后在每个区间上应用Simpson公式5）Romberg积分法逼近的阶为 6）求积节点为n+1的机械求积公式的代数精度=2n+1；7）Gauss求积公式在a，b上与所有次数m1时，f(n+1) (x)=0，Rn(x)=0，所以 注意：对于次数不超过n的多项式,利用上结果，有 = =上式正是Qn(x)的拉格朗日插值多项式。可见，Qn(x)的拉格朗日插值多项式

7、就是它自身，即次数不超过n的多项式在n+1个互异节点处的拉格朗日插值多项式就是它自身。例5 已知数据如表的第2，3列，试用直线拟合这组数据。 解 计算列入表中。n=5。a0,a1满足的法方程组是kxkykxkyk11414224.5493369184481632558.52542.5S153155105.5 解得a0=2.45, a1=1.25。所求拟合直线方程为 y=2.45+1.25x例6选择填空题1. 设y=f(x), 只要x0,x1,x2是互不相同的3个值，那么满足P(xk)=yk(k=0,1,2)的f(x)的插值多项式P(x)是 (就唯一性回答问题)答案：唯一的3. 拉格朗日插值多项

8、式的余项是( ),牛顿插值多项式的余项是( ) (A) (B) f(x,x0,x1,x2,xn)(xx1)(xx2)(xxn1)(xxn) (C) (D) f(x,x0,x1,x2,xn)(xx0)(xx1)(xx2)(xxn1)(xxn)答案：(A)，(D)。见教材有关公式。 第四章典型例题例1 试确定求积公式的代数精度。依定义，对xk(k=0,1,2,3,)，找公式精确成立的k数值解 当f(x)取1,x,x2,时，计算求积公式何时精确成立。(1) 取f(x)=1,有左边, 右边(2) 取f(x)=x,有 左边, 右边(3) 取f(x)=x2,有 左边=, 右边=(4) 取f(x)=x3,有

9、 左边=, 右边=(5) 取f(x)=x4,有 左边=, 右边=当k3求积公式精确成立，而x4公式不成立，可见该求积公式具有3次代数。例5 试确定求积公式中的参数a，并证明该求积公式具有三次代数精度。解 公式中只有一个待定参数a。当f(x)=1,x时，有,即h=h ,不能确定a，再令f(x)=x2, 代入求积公式，得到,即 得. 求积公式为将f(x)=x3代入上求积公式，有 可见，该求积公式至少具有三次代数精度。再将f(x)=x4代入上公式中，有 所以该求积公式具有三次代数精度。例6 选择填空题1. 牛顿科茨求积公式与高斯型求积公式的关键不同点是 。解答：牛顿科茨求积公式的节点和求积系数确定后

10、，再估计其精度；高斯型求积公式是由精度确定其节点和求积系数。第五章典型例题例1 证明方程1xsinx0在区间0,1内有一个根，使用二分法求误差不超过0.5104的根要迭代多少次？证明 令f(x)1xsinx f(0)=10，f(1)=sin10(x0，1),故f(x)0在区间0，1内有唯一实根。给定误差限e0.5104，有只要取n14。例2 用迭代法求方程x54x20的最小正根。计算过程保留4位小数。分析 容易判断1，2是方程的有根区间。若建立迭代格式，此时迭代发散。建立迭代格式，此时迭代收敛。解 建立迭代格式 (可任取1，2之间的值)1.431 0 1.505 1 1.516 5 1.518

11、 2 1.5185 取1.5185例3 试建立计算的牛顿迭代格式，并求的近似值，要求迭代误差不超过105分析首先建立迭代格式。确定取几位小数,求到两个近似解之差的绝对值不超过105。解 令，求x的值。牛顿迭代格式为迭代误差不超过105，计算结果应保留小数点后6位。当x=7或8时，x3=343或512，,取x0=8,有 7.478 0787.439 956 7.4397607.439760于是，取7.439760例4 用弦截法求方程x3x210，在x=1.5附近的根。计算中保留5位小数点。分析 先确定有根区间。再代公式。解 f(x)= x3x21，f(1)=1，f(2)=3，有根区间取1,2取x

12、1=1, 迭代公式为(n=1,2,) 1.37662 1.48881 1.46348 1.46553取1.46553，f(1.46553)0.000145例4 选择填空题1. 设函数f(x)在区间a,b上连续，若满足 ，则方程f(x)=0在区间a,b一定有实根。答案：f(a)f(b)04牛顿切线法是用曲线f(x)上的 与x轴的交点的横坐标逐步逼近f(x)0的解；而弦截法是用曲线f(x)上的 与x轴的交点的横坐标逐步逼近f(x)0的解。答案：点的切线；两点的连线解答：见它们的公式推导.三、复习资料一、 用牛顿法求解方程的解，收敛精度二、 应用列主元素消元法求解的解，保留4位有效数字三、 应用雅可

13、比迭代法求解下列方程组，收敛精度求解过程如表：012345600.720.9711.0571.08531.09511.098300.831.0701.15711.18531.19511.198300.841.1501.24821.28281.29411.29800.840.310.09820.03460.01130.0039四、 设，应用幂法求解其最大特征值，以及对应的特征向量五、 用拉格朗日差值法构造三次多项式，求解处的函数值x0.460.470.480.49y0.4846550.4937450.502750.511668要求小数点后4位六、 已知如下表的函数，试用最小二乘法求二次多项式来拟

14、合这组数据x-1.00-0.75-0.5-0.2500.250.50.751y-0.22090.32950.88261.43922.00032.56453.13343.70614.2836七、 用龙贝格求积公式计算积分的近似值，要求收敛精度八、 用预测-校正的改进欧拉法求解如下常微分方程，取步长为0.1，计算到1.0各阶段计算结果如下：00.10.20.30.40.5000.01600.04370.08410.137800.00550.02190.05010.09090.14500.60.70.80.91.00.20550.28770.38890.49760.62630.21290.29540

15、.39290.50590.6348四、复习总结 第二章 数值分析基本概念教学内容： 1. 误差与有效数字误差、误差限、相对误差、相对误差限和有效数字的定义及相互关系；误差的来源和误差的基本特性；误差的计算（估计）的基本方法。 2. 算法的适定性问题数值分析中的病态和不稳定性问题；病态问题和不稳定算法的实例分析。3. 数值计算的几个注意问题 数值计算的基本概念l 误差概念和分析误差的定义：设x是精确值，p是近似值，则定义两者之差是绝对误差: 由于精确值一般是未知的,因而不能求出来,但可以根据测量误差或计算情况估计它的上限相对误差定义为绝对误差与精确值之比l 误差的来源：舍入误差将无限位字长的精确

16、数处理成有限位字长近似数的处理方法称为舍入方法。带来舍人误差。截断误差用数值法求解数学模型时，往往用简单代替复杂,或者用有限过程代替无限过程所引起的误差。 l 有效数字对于a=a0 a1 am . am+1 am+n(a00)的近似数， 若|0.5x10-n，则称a为具有m+n+1位有效数字的有效数，其中每一位数字都叫做a的有效数字。有效数和可靠数的最末位数字称为可疑数字有效数位的多少直接影响到近似值的绝对误差与相对误差的大小。推论1 对于给出的有效数，其绝对误差限不大于其最末数字的半个单位。推论2 对于给出的一个有效数，其相对误差限可估计如下：例:计算y = ln x。若x 20，则取x的几

17、位有效数字可保证y的相对误差 0.1% ?l 数值计算的算法问题“良态”问题和“病态”问题在适定的情况下，若对于原始数据很小的变化X，对应的参数误差y也很小，则称该数学问题是良态问题；若y很大，则称为病态问题。病态问题中解对于数据的变化率都很大，因此数据微小变化必将导致参数模型精确解的很大变化。数学问题的性态完全取决于该数学问题本身的属性，在采用数值方法求解之前就存在，与数值方法无关。稳定算法和不稳定算法如果用数值方法计算时,误差在计算过程中不扩散的算法称为稳定算法。否则称为不稳定算法。l 数值计算应注意的问题避免相近二数相减；避免小分母；避免大数吃小数；选用稳定的算法。绝对误差的运算：第三章

18、 线性方程组求解的数值方法教学内容：1. 高斯消元法消元法的实现过程；主元问题。2. 矩阵分解矩阵LU分解的一般计算公式；利用LU分解的线性方程组求解方法；Cholesky分解；Matlab的Cholesky分解函数。3. 向量范数与矩阵范数向量范数及其性质；矩阵函数及其性质；常用范数形式。4. 线性方程组的迭代法求解迭代求解的思路；Jacobi迭代法；高斯_赛德尔迭代法；迭代法的收敛性。 5. 方程组的病态问题与误差分析线性方程组解的误差分析；条件数和方程组的病态性。消元法：问题：消去法是按照系数矩阵的主对角线上的元素（主元）进行消元。从而可能出现：（1）某个主元为零，导致消元过程无法进行。

19、（2）当某个主元的绝对值很小时，计算结果误差很大。定理：若A的所有顺序主子式 均不为0，则高斯消元无需换行即可进行到底，得到唯一解。全主元消去法每一步选绝对值最大的元素为主元素。列主元消去法省去换列的步骤，每次仅选一列中最大的元。矩阵三角分解法 计算公式：算法：Cholesky分解：定理：设矩阵A对称正定，则存在非奇异下三角阵L使得 。若限定L对角元为正，则分解唯一。Matlab中的Cholesky分解函数：chol（）向量和矩阵的范数为了研究线性方程组近似解的误差估计和迭代法的收敛性，引进向量（矩阵）的范数的概念。向量范数定义：空间的向量范数 | | 对任意 满足下列条件：(齐次性) （三角

20、不等式)（正定性)常用范数：矩阵范数定义：空间的向量范数 | | 对任意 满足下列条件： (4)* | AB | | A | | B |常用矩阵范数：Frobenius范数：由向量范数 | |p导出关于矩阵 A Rnn的p范数:特别有：（行和范数）（列和范数）（谱范数 ）谱半径：矩阵A的谱半径记为r (A) =，其中li为A的特征根。定理：对任意算子范数 | | 有定理：若A对称，则有定理：若矩阵B对某个算子范数满足 |B| 1，则必有解线性方程组的迭代法 研究内容：l 如何建立迭代公式？ l 收敛速度？l 向量序列的收敛条件？l 误差估计？思路：收敛问题：雅可比(Jacobi)迭代法高斯塞德

21、尔迭代法迭代法的收敛性谱半径小于1. 迭代法的误差估计：误差分析：问题的提出：l b有扰动，A无扰动l A有扰动，b有扰动l A有扰动，b有扰动定义：条件数: cond(A) = |A-1 | |A|条件数的性质： 结论：当条件数很大时,方程组 Ax = b是病态问题;当条件数较小时,方程组 Ax = b是良态问题。第四章 函数的数值逼近1. 代数多项式插值问题插值多项式的存在唯一性；插值基函数和插值多项式的一般形式；插值的误差分析；多项式插值的Runge现象。2. 分段低次插值分段线性插值；Hermite插值和分段Hermite插值。3. 三次样条插值样条插值的定义；三次样条函数的计算；Ma

22、tlab中的插值函数。4. 曲线拟合的最小二乘法曲线拟合的最小二乘法法；多项式拟合方法；Matlab中的多项式拟合函数；5. 最佳平方逼近权内积；正交多项式的最佳平方逼近。插值问题：函数解析式未知，或计算复杂，用函数p(x)去近似代替它，使得p(xi)= yi (i=0,1,2,n) 函数p(xi)称为插值函数。 x0,x1, xn称为插值节点或简称节点。 插值节点所界的区间称为插值区间。 p(xi)= yi 称为插值条件。多项式的插值问题构造n次多项式Pn(x)= a0 + a1x + a2x2+ anxn使满足 Pn(xi)= yi (i=0,1,2,n)，讨论的主要内容：l 如何求出插值函数；l 插值函数的存在性；l 收敛性和误差估计。拉格朗日插值插值多项式的存在唯一性：结论通过n+1个节点的n阶插值多项式唯一存在。一次基函数 二