浅析二分法及其Matlab和C程序实现

1、1浅析二分法及其 Matlab 和 C 程序实现第一部分：二分法浅析用二分法求方程的近似解是紧跟在“函数的零点”之后的教学内容。从联系的角度看， 前面一节，学生已经学习了方程的根与函数的零点之间存在着对立统一的关系，这一节 则是介绍一种具体的方法来运用这一关系解决问题。从整个教材来分析，这一部分的内容是在“函数的应用”这一大章节之下。新课程标 准中强调函数的应用性，这里包括两个方面：一方面是函数在生活实践中的应用，函数 建模等内容属于这个范畴；另一方面则是函数在数学自身范围内的应用，“二分法”即是其中的代表。基于以上的分析，笔者给出了以下的一些教学建议，与读者朋友们分享。一、为什么要用二分法就

2、通过试验缩小搜索区间来讲，试验点不一定取中点，取其他的点也可以，那么为什么取中点呢？下面以搜索区间为0，1的情况作讨论。一种对所有搜索区间为0，1的方程f(x) = 0都适用的方法，即对集合G=f(x)=0，f(X)连续，且f(0)f(1) 0中的所有方程都适用的方法.一个合理的假设是：G中所 有方程f(x)= 0的根在0，1上均匀分布.设试验点是c,那么c将0，1分成0，c和c,1两部分，它们的长度分别是c和1-c.由假设，通过试验保留的搜索区间是0，c(即 方程f(x)=0的根在0 ,C中)的概率是c，通过试验保留的搜索区间是C，1的概率是1-c.因此，通过一次试验保留的搜索区间的期望长度

3、为c2+ (1 -c)2= 2c2- 2c+ 1=2(c- )2+，2 2容易看出，当c=的时候，通过一次试验保留的搜索区间的期望长度最小。这就是取中2点作为试验点的原因。二、引入方法方法1:已知商店里一件商品的利润y与它的价格x之间满足函数关系y=x2-4x+3， 请画出这个函数的图像，并思考当价格为多少元的时候商店不盈也不亏方法2:创设问题情景： 蹦极运动.设下落的时间t秒.人离开参照点“礁石尖端”的 位移为S(S=0表示人在礁石点处，向下取负，向上取正)，开始下落时，时间t=0，在t4，6时的变化如2下表：3t4. 04. 34. 64. 95. 25. 56. 0S-5101-1018

4、-3问：这段时间内人有几次通过礁石尖端处？方法3:使用“幸运52”猜测商品价格的游戏作情景.方法4：请同学们思考下面的问题：能否解下列的方程1X2-2X-1=02lgx=3-x3X4-3X-1=0特殊入手：不解方程求方程X2-2X-1=0的近似解(精确到0. 01).方法1、2、3都是以“实际问题”为情境引入.方法4以学生已有的认知水平：会求一元二次方程的实数解，对应二次函数的图像与二轴的交点坐标.让学生探究具体的一元二次方程的根与其对应的一元二次函数的 图像与二轴的交点的横坐标的关系，再探究一般的一元二次方程的根与其对应的一元 二次函数的图像与X轴的交点的横坐标的关系.三、函数零点的处理用二

5、分法求方程近似解的理论基础是零点存在定理.下面我们来看看教材上描述的零点存在定理.如果函数y =f (X)在区间a，b上的图像是连续不断的一条曲线，并且有f (a ) f (b)0，那么函数y =f (x)在区间a，b内有零点即存在c(a，b)，使f (C)=0.由此可见，定理的题设部分有两个条件：(1) y =f (X)在区间a，b上的图象是连续不断的一条曲线；f (a ) f (b)0，则y =f (X)在(a，b)是否就不存在零点呢？对于问题，我们可以告诉学生我们前面所学的一次函数、二次函数、指数函数、对数函数、幕函数在它们各自的定义域内图像都是连续的.这些函数经过加减乘除或经过复合而成

6、的新的函数在各自的定义域内图像仍然是连续的.对于问题，主要通过观察函数图像来总结.4567A 3211人必恥0有四个零点/Xf/ / 1K /1f* /t1/ 1:A 1i1./ /1140-11 2|4写*- / /t1*11 /-2图61 /(1)对全部零点为单重零点既对应方程无重根的情况.y=f (x)在区间a,b上的图象是连续不断的一条曲线且f(a) f (b)0，则y =f (x)在(a,b)上有偶数个零点若y =f (x)在区间a,b上单调则y =f (x)在(a,b)上有无零点.可以看出连续函数的零点具有一个很重要的性质：函数的图象如果是连续的，当它通 过零点时，函数的值变号，也

7、就是图象要经过该点要穿越x轴.(2)对多重零点的情况8/(X)=(A*-1)2(X-2)CV3)入=1是二重零点“-图X从图7、图8可以看出偶数重零点不穿过x轴；奇数重零点穿过x轴函数若有一零 点为多重零点，当该零点为偶重零点时，图象通过该零点时，函数值不变号，也就是图 象经过该零点而不穿越x轴.当该零点为奇重零点时，图象经过该点时函数值要变号， 也就是图象经过该零点且穿越x轴.处理好这个问题是本节课的关键.四、精度/U戶何1他2) )x=l是三重零点9精确度的说明是一个无法避免的问题， 而且需要和初中学习的 “精确到”有所区分.教 学中不可能让学生掌握严格的、形式化的定义，而且教科书对此也作

