二分法matlab程序.doc

二分法二分法基本思路一般地，对于函数f(x),如果存在实数c,当x=c时，若f(c)=0,那么把x=c叫做函数f(x)的零点。 解方程即要求f(x)的所有零点。 假定f(x)在区间（x，y）上连续 先找到a、b属于区间（x，y），使f(a)，f(b)异号，说明在区间(a,b)内一定有零点，然后求f(a+b)/2, 现在假设f(a)0,ab 如果f(a+b)/2=0，该点就是零点， 如果f(a+b)/2=a，从开始继续使用 中点函数值判断。 如果f(a+b)/20，则在区间(a,(a+b)/2)内有零点，(a+b)/2d c=(a+b)/2; if 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,2）内的实根的近似值，要求精度为0.0001解：(输入函数形式)fx=x4-2*x3+4*x+10（输入二分法下限）a=-2（输入二分法上限）b=2输入误差限 d=0.0001得到结果d = 1.0000e-004x = 2.0000k = 16