数值计算及误差分析.ppt

第一章 绪 论 第一章 绪 论 第一章 绪 论 基本要求 第二节 数值问题与数值算法 一、误差的来源与分类 一、误差的来源与分类 一、误差的来源与分类 二、截断误差分析 二、截断误差分析 三、舍入误差分析 1、计算机数系 3、数值运算中的误差估计1）数值运算的绝对误差和相对误差 3、数值运算中的误差估计2）和、差、积、商的误差估计式 3、数值运算中的误差估计2）和、差、积、商的误差估计式 4、数值计算中应注意的问题 计算机中数的计算特点: 4、数值计算中应注意的问题 四、误差的定性分析 四、误差的定性分析 病态问题与条件数 习题一  P27-----3,4, 7, 11,12,13,14,15  对数学问题而言，如果输入数据有微小扰动，引起输出数据（即数学问题的解）有很大扰动，则称数学问题是病态问题，否则称为良态问题。 条件数 例4：考虑Matlab程序 x=0.988:0.0001:1.012; y=x.^7-7*x.^6+21* x.^5-35* x.^4+35* x.^3-21* x.^2+7* x -1; plot(x, y) syms x; y=x^7-7*x^6+21* x^5-35* x^4+ 35* x^3-21* x^2+7* x -1; ezplot(y,[0.988,1.012]) 在实数系中，每一个实数可以有无穷位，不同的 实数代表数轴上不同的点； 在计算机数系中，每一个数只有有限位，只有部分 有理数能被计算机数系中的数精确表示。 浮点数： 36.83=0.3683×102=0.03683×103 这种允许小数点位置浮动的表示法称为数的 浮点形式。 实数x的十进制浮点形式为 x= ? 0.a1 a2... ak …?10c, 尾数 阶码 ai?{0,1,2,…,9}, c∈Z 基数 a1≠0, (1)称为x 的 规格化的浮点形式 (1) x的k位十进制机器数fl(x)可用两种方法定义： （1）截断式 fl(x)=? 0. a1 a2... ak ?10c x= ? 0. a1a2…akak+1 …?10c (2)四舍五入式 x 的k位规格化十进制机器数 y= ? 0. a1 a2... ak?10c, y=fl(x) ai?{0,1,2,…,9}, a1?0，L?c ?U, k是机器数的字长； L、U 是常数。 一般数制情况： k位规格化机器数 y= ? 0. a1 a2... ak??c , ?=2,8,10,16, ai?{0,1,2,…, ?-1}, L?c ?U,a1≠0 F(?,k,L,U)表示以上数集全体加数0，它是计算机中使用的有限离散数集(机器数系）。 F(?,k,L,U)中的数称为机器数。 F(10,4,-33,33)， y= ? 0.a1a2a3a4?10c 在此机器数系 中表示 fl( ). p 例5 在机器数系 F(10,4,-33,33)中表示 fl(п). 若浮点数的阶码不在[L,U]内，则出现上溢或下溢。 采用截断式 fl(п )=0.3141×10 采用四舍五入式 fl(п )=0.3142×10 但是 例如 在4位机器数系 F(10,4,-33,33)中输入 出现下溢，输入 出现上溢。 Floatgui.m 2、数的误差和有效数字 定义1-1： 设数a是精确值，x是a的一个近似值，记 例6 已知准确值a=3.1415926…是一个无限不循环小数，求截取不同位数后的近似值和误差界。 解： 例7: 以下近似数是经四舍五入得到的 x1=0.12, x2=0.120, x3=1.73, x4=0.00073, x5=73000 x6=0.73*103 问它们分别有几位有效数字? 解： x1=0.12有2位有效数字, x2=0.120, 有3位有效数字。 x3=1.73有3位有效数字, x4=0.00073有2位有效数字。 x5=73000有5位有效数字， x6=0.73*103有2位有效数字。 例8:设x=219.15456为真值x*=219.15123的近似， 则x有几位有效数字? 答案:5位 例:考虑 1.x =101, a=100 , e=1, er=0.01 2. x =42196, a=42195, e=1, er=0.000023699 当某个量的准确值很小或很大时，相对误差比绝对 误差更能反映准确数与近似数的差异。