C语言实现导数近似计算：代码运行异常，请求问题排查

代码中的核心问题 #

• 类型错误：hn函数中把hn声明为int类型，pow(2,-n)*h0的结果是浮点数，强制转换为整数会直接截断小数部分，导致hₙ计算完全错误，应该改为double类型。
• 参数传递错误：func函数的参数定义是(double h, double x)，但调用时却传入(h0, n)（比如func(h0,n)），参数类型和含义完全不匹配，正确的做法是先计算出hₙ，再传入函数计算差商。
• 未使用函数指针：func函数直接调用全局的f，忽略了diff函数传入的函数指针，导致diff失去通用性，无法适配其他函数的导数计算。
• 终止逻辑错误：循环结束后返回的是Φ(hₙ₊₁)，但题目要求返回满足条件的Φ(hₙ)；同时初始的分支判断逻辑冗余，应该在每次迭代中动态判断终止条件，而不是提前分支。
• 缺少参数校验：未按照题目要求用assert检查h0 > 0和tau > 0，无法保证输入合法性。

修正后的代码 #

#include <stdio.h>
#include <assert.h>
#include <math.h>
#include <stdbool.h>

// 目标函数示例：f(x) = x² + 2
double f(double x) {
    return x*x + 2;
}

// 计算hₙ = 2^(-n) * h0
double calc_hn(double h0, int n) {
    return pow(2, -n) * h0;
}

// 计算差商Φ(h) = (f(x+h) - f(x))/h，使用传入的函数指针
double calc_phi(double (*f)(double), double x, double h) {
    return (f(x + h) - f(x)) / h;
}

double diff(double (*f)(double), double x, double h0, double tau) {
    // 校验输入合法性
    assert(h0 > 0 && tau > 0);
    
    int n = 0;
    double phi_n = calc_phi(f, x, calc_hn(h0, n));
    double phi_n1 = calc_phi(f, x, calc_hn(h0, n+1));
    
    // 循环直到满足终止条件
    while (1) {
        double diff_val = fabs(phi_n - phi_n1);
        bool condition_met = false;
        
        if (fabs(phi_n) <= tau) {
            condition_met = (diff_val <= tau);
        } else {
            condition_met = (diff_val <= tau * fabs(phi_n));
        }
        
        if (condition_met) {
            break;
        }
        
        // 更新n和对应的差商值
        n++;
        phi_n = phi_n1;
        phi_n1 = calc_phi(f, x, calc_hn(h0, n+1));
    }
    
    return phi_n;
}

int main(void) {
    double x, h0, tau;
    printf("请输入x、h0、tau的值（用空格分隔）：");
    // 读取输入，处理可能的输入错误
    if (scanf("%lf %lf %lf", &x, &h0, &tau) != 3) {
        printf("输入错误！\n");
        return 1;
    }
    
    double result = diff(f, x, h0, tau);
    printf("导数近似值：%.6f\n", result);
    return 0;
}

修正说明 #

1. 修复了hₙ的类型错误，确保计算精度正确。
2. 重构了差商计算函数，使用传入的函数指针，保证diff函数的通用性。
3. 修正了终止逻辑：每次迭代判断当前Φ(hₙ)和Φ(hₙ₊₁)是否满足终止条件，满足时返回Φ(hₙ)，符合题目要求。
4. 添加了assert校验h0和tau的合法性。
5. 优化了主程序的输入提示和错误处理，提升用户体验。

内容的提问来源于stack exchange，提问作者annnna