计算方法3new.pptVIP

第一章总结 一、误差来源 模型误差/描述误差 观测误差/测量误差 截断误差：把无穷数用有限数表示。 舍入误差：受字长限制，取有限位数字运算。 作业：[数]P16:1.1,1.3,1.4 [科]P26:1,2,3 第二章 方程求根 可归结为f(x)=0求根或者求f(x)零点。 求n次代数方程 求超越方程 问题提出： 一、 根的分类 （一） 实根和复根 1、 实根—— 2、 复根——复数根 实数根 （二） 单根和重根 1、 若f (x)可分解因式 则α是f(x)的m重零点，或是方程 f(x) =0的m重根。 m=1 α——单根 m1 α——重根 2、 函数f (x) 对x=α有 f(α)=0且f′(α)≠0 。 α— 单根 α— m重根 二、 定理 设f(x)有m阶连续导数，则x*是f (x) = 0的m重 根的充要条件是 1、 已知大概范围，位置求单根。 三、 求根的情况 2、 未知根的数目、范围和位置，求全部根。 四、 举例 例1：已知f(x)=3(x-α) ,问x=α是几重根？ 解： x=α是单根 例2： 解： x=0是方程f (x)=0的三重根 2.1 二分法 2.1.1 初始值的搜索 2.1.2 区间二分法 2.2 迭代法的一般理论 2.2.1 不动点迭代 2.2.2 迭代法的收敛性 2.2.3 迭代法的收敛阶 2.2.4 迭代过程的加速 2.3 牛顿迭代法 2.3.1 牛顿迭代公式 2.3.2 局部收敛和收敛阶 第二章 方程求根 2.1.1 初始值的搜索 2.1 二分法 一、 有根区间 定理：若函数f(x)在[a,b]连续且两端异号，即 f(a)f(b)0,则方程f(x)=0在(a,b)内至少 有一个实根。将[a,b]称为有根区间。 二、 逐步搜索法 y x f(x) a b x1 x2 x3 x4 例1：已知方程f (x) =x3-x-1=0,利用逐步搜索法确 定有根区间。 解： f (x) 为连续函数， f (0)0, f(2)0, 有根区间[0,2] 取h=0.5，从x=0出发，逐步搜索。 x f(x) 0 0.5 1 1.5 2 _ _ _ + + 确定缩小的有根区间为 [1，1.5] 例2：求方程f (x) =x3-1.8x2+0.15x+0.65=0的有 根区间。 解： 设 f (x) =x3-1.8x2+0.15x+0.65 取 a =-1，b=2。 f (-1) 0, f(2)0 所以方程在[1，2]内至少有一个实根 取N=4，h=(2+1)/4=0.75 x f(x) -1 -0.25 0.5 1.25 -2.3 0.484 0.55 -0.0219 2 1.75 可以看出在区间[-1，-0.25],[0.5,1.25],[1.25,2]内各有一根。 * 1.3.4 有效数字与相对误差 一、 定理1及说明 1、定理1： 若近似数x* =± 0.x1x2…xn×10m具 有n位有效数字，则其相对误差 例： 取3.14作为π的近似值，求er* 解： 3.14作为近似值，n=3,x1=3 二、 定理2及说明 1、定理2： 若近似数x* =± 0.x1x2…xn×10m的相对误差 则该近似数具有n位有效数字。 例： 已知近似数的相对误差为0.25%，问至少 有几位有效数字？ 解： x1=1，n=3 x1=9，n=2.3 则x*至少有2位有效数字。 1.4.1 函数运算误差 一、 一元函数f(x)运算误差 1.4 运算误差分析 二、多元函数y=f(x1,x2, …,xn)运算误差 1.4.2 算数运算误差 二、 乘除运算误差 1、 乘法运算误差 2、 除法运算误差 1.5.1 数值稳定性 一、 运算的数误差来源 1.5 数值稳定性和减小运算误差 1、 初值本身误差： 观测误差、 估算误差 2、 受计算机字长限制所造成的舍入误差 二、 数值稳定性 1、 误差的积累和传播受控 2、 对初始误差和舍入误差不敏感 三、 应选用数值稳定性好的计算公式 0.1008 0.0760 -0.7280 I8 0.1124 0.1155 0.2160 I7 0.1268 0.1264 0.1120 I6 0.1455 0.1456 0.1408 I5 0.1709 0.1709 0.1704 I4 0.2073 0.2073 0.2074 I3 0.2642 0.2642 0.2642 I2 0.3