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

【SICP练习】31 练习1.37

时间:2015-02-17 22:27:53      阅读:196      评论:0      收藏:0      [点我收藏+]

标签:sicp



练习1.37

根据题目中的意思通过观察得到k项有项连分式的一种表达方式:

f=N1/(D1+(N2/(…+Nk/Dk)))

这个式子可以不断展开,但如果我们把每一个”+”后面的式子记作T(i)。不对,我们应该将每一个N/D记作T(i),因为这组式起始于N/D,且中止与N/D。计N1/D1T(1)N2/D2T(2),Nk/DkT(k)。在数学上可能不会联想到递归,而是联想到一个表达式,以谋求能够化简。但在这里大家应该都发现了这是一个递归过程。

接着我们来写出题目要求的(cont-fracn d k)

(define (cont-frac n d k)

    (define (cont-frac-t i)

        (if (= k i)

           (/ (n k) (d k))

           (/ (n i) (+ (d i) (cont-frac-t (+ i1))))))

    (cont-frac-t 1))

下面我们利用题中给出的过程来写一个黄金分割率的定义。

(define(cont-frac-golden-ratio k)

    (+ 1 (cont-frac (lambda (i) 1.0)

                   (lambda (i) 1.0)

                   k)))

接下来我们按照题中的要去来测试一下刚才所写的过程。

通过几次测试之后得出的结论是当k11时即有十进制的4位精度,而且k值越大精度越大。

(cont-frac-golden-ratio 11)

;Value: 1.6180555555555556

好了,现在我们还要来考虑迭代版本。我们只要记住迭代的空间需求是不变的,在题中的式子中,我们可以发现如下规律:

T(1)=N1/(D1+T(2))

T(2)=N2/(D2+T(3))

T(3)=N3/(D3+T(4))

……

T(k-2)=N(k-2)/(D(k-2)+T(k-1))

T(k-1)=N(k-1)/(D(k-1)+T(k))

Tk=Nk/Dk

根据如下的表达式的规律,我们可以写出迭代版本的cont-frac函数,像前面博文中所讲,我们依旧应该添加一个存储器,比如other。而不断变化的other就是前面表达式从依次往下的最右边的T

(define (cont-frac n d k)

    (define (cont-frac-t i other)

        (if (= i 0)

           other

           (cont-frac-t (- i 1) (/ (n i)

                             (+ (d i) other)))))

    (cont-frac-t (- k 1) (/ (n k)

                        (d k))))

这道习题我们就此完成了,再接再厉完成这一节的最后一题。不对,反面还有一题。

【SICP练习】31 练习1.37

标签:sicp

原文地址:http://blog.csdn.net/nomasp/article/details/43867397

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