8、了简化处理：对于 达到精确度&的界定是只要精确值所在区间的长度小于,那么这个区间的所有的值就都是满足精确度&的近似值.那么，如何让学生明白这个含义呢？一个可行的方法就是通 过简单例子来说明问题最后，在学生思考、讨论及进一步分析的基础上给出精确度的 含义：“一般地，对于数值X,如果要获得它的满足精确度0. 01的近似值，就是找到一个包含X的区间a,b,只要I a-b|0,f(1)0,求证:a(I)a0且-2v v-1；b(n)方程f(x)=0在(0,1)内有两个实根.证明：(略)(Ji)解法1：抛物线的顶点坐标为:I,3a 3a利用二分法思想L1在- 的两边乘以-，得a31b2-

9、 -.33a3又因为:-古“b、(?+八血f、而 |：二 0抛物线开口方向向上)综合、，由函数的图象 可知：方程f(x) = 0在(0,1)内必有两个不同实根.1.5001221.5001831.5002440.0004026.75E -05-0.00027-1.07627E -072.71495E -08-1.8E-080.00012207 141.5001831.5002141.5002446.75E-05-1E-04-0.00027-1.80465E-08 -6.74688E-092.67E-086.10352E-05151 5001831 50019815002146 75E-05-1

10、 6E-05-1E-04-6 74688E -09-1 09724E-091 63E-093 05176E-05161.5001831.500191 1.5001986.75E-05 2.56E-05-1.6E-05-1.09724E-091.72755E -09-4.2E-101.52588E-05171.5001911.5001951.5001982.56E-05 4.67E-06-1.6E-05-4.16389E-101.19598E-10-7.6E-117.62939E-06181 5001951 50019615001984 67E-06-5 8E-06-1 6E-05-759619

11、E -11-2 70718E-119 43E-113 8147E-06 191.5001951.5001961.5001964.67E-06 -5.6E-07-5.8E-06-2.70718E-11-2.62675E-123.26E-121.90735E-06 201.5001951.5001951.5001964.67E-06 2.05E-06-5.6E-07-2.62675E-12 9. 59574E-12-1.2E-129.53674E-07 211.5001951.5001951.5001962.05E-067.46E -07-5.6E-07-1.15526E-12 1.53249E-

12、12-4.2E-134.76837E-07 221211?1在区间(0，1)内选取二等分点一，因f)= a+b+c= a+ (-a)=a 0，22444113所以结论得证.(若f( ) 0不成立，可看f()是否为负，若还不成立，再看f()是否244为负，总之，在区间(0，1)内存在一个分点，使它对应的函数值为负即可.)例2( 2006全国1112题)函数f(x)=X-1| + |x-191的最小值为()(A)190(B)171(C)90(D)45分析 因x-n|表示数轴上的动点x到点n之间的距离.当x-1| + x-19|最小时，x为 区间1,19内的任意一个分点;,以此类推，当|x- 9|

13、+ |x- 11|最小时，x为区间9，11内的任意一个分点；当x-10最小时，x=10.利用“二分法”的思想方法，当x是区间1，19，2，18，3，17，,，9，11共同二等分点，即x=10时，f(x)取得最小值，所以f(x)min= 10 - 1 + 10- 2 + 10-3 + 1Q9 + 10-10 + 10 - 11 + 10- 19 = 90.第一部分：二分法的Matlab和C程序实现第二部分主要通过一个实例来研究单变量非线性方程f(x)=0的二分法求解及此方 法的收敛性，根据误差估计确定二分次数并进行求解。同时实现Matlab和C语言程序编写。从而掌握过程的基本形式和二分法的基本思

14、想，在以后的学习过程中得以应用。一、弓丨言在科学研究与工程技术中常会遇到求解非线性方程f(x)=0的问题。而方程f(x)是多项式或超越函数又分为代数方程或超越方程。对于不高于四次的代数方程已有求根公式， 而高于四次的代数方程则无精确的求根公式，至于超越方程就更无法求其精确解了。因 此，如何求得满足一定精度要求的方程的近似根也就成为了我们迫切需要解决的问题。 近年来，随着数学科学研究的不断进展，又更新了许多方程求解的方法。我们知道，对 于单变量非线性方程f(x)=0, 一般都可采用迭代法求根，由此产生了二分法。二、Matlab实现function c,err,yc=bisect(f,a,b,de

15、lta)%f是要求解的函数%a和b分别是有根区间的左右限%delta是允许的误差界%c为所求近似解13%yc为函数f在c的误差估计if nargin0disp(a,b)不是有根区间)；return,end max1=1+round(log(b-a)-log(delta)/log(2);for k=1:max1 c=(a+b)/2; ye=feval(f,c);if yc=0 a=c;b=c;break, elseif yb*yc0 b=c;yb=yc;else a=c;ya=c;endif(b-a)endk,c=(a+b)/2,err=abs(b-a),yc=feval(f,c)例如，要求f(

16、x)=x3-x-1=0在区间1，2内的根 先在命令窗口中输入：fplot(xX3,0,1 2);grid回车后输出曲线图现在再编写所求的非线性函数，输入：f=inline(xA3-x-1);bisect(f,1,2)敲回车就能得到你想要的了。、C实现#i ncludestdio.h#i ncludemath.hmai n()double x1,x2;double x12;double y1,y12;x1=1;x2=3;14doy1=exp(x1)-si n(x1)-sqrt(x1)-2; x12=(x1+x2)/2.0;y12=exp(x12)-si n(x12)-sqrt(x12)-2; if(y1*y121e-12);prin tf(x is:%.6gn,x12);运行结果：编程思路:二分法的思是，取自变量两个边界的中间值，并求此自变量值下的函数值，与其中的 一个边界值(假若为第一个边界值且为负)相乘，如结果为负，则表明该中间值包含函 数值为0的x值，再减小后边界继续循环即可；如结果为正，则表明该中间