黄金连分数
题目
黄金分割数0.61803… 是个无理数,这个常数十分重要,在许多工程问题中会出现。有时需要把这个数字求得很精确。
对于某些精密工程,常数的精度很重要。也许你听说过哈勃太空望远镜,它首次升空后就发现了一处人工加工错误,对那样一个庞然大物,其实只是镜面加工时有比头发丝还细许多倍的一处错误而已,却使它成了“近视眼”!!
言归正传,我们如何求得黄金分割数的尽可能精确的值呢?有许多方法。
比较简单的一种是用连分数:
1
-----------------------------
1
1 + ------------------------
1
1 + -----------------
1
1 + ---------
1 + ...
这个连分数计算的“层数”越多,它的值越接近黄金分割数。
请你利用这一特性,求出黄金分割数的足够精确值,要求四舍五入到小数点后100位。
小数点后3位的值为:0.618
小数点后4位的值为:0.6180
小数点后5位的值为:0.61803
小数点后7位的值为:0.6180340
(注意尾部的0,不能忽略)
你的任务是:写出精确到小数点后100位精度的黄金分割值。
注意:尾数的四舍五入! 尾数是0也要保留!
显然答案是一个小数,其小数点后有100位数字,请通过浏览器直接提交该数字。
注意:不要提交解答过程,或其它辅助说明类的内容。
题目分析
第一项: 1
-----------------
2
第二项: 2
-----------------
3
第三项: 3
-----------------
5
第四项: 5
-----------------
8
第三项: 8
-----------------
13
… \ldots … … \ldots … … \ldots … … \ldots … … \ldots …
第n项: f(n-2)
-----------------
f(n-1)
典型的斐波那契数列的运用
第一项:1
第二项:1
第三项:2
第四项:3
第五项:5
第六项:8
… \ldots … … \ldots … … \ldots … … \ldots … … \ldots …
第n项:f(n-1) + f(n-2)
解题思路
第一:我们知道了,该题需要使用斐波那契来求解,但是我们不能使用递归,而是应该使用迭代来实现代码。
高精度使用的静态变量:
static BigInteger ONE BigInteger常数。 表示的是1
static BigInteger TEN BigInteger常数十。 表示的是10
static BigInteger ZERO BigInteger常数零。 表示的是0
BigInteger add(BigInteger val) 返回值为 (this + val) 。
BigInteger divide(BigInteger val) 返回值为 (this / val) 。
String toPlainString() 返回没有指数字段的此 BigDecimal的字符串表示形式 。
String substring(int intIndex, int lastIndex) 从下标intIndex到lastIndex开始截取。
误区分析
第一:在这里double为浮点数,在Java中占用空间64位,也就是八个字节。在具体的数中能够精确到小数点的后16位。
因此,我们需要使用高精度运算,即BigInteger和BigDecimal,大数运算。
代码如下
public class Analyses4_黄金连分数 {
public static void main(String[] args) {
// TODO Auto-generated method stub
System.out.println(goldNumber(100));
}
//这里需要的是小数,故此使用浮点类型double
public static double goldNumber(int n) {
double a = 1.0;
double b = 1.0;
//使用for循环进行100次迭代,注意,此处需要从3开始
for(int i = 3; i < 100; i++) {
double t = b;
b = a + b;
a = t;
}
double result = a / b;
return result;
}
}
程序的答案:0.6180339887498948
由于程序需要的是需要精确到小数点后100位,可以在100位采用四舍五入的方法,但是由于黄金连分数是一个无理数,精确到小数点后100位,则黄金连分数逐渐趋于稳定。
随着for循环范围的不同,所得到的黄金连分数也有差异:
循环到100:
0.618033988749894848204586834365638117720299848718891653939793515920494891335798525795260376280960598
循环到300:
0.6180339887498948482045868343656381177203091798057628621354486227052604628189024497072072041893911374
循环到500:
0.6180339887498948482045868343656381177203091798057628621354486227052604628189024497072072041893911374
循环到800:
0.6180339887498948482045868343656381177203091798057628621354486227052604628189024497072072041893911374
正确的代码如下:
import java.math.BigDecimal;
import java.math.BigInteger;
public class Analyses4_黄金连分数 {
public static void main(String[] args) {
// TODO Auto-generated method stub
goldNumber();
}
public static void goldNumber() {
BigInteger a = BigInteger.ONE;//表示的是整数1
BigInteger b = BigInteger.ONE;//表示的是整数1
for(int i = 3; i < 800; i++) {
BigInteger t = b;
b = a.add(b);//超高精度BigInteger方法,将a与b相加
a = t;
}
//采取超高精度的BigDecimal的面向对象的方式来为黄金连分数进行除运算,并且指定a和b的范围,将结果四舍五入处理
BigDecimal result = new BigDecimal(a,102).divide(new BigDecimal(b,102), BigDecimal.ROUND_HALF_DOWN );
System.out.println(result.toPlainString().substring(0, 102));//以字符串BigDecimal开始获取,并且截取0到101位的值
}
}
正确的结果:题目需要的是小数点后一百位,故此采用四舍五入的方法,其结果如下:
0.6180339887498948482045868343656381177203091798057628621354486227052604628189024497072072041893911375
借助的资料
链接:https://pan.baidu.com/s/1h7xKHDlBhpmLhmMDEyCJqA
提取码:iqsp
在这里,建议大家坚持阅读Java API帮助文档,从上面你能学到课本上学习不到的知识。善用工具是我们学习知识的前提,希望能够帮助到你。
资料推荐:《Java 基础入门》黑马程序员 / 编著