基于Verilog的正弦波相位差计算算法设计问询

已知正弦波频率为1kHz，最大幅值为4.97V（采用IEEE 754浮点格式），波形其余采样值已转换为16位二进制，存储在包含60个地址的RAM中。现需要基于这些存储值，推导计算正弦波相位差的算法思路。当前已实现用于读取RAM值的查找表（LUT），Verilog代码如下：

// 存储器实现
module ram_memory_2(clk,add,data);

    input [5:0]add;
    input clk;
    output reg [15:0]data;
   

    always @ (*)
    begin
       case(add)
            0  : data <= 16'b0011011010010011;
            1  : data <= 16'b0011011011111010;
            2  : data <= 16'b0011011011111010;
            3  : data <= 16'b0011011101110001;
            4  : data <= 16'b0011011110101010;
            5  : data <= 16'b0011011111010111;
            6  : data <= 16'b0011011111100111;
            7  : data <= 16'b0011011111110100;
            8  : data <= 16'b0011011111011111;
            9  : data <= 16'b0011011111001111;
            10 : data <= 16'b0011011110011110;
            11 : data <= 16'b0011011101101101;
            12 : data <= 16'b0011011100110011;
            13 : data <= 16'b0011011011100001;
            14 : data <= 16'b0011011001111111;
            15 : data <= 16'b0011011000100001;
            16 : data <= 16'b0011010110100010;
            17 : data <= 16'b0011010100100011;
            18 : data <= 16'b0011010010011100;
            19 : data <= 16'b0011010000001100;
            20 : data <= 16'b0011001011100001;
            21 : data <= 16'b0011000110100010;
            22 : data <= 16'b0011000001010010;
            23 : data <= 16'b0010111000101011;
            24 : data <= 16'b0010101010101110;
            25 : data <= 16'b0010000110100010;
            26 : data <= 16'b1010100000001011;
            27 : data <= 16'b1010110010011000;
            28 : data <= 16'b1010111100111011;
            29 : data <= 16'b1011000011101110;
            30 : data <= 16'b1011001000111101;
            31 : data <= 16'b1011001101011100;
            32 : data <= 16'b1011010001001010;
            33 : data <= 16'b1011010011011101;
            34 : data <= 16'b1011010101100000;
            35 : data <= 16'b1011010111010011;
            36 : data <= 16'b1011011001000010;
            37 : data <= 16'b1011011010101000;
            38 : data <= 16'b1011011011111010;
            39 : data <= 16'b1011011101001100;
            40 : data <= 16'b1011011110001101;
            41 : data <= 16'b1011011110111010;
            42 : data <= 16'b1011011111010011;
            43 : data <= 16'b1011011111011011;
            44 : data <= 16'b1011011111101100;
            45 : data <= 16'b1011011111100011;
            46 : data <= 16'b1011011111000011;
            47 : data <= 16'b1011011110011110;
            48 : data <= 16'b1011011101011100;
            49 : data <= 16'b1011011100100011;
            50 : data <= 16'b1011011011001001;
            51 : data <= 16'b1011011001101010;
            52 : data <= 16'b1011011000000100;
            53 : data <= 16'b1011010110001101;
            54 : data <= 16'b1011010100001010;
            55 : data <= 16'b1011010001111011;
            56 : data <= 16'b1011001111010111;
            57 : data <= 16'b1011001010011000;
            58 : data <= 16'b1011000101011000;
            59 : data <= 16'b1011000000011001;
            default: data <= 16'b0000000000000000; // 地址未匹配时的默认输出
        endcase
    end
endmodule

相位差计算算法思路 #

前提说明 #

1kHz正弦波的周期为1ms，RAM中存储了60个采样点，因此采样间隔约为1ms / 60 ≈ 16.67μs，每个采样点对应的相位增量为360° / 60 = 6°（或π/30弧度）。

核心算法方案 #

方案1：采样点索引匹配法

• 假设需计算两个正弦波（A和B）的相位差，先分别读取两个波形的采样序列。
• 找到序列中相同幅值（或最接近幅值）且处于相同上升/下降沿的采样点，记录它们的地址索引idx_A和idx_B。
• 相位差计算公式：Δθ = |idx_A - idx_B| × 6°，若差值超过30个索引（对应180°），则取360° - Δθ作为最终相位差（保证结果在0~180°范围内）。
• 注意：需先将16位二进制值解码为实际电压值（根据IEEE 754格式或自定义量化规则），再进行幅值匹配。

方案2：反正弦查表还原相位法

• 对每个采样值解码为实际电压值V，计算归一化幅值V_norm = V / 4.97V。
• 通过反正弦函数计算单个采样点的相位：θ = arcsin(V_norm)，结合波形的上升/下降沿判断相位象限（比如当前采样点前一个值小于当前值，说明处于上升沿，相位在0_{180°；反之则在180°}360°）。
• 分别得到两个波形对应时刻的相位θ_A和θ_B，相位差为|θ_A - θ_B|（取0~180°范围内的最小值）。
• 硬件实现提示：可预先将反正弦函数的映射表存储在另一块RAM中，用归一化幅值的量化值作为地址查表得到相位值。

方案3：互相关运算定位法

• 对两个正弦波的采样序列做互相关运算，找到互相关值最大的位置偏移量k。
• 相位差计算公式：Δθ = k × 6°，同样根据偏移量大小取最小相位差。
• 优势：抗噪声能力强，适合存在干扰的场景；缺点是硬件实现复杂度较高，需要乘法和累加运算单元。

实现注意事项 #

• 必须明确16位二进制值的编码规则：是IEEE 754半精度浮点，还是定点量化（偏移二进制、补码等），这直接影响幅值解码的准确性。
• 若RAM中存储的是单个周期的正弦波，需确保两个待比较波形均基于该周期采样，否则需先做周期对齐。
• 硬件实现时优先选择索引匹配或查表法，避免复杂浮点运算，降低资源占用。

内容的提问来源于stack exchange，提问作者Vardan Chetnani