二分法,牛顿法,梯形法原理及流程图..docVIP

开始 输入区间 [a,b] ，精度 x=(a+b)/2 f(x)= Y f （ x）=0 N N Y f(x)f(a)0 a=x b=x N /x 1-x 2/ ＜ Y 输出 x 结束 二分法基本思路： 一般地，对于函数函数 f(x) 的零点。  f(x),  如果存在实数解方程即要求  c, 当 f(x)  x=c 时，若 f(c)=0, 的所有零点。  那么把  x=c  叫做 假定 f(x) 在区间（ x，y）上连续 先找到 a、 b 属于区间（ x，y），使 f(a) ，f(b) 异号，说明在区间 (a,b) 内一定有 零点，然后求 f[(a+b)/2], 现在假设 f(a)0,f(b)0,ab ① 如果 f[(a+b)/2]=0 ，该点就是零点， 如果 f[(a+b)/2]0, 则在区间（ (a+b)/2 ， b) 内有零点， (a+b)/2=a ，从①开始 继续使用 ② 中点函数值判断。 如果 f[(a+b)/2]0 ，则在区间  (a,(a+b)/2)  内有零点， (a+b)/2=b  ，从①开始继 续使用 中点函数值判断。 这样就可以不断接近零点。 通过每次把 f(x) 的零点所在小区间收缩一半的方法，使区间的两个端点逐步迫近函数的零点，以求得零点的近似值，这种方法叫做二分法。 从以上可以看出，每次运算后，区间长度减少一半，是线形收敛。另外，二分法不能计算复根和重根。 二分法步骤： 用二分法求方程 f (x) 0 的根 x* 的近似值 xk 的步骤 ① 若对于 a b 有 f ( a) f (b) 0，则在 (a, b) 内 f ( x) 0 至少有一个根。 取 a, b 的中点 x1 a b ② 计算 f ( x1 ) 2 ③ 若 f ( x1 ) 0 则 x1 是 f (x) 0 的根，停止计算， 运行后输出结果 x* x1 若 f (a) f ( x1 ) 0 则在 (a, x1 ) 内 f (x) 0 至少有一个根。取 a1 a,b1 x1 ； 若 f (a) f ( x1 ) 0 ，则取 a1 x1, b1 b ； ④ 1 bk ak （ 为预先给定的要求精度）退出计算，运行后输出结果 若 2 x* ak bk 2 ，反之，返回步骤 1，重复步骤 1,2,3 二分法 Mtalab 程序 syms x; fun=input(( 输入函数形式 )fx=); a=input( （输入二分法下限） a=); b=input( （输入二分法上限） b=); d=input( 输入误差限 d=)% 二分法求根 %f=inline(x^2-4*x+4); %修改需要求解的 inline 函数的函数体 f=inline(fun);% 修改需要求解的 inline 函数的函数体 e=b-a; k=0 ; while ed c=(a+b)/2; if f(a)*f(c)0 b=c; elseif f(a)*f(c)0 a=c; else a=c;b=c end e=e/2; k=k+1; end x=(a+b)/2; x%x 为答案 k%k 为次数 2，牛顿法及流程图： 方程 f(x)=0 的根就是曲线 y=f(x) 与 x 轴交点的横坐标 x* ，当初始近似值 x0 选取后， 过( x0,f(x0))作切线，其切线方程为： y- f(x0)=f ′(x0)(x -x0) 它与 x 轴交点的横坐标为 x 一般地，设 是 x* 的第 n 次近似值，过 ( x,f(x)) 作 y=f(x) 的横坐标为： x = - 即用切线与 x 轴交点的横坐标近似代曲线与 x 轴交点的横坐标，如图  的切线，其切线与  x 轴交点 牛顿法正因为有此明显的几何意义，所以也叫切线法。 流程图如下： 开始 输入 x0 ， ， N 1=k f =0? (x0 )  Y N f (x0 ) =x1 x0 f ( x0 ) ∣x1-xo ∣ ? Y N k+1=k N K=N ? x1=x0 Y 输出迭代失败标志 输出 x1 输出奇异标志 结束 3，梯形法及流程图： 梯形法就是将该积分约等于若干个小梯形面积之和，第一个小梯形的面积等 为 s1 = h( f (a) + f (a +h)) / 2 , 第二个小梯形的面积为 s2 = h( f (a + h) + f (a + 2h)) / 2 ， ， 第 i 个小梯形的面积为 si = h( f (a + (i -1)h) + f (a +ih )) / 2 故有 b f ( x) = n si = h[ 1 ( f (a) + f (b)) + n -1 f (a +ih )] a i 1 2