Scheme(DrRacket)中实现指定公式生成α值列表的技术求助

作为刚接触Scheme/Racket的新手，能尝试上手数值计算的公式实现已经超棒啦！咱们一步步拆解这个问题，把你的alpha_values函数写出来～

首先先把你给出的公式简化一下，这样代码写起来更清晰：

αᵢ = (3/h)((aᵢ₊₁ – aᵢ ) – (aᵢ – aᵢ₋₁ ))
整理后可以得到：αᵢ = (3/h) * (aᵢ₊₁ - 2aᵢ + aᵢ₋₁)

函数实现思路 #

你的函数需要接收三个参数：

• y_values：对应公式里的a序列（列表形式，比如(a₀ a₁ a₂ ... aₙ)）
• h：公式中的常数分母
• n：指定计算i从1到n-1的α值（注意：这要求y_values至少包含a₀到aₙ，否则会出现索引越界）

Racket里用for/list循环来生成结果列表非常直观，适合新手理解，下面是两种实现方式：

方式一：严格对应原公式

(define alpha_values
  (lambda (y_values h n)
    (for/list ([i (in-range 1 n)])  ; 生成1到n-1的整数序列，对应公式里的i
      (/ (* 3 (- (- (list-ref y_values (+ i 1))  ; 取aᵢ₊₁
                    (list-ref y_values i))       ; 减去aᵢ
                 (- (list-ref y_values i)        ; 减去(aᵢ - aᵢ₋₁)
                    (list-ref y_values (- i 1)))))
         h))))

方式二：用简化后的公式（更简洁）

(define alpha_values
  (lambda (y_values h n)
    (for/list ([i (in-range 1 n)])
      (/ (* 3 (+ (list-ref y_values (+ i 1))       ; aᵢ₊₁
                 (* -2 (list-ref y_values i))      ; -2aᵢ
                 (list-ref y_values (- i 1))))     ; +aᵢ₋₁
         h))))

测试例子 #

咱们用一个简单的序列来验证：
假设y_values = '(0 1 3 6 10)（对应a₀到a₄），h=1，n=4（计算i=1、2、3）。
按照公式计算：

• α₁ = (3/1)(3 - 21 + 0) = 3
• α₂ = (3/1)(6 - 23 + 1) = 3
• α₃ = (3/1)(10 - 26 + 3) = 3

调用函数：

(alpha_values '(0 1 3 6 10) 1 4)

会返回(3 3 3)，和预期结果一致。

递归版本（可选） #

如果你想练习递归写法，这里也提供一个递归实现的版本，逻辑和循环版本完全一致：

(define alpha_values
  (lambda (y_values h n)
    ; 定义内部递归辅助函数
    (define (calculate-alpha i)
      (if (>= i n)
          '()  ; 递归终止条件：i达到n时返回空列表
          (cons (/ (* 3 (- (- (list-ref y_values (+ i 1))
                              (list-ref y_values i))
                           (- (list-ref y_values i)
                              (list-ref y_values (- i 1)))))
                   h)
                (calculate-alpha (+ i 1)))))  ; 递归计算下一个i
    (calculate-alpha 1)))  ; 从i=1开始计算

新手注意事项 #

• 确保y_values的长度足够：因为i最大是n-1，需要访问a_{n}，所以y_values的长度至少是n+1，否则list-ref会抛出索引越界的错误。
• Racket的列表索引是从0开始的，这正好和公式里的a₀、a₁下标对应，不用额外调整。

内容的提问来源于stack exchange，提问作者user-2147482428