数值计 课件.pptx

何光辉 董海云 魏曙光 编 王开荣 主审 ;1.1 算 法 (1)有穷性 (2)确定性 (3)可行性 (4)有输入 (5)有输出;1.1.1 算法常用描述形式 (1)数学公式和文字说明描述 (2)框图描述 (3)伪代码描述 (4)程序语言描述 1.1.2 数值型算法的基本特点 1)无穷过程的截断 2)连续过程的离散化 3)迭代计算;图1.1;1.2 误 差 1.2.1 误差的来源 ①截断误差 ②舍入误差 1.2.2 误差的基本概念 1)误差与误差限 定义1.1 设x*是精确值，x是它的一个近似值，称 e=x-x*;为近似值x的绝对误差，简称误差。 2)相对误差与相对误差限 定义1.2 绝对误差与精确值的比值 称为x 的相对误差。相对误差是无量纲的量，常用百分比表示。 1.2.3 有效数字;定义1.3 若x的误差限ε是它某一数位的半个单位，就称x准确到该位，从这一位起直到前面第一个非零数字为止的所有数字称为x的有效数字。 定理1.1 设近似值x=±0.a1a2…an×10m, 有n位有效数字,则其相对误差限 定理1.2 设近似值x=±0.a1a2…an×10m的相对误差限不大于;则它至少有n位有效数字。 1.3 设计数值型算法的基本原则 1.3.1 函数计算的误差传播 1.3.2 四则运算中误差的传播 1)避免两个相近数相减 2)避免除法中除数的数量级远小于被除数 3)合理安排运算顺序 4)注意运算步骤的简化;2.1 消元法 2.1.1 Gauss消元法 1.消元 2.回代 定理2.1 Ax=b能用Gauss消元法求解的充要条件是A的各阶顺序主子式不为零。;2.1.2 列主元Gauss消元法 2.1.3 Gauss—Jordan消元法 Gauss—Jordan消元法是Gauss消元法的一种变形。;2.2 矩阵三角分解 设A∈Rn×n，若存在下三角矩阵L和一个上三角阵U使得A=LU，则称LU为矩阵A的三角分解。 定理2.2 若矩阵A的各阶顺序主子式非零，则A可进行Doolittle分解和Crout分解，且这种分解是唯一的。 2.2.1 Doolittle分解 Doolittle算法实际上就是Gauss消元法的另一种形式。;2.2.2 Crout分解 2.2.3 平方根法 当A是对称正定矩阵时，矩阵的三角分解就具有其特殊形式。 定理2.3 设A对称正定，则存在三角分解A=LLT，L是非奇异下三角矩阵，且当限定L的对角元为正时，这种分解是唯一的，称为对称正定矩阵的Cholesky分解。 2.2.4 追赶法 ;(1)平方根法 (2)追赶法 2.3 向量和矩阵的范数 2.3.1 向量范数 定义2.1 x和y是Rn中的任意向量，向量范数‖·‖是定义在Rn上的实值函数 2.3.2 矩阵范数 定义2.2 设A是n×n矩阵，x∈Rn，定义;为矩阵A的范数。 定理2.4(范数等价) 对任意给定的两种范数‖A‖α， ‖A‖β有下列关系： 其中的m,M是正的常数， ‖A‖α表示向量(或矩阵)的α范数。 定义2.3 如果向量x是精确值,x(k)是它的一个 近似值,‖x(k)-x‖是x(k)对x的误差, 是x(k)对x的相对误差。;定义2.4 如果 ‖x(k)-x‖=0，称Rn中的向量序列{x(k)}收敛于Rn中的向量x。 定义2.5 如果 ‖A(k)-A‖=0，称n×n矩阵序列{A(k)}收敛于n×n矩阵A。 定理2.5 Rn中的向量序列{x(k)} 收敛于Rn中的向量x的充要条件是： x(k)j=xj, j=1,2,…,n，其中x(k)j和xj是x(k) 和x中的第j个分量。 定理2.6 n×n矩阵序列{A(k)} 收敛于n×n矩阵A的充要条件是：; a(k)ij=aij, i,j=1,2,…,n，其中a(k)ij和aij 分别是A(k) 和A在(i,j)位置上的元素。 2.3.3 谱半径 定义2.6 设n×n矩阵A的特征值为λi,i=1,2,…,n，则称ρ(A)= |λi|为A的谱半径。 定理2.7 设矩阵A谱半径为ρ(A)，则矩阵谱半径与范数满足：ρ(A)≤‖A‖。;定理2.8 设A是n×n阶矩阵， A的各次幂组成的矩阵序列I，A,A2,…,Ak,…收敛于零，即 Ak=0的充要条件是ρ(A)1。 2.3.4 条件数及病态方程组 定义2.7 若n×n方阵A非奇异,则称‖A-1‖ ×‖A‖为A的条件数，记为 2.4 线性方程组的迭代法解法 2.4.1 Jacobi迭代法;2.4.2 Gauss—Seidel迭代法 2.4.3