第二节误差分析.ppt

工程数学 杨圣洪 (手机) 2264454709(QQ) ysh@hnu.edu.cn 工程数学：复变函数、线性代数、概率论、场 论、数值计算或计算数学 ，等工程方面数学 计算过程 建模。用一个方程表示，还一个图。 设计求解数模的算法。伪代码，自然语言 编码实现，语言不限，matlab.m也可 用观测数据或原始公知数据进行运算。 可视化运算结果。图、表化，画很多的曲线 分析曲线，引出结果，提升自己算法的江湖地位。 工程数学—引言 问题可算化： 有限维--代—无限维，多项式?连续函数 有限过程—代-无限过程 定积分-数值积分 无穷级数取有限项 简单—代—复杂 ，尤其简单不断重复 如 简单重复xk+1=(xk+c/xk)/2求sqrt(c) 算法简单了，但在一定误差范围内与原问题等价。 设计求解数模的算法—核心 设计求解数模的算法—核心 秦九韶(中国 1202-1261)或Horner算法 p(x)=5+6x+7x2+8x3+9x4. 乘法：1+2+3+4 加法：4 p(x)=5+x(6+7(x+8(x+9x))) 乘法：4---n 加法：4 对于一般式p(x)=a0+a1x1+a2x2…anxn. b0=an, bk=an-k+bk-1x, (算例写黑板) 乘法：n(n+1)/2?n，效率提高了！ 加法：n不变 设计求解数模的算法—核心 求行列式的值Dn=|…|(板) Dn=?(-1)t(i1,i2,…,in) ai1,1ai2,2…ain,n. 乘法次数 n!(n-1) 4阶实例 n=20, 乘法：20*19=4.6*1019， 3.3亿亿=3.3万万亿=3.3*108*108. 4.6*1019/3.3*1016=4.6*1000/3.3=1394=23.2 4G=4*1024*1024*1024=4294967296 4.6*1019/4294967296/3600/24/365=339年 转换三角形 (n-1)2+(n-2)2+…+1+(n-1) =n(n-1)(2n-1)/6+(n-1) 4阶实例 误差分类 模型误差：少次要因素如飞行器， 观测(数据)误差：参数值测不准，预处理 公式(截断)误差：公式的近似式。sin(x) 舍入误差：能表示数有限，故存贮与计算总要取舍,3.1415926897,8位机要舍 例：地球表面积A 第二节：误差分析 绝对误差：?(x)=x-x , x=x+?(x) 由于x不知求不出?, |?(x)|??，为最后位半个单位 如：1800mm?0.5mm, 0.5*10-n.到10-n是精确的，精确点后n位 10-(n+1)就看不准了，画图 相对误差：?(x)=?(x)/x ? ?(x)/x’.(2者差去推导) 第二节：误差分析 浮点数与有效数字 x=3 x’=3.14 y=314.159265 y’=314.16 精度高！ z=0.0314159265 z’=0.03 有效数字多，精度高！浮点数标准式 fl(x)=?0.a1a2…at?m，t字长 ?进制 如t=3, ? =2，m=0,1,2 第一个数字非0 3位2进制：6组+0，其他以近似表示 第二节：误差分析 误差分类 模型误差：少次要因素如飞行器， 观测(数据)误差：参数值测不准，预处理 公式(截断)误差：公式的近似式。sin(x) 舍入误差：能表示数有限，故存贮与计算总要取舍,3.1415926897,8位机要舍 例：地球表面积A 第二节：误差分析 如t=3, ? =2，m=0,1,2 a1非0,浮点到非0 fl(x)=?0.a1a2a32m m=0,1,2 a1=1,a2=0,1,a3=0,1 m=0,1,2共2(t-1)*3 0.100,0.101,0.110, 0.111 m=0 ,乘20 0.5, 0.625, 0.750, 0.875 0.100,0.101, 0.110, 0.111 m=1 ,乘21点右移1位 1.0, 1.25, 1.50, 1. 75 0.100,0.101,0.110,0.111 m=1 ,乘22点右移2位 2.0, 2.5, 3. 0, 3. 5 12个与0，3位实数很多 如0.51,0.53 表示为0.50,舍入误差:0.01,0.03,少于 0.125(1/8)的一半0.0625, 又如0.57,近似为0.625 第二节：误差分析—舍入误差 fl(x)=0.a1a2…at?m进制，t字长 x=?a1a2…am.b1b