Защо звуковият сигнал не приближава правилно моята непрекъсната дроб?

Чета още SICP и съм заседнал на упражнение 1.3.8. Моят код работи правилно за приблизително 1/phi, но не работи за приблизително e - 2.

(define (cont-frac n d k)
  (define (frac n d k)
        (if (= k 0)
            1.0
            (+ (d k) (/ (n (+ k 1)) (frac n d (- k 1))))))
    (/ (n 1) (frac n d k)))

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

(define (1-over-phi)
  (cont-frac (lambda (i) 1.0)
             (lambda (i) 1.0)
             100))

Вместо да получавам .7 бла бла бла за e-2, получавам .5 бла бла нещо. Не мога да разбера защо. Почти съм сигурен, че съм дефинирал правилно "d" във функцията "eulers-e-2".

Редактиране: Благодаря, момчета, изчислявах го наобратно. Ето фиксирания код.

(define (cont-frac n d k)
  (define (frac n d i)
        (if (= k i)
            (d i)
            (+ (d i) (/ (n (+ i 1)) (frac n d (+ i 1))))))
    (/ (n 1) (frac n d 1)))

person hraesvelgr    schedule 06.10.2010    source източник
comment
@Seth: това би било „толкова много скоби“.   -  person High Performance Mark    schedule 06.10.2010


Отговори (1)


Изглежда, че изчислявате следното:

N1/(D100 + (N101/ D99 + N100/(D98 + N99/(..))))

Вместо

N1/(D1 + N2/(D2 + ...))

Тъй като N и D са еднакви (всички единици) за 1/phi, вие получавате правилния отговор там.

person Community    schedule 06.10.2010