数值分析方程求根二分法.pptVIP

数值分析方程求根二分法第1页，共17页，星期日，2025年，2月5日

数值计算中的基本原则(1)避免绝对值小的数做除数;(2)避免两相近数相减;(3)防止大数“吃”小数现象a=109，b=9，在8位浮点数系统中做加法a+b=1.0000000×109+0.000000009×109由于只保留8位有效数，处于第九、十位的数09被舍去,实际操作是:将a的数据作为加法计算的最终结果.2/16第2页，共17页，星期日，2025年，2月5日

(4)尽量减少计算工作量(乘、除法次数)例计算P(x)=1+2x+3x2+4x3+5x4的值秦九韶算法P(x)=1+x(2+x(3+x(4+5x)))应用:2进制数转换为10进制数算法(11101110)2=27+26+25+0+23+22+2+0=((((((1·2+1)2+1)2+0)2+1)2+1)2+1)2+0=2383/16第3页，共17页，星期日，2025年，2月5日

例1计算(n=0,1,···,20)4/16(分部积分法)第4页，共17页，星期日，2025年，2月5日

初值：I0=1–e–1≈0.63212055882856n=20时，S20=-30.19239488558378递推公式:In=1–nIn-1(I0=1-e-1)S0=1-exp(-1);S(1)=1-S0;forn=2:20S(n)=1-n*S(n-1)end实际计算:Sn=1-nSn-1,S0=0.63212055882856|e(S0)|=|S0–I0|0.5·10-14error5/16第5页，共17页，星期日，2025年，2月5日

由In=1-nIn-1Sn-In=-n(Sn-1-In-1)则e(Sn)=–ne(Sn-1)=······=(n!)(–1)ne(S0)新算法:In-1=(1-In)/nS(30)=1/31forn=30:-1:2S(n-1)=(1-S(n))/n;endS0=1-S(1),S(1:21)初值误差在算法执行过程中不断增大,这种算法称为数值不稳定算法。6/16Sn-1-In-1=–(Sn-In)/n和Sn=1-nSn-1得Sn-1=(1-Sn)/n第6页，共17页，星期日，2025年，2月5日

在算法执行过程中,舍入误差对计算结果影响不大的一类算法被称为数值稳定算法;否则称为不稳定算法.初始误差在算法执行过程中不断减小,这种算法称为数值稳定算法。|e(S20)|=|S20-I20|=|(1-S21/21)-(1-I21/21)|=|S21-I21|/21=·······=|S30-I30|/(21·22·23·····30)7/16第7页，共17页，星期日，2025年，2月5日

rd例2.水中浮球问题有一半径r=10cm的球体,密度?=0.638.球体浸入水中后,浸入水中的深度d是多少？根据阿基米德定律,物体排开水的质量就是水对物体的浮力。整理得:d3–3rd2+4r3?=0注意：上式已知物体密度?=0.638,r=10,水的密度为1.第8页，共17页，星期日，2025年，2月5日

由?=0.638,r=10.代入,得d3–30d2+2552=0令f(x)=x3–30x2+2552,函数图形如下所示求解方程f(x)=0,即是求函数f(x)的零点.f(x)的零点所在区间为:[0,20]9/16roots([l-3002552])ans=26.314611.8615-8.1761第9页，共17页，星期日，2025年，2月5日

第一步:对根进行隔离,找出隔根区间,或在隔根区间内确定一个解的近似