发明名称 一种宽范围的码辅助载波同步实现方法
摘要 本发明涉及一种宽范围的码辅助载波同步实现方法,将码辅助载波同步分为粗同步和细同步两部分。首先以最大化评价函数为准则,利用二维搜索窗算法将载波偏移限制在一定范围内,然后再选择具有较大参数估计范围且复杂度较低的精估计算法来准确实现载波恢复。本发明的优势在于能够大幅提高同步范围,具有较高的估计精度,可以有效实现载波同步,且获得接近理想同步的译码性能。
申请公布号 CN103298101A 申请公布日期 2013.09.11
申请号 CN201310245034.2 申请日期 2013.06.19
申请人 中国人民解放军空军工程大学 发明人 许华;沈海鸥;李晓东
分类号 H04W56/00(2009.01)I;H04L1/00(2006.01)I 主分类号 H04W56/00(2009.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种宽范围的码辅助载波同步实现方法,其特征在于步骤如下:步骤1:将接收信号序列r<sub>k</sub>可偏移的最大相位、频率值作为相位-频率二维搜索空间的边界,在此空间内等间隔的选择n个离散的相位-频率点(θ<sub>i</sub>,Δf<sub>j</sub>),并分别利用下式对接收信号进行载波偏差的补偿:<maths num="0001"><![CDATA[<math><mrow><msup><msub><mi>r</mi><mi>k</mi></msub><mo>&prime;</mo></msup><mo>=</mo><msub><mi>r</mi><mi>k</mi></msub><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>+</mo><mn>2</mn><mi>&pi;k</mi><msub><mi>&Delta;f</mi><mi>j</mi></msub><mi>T</mi><mo>)</mo></mrow></mrow></msup></mrow></math>]]></maths>式中,k=0,1,2…K-1,K是信道编码的码元长度,T是码元间隔周期;对n个补偿后的信号r<sub>k</sub>′依次进行解调、译码;步骤2:将各组译码输出的K个软信息L(c<sub>k</sub>)取绝对值后进行门限判决,选择大于门限值的|L(c<sub>k</sub>)|进行累加并求平均得<img file="FDA0000336735320000011.GIF" wi="150" he="68" />然后比较n个<img file="FDA0000336735320000012.GIF" wi="121" he="77" />值的大小,选择<img file="FDA0000336735320000013.GIF" wi="126" he="76" />最大时对应的那个相位-频率点作为粗同步的估计值<img file="FDA0000336735320000014.GIF" wi="201" he="85" />而其对应的已补偿信号r<sub>k</sub>′即为载波细同步的输入信号y<sub>k</sub>;步骤3:对载波细同步进行初始化,令迭代次数i=0,计数器<img file="FDA0000336735320000015.GIF" wi="551" he="84" />设细同步的最大迭代次数为I;步骤4:利用<img file="FDA0000336735320000016.GIF" wi="79" he="85" />和<img file="FDA0000336735320000017.GIF" wi="101" he="85" />更新载波细同步的输入信号向量:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>y</mi><mi>k</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msub><mi>y</mi><mi>k</mi></msub><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><msubsup><mover><mi>&theta;</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>+</mo><mn>2</mn><mi>&pi;k&Delta;</mi><msubsup><mover><mi>f</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mi>T</mi><mo>)</mo></mrow></mrow></msup></mrow></math>]]></maths>对信号<img file="FDA0000336735320000019.GIF" wi="74" he="78" />进行解调、译码得输出软信息L<sub>i</sub>(c<sub>k</sub>),然后采用与步骤2同样的方法得<img file="FDA00003367353200000110.GIF" wi="166" he="67" />步骤5:当i&gt;1时,计算第i次和第i-1次迭代所得<img file="FDA00003367353200000111.GIF" wi="133" he="77" />的差值,如果<img file="FDA00003367353200000112.GIF" wi="475" he="77" />计数器count加1,否则计数器清零;若i=1则直接跳过此步,直接进行步骤6;步骤6:利用步骤4中大于门限值的|L<sub>i</sub>(c<sub>k</sub>)|所对应的L<sub>i</sub>(c<sub>k</sub>)通过下式求调制星座图上所有可能符号对后验概率的均值<img file="FDA00003367353200000113.GIF" wi="201" he="71" /><img file="FDA00003367353200000114.GIF" wi="979" he="141" />式中,M是星座点的个数,a<sub>k</sub>是接收到的调制码元序列,A<sub>m</sub>是调制星座图上第m个点的值,tanh是双曲正切运算;步骤7:利用<img file="FDA0000336735320000021.GIF" wi="206" he="78" />进行频偏估计,利用简化的EM估计方法<img file="FDA0000336735320000022.GIF" wi="886" he="142" />通过最大值搜索即可得到频率偏移的估计值,搜索区间利用EM算法的恒收敛特性根据前两次的迭代结果自适应地减小下一次频率搜索区间的大小,第i+2次迭代的参数的搜索区间为:<maths num="0003"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mover><mi>f</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>2</mn><mo>)</mo></mrow></msubsup><mo>&Element;</mo><mrow><mo>(</mo><msubsup><mrow><mi>&Delta;</mi><mover><mi>f</mi><mo>^</mo></mover></mrow><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>-</mo><msubsup><mi>f</mi><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><mi>&Delta;</mi><msubsup><mover><mi>f</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>+</mo><msubsup><mi>f</mi><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>其中,<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>f</mi><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>=</mo><mo>|</mo><mi>&Delta;</mi><msubsup><mover><mi>f</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>-</mo><mi>&Delta;</mi><msubsup><mover><mi>f</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>|</mo><mo>;</mo></mrow></math>]]></maths>步骤8:利用<img file="FDA0000336735320000025.GIF" wi="204" he="78" />进行相偏估计,<img file="FDA0000336735320000026.GIF" wi="600" he="142" />且令n=n+1;步骤9:重复上述步骤(4)~(8),直到count=P或i=I时,得到载波偏差精估计值<img file="FDA0000336735320000027.GIF" wi="58" he="84" />和<img file="FDA0000336735320000028.GIF" wi="108" he="85" />退出迭代循环,执行下一步;步骤10:利用最终得到的载波偏差精估计值<img file="FDA0000336735320000029.GIF" wi="50" he="87" />和<img file="FDA00003367353200000210.GIF" wi="78" he="84" />对细同步输入信号y<sub>k</sub>进行补偿得z<sub>k</sub>:<img file="FDA00003367353200000211.GIF" wi="446" he="86" />再经过解调、译码判决之后即可恢复出发送端的原始信息。
地址 710038 陕西省西安市霸桥区长乐东路甲字1号