王朝百科
分享
 
 
 

布伦特法

王朝百科·作者佚名  2010-06-11  
宽屏版  字体: |||超大  

the method of Brent

布伦特(Brent)方法是在二分法或试位法的基础上,借助二次插值方法进行加速,有利用反插值方法来简化计算而形成的一种方法,初步感觉有些小题大做,仔细分析又不无道理。

假如知道f(x)的零点x’在一个不太大的区间[x0,x1]内,而且已知f(x)在区间的端点a,b处的函数值y0=f(x0),y1=f(x1),f(x0) f(x1)<0,以及f(x)=0在(x0,x1)内的近似解x=x2和f(x2),接下来利用这已知的三个点以及它们所对应的函数值作插值抛物线。与Muller方法不同的是,把

x0 ,x1 ,x2 与f(x0), f(x1), f(x2)的对应关系反过来用,相当于用y0 ,y1 ,y2 替代方程组

f(x0)= a(x0-x2) ^2+b(x0-x2)+c **************** (1)

f(x1)= a(x1-x2) ^2+b(x1-x2)+c **************** (2)

f(x2)= a(x2-x2) ^2+b(x2-x2)+c **************** (3)

中的x0 ,x1 ,x2 ,而方程组的常数项f(x0), f(x1), f(x2) 则用这里的x0 ,x1 ,x2替换。所以对应的插值抛物线的一般形式为

x= a(y-y2) ^2+b(y0-y2)+c **************** (4)

x0= a(y0-y2) ^2+b(y0-y2)+c **************** (5)

x1= a(y1-y2) ^2+b(y1-y2)+c **************** (6)

x2= a(y2-y2) ^2+b(y2-y2)+c **************** (7)

由(5)(6)(7)式可以得c=x2,利用Muller方法,可以写成a,b 的表达式。

在(4)中令y=0,可以得到下一个近似点为

X=x2+(ay2 ^2+by2) **************** (8)

其中ay2 ^2+by2相当于校正项。

Brent方法的下一个规则是,如果得到的x仍然在区间[x0,x1]内,则用x2根据f(x2)

的符号替换x0或x1 ,用x替换x2;如果所得到的x不在区间[x0,x1]内,则暂时放弃反抛物线插值法,继续用二分法或试位法。

实际上,我们可以利用二分法所得到函数顺便采用反插值方法试探一下,看能否捡到一些便宜。

在这种平常心态下,我们也可以得到类似的效果,而且方法还简单得多。

 
 
免责声明:本文为网络用户发布,其观点仅代表作者个人观点,与本站无关,本站仅提供信息存储服务。文中陈述内容未经本站证实,其真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
如何用java替换看不见的字符比如零宽空格&#8203;十六进制U+200B
 干货   2023-09-10
网页字号不能单数吗,网页字体大小为什么一般都是偶数
 干货   2023-09-06
java.lang.ArrayIndexOutOfBoundsException: 4096
 干货   2023-09-06
Noto Sans CJK SC字体下载地址
 干货   2023-08-30
window.navigator和navigator的区别是什么?
 干货   2023-08-23
js获取referer、useragent、浏览器语言
 干货   2023-08-23
oscache遇到404时会不会缓存?
 干货   2023-08-23
linux下用rm -rf *删除大量文件太慢怎么解决?
 干货   2023-08-08
刀郎新歌破世界纪录!
 娱乐   2023-08-01
js实现放大缩小页面
 干货   2023-07-31
生成式人工智能服务管理暂行办法
 百态   2023-07-31
英语学习:过去完成时The Past Perfect Tense举例说明
 干货   2023-07-31
Mysql常用sql命令语句整理
 干货   2023-07-30
科学家复活了46000年前的虫子
 探索   2023-07-29
英语学习:过去进行时The Past Continuous Tense举例说明
 干货   2023-07-28
meta name="applicable-device"告知页面适合哪种终端设备:PC端、移动端还是自适应
 干货   2023-07-28
只用css如何实现打字机特效?
 百态   2023-07-15
css怎么实现上下滚动
 干货   2023-06-28
canvas怎么画一个三角形?
 干货   2023-06-28
canvas怎么画一个椭圆形?
 干货   2023-06-28
canvas怎么画一个圆形?
 干货   2023-06-28
canvas怎么画一个正方形?
 干货   2023-06-28
中国河南省郑州市金水区蜘蛛爬虫ip大全
 干货   2023-06-22
javascript简易动态时间代码
 干货   2023-06-20
感谢员工的付出和激励的话怎么说?
 干货   2023-06-18
 
>>返回首页<<
 
 
静静地坐在废墟上,四周的荒凉一望无际,忽然觉得,凄凉也很美
© 2005- 王朝网络 版权所有