oneChirp_fft_mat_512=fft(oneChirp_MCU,512); figure() plot(abs(oneChirp_fft_mat_512(1:nSample/2)));title('512fft') oneChirp_fft_mat_512_int = round(real(oneChirp_fft_mat_512));用c语言实现此matlab程序

以下是用 C 语言实现的代码：

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

#define nSample 512

float oneChirp_MCU[nSample];
float oneChirp_fft_mat_512[nSample];

int main()
{
    // 初始化 oneChirp_MCU 数组
    for (int i = 0; i < nSample; i++) {
        oneChirp_MCU[i] = sin(i * M_PI / 16);
    }

    // 计算 512 点 FFT
    for (int k = 0; k < nSample; k++) {
        float re = 0;
        float im = 0;
        for (int n = 0; n < nSample; n++) {
            float cos_val = cos(2 * M_PI * k * n / nSample);
            float sin_val = sin(2 * M_PI * k * n / nSample);
            re += oneChirp_MCU[n] * cos_val;
            im += oneChirp_MCU[n] * sin_val;
        }
        oneChirp_fft_mat_512[k] = sqrt(re * re + im * im);
    }

    // 输出结果
    for (int i = 0; i < nSample / 2; i++) {
        printf("%f\n", oneChirp_fft_mat_512[i]);
    }

    return 0;
}

需要注意的是，MATLAB 中的 `fft` 函数计算的是单边频谱，而 C 语言实现的代码计算的是双边频谱，因此在输出结果时只输出了前一半的结果。另外，MATLAB 中的 `round` 函数在 C 语言中可以使用强制类型转换实现，即将 `float` 转换为 `int`，这里省略了这一步转换。