布伦特方法(Brent‘s method)---结合二分法、割线法和逆二次插值法的求根方法

基础介绍:

给定给定区间\left [ a,b \right ],函数连续且f(a)\cdot f(b)<0,那么根据介值定理,函数必然在区间内有根。

  1. 二分法:将区间不断二分,使端点不断逼近零点。下一次迭代的区间为\left [ a,c \right ]\left [ c,b \right ],其中c=\frac{a+b}{2}
  2. 割线法(线性插值):基本思想是用弦的斜率近似代替目标函数的切线斜率,并用割线与横轴交点的横坐标作为方程式的根的近似。即给定两个点\left ( a,f(a) \right ),\left ( b,f(b) \right )。其割线方程为y=\frac{f(b)-f(a)}{b-a}\cdot (x-b)+f(b),那么令y=0,x的值即为下一次迭代的结果c=b-\frac{f(b)\cdot (b-a)}{f(b)-f(a)}
  3. 逆二次插值法:为割线法的进化版本。使用三个点确定一个二次函数,二次函数与横轴交错的点即为下次迭代的值。但是,其二次函数可能不会和横轴相交,因此做出一点改变,以y值作为自变量。给定三个点\left ( x_{i-2},f(x_{i-2}) \right ),(x_{i-1},f(x_{i-1})),(x_{i},f(x_{i})),则通过这三个点确定的二次函数为x=\frac{(y-f(x_{i-1}))(y-f(x_{i}))}{(f(x_{i-2})-f(x_{i-1}))(f(x_{i-2})-f(x_{i}))}\cdot x_{i-2}+\frac{(y-f(x_{i-2}))(y-f(x_{i-1}))}{(f({x_{i}})-f(x_{i-2}))(f(x_{i})-f(x_{i-1}))}\cdot x_{i}+\frac{(y-f(x_{i-2}))(y-f(x_{i}))}{(f(x_{i-1})-f(x_{i-2}))(f(x_{i-1})-f(x_{i}))}\cdot x_{i-1},令y=0,求得x_{i+1}=\frac{f(x_{i-1})f(x_{i})}{(f(x_{i-2})-f(x_{i-1}))(f(x_{i-2})-f(x_{i}))}\cdot x_{i-2}+\frac{f(x_{i-2})f(x_{i-1})}{(f({x_{i}})-f(x_{i-2}))(f(x_{i})-f(x_{i-1}))}\cdot x_{i}+\frac{f(x_{i-2})f(x_{i})}{(f(x_{i-1})-f(x_{i-2}))(f(x_{i-1})-f(x_{i}))}\cdot x_{i-1}

布伦特方法:

初始化区间(a_{0},b_{0})使得f(a_{0})\cdot f(b_{0})<0。其中b_{k}是上次迭代中的根估计值。如果\left | f(a_{0}) \right |<\left | f(b_{0}) \right |,那么赋值互换(我们认为对应函数值的绝对值较小的点更接近真正的根值)。

每次迭代包含四个点:

  1. b_{k}:为当前迭代的根估算值;
  2. a_{k}:对位点,即满足\left | f(a_{k}) \right |<\left | f(b_{k}) \right |f(a_{k})\cdot f(b_{k})<0的值。
  3. b_{k-1}:上一次迭代的根估算值,第一次迭代设置为b_{k-1}=a_{0}
  4. b_{k-2}:上上此迭代的根估算值(不用初始化,在首次迭代过程中,不会用到他来进行判断,结尾进行赋值)。

有以下四个不等式:

\left | \delta \right |<\left | b_{k}-b_{k-1} \right |  ①

\left | \delta \right |<\left | b_{k-1}-b_{k-2} \right |  ②

\left | s-b_{k} \right |<\frac{1}{2}\left | b_{k}-b_{k-1} \right |  ③

\left | s-b_{k} \right |<\frac{1}{2}\left | b_{k-1}-b_{k-2} \right | ④

上次迭代为二分法且①为假;上次迭代为二分法且③为假;上次迭代为插值法且②为假;上次迭代为插值法且④为假;以插值法计算的临时值不在\frac{3a_{k}+b{k}}{4}b_{k} 中间,以上五个条件满足一个,那么本次迭代的值采用二分法,否则采用插值法。

而插值法的选择如下:如果三点各不同,则用二次插值;否则用线性插值。

本次迭代的临时值s作为区间的一个端点,另一个端点在a_{k}b_{k}中选择,二者作为a_{k+1},b_{k+1},且满足f(a_{k+1})\cdot f(b_{k+1})<0,|f(a_{k+1})|>|f(b_{k+1})|

 

 

 

 

 

 

 

 

 

布伦特方法Brent's Method)是一种结合二分法割线法二次插值法的数值求解根算法。它利用这三种技术的优点,在保证收敛的同时提高效率。下面通过一个具体的例子来说明如何应用布伦特方法。 ### 示例题目 假设我们有函数 f(x) = x^3 - 2x + 1,并需要找到其在区间 [0, 2] 内的一个实数根,精度达到小数点后6位。 #### 步骤解析: **初始条件** 已知 f(0) * f(2) < 0 (即两端点处f(a)f(b)<0),所以按照介值定理可以确认该区间内存在至少一个零点。 **迭代步骤** 1. **初始化阶段**: 设置最大允许误差ε=1e−6,以及a,b初值分别为02。 2. **选择合适的策略**: 每次迭代会依据当前状态判断采用哪种方式更合适: - 如果能显著提升逼近速度,则优先考虑“二次插值”; - 若无法满足安全条件或失败,则退而求其次使用“割线法”; - 最保守的方式就是“二分法”,当其他两种都不可行时必选此方案。 3. **逐步计算直至满足终止准则**: - 计算新的试探点c及其对应的y=f(c) - 更新端点[a,b] 根据新得到的信息调整左右边界位置, 并重新评估是否继续循环还是结束程序返回结果. 最终经过多次迭代运算得出的结果应该是接近于黄金分割比φ≈1.618...附近某一点作为近似解答。 --- 以下是伪代码形式表达整个过程: ``` function brent(f,a,b,tol): fa,fm,fb=f(a),f((a+b)/2.),f(b); if abs(fm)<=tol: return m; while b-a > tol: ... return root_value ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值