C#,数值计算,求平方根之巴比伦算法(Babylonian algorithm)的源代码

平方根的巴比伦算法。

1 巴比伦算法介绍一

巴比伦算法可能算是最早的用于计算sqrt{S}的算法之一,因为其可以用牛顿法导出,因此在很多地方也被成为牛顿法。其核心思想在于为了计算x的平方根,可以从某个任意的猜测值g开始计算。在真实的运算中,我们往往将g直接设置为x,不过也可以选择其他任何的正数值。那么其计算的迭代过程为:

1.如果猜测值g已经足够接近于正确的平方根,算法结束,函数将g作为结果返回。

2.如果猜测值g不够精确,那么使用g和x/g的平均值作为新的猜测值。因为这两个值中的一个小于确切的平方根,另一个则大于确切的平方根,选择平均值有助于你得到一个更接近于正确答案的值。

3.把新的猜测值赋予给变量g,重复第一步的判断。

2 巴比伦算法介绍二

  1. 假设要求a的平方根,先找到一个近似值x0,如a/2;

  2. 被开方数a除以近似值x0;

  3. 求第一步和第二步得到的值的平均数;

  4. 令平均数为新的近似值,返回第二步,循环计算(根据精度要求,进行循环控制)

经过前面四个步骤,可以得到如上公式;

对公式的理解:可以把xn理解为x0,xn+1理解为x1,然后套用公式即可。

3 核心代码

/// <summary>

/// 计算 v 的平方根

/// 据说来自于远古巴比伦时代的算法

/// </summary>

/// <param name="v"></param>

/// <returns></returns>

public static double Sqrt(double v)

{

const double e = 0.000001;

double x = v;

double y = 1.0;

while ((x - y) > e)

{

x = (x + y) / 2;

y = v / x;

}

return x;

}


POWER BY TRUFFER.CN

相关推荐
惊鸿一博4 分钟前
Plücker直线_ 用两个三维向量描述一条空间直线
线性代数·算法·直线·几何表示
估值探索者8 分钟前
【Python量化技术指标实战 02】MACD金叉也分强弱?手工算一次 DIF/DEA 就懂了
java·开发语言·python·股票数据·股票数据api接口
讳疾忌医丶9 分钟前
编译器眼中的 C++ 循环:从指针别名、IEEE 754 结合律到 SIMD 语义证明
开发语言·c++
小道士写程序14 分钟前
Rust 中的字符串与切片(String 与 &str)
算法
重生六道20 分钟前
B站视频时长统计:一段JS脚本搞定分P/合集总时长(附倍速计算)
开发语言·javascript·音视频
Tisfy30 分钟前
LeetCode 3863.将一个字符串排序的最小操作次数:分类讨论(一次遍历)——简单易懂版
算法·leetcode·字符串·题解·分类讨论
XiYang-DING1 小时前
P5528.带过期时间的 LRU 缓存
java·开发语言·缓存
小徐不徐说1 小时前
C++ 机器视觉入门:从零搭建 AOI 开发环境(OpenCV + CMake)
开发语言·c++·opencv
Young Doro1 小时前
从平衡表征到最小充分因果表征:端到端神经网络因果效应估计
算法
0+1111 小时前
算法 --递归
算法·leetcode