码迷,mamicode.com
首页 > 其他好文 > 详细

SICP 习题 (1.38)解题总结

时间:2014-06-22 21:55:22      阅读:222      评论:0      收藏:0      [点我收藏+]

标签:style   class   blog   code   http   ext   

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


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

bubuko.com,布布扣


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


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

(lambda (i) 1.0)

而调用cont-fract时给的D复杂一点点,D是一个lambda过程,能根据下标生成1 ,  2 , 1 , 1 , 4 , 1 , 1 , 6 , 1 , 1 , 8 …这样的数列。

这个数列规律还是比较明显,稍微费点脑筋可以做个lambda过程来生成,我做的过程如下:

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


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

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

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


(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)))


结束!!


SICP 习题 (1.38)解题总结,布布扣,bubuko.com

SICP 习题 (1.38)解题总结

标签:style   class   blog   code   http   ext   

原文地址:http://blog.csdn.net/keyboardota/article/details/23475541

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!