Question

我写了一个过程golden-ratio来计算黄金分割率。第一个参数是我们开始计算黄金分割率的数字，第二个参数是我们将重复递归的次数。

(defun golden-ratio (start-num times)
  (if (eq times 0)
      start-num
    (golden-ratio (+ (/ 1.0 start-num) 1) (- times 1))))

(golden-ratio 30 250)

得到结果1.618033988749895。但是当我尝试

(golden-ratio 30 300)

我明白了 Lisp nesting exceeds `max-lisp-eval-depth' 我的问题是，还有另一种算法能让我深入到250左右吗？

Answer 1

也许更好的方法是将注意力集中在所需的精度上，而不是重复次数。这样，您可以迭代地改善部分结果，直到改善率低于准确性为止。我发现do宏特别适合这种类型的计算。您可以按以下方式解决该问题：

(defun golden-ratio (start delta)
  "Calculate the golden ratio starting at START with accuracy
DELTA."
  (do ((old (float start))
       (new (1+ (/ (float start))) (1+ (/ new))))
      ((< (abs (- old new)) delta) new)
    (setq old new)))

您可以这样做：

ELISP> (golden-ratio 1 1.0e-40)
1.618033988749895

Elisp中的高效算法以计算黄金分割率

1 个答案: