<p><span style="font-size:14px;">sqrt()函数,是绝大部分语言支持的常用函数,它实现的是开方运算;开方运算最早是在我国<span style="font-family:arial, '宋体', sans-serif;line-height:24px;">魏晋时数学家刘徽所著的</span><span style="font-family:arial, '宋体', sans-serif;"><span style="line-height:24px;">《九章算术》被提及。今天写了几个函数加上国外大神的几个神级程序带大家领略sqrt的神奇之处。</span></span></span></p>
<p><span style="font-size:14px;"><br></span></p>
<h1><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;">1.古人算法(暴力法) </span></span></h1>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;">原理:从0开始0.00001,000002...一个一个试,直到找到x的平方根,代码如下:</span></span></p>
<p></p>
<pre class="blockcode"><code class="language-java">public class APIsqrt {
static double baoliSqrt(double x) {
final double _JINGDU = 1e-6;
double i;
for (i = 0; Math.abs(x - i * i) > _JINGDU; i += _JINGDU)
;
return i;
}
public static void main(String[] args) {
double x = 3;
double root = baoliSqrt(x);
System.out.println(root);
}
}
</code></pre><br>
测试结果:
<p></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;">1.7320509999476947<br></span></span></p>
<h1><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;">2.牛顿迭代法</span></span></h1>
<p></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;">计算机科班出身的童鞋可能首先会想到的是《数值分析》中的牛顿迭代法求平方根。原理是:随意选一个数比如说8,要求根号3,我们可以这么算:</span></span></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;">(8 + 3/8) =4.1875</span></span></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;">(<span style="font-family:arial, '宋体', sans-serif;line-height:24px;">4.1875</span> + 3/<span style="font-family:arial, '宋体', sans-serif;line-height:24px;">4.1875</span>)
=</span><span style="font-family:arial, '宋体', sans-serif;line-height:24px;">2.4519</span><br></span></span></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;">(<span style="font-family:arial, '宋体', sans-serif;line-height:24px;">2.4519</span>
+ 3/<span style="font-family:arial, '宋体', sans-serif;line-height:24px;">2.4519</span>) =</span><span style="font-family:arial, '宋体', sans-serif;line-height:24px;">1.837</span><br></span></span></span></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;">(<span style="font-family:arial, '宋体', sans-serif;line-height:24px;">1.837</span>+
3/<span style="font-family:arial, '宋体', sans-serif;line-height:24px;">1.837</span>) =</span><span style="font-family:arial, '宋体', sans-serif;line-height:24px;">1.735</span><br></span></span></span></span></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;">做了4步基本算出了近似值了,这种迭代的方式就是传说中的牛顿迭代法了,代码如下:</span></span></span></span></span></p>
<p><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"><span style="font-family:arial, '宋体', sans-serif;line-height:24px;"></span></span></span></span></span></p>
<pre class="blockcode"><code class="language-java">public class APIsqrt {
static double newtonSqrt(double x) {
if (x < 0) {
System.out.println("负数没事开什么方");
return -1;
}
if (x == 0)
return 0;
double _avg = x;
double last_avg = Double.MAX_VALUE;
final double _JINGDU = 1e-6;
while (Math.abs(_avg - last_avg) > _JINGDU) {
last_avg = _avg;
_avg = (_avg + x / _avg) / 2;
}
return _avg;
}
public static void main(String[] args) {
double x = 3;
double root = newtonSqrt(x);
System.out.println(root);
}
}</code></pre>
<p></p>
<p>测试结果:</p>
<p>1.7320508075688772</p>
<br><p></p>
<h1><span style="font-family:arial, '宋体', sans-serif;font-size:14px;"><span style="line-height:24px;"><span style="line-height:24px;"><span style="line-height:24px;"><span style="line-height:2 |
|