matlab在二分法中的应用张扬.docVIP

MATLAB在二分法中的应用张扬 ———————————————————————————————— 作者： ———————————————————————————————— 日期： ?MATLAB语言?课程论文 MATLAB在二分法中的应用 姓名：张扬 专业：电子信息工程 班级：10级一班 指导教师：汤全武 学院：物理电气信息学院 MATLAB在二分法中求方程的近似解的应用 【摘要】计算一个函数f的零点(即求方程的根)是在科学计算中经常遇到的问题,二分法是解决此类问题的方法之一,MATLAB在处理用二分法求方程的近似解问题上非常方便,而且可以为学生创造出图文并茂、丰富多彩人机交互、即时反应的学习环境,为学生更好的学习数学提供了一个有力的辅助工具. 【关键词】二分法 精度 近似值 迭代次数 一、二分法 二分法也称逐次分半法。它的根本思想是：先确定方程含根的区间(a,b),再把区间逐次二等分。 〔一〕判别有根区间 设函数在区间[a,b]上连续，且,那么根据零点定理知方程在开区间〔a,b)内至少有一个实根，同时假设在[a,b]上是单调的，那么方程在开区间(a,b)有唯一一个实根。这是使用二分法的前提。 〔二〕用二分法求方程的根的近似值的步骤 步骤1 假设对于ab,有,那么在(a,b)内至少有一个实根 步骤2 取(a,b)的中点,计算. 步骤3 假设,那么是的根，停顿计算，运行后输出结果=；假设,那么在(a,〕内至少有一个根。取,;假设那么取,。 步骤4 假设 ( 为预先给定的要求精度〕，退出计算，运行后输出结果；；反之，返回步骤1，重复步骤1，步骤2和步骤3. 〔三〕二分法的收敛性 以上方法可得到每次缩小一半的区间序列｛ [,]},(,)中总有方程的根。当区间长-很小时，取其中点= (+)/2为根的近似值，显然有 (k=1,2,……), (1) 由式(1)可知，总有，且近似根可以到达任意精度。 事实上，对于给定的精度要求 ，要使 只要按次数 〔2〕 逐次积分即可。 〔四〕 求迭代次数的MATLAB命令 闭区间[a,b]和误差 。用二分法求方程误差根的迭代次数的MATLAB命令 k=-1+ceil((log(b-a)-log(abtol))/log(2)) %ceil是向无穷大方向取整，abtol是误差 。 〔一〕二分法的MATLAB主程序 二分法需自行编制程序，根据用二分法求方程的根的近似值的步骤和式〔1)编写一个名为的二分法的MATLAB主程序如下 function[k,x,wuca,yx]=erfen(a,b,abtol) %a和b是闭区间[a,b]的左右端点； %k是是用二分法的次数； %x是实根x*的近似值； %wuca是x的绝对精度限； %yx是方程在f(x)=0在x处的函数值； a(1)=a;b(1)=b; ya=fun(a(1));yb=fun(b(1)); %程序中调用的为函数 if ya*yb0, disp(注意：ya*yb0,请重新调整区间端点a和b.),return end max1=-1+ceil((log(b-a)-log(abtol))/log(2)); %ceil是向正无穷大方向取整 for k=1:max1+1 a;ya=fun(a);b;yb=fun(b);x=(a+b)/2; yx=fun(x);wuca=abs(b-a)/2;k=k-1; [k,a,b,x,wuca,ya,yb,yx] if yx==0 a=x;b=x; elseif yb*yx0 b=x;yb=yx; else a=x;ya=yx; end if b-aabtol,return,end end k=max1;x;wuca;yx=fun(x); 〔二〕使用二分法的主程序求解方程在开区间(a,b)内的一个近似根的值的步骤 步骤1 建立名为 的M文件如下 function y1=fun(x) y1=(