GnuRadio 3.8中Costas Loop模块无法在理想信号上正确锁相

我正在使用GnuRadio 3.8.x版本测试Costas Loop模块，输入为理想无移相的1kHz正弦波，预期输出相位稳定在0，但实际输出相位波动于非零值（如带宽0.05时为2.35619，带宽0.1时为0.785398）。

测试设置：

• 信号源：1kHz正弦波（相位设为0）
• Costas环带宽：0.05
• Costas环阶数：4
• 测试目标：理想信号下的相位锁定
• 预期输出：相位稳定在0

测试代码：

int main() {
auto tb = gr::make_top_block("costas_loop_test");

double samp_rate = 32000;
double freq = 1000;
float loop_bw = 0.05;     // Costas loop bandwidth 
int order = 4;           

// Signal source: complex sine wave
auto src = gr::analog::sig_source_c::make(samp_rate, 
                                          gr::analog::GR_SIN_WAVE, 
                                          freq, 
                                          1.0, 
                                          0.0);  // Phase set to 0


auto costas = gr::digital::costas_loop_cc::make(loop_bw, order); 


auto mag_block = gr::blocks::complex_to_mag::make();
auto phase_block = gr::blocks::complex_to_arg::make();

auto mag_probe = gr::blocks::probe_signal_f::make();
auto phase_probe = gr::blocks::probe_signal_f::make();

auto sink = gr::blocks::null_sink::make(sizeof(float));

// Connect the flowgraph
tb->connect(src, 0, costas, 0);               // Source -> Costas Loop
tb->connect(costas, 0, mag_block, 0);         // Costas -> Magnitude Block
tb->connect(costas, 0, phase_block, 0);       // Costas -> Phase Block
tb->connect(mag_block, 0, mag_probe, 0);      // Magnitude -> Probe
tb->connect(phase_block, 0, phase_probe, 0);  // Phase -> Probe
tb->connect(mag_block, 0, sink, 0);           // Magnitude -> Sink (output)


tb->start();
std::cout << "Running Costas loop with perfect sine wave input...\n";


for (int i = 0; i < 100; ++i) {
    float amplitude = mag_probe->level();  
    float phase = phase_probe->level(); // Phase (should stabilize near 0 for perfect signal)

    std::cout << "Amplitude: " << amplitude << " | Phase: " << phase << std::endl;
    std::this_thread::sleep_for(std::chrono::milliseconds(100));

    // Check for phase lock (near zero)
    if (std::abs(phase) < 0.01) {
        std::cout << "Phase detected at sample " << i << std::endl;
        
    }
}

// Stop the flowgraph
tb->stop();
tb->wait();
return 0;
}

疑问：是GnuRadio中Costas Loop的实现存在bug，还是我对相位锁定的理解有误？为何理想输入信号下输出相位无法稳定在0？

这不是GnuRadio的实现bug，而是你对高阶Costas环的相位锁定特性理解有误。

对于4阶Costas环，它的锁定相位存在固有模糊性——环会稳定在输入信号相位加上π/2的整数倍的位置，也就是0、π/2（≈0.785）、π（≈3.1416）、3π/2（≈2.356）这四个相位中的任意一个，具体锁定到哪个取决于环路的初始状态和瞬态过程（哪怕是理想信号，环路启动时的微小波动也会触发选择其中一个稳定点）。

你观察到的数值正好对应这些稳定相位：

• 带宽0.05时的2.35619就是3π/2
• 带宽0.1时的0.785398就是π/2

这是4阶Costas环鉴相器的特性决定的：它的误差信号在上述四个相位点处都会归零，因此环路会锁定在这些点上，而非唯一的0相位。如果需要锁定到特定相位，你需要在Costas环之后添加额外的相位校准机制，或者使用差分编码来消除相位模糊。

内容的提问来源于stack exchange，提问作者James S