Производная схема - PullRequest
Новая
       53

Производная схема

1 голос
/ 14 октября 2019

Функция d должна возвращать производную введенного выражения в упрощенной форме. Я понимаю, что есть функция производной схемы, но вместо этого я использую списки, чтобы бросить вызов себе. Я новичок в языке. Я написал контрольные примеры для проверки своего решения, но ожидаемый результат неправильный и не упрощенный. 

(define d
  (λ (e)
    (cond ((number? e) 0)
      ((equal? e 'x) 1)
      (else
       (let ((op (car e)) (args (cdr e)))
         (apply (lookup op d-op-table) args))))))


(define d-op-table
  (list(list '+   (λ (u v)(list '+ (d u) (d v))))
       (list '-   (λ (u v)(list '- (d u) (d v))))
       (list '*   (λ (u v)(list '+ (list '* u (d v))(list '* v (d u)))))))
       (list 'sin (λ (u)(list '*( list 'cos (d u)))))
       (list 'cos (λ (u)(list '*( list '-sin (d u)))))
       (list 'log (λ (u)(list '* (list '/1 u) (d u))))
       (list 'exp (λ (u)(list '* (d u)(list 'exp u))))
       (list 'expt (λ (u v) (list 'expt (list '* v  u) (- v 1))))

(define lookup
  (λ (op
      table)
    (if (equal? op (caar table))
    (cadar table)
    (lookup op (cdr table)))))



;; Test cases

;; (d '(* (+ x 4) (+ x -7)))
;; '(+ (* (+ x 4) (+ 1 0)) (* (+ x -7) (+ 1 0)))

;; (d '(* x (* x (* x (* x x)))))
;; '(+ (* x (+ (* x (+ (* x (+ (* x 1) (* x 1))) (* (* x x) 1))) (* (* x (* x x)) 1))) (* (* x (* x (* x x))) 1))

Пример:

Путем ввода функции 

(d '(* (+ x 1) (+ x -1)))

Я ожидаю 

(+ (* (+ x 1) (+ 1 0)) (* (+ 1 0) (+ x -1)))

но я получаю 

(+ (* (+ x 1) (+ 1 0)) (* (+ x -1) (+ 1 0)))

1 Ответ

1 голос
/ 16 октября 2019

Это потому, что случай умножения говорит так: 

(list '+ (list '* u (d v)) (list '* v (d u)))

То есть d(u*v) = u * dv + v * du.

Второй термин «перевернут» по сравнению с тем, что вы говорите, хотите,что d(u*v) = u * dv + du * v: 

(list '+ (list '* u (d v)) (list '* (d u) v))
...