matlab fir1 C实现

Matlab 的 `fir1` 函数用于设计一线性相位有限冲积（FIR）低通滤波器，它在 MATLAB 中是一种方便的工具。然而，如果你想要在 C 语言环境下实现类似的功能，你需要手动计算滤波器系数。C 语言没有内置的函数直接对应于 `fir1`，但你可以参考以下步骤来完成：

1. 确定滤波器参数，如截止频率、采样率、滤波器阶数等。
2. 使用滤波器设计公式，例如 Parks-McClellan 或窗函数法（如汉明窗、海明窗等）来计算滤波器系数。
3. 编程实现一维数组，存储滤波器系数。

下面是一个简单的示例代码片段，展示如何使用窗口函数法在 C 语言中实现 FIR 滤波器：

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

// 定义滤波器长度和窗口类型
#define FILTER_ORDER 10
#define WINDOW_HAMMING

double window[FILTER_ORDER + 1]; // 存储滤波器系数
double fs; // 采样率
double cutoff_freq; // 截止频率

void fir_design(double fs, double cutoff_freq) {
    double alpha = cutoff_freq / (fs / 2);
    if (WINDOW_HAMMING) {
        for (int i = 0; i <= FILTER_ORDER; ++i) {
            window[i] = 0.54 - 0.46 * cos(2 * M_PI * i / FILTER_ORDER);
        }
    } else { // 也可以选择其他窗口函数
        // ... 其他窗口函数实现 ...
    }
}

// 更多功能实现
// ...

int main() {
    fs = 44100; // 示例采样率
    cutoff_freq = 1000; // 示例截止频率
    fir_design(fs, cutoff_freq);

    // ... 使用window数组计算滤波...

    return 0;
}