温馨提示×

c语言傅里叶变换程序怎么实现

小亿
106
2024-08-05 18:41:13
栏目: 编程语言

实现傅里叶变换的C语言程序可以使用库函数或者自己编写代码来实现。以下是一个使用库函数的示例代码:

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

int main() {
    int N = 8; // 采样点数
    double *in = (double*) fftw_malloc(sizeof(double) * N);
    fftw_complex *out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * N);

    fftw_plan p;
    p = fftw_plan_dft_r2c_1d(N, in, out, FFTW_ESTIMATE);

    // 初始化输入信号
    for (int i = 0; i < N; i++) {
        in[i] = sin(2 * M_PI * i / N);
    }

    // 执行傅里叶变换
    fftw_execute(p);

    // 输出变换结果
    for (int i = 0; i < N; i++) {
        printf("out[%d] = %f + %f*I\n", i, out[i][0], out[i][1]);
    }

    fftw_destroy_plan(p);
    fftw_free(in);
    fftw_free(out);

    return 0;
}

在上面的示例代码中,首先使用fftw库函数进行了傅里叶变换的初始化,然后生成了一个正弦波信号作为输入信号,在执行傅里叶变换后输出了变换结果。通过这种方式可以很方便地实现傅里叶变换的计算。

0