[实战]带限高斯白噪声的C语言实现（完整代码）

带限高斯白噪声的C语言实现（完整代码）

不依赖任何第三方库，用C语言生成带限高斯白噪声，并使用python进行数据频谱分析，确认C语言实现的正确性。
话不多说，实战开始。

---

C语言实现

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>

#define PI 3.14159265358979323846

// 生成高斯随机数
double gauss_rand(double sigma) {
   
   
    double u1 = rand() / (RAND_MAX + 1.0);
    double u2 = rand() / (RAND_MAX + 1.0);
    return sigma * sqrt(-2.0 * log(u1)) * cos(2 * PI * u2);
}

// 计算sinc函数
double sinc(double x) {
   
   
    if (fabs(x) < 1e-6) return 1.0;
    return sin(PI * x) / (PI * x);
}

int main() {
   
   
    // 参数设置
    const double Fs = 1000.0;       // 采样率 (Hz)
    const double BW = 100.0;        // 带宽 (Hz)
    const int N = 4096;             // 样本点数
    const double sigma = 1.0;       // 高斯分布标准差
    const char* filename = "noise_signal.bin";
    
    // 滤波器参数
    const int filter_length = 101;  // 滤波器长度 (奇数)
    const double cutoff = BW / Fs;  // 归一化截止频率

    srand(time(NULL));  // 初始化随机种子
    
    // 分配内存
    double *real = (double*)malloc