快速傅里叶变换载波频偏估计算法 matlab,快速傅里叶变换载波频偏估计算法

第11 卷 第2期 电路与系统学报 Vol.11, No.2 2006 年 4 月 JOURNAL OF CIRCUITS AND SYSTEMS April, 2006 文章编号：1007-0249 (2006) 02-0128-05 快速傅里叶变换载波频偏估计算法* 陈大夫， 张尔扬， 朱江 (国防科技大学 电子科学与工程学院，湖南 长沙 410073) 摘要：基于最大似然法参数估计，本文提出一种快速傅里叶变换(FFT)的频偏估计算法，它能直接检测得到载 波频偏绝对值大小，且精度与输入端信噪比无关。它可应用于中继卫星通信中高速数传且频带有效 QPSK 的同步解调。 关键词：卫星通信；QPSK；载波频偏；FFT 中图分类号：TN927 文献标识码：A 1 导言 载波恢复是数字通信调制解调器中一个关键组成部分。特别是在中继卫星通信高速数据传输系统中多卜勒频率参数使得接收机信号中存在高动态范围的载波频率偏差，它需要有一个频率跟踪控制环路辅助载波相位跟踪环路。为实现载波相位的快速准确捕获和低相位误差抖动性能的锁相环路跟踪，希望载波频偏越小越好，这就要求设计一个高分辨率的载波频偏估计器算法。例如对于频带有效的QPSK 调制，高速数据的宽带接收机中载波相位误差估计器除了采用 Costas 环外，也可采用 Viterbi 非线性相位误差估计算法[1]，当载波频偏远远小于符号速率时，它通过对 IQ 观测采样点的非线性坐标变 换和累加平均，可以在观测间隔内得到一个无偏的载波相位误差估计值；但载波频偏大小对相位跟踪环路的性能或者对载波相位估计值是很敏感的，它将会导致相位误差环路出现较大的相位偏差和相位抖动方差，从而使得系统信噪比实现损耗增大。解决问题的最有效方法就是预先利用载波频偏估计进 行检测与校正，实现载波频率的捕获与跟踪。 在以往的载波频率跟踪环路[2]或者载波频偏参数估计算法[3,4]中，交叉噪声引入频率跟踪环路，使得频偏误差控制精度与输入信噪比之间存在很强的依赖性，其性能受到限制。本文提出一种新的快速傅里叶变换(FFT)载波频偏估计器算法，它无需辅助数据(NDA)，能直接检测得到载波频偏绝对值大小。当输入端信噪比大于一个门限时，其频偏估计值的误差精度只与频偏检测器设计参数有关，而与输入端信噪比无关，并且非常适合于可编程数字电路 FPGA 实现。本文第 2 节基于 QPSK 信号和最大似然估计方法推导出载波频偏检测器 FFT 算法；第 3 节给出 FFT 频偏检测估计器设计与计算机仿真结果；第 4 节简述 FPGA 实现；最后第 5 节为结束语。 2 载波频偏估计 FFT 算法 多卜勒频率是空间信息源相对于中继卫星基站的径向运动速度所导致，而这种速度变化通常是慢变化的，因此在一段观测时间内这种载波频偏可以认为是一常数。对于 QPSK 调制信号，经过定时恢复后得到符号的峰值采样点。设载波频偏为 d F ，则 IQ 正交混频后输出信号可表示为： QntdFMQ I nt d F M I +++= +++= ) 0 2sin( ) 0 2cos( ϕπ φ ϕπ φ (1) 式中 QI nn 与 为 IQ 支路上加性的高斯噪声分量， M φ 、 0 ϕ 分别为 QPSK 调制信息相位和初始相位， M φ 具 体取值为 3,2,1,0,4/4/ =+= mm M ππ φ 。利用 Viterbio 非线性相位估计算法中坐标变换[1]，可以计算得 到 IQ 极坐标相位角，并且对其进行 4 倍放大处理，然后再对其进行正弦和余弦三角函数运算，即分 * 收稿日期：200