载波同步算法程序.doc

5.2.1 载波频偏的捕获 A．前导序列估计 由于发送端和接收端载波频率的不同，每一个采样信号在时间时包含一个未知的相位因素，这里指的是未知的载波频偏。这个未知的相位因素在接收端必须被估计和补偿，否则子载波的正交性将会被破坏。例如，当载波频率为5GHz时，那么100ppm的晶振偏移相对应的频率偏移为50kHz。若符号的间隔时期为。 前面对IEEE802.11a的分析，我们知道在它的帧结构中包含10个完全相同的短前导序列和两个相同的背靠背的长前导序列。其中短前导序列主要用于自动增益控制、分集选择、定时估计以及粗频率估计，而长前导序列主要用于信道估计和精确的频率估计。故结合这两个序列可以较精确的估计出载波的频偏，其中具体算法主要是利用它们良好的相关性[21]。首先设，则两个长前导字的相关值为： （39） 因此我们可以估计出，这里的指的是接收信号。 然而我们知道实际的值会比1大（如前面提到的100ppm的晶振偏移对应的为1.6），而长前导对其估计只能限制在内，故必须使用短前导字对其进行粗频率估计。短前导字的相关值为： （40） 故可以得到，可见短前导字的估计范围扩大到长前导字估计范围的4倍，也就是说精频偏估计的精度为粗频偏估计的4倍。结合上面提到的粗估计和精确估计，可以得到 （41） 这里指的是向下取整的意思。通过仿真结果表明此算法在理想状态下可估计的最大范围为（归一化值）。 B．载波频偏综合估计算法[12] 对于复杂的无线信道来说，上述算法的频偏估计范围偏小了一些。因此在文献[12]中提出了载波频偏估计算法，它的主要思想上在捕获阶段利用报头中长短训练符号的特殊结构，在频域和时域上分别实现整数载波频偏和分数载波频偏的捕获。 载波频差Δf 可以分为整数载波频差ΔfI 部分和分数载波频差ΔfF 部分，见下式 (42) (43) 整数载波频差ΔfI 是子载波频宽1/Tu 的整数倍，、和分别为ΔfI 、ΔfF和Δf相对于子载波频宽1/Tu的归一化值。频偏捕获阶段的任务是在尽可能短的时间内准确地建立子载波间的正交关系，即要求能够同时较快地估计出和。本文提出的捕获算法同样是基于IEEE802.11aPreamble 结构分两步，首先利用报头Preamble 中长训练符号已知和FFT 变换的频域移位特性在频域上捕获整数载波频偏 ，然后利用长短训练符号的结构特性在时域内捕获分数载波频偏。 整数频偏估计 设 和 分别表示接收端FFT变换之前preamble中第m个短（长）符号的第n个采样点，假设经过高斯白噪声信道，由于存在频偏，对于接收到的2个长训练符号中的任何一个都可以写成： (44) 假设小数部分频偏已经得到补偿，(5-12)式可以改写成： (45) 式中n(m,n)为高斯白噪声项。对上式两边同时做FFT变换，由于FFT变换频域移位特性可知，时域上的整数频率偏差经过FFT变换以后成为频域上的圆周移位。 (46) 其中是经过FFT变换解调子载波后在频域上的表示， 是发射端preamble中长训练符号在FFT调制子载波前在频域的表达式，表示长度为的圆周移位，是噪声项在频域的表达式。由于对接收端来说是已经信号，对任何一个经过FFT解调子载波的长训练符号，我们都用一个带偏移量i的本地长训练符号与之相关，并逐步改变i值，由(46)性质可知，最大相关时刻对应的i值就是所求的整数载波频偏。 (47) 文献指出此算法估计范围为(-32,31)(归一化值)，下面将进行仿真来验证，并将两者进行比较。 C．两种算法的仿真比较 在给定相同条件下，即相同的一帧数据、信道和频偏，改变信道和频偏参数对两种算进行比较。 通常设信噪比（SNR）为25～30，这里给定为25： 5.2.2 载波频偏的跟踪 前面对频偏捕获之后，我们认为残余的频偏足够小，可以使用跟踪模式进行补偿，而为了使得跟踪算法的可靠性，捕获后残余的频偏精度必须在0.5个子载波的空间范围内。到目前为止，可使用的频偏跟踪算法总共可以分为以下三种[36]：基于导频的算法（pilot tone-aided PTA），基于循环前缀的算法（cyclic prefix-based CPB）以及直接判决法(decision-directed DD)。PTA算法估计频偏主要通过在特定的子载波位置插入导频以及把接收到的符号和已知的导频符号相关，而对于IEEE80