FPGA 快速傅里叶变换（FFT）IP 核配置与实现

2. 数据准备

关键在于起始符（SOP）和结束符（EOP）的时序控制。SOP 上升沿标志数据开始，EOP 上升沿标志结束。建议在仿真环境中验证时序，避免标志符错位导致数据错乱。

3. 数据处理

获取输出数据的实部、虚部及指数后，首先将补码转换为原码，再进行能量谱转换。注意幅频特性线的对称性，N 点 FFT 只需取前 N/2 部分频谱，后半部分与前半部分对称。

4. 代码实现

以下为中间截取的关键逻辑代码，展示了 1024 点 FFT 的峰值搜索与幅值计算流程。

// 1024 points FFT and search peak
// Use 1024kHz sampling frequency, for the resolution frequency will be about 1kHz
// The two signal frequencies in the FFT will be saved in Var fft_out[0] and fft_out[1]

wire signed [11:0] Data_signed;
assign Data_signed = Din[11] ? Din - 12'h800 : Din + 12'h800; // 偏移至有符号范围（-2048~2047）

reg [9:0] sample_count; // 0~1023
always @(posedge clk1024k) begin
    if (sample_count == 10'd1023)
        sample_count <= 0;
    else
        sample_count <= sample_count + 1;
end

// Gain FFT frame
wire sink_sop;
wire sink_eop;
assign sink_sop = (sample_count == 10'd0);
assign sink_eop = (sample_count == 10'd1023);

wire [11:0] real_out;
wire [11:0] imag_out;
wire source_valid;
wire freq_start;
wire freq_end;
wire [5:0] exp;

fft1024 u_fft (
    .clk(clk1024k),
    .reset_n(1),
    .inverse(0),
    .sink_valid(1),
    .sink_sop(sink_sop),
    .sink_eop(sink_eop),
    .sink_real(Data_signed),
    .sink_imag(12'b0),
    .sink_error(2'b00),
    .source_ready(1),
    .source_valid(fft_valid),
    .source_real(real_out),
    .source_imag(imag_out),
    .source_error(),
    .source_sop(freq_start),
    .source_eop(freq_end),
    .source_exp(exp)
);

// Gain magnitude
wire [11:0] real_num;
wire [11:0] imag_num;
wire [23:0] power;
assign real_num = real_out[11] ? (~real_out + 1) : real_out;
assign imag_num = imag_out[11] ? (~imag_out + 1) : imag_out;
assign power = (real_num * real_num) + (imag_num * imag_num);

wire [11:0] magnitude_temp;
wire [17:0] magnitude;
sqrt u_sqrt(
    .radical(power),
    .q(magnitude_temp)
);
assign magnitude = magnitude_temp << exp;

// Get spectrum
reg [9:0] cnt;
reg freq_end_d;
always @(posedge clk1024k) begin
    if (freq_start)
        cnt <= 0;
    else if (fft_valid)
        cnt <= cnt + 1;
    
    if (fft_valid && cnt <= 10'd512) begin
        // magnitude 即为当前频率下的幅频特性
    end
end

四、注意事项

在实际调试过程中，常见的问题集中在数据输入模式的选择与时序控制上。

1. 数据流模式：如果不熟悉内存管理，强烈建议使用 Streaming 模式。其他模式如 Burst 可能涉及复杂的内存读写逻辑，容易导致输出错乱。
2. 标志符时序：使用 Streaming 模式时，务必仔细推演 SOP 和 EOP 的时序。标志符错位会导致整个帧的数据解析错误，建议配合仿真工具验证。

总结

对于满足工程需求的场景，推荐使用系统提供的 IP 核进行配置，效率更高且稳定性好。虽然基于原理手写代码（如 8 点 FFT）有助于理解底层机制，但资源消耗较大，难以直接应用于复杂系统。若需深入练习算法原理，手写实现不失为一种好的方法，但在实际工程中应优先考虑成熟 IP 核方案。

以上方案基于实际调试经验整理，难免存在细节疏漏，欢迎指正讨论。