最优化第二讲—一维搜索算法(二分法、等区间法)

本讲主要列一下单峰搜索算法

问题:f(x)在区间[a,b]内只有一个极小值点,要找到这个极小值点或者这个极小值点所在的区间[x1,x2],其中[x1, x2]要远远小于[a, b]

方法:

一个通用的结论

要缩小区间,必须计算两个点,如果所示,必须计算x1和x2,然后对函数值进行比较。如果f(x1)小于f(x2),那么就说明极小值点一定在a到x2之间,反之也是这个思路。途中后面的两个式子将在以后的方法中反复使用。注意这了仅仅指的是单峰函数

具体的方法基本都是基于以上的思路,不同的是怎么确定x1和x2

具体的方法有

1. 二分搜索法(dichotomous search)

步骤一:首先找到[a,b]的中间点c,c=(a+b)/2,这就是“二分”的意思

步骤二:事先确定一个值sigma,在c左右各找二分之sigma,产生x1、x2。这个就是通用结论中的x1、x2

步骤三:按照通用结论中的方法来做

步骤四:循环确定区间[x1, x2],直到满足要求为止

查找速度:新区间的长度L(n+1),上一个区间的长度L(n),他们的关系是:L(n+1) 

  • 1
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
二分法是一种用于求解函数零点的常见算,它的基本思想是将一个区间逐步分割为两个子区间,直到最后找到一个足够小的区间,使得函数在该区间内的函数值与零点的误差小于给定的容忍度。在Matlab中,我们可以用递归的方实现二分法分割区间,具体步骤如下: 1. 定义函数。首先要定义一个用于求解函数零点的函数,可以用Matlab自带的函数fzero(),也可以自己编写一个函数。 2. 定义区间。在主函数中,定义需要求解零点的区间[a,b]。 3. 分割区间。用while循环依次对区间进行分割,直到满足给定的精度要求。 4. 判断根的位置。在每个子区间中,判断根的位置,如果根在左边区间,则进入左边区间继续分割,否则进入右边区间分割。 5. 输出结果。当找到根之后,输出求解结果。 二分法分割区间的Matlab代码如下: function [x,fval,exitflag,output] = bisection(funfcn,a,b,options,varargin) if nargin < 4 || isempty(options), options = optimset; end if any(strcmp(fieldnames(options),'TolFun')) && strcmp(options.TolFun,'eps') options = optimset(options,'TolFun',1000*eps); end if (b-a)<options.TolX, x=a; fval=feval(funfcn,x,varargin{:}); exitflag=1; output=[]; return; end fa = feval(funfcn,a,varargin{:}); fb = feval(funfcn,b,varargin{:}); while (b-a)/2 > options.TolX c = (a+b)/2; fc = feval(funfcn,c,varargin{:}); if fc==0, x=c; fval=fc; exitflag=1; output=[]; return; end if sign(fc)==sign(fa), a = c; fa = fc; else b = c; fb = fc; end end x = (a+b)/2; fval = feval(funfcn,x,varargin{:}); exitflag = 1; output.iterations = NaN; output.funcCount = NaN; output.algorithm = 'bisection';

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值