- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
二分法跟其matlab程序-经典
* 三. 二分法及其MATLAB程序 1. 用二分法求方程 f(x)=0的根 x*的近似值 xk 的步骤 步骤1. 若对于ab, 有f(a)f(b)0, 则在(a, b)内f(x)=0至少有一个根. 步骤2. 取a, b的中点 计算 步骤3. 若 则x1是f(x)=0的根, 停止计算, 若 则在(a, x1 )内f(x)=0至少有一个根. 取a1=a, b1=x1; 若 则取a1=x1, b1=b; 运行后输出结果x*=x1. 1. 用二分法求方程 f(x)=0的根 x*的近似值 xk 的步骤 步骤1. 若对于ab, 有f(a)f(b)0, 则在(a, b)内f(x)=0至少有一个根. 步骤2. 取a, b的中点 计算 步骤3. 若 则x1是f(x)=0的根, 停止计算, 若 则在(a, x1 )内f(x)=0至少有一个根. 取a1=a, b1=x1; 若 则取a1=x1, b1=b; 步骤4. 若 退出计算, 运行后输出结果x*=x1. 运行后输出结果 反之, 返回步骤1, 重复步骤1,2,3. 2. 二分法的matlab主程序 求解方程f(x)=0在开区间(a,b)内的一个根的前提条件是f(x)在闭区间[a,b]上连续, 且f(a)f(b)0. 输入的量: a和b是闭区间[a,b]的左右端点, abtol是预先给定的精度. 运行后输出的量: k是使用二分法的次数. x是方程在(a,b)内的实根x*的近似值, 其精度是abtol. wuca=|bk-ak|/2是使用k次二分法所得到的小区间的长度的一半, 即实根x*的近似值x的绝对精度限, 满足wuca≤abtol. yx=f(xk), 即方程f(x)=0在实根x*的近似值x处的函数值. 2. 二分法的matlab主程序 function [k,x,wuca,yx]=erfen(a,b,abtol) 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)); 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-a abtol , return, end end k=max1; x; wuca; yx=fun(x); % ceil是上取整 例: 确定方程x3-x+4=0的实根的分布情况,并用二分法求在开区间 (-2,-1)内的实根的近似值,要求精度为0.001. 0.003 7 0.012 1 -0.004 8 0.001 0 -1.795 9 -1.794 9 -1.796 9 9 0.012 1 0.029 0 -0.004 8 0.002 0 -1.794 9 -1.793 0 -1.796 9 8 0.029 0 0.062 7 -0.004 8 0.003 9 -1.793 0 -1.789 1 -1.796 9 7 0.062 7 0.129 6 -0.004 8 0.007 8 -1.789 1 -1.781 3 -1.796 9 6 -0.004 8 0.129 6 -0.141 8 0.015 6 -1.796 9 -1.781 3 -1.812 5 5 0.129 6 0.390 6 -0.141 8 0.031 3 -1.781 3 -1.750 0 -1.812 5 4 -0.141 8 0.390 6 -0.716 8 0.062 5 -1.812 5 -1.750 0 -1.875 0 3 -0.716 8 0.390 6 -2.000 0 0.125 0 -1.875 0 -1.750 0 -2.000 0 2 0.390 6 2.125 0 -2.000 0 0.250 0 -1.750 0 -1.500 0 -2.000 0 1 2.125 0 4.000 0 -2.000 0 0.500 0 -1.500 0 -1.000 0 -2.000 0 0 函数值f(xk) 函数值f(bk) 函数值f(ak) 中点xk 右端点bk 左端点ak 次数k 3. 用二分法求解方程f(x)
您可能关注的文档
- 电气控制跟可编程控制器技术 第三版 史国生课后习题答案.pdf
- 电影视听语言经典分析解析.ppt
- 电子商务培训章节稿-农村电子商务基础培训.ppt
- 顶岗实习校本科学研究探究基金项目申请书.pdf
- 定理跟证明.ppt
- 定义跟命题课件课本.ppt
- 动词不定式的用法回顾跟复习.ppt
- 动物的运动跟行为.ppt
- 动物的运动跟行为复习课课件课本.ppt
- 动物跟仿生学-【自制ppt】.ppt
- 四川省德阳市罗江中学2025届高三考前热身化学试卷含解析.doc
- 山东省枣庄现代实验学校2025届高三下学期第五次调研考试化学试题含解析.doc
- 吉林省长春市十一高中等九校教育联盟2025届高三一诊考试生物试卷含解析.doc
- 2025届江苏省盐城市伍佑中学高考仿真模拟化学试卷含解析.doc
- 2025届广西贺州中学高考冲刺押题(最后一卷)生物试卷含解析.doc
- 安徽省池州市贵池区2025届高三第一次模拟考试生物试卷含解析.doc
- 宁夏银川一中2025届高三(最后冲刺)化学试卷含解析.doc
- 广东省广州市增城区四校联考2025届高考压轴卷化学试卷含解析.doc
- 2025届邯郸市第一中学高考生物必刷试卷含解析.doc
- 2025届安徽省安庆市石化第一中学高考仿真卷化学试卷含解析.doc
文档评论(0)