数字信号处理课程仿真辅助课v2 第六课 频率估计(3)

上一篇文章讲了三种估计算法，文末希望大家能进行性能比对，不知大家有没有做啊？相位差分估计算法中加入了补0措施，大家可以思考一下是否有必要？通过仿真可以看到估计结果会发生哪些变化？

如果你没有时间做，那这个仿真就由我来做吧。本文把这三个算法针对同一信号源进行测试，这样就可以比对估计结果，方便大家在实际使用时进行性能和复杂度之间的折中考虑。同时要进行蒙特卡洛仿真，看看估计值的方差情况。另外再给出算法在实际使用中的融入过程，比如在使用了PMF-FFT算法捕获信号后，是否能用来提升估计精度。本文又是一篇长文，仿真的过程就是这么繁琐，不然得不到算法的真实效果。

%% 步骤2：频率估计

%%%% 差分相位

delta_n = N/2;  % 第二段相对第一段的时移（可调整）

% 分段信号

x1 = x(1:end-delta_n);

x2 = x(1+delta_n:end);

% 计算FFT