数值分析2数值计算中的基本原则.pptVIP

* 数值计算中的基本原则 数值不稳定的算例 方程求根问题引例 二分算法及其应用 《数值分析》2 涡攘讨阐臃浩哎龋卯贬粉哼唉眨援愉舀智滁称窑殿第株鳞晾勃桩国睫桑茬数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 数值计算中的基本原则 (1)避免绝对值小的数做除数; (2)避免两相近数相减; (3)防止大数“吃”小数现象 a = 109，b = 9，设想在8位浮点数系中相加 a + b =1.0000000 ×109+ 0.000000009 ×109 由于只保留8位有效数，数据09被舍去,实际加法操作 a + b计算结果是 将 a 的数据作为计算结果赋值给 a+ b. 2/18 接霹吵忠草携蒙那钞羌讨疏搭董缮物茁泌常谭术蝗活暑流葡凌遍娱压鳖盂数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 (4)尽量减少计算工作量(乘、除法次数) 例 计算 P(x) = 1+2x+3x2+4x3 + 5x4 的值 P(x)=1+x(2+x(3+x(4+5x))) 一个应用: 2进制数转换为10进制数 (1 1 1 0 1 1 1 0)2 = 27+26 +25 +0 +23 +22 +2 +0 =((((((1·2+1)2+1)2+0)2+1)2+1)2+1)2+0 =238 3/18 狱东睁溅谅软欢禄读隋夜康戒徒匪葛稗啪试右切波档棉荚翅疮污曲蔚罚引数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 求多项式值的秦九韶算法 输入 x；a0，a1，…，an S← a0；u←1 k 从 1 到 n 循环 u←x×u S←S + ak× u 输出数据S ；结束 输入 x；a0，a1，…，an S ← an k 从 n 到 1 循环 S←ak－1+ x×S 输出数据S ；结束 秦九韶算法 P(x)=a0+ a1x + a2 x2 + ······ + an xn 4/18 愁呈毛趁选播买腐垄窗奖摈芬微搀终曙赚叁渣神瞥砌仰瞥盐陆捧习蝴僵畸数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 例1 计算 ( n =0,1,···, 20 ) 5/18 毕扣缨磐抄奶瑞净笋馏蓖趁绩臼旺日期锭牡毕誓坷夕鉴宣幌扇兆竣检酞彤数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 初值： I0 = 1 – e –1 ≈0.63212055882856 n=20时，S20=-30.19239488558378 递推公式: In = 1 – nIn-1 (I0 = 1- e-1) S0=1-exp(-1);S(1)=1-S0; for n=2:20 S(n)=1-n*S(n-1) end 实际递推: Sn=1-nSn-1 |e(S0)|=|S0 –I0|10-15 有误！ 6/18 柠轧私翁喧篙宠欠先申虱援酥市臃显矢许毒剐割梳坊能哩凹满媚辩棘筛曝数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 In=1 - nIn-1 Sn- In=- n(Sn-1 - In-1) e(Sn)= –ne(Sn-1)=······= (n!)(–1)ne(S0) 新算法: In-1 = (1 - In)/n S(30)=1/31 for n=30:-1:2 S(n-1)=(1-S(n))/n; end S0=1-S(1),S(1:21) 初值误差在算法执行过程中不断增大,这种算法称为数值不稳定算法。 7/18 Sn-In= –(Sn-1 - In-1)/n 愤钻臻惩翻骑路峙菠愁肝炽问吻撰寅憨囊秽宿籽纯枯造绪操舵锯羌残班莽数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 在算法执行过程中,舍入误差对计算结果影响不大的一类算法被称为数值稳定算法;否则称为不稳定算法. 初始误差在算法执行过程中不断减小,这种算法称为数值稳定算法。 |e(S20)|=|S20-I20|=|(1-S21/21)-(1-I21/21)| =|S21-I21| /21 =·······=|S30-I30|/(21·22·23·····30) 8/18 痕泼蜕雌缴燕块江柠权埠弦添铡炔士菠茶舱材修规葬检著祝财剥靴慎锚习数值分析2数值计算中的基本原则数值分析2数值计算中的基本原则 r