基于FPGA的信号处理算法之FFT相位差检测用Verilog实现（采用Altera芯片，Qua...

基于FPGA的信号处理算法，FFT法相差检测verilog实现 1.硬件平台：altera芯片 2.软件平台：Quartusii 13.1 Verilog 3.实现功能：检测两路正弦波的相位差

打开QuartusII工程就看到满屏的always块在跳动，这年头搞FPGA信号处理就像在钢丝绳上跳华尔兹。今天咱们要折腾的是两路正弦波相位差检测，用FFT暴力破解相位信息，Altera家的Cyclone IV芯片此刻正在实验台上瑟瑟发抖。

先甩个干货——FFT核的例化参数得调教到位。在Altera的FFT IP核配置界面，选个1024点浮点运算模式（虽然定点数更省资源，但相位精度要求必须上浮点）。关键参数得这么玩：

fft_core fft_inst ( .clk(clk_100M), .reset_n(!sys_rst), .inverse(1'b0), .sink_valid(adc_valid), .sink_sop(adc_sop), .sink_eop(adc_eop), .sink_real(adc_data_real), .sink_imag(12'd0), // 实信号输入时虚部接地 .source_ready(1'b1) );

注意第9行这个骚操作，处理实信号时直接把虚部接零，相当于白嫖了FFT核的实信号优化模式。但有个坑爹的地方：Altera的FFT核输出相位是弧度制，想转角度得乘以57.3，这时候千万别手贱用浮点乘法器，直接上查找表才是老司机的选择。

相位提取部分得搞个CORDIC核来算反正切，Verilog代码得这么写才不会被时序打脸：

// 相位计算流水线 always @(posedge clk) begin if(fft_valid) begin quad_judge <= {fft_imag[31], fft_real[31]}; //符号位捕获 abs_real <= fft_real[30:0]; abs_imag <= fft_imag[30:0]; end end cordic_atan cordic_inst( .clk(clk), .areset(sys_rst), .x_in(abs_real), .y_in(abs_imag), .phase_out(raw_phase) );

看到第2行的符号位拼接没？这波操作能准确判断相位所在象限，比直接算绝对值再判断象限省了三级流水线。CORDIC核输出的相位范围是-π到+π，处理两路信号的相位差时记得做循环补偿，否则当相位差超过180度时会突然反向跳变。

两路相位做差时要玩点环形缓冲区的把戏，这里给出核心差值计算代码：

// 相位差环形修正 always @(posedge clk) begin phase_diff_raw <= phase_a - phase_b; phase_diff <= (phase_diff_raw > 32'h40000000) ? phase_diff_raw - 32'h80000000 : (phase_diff_raw < 32'hC0000000) ? phase_diff_raw + 32'h80000000 : phase_diff_raw; end

这个三目运算符套娃看着头大？其实就干一件事：当原始相位差超过±π时，通过±2π调整使其落在[-π, +π]范围内。用32'h40000000表示π值是因为用了32位定点数，最高位符号位，后面31位用Q31格式表示。

最后在SignalTap里抓波形时，发现相位差结果偶尔抽风。排查发现是FFT帧同步信号没对齐，解决办法是在两路信号输入前加个同步FIFO，代码加这么几行：

sync_fifo fifo_dual( .aclr(sys_rst), .data({adc_a, adc_b}), .wrreq(adc_valid), .rdreq(fft_sink_ready), .q({sync_a, sync_b}) );

这个双通道同步FIFO确保了两路信号采样时刻严格对齐，否则测出来的相位差会有±1个采样周期的抖动。实测在50MHz采样率下，1024点FFT能做到0.35度的相位分辨率，足够检测电机控制中的转子位置偏差。

当把工程全编译完，发现用了8923个LE，36个18位乘法器，这资源消耗在Cyclone IV上还算能接受。不过要是想更省资源，可以把FFT点数砍到512，但相位分辨率会降到0.7度——做项目嘛，总是在性能和资源之间走钢丝。