C89下多类型滤波器通用幅值响应计算的代码去重方案问询

针对你在C89环境中面临的重复编写滤波器幅值响应计算函数的问题，以下几种方案可以有效减少代码冗余，同时兼顾兼容性和可读性：

1. 预处理器宏生成专属函数 #

利用C预处理器的字符串拼接与宏替换能力，编写一个通用宏来生成不同滤波器对应的幅值响应计算函数。这种方式无需额外工具，且生成的代码类型安全，避免void指针的强制转换问题。

示例实现：

// 通用宏：生成指定滤波器的幅值响应计算函数
#define GENERATE_MAG_RESPONSE(FILTER_TYPE, PROCESS_FUNC) \
void compute_##FILTER_TYPE##_mag_response(FILTER_TYPE *filter, int n_points, double *freq, double *mag) { \
    /* 通用逻辑：生成测试输入 */ \
    double input[n_points]; \
    generate_sweep_input(input, n_points, SAMPLE_RATE); \
    \
    /* 调用对应滤波器的处理函数 */ \
    PROCESS_FUNC(filter, input, input, n_points); \
    \
    /* 通用FFT与幅值计算逻辑 */ \
    fftw_plan plan = fftw_create_plan(n_points, FFTW_FORWARD, FFTW_ESTIMATE); \
    fftw_execute(plan, input, mag); \
    compute_normalized_magnitude(mag, n_points); \
    generate_freq_axis(freq, n_points, SAMPLE_RATE); \
    \
    fftw_destroy_plan(plan); \
}

// 生成不同滤波器的函数
GENERATE_MAG_RESPONSE(OnePoleFilter, onepole_process)
GENERATE_MAG_RESPONSE(BiquadFilter, biquad_process)
GENERATE_MAG_RESPONSE(AllPassFilter, allpass_process)

宏展开后会生成独立的类型安全函数，通用逻辑仅需维护一次，且调用方式与手写函数完全一致。

2. 通用核心函数+类型安全封装层 #

提取所有滤波器共有的FFT、幅值计算逻辑到一个通用核心函数，通过函数指针传递不同滤波器的处理逻辑，再为每个滤波器编写极简的类型安全封装入口。

示例实现：

// 通用核心函数：处理公共逻辑
static void compute_mag_response_core(void *filter, 
                                      void (*process)(void*, double*, double*, int),
                                      int n_points, double *freq, double *mag) {
    double input[n_points];
    generate_impulse_input(input, n_points);
    
    // 调用传入的滤波器处理函数
    process(filter, input, input, n_points);
    
    // 公共FFT与幅值计算
    fftw_plan plan = fftw_create_plan(n_points, FFTW_FORWARD, FFTW_ESTIMATE);
    fftw_execute(plan, input, mag);
    compute_magnitude_db(mag, n_points);
    generate_freq_axis(freq, n_points, SAMPLE_RATE);
    
    fftw_destroy_plan(plan);
}

// 类型安全的封装入口（每个滤波器仅需一行）
void compute_onepole_mag_response(OnePoleFilter *filter, int n_points, double *freq, double *mag) {
    compute_mag_response_core(filter, (void(*)(void*, double*, double*, int))onepole_process, n_points, freq, mag);
}

void compute_biquad_mag_response(BiquadFilter *filter, int n_points, double *freq, double *mag) {
    compute_mag_response_core(filter, (void(*)(void*, double*, double*, int))biquad_process, n_points, freq, mag);
}

这种方案将通用逻辑集中维护，封装层仅负责类型转换，既保证了调用时的类型安全，又避免了大量代码重复。

3. 外部脚本自动生成代码 #

如果滤波器类型较多，可使用简单的脚本（如Python、Perl）读取滤波器类型定义，通过模板生成对应的幅值响应函数代码。生成的代码为纯C89格式，可读性强，且无需在C代码中编写复杂宏。

示例流程：

1. 编写代码模板文件mag_response_template.c：

void compute_@FILTER@_mag_response(@FILTER@ *filter, int n_points, double *freq, double *mag) {
    double input[n_points];
    generate_sweep_input(input, n_points, SAMPLE_RATE);
    
    @PROCESS_FUNC@(filter, input, input, n_points);
    
    fftw_plan plan = fftw_create_plan(n_points, FFTW_FORWARD, FFTW_ESTIMATE);
    fftw_execute(plan, input, mag);
    compute_normalized_magnitude(mag, n_points);
    generate_freq_axis(freq, n_points, SAMPLE_RATE);
    
    fftw_destroy_plan(plan);
}

2. 用脚本遍历滤波器类型列表，替换模板中的@FILTER@和@PROCESS_FUNC@占位符，生成generated_mag_responses.c文件。
3. 将生成的文件纳入项目编译流程。

这种方式适合大规模扩展，且生成的代码完全符合C89标准，调试与维护成本低。

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