您的位置:首页 > 其它

SICP 习题 (1.38)解题总结

2014-06-18 08:16 281 查看
SICP 习题1.38 紧跟着习题1.37的方向,要求我们用习题1.37中定义的cont-frac过程计算数学家欧拉大师在论文De Fractionibus Continuis 中提到的e-2的连分式。说实话,我不知道论文De Franctionibus Continuis讲的是什么,我甚至不知道论文的题目是什么意思。不过,这一切都不能阻止我这个数学盲去解答这道SICP习题。

仔细阅读题目,我们可以发现题目要求我们计算的是下面这样的无穷连分式:



其中N永远等于1, D等于1 ,  2 , 1 , 1 , 4 , 1 , 1 , 6 , 1 , 1 , 8 …

所以我们调用cont-frac时给的N比较简单,就是一个永远返回1的lambda过程,过程如下:

[plain]
view plaincopy





(lambda (i) 1.0)  

而调用cont-fract时给的D复杂一点点,D是一个lambda过程,能根据下标生成1 ,  2 , 1 , 1 , 4 , 1 , 1 , 6 , 1 , 1 , 8 …这样的数列。
这个数列规律还是比较明显,稍微费点脑筋可以做个lambda过程来生成,我做的过程如下:

[plain]
view plaincopy





(lambda (i)   
     (if (= (remainder (+ i 1) 3) 0)  
     (* (/ (+ i 1) 3) 2)  
     1))  

所以,下面的调用就可以得出e-2的值了:

[plain]
view plaincopy





(cont-frac   
   (lambda (i) 1.0)  
   (lambda (i)   
     (if (= (remainder (+ i 1) 3) 0)  
     (* (/ (+ i 1) 3) 2)  
     1))  
   k)  

以上结果再+2就等于e了,我的完整测试过程如下:

[plain]
view plaincopy





(define (e-test k)  
 (+ 2 (cont-frac   
   (lambda (i) 1.0)  
   (lambda (i)   
     (if (= (remainder (+ i 1) 3) 0)  
     (* (/ (+ i 1) 3) 2)  
     1))  
   k)))  

结束!!
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: