MATLAB在分法中的应用 张扬.docVIP

《MATLAB语言》课程论文 MATLAB在二分法中的应用 姓名：张扬 学号：12009243762 专业：电子信息工程 班级：10级一班 指导老师：汤全武 学院：物理电气信息学院 完成日期：2011.11.31 MATLAB在二分法中求方程的近似解的应用 （张扬 12009243762 物理电气信息学院 10级电子1班） 【摘要】计算一个函数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程序 （一）二分法的MATLAB主程序 二分法需自行编制程序，根据用二分法求方程的根的近似值的步骤和式（1)编写一个名为erfen.m的二分法的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)); %程序中调用的fun.m为函数 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 建立名为 fun.m的M文件如下 function y1=fun(x)