SICP 习题 (1.38)解题总结

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 plain copy 在CODE上查看代码片 派生到我的代码片
  1. (lambda (i) 1.0)  

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

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

[plain] view plain copy 在CODE上查看代码片 派生到我的代码片
  1. (lambda (i)   
  2.      (if (= (remainder (+ i 1) 3) 0)  
  3.      (* (/ (+ i 1) 3) 2)  
  4.      1))  


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

[plain] view plain copy 在CODE上查看代码片 派生到我的代码片
  1. (cont-frac   
  2.    (lambda (i) 1.0)  
  3.    (lambda (i)   
  4.      (if (= (remainder (+ i 1) 3) 0)  
  5.      (* (/ (+ i 1) 3) 2)  
  6.      1))  
  7.    k)  

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


[plain] view plain copy 在CODE上查看代码片 派生到我的代码片
  1. (define (e-test k)  
  2.  (+ 2 (cont-frac   
  3.    (lambda (i) 1.0)  
  4.    (lambda (i)   
  5.      (if (= (remainder (+ i 1) 3) 0)  
  6.      (* (/ (+ i 1) 3) 2)  
  7.      1))  
  8.    k)))  


结束!!

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值