基于中国余数定理的欠采样下余弦信号的频率估计-物理学报.pdfVIP

物理学报 Acta Phys. Sin. Vol. 63, No. 19 (2014) 198403 基于中国余数定理的欠采样下余弦 信号的频率估计 黄翔东 丁道贤 南楠 王兆华 (天津大学电子信息工程学院, 天津 300072) ( 2014 年4 月25 日收到; 2014 年5 月22 日收到修改稿) 基于中国余数定理的重构算法的信号频率估计是近年来信号处理、电磁学以及光学等领域的前沿问题, 但目前这些研究仅限于对复指数信号做粗略频率估计. 因而, 本文把基于中国余数定理的频率估计从复指数 信号粗估计拓展到实余弦信号精细估计领域, 其所提出的估计方案处理过程如下: 1) 对高频余弦波形进行过 零点检测, 确定信号的相位信息; 2) 对各路欠采样信号做快速傅里叶变换, 并借助Candan 估计器对各路谱峰 值做频率校正以获取高精度余数估计, 基于此算出频偏值以做相位校正; 3) 用提出的基于相位特征分类方法 对校正得到的余数做筛选; 4) 将筛选出的频率余数代入闭合形式的中国余数定理得到原信号频率的高精度估 计. 此外, 本文还推导出了频率估计方差的理论表达式. 数据模拟实验验证了该表达式的正确性, 实验结果还 反映了本文提出的方案具有高精度和高抗噪性能. 关键词: 中国余数定理, 欠采样, 频率估计, 余数筛选 PACS: 84.40.Ua, 43.58.Kr, 02.10.De, 43.60.Dh DOI: 10.7498/aps.63.198403 足互素关系的整数模值: M , , M , 及相应的 余数r (即满足r N modM , i , , L, ‘mod’ 1 引 言 表示模除操作) 来估计N . 频率估计19 在信号处理, 电磁学以及光学等 CRT 重构算法主要分为两类: 余数数目冗 领域一直是关键问题, 并因此产生了多种频率估计 余方法11213 和余数冗余方法561416 . 其中后 算法. 然而欠采样(即采样速率fs 远小于2 倍信号 者非常适用于欠采样速率估计, 特别是文献[16] 频率f ) 下的高精度频率估计一直是个难题, 因为 提出新的余数冗余方法—– 闭合形式的CRT 重构 奈奎斯特定理要求一个信号周期内采到2 个以上样 法是目前实现欠采样速率估计的最佳CRT 算法, 点, 随着信号频率升高, 必然要求其采样速率相应 这是因为该方法放宽了经典CRT 要求L 个模值 升高, 这就会对模数转换器(analog to digital con- M —M 互素的限制, 允许其最大公约数(greatest verter, ADC) 的转换速率、功耗以及硬件成本提出 common divider, gcd)M =gcd{M , , M } 为不 更高要求, 在某些特定场合(如fs samples/s 等于1 的正整数, 且文献[16] 将CRT 重构对象从整 的