发明名称 基于导频联合编码辅助的SOQPSK载波同步方法
摘要 本发明公开一种基于导频联合编码辅助的SOQPSK载波估计方法,主要解决现有技术估计复杂度高和估计范围小的问题。其实现步骤是:(1)确定粗估计参数和细估计参数;(2)利用粗估计参数进行粗估计,得到频偏粗估值f<sub>co</sub>和相偏粗估值φ<sub>co</sub>;(3)利用细估计参数进行细估计得到频偏细估值f<sub>in</sub>和相偏细估值φ<sub>in</sub>;(4)由频偏粗估计值f<sub>co</sub>和频偏细估值f<sub>in</sub>,得到最终的频偏估计值f=f<sub>co</sub>+f<sub>in</sub>;(5)由相偏粗估值φ<sub>co</sub>和相偏细估值φ<sub>in</sub>,得到最终的相偏估计值φ=φ<sub>co</sub>+φ<sub>in</sub>。本发明可将估计的归一化频偏范围提高到-0.5~0.5,且估计复杂度小,可用于无线通信系统的载波同步和无线信道估计。
申请公布号 CN103916357A 申请公布日期 2014.07.09
申请号 CN201410177855.1 申请日期 2014.04.29
申请人 西安电子科技大学 发明人 孙锦华;韩会梅;余忠洋;丁玉杰;王雪梅;刘鹏
分类号 H04L27/26(2006.01)I;H04L25/02(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于导频联合编码辅助的SOQPSK载波同步方法,包括:(1)细估计参数确定步骤:(1a)产生数据长度为N<sub>0</sub>=N×t的二进制比特序列Y,其中N为正整数,t为一个Turbo乘积码TPC块编码前的数据长度;(1b)利用二进制比特序列Y,通过Turbo乘积码TPC编码,得到N个Turbo乘积码TPC块C<sub>N</sub>;(1c)利用N个Turbo乘积码TPC块C<sub>N</sub>,通过成形偏移正交相移键控SOQPSK调制,得到调制后的信号S<sub>N</sub>;(1d)根据不同频偏对调制后的信号S<sub>N</sub>的影响,确定满足0.2dB信噪比损失条件下的归一化频率偏移值Δf<sub>1</sub>T<sub>s</sub>和满足0.6dB信噪比损失条件下的归一化频率偏移值Δf<sub>2</sub>T<sub>s</sub>,其中,Δf<sub>1</sub>、Δf<sub>2</sub>指频率偏移值,T<sub>s</sub>指码元周期;(1e)利用调制后的信号S<sub>N</sub>,通过细估计算法,得到N<sub>1</sub>个归一化剩余频偏<img file="FDA0000498677330000011.GIF" wi="92" he="78" />(1f)根据归一化剩余频偏<img file="FDA0000498677330000012.GIF" wi="62" he="77" />及归一化频率偏移值Δf<sub>1</sub>T<sub>s</sub>,确定细估计所需的Turbo乘积码TPC块数B<sub>in</sub>;(2)粗估计参数确定步骤:(2a)产生数据长度都为L<sub>P</sub>的前段导频序列R1和后段导频序列R2,其中L<sub>P</sub>为200;(2b)利用前段导频序列R1、后段导频序列R2和N个Turbo乘积码TPC块C<sub>N</sub>,通过时域相关算法,得到N<sub>2</sub>个归一化剩余频偏<img file="FDA0000498677330000013.GIF" wi="103" he="70" />(2c)根据归一化剩余频偏<img file="FDA0000498677330000014.GIF" wi="66" he="74" />及步骤(1d)中确定的频率偏移值Δf<sub>2</sub>T<sub>s</sub>,确定粗估计所需的Turbo乘积码TPC块数B<sub>co</sub>;(2d)利用粗估计所需的Turbo乘积码TPC块数B<sub>co</sub>,得到<img file="FDA0000498677330000015.GIF" wi="90" he="76" />个最大归一化频偏<img file="FDA0000498677330000016.GIF" wi="582" he="86" />其中L<sub>C</sub>为一个Turbo乘积码TPC块编码后的码元长度,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>L</mi><mi>P</mi><mo>&prime;</mo></msubsup><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>L</mi><mi>P</mi><mi>max</mi></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000498677330000017.GIF" wi="351" he="79" /></maths>(2e)产生长度为L′<sub>P</sub>的导频信号集<img file="FDA0000498677330000018.GIF" wi="98" he="75" />利用导频信号集<img file="FDA0000498677330000019.GIF" wi="74" he="78" />和最大归一化频偏<img file="FDA00004986773300000110.GIF" wi="173" he="82" />通过旋转平均周期图RPA算法,得到<img file="FDA00004986773300000111.GIF" wi="86" he="74" />个标志信息<img file="FDA00004986773300000112.GIF" wi="95" he="84" />(2f)利用<img file="FDA00004986773300000113.GIF" wi="94" he="77" />个标志信息<img file="FDA00004986773300000114.GIF" wi="100" he="84" />得到粗估计所需的导频符号长度<img file="FDA00004986773300000115.GIF" wi="83" he="82" />(3)将数据帧中的数据结构设置为前段导频、第一段数据、后段导频、第二段数据,并根据步骤(1b)、(2c)、(2f)得到的结果,设前段导频和后段导频的符号长度均为<img file="FDA0000498677330000021.GIF" wi="102" he="79" />第一段数据的符号长度为d<sub>1</sub>=B<sub>co</sub>×L<sub>C</sub>、第二段数据的符号长度为d<sub>2</sub>=(N‑B<sub>co</sub>)×L<sub>C</sub>;(4)粗估计步骤(4a)在接收端,将接收到的基带信号Z,送到载波估计器,取出基带信号Z前<img file="FDA0000498677330000022.GIF" wi="60" he="80" />个数据,作为基带信号Z的前段导频P1;(4b)利用前段导频P1,通过旋转平均周期图RPA算法,得到频偏估计值f<sub>R</sub>;(4c)利用频偏估计值f<sub>R</sub>,通过复共轭相乘的方法,对接收信号Z进行校正,得到频偏校正信号Z′;(4d)取出频偏校正信号Z′的前<img file="FDA0000498677330000023.GIF" wi="60" he="80" />个数据作为Z′的前段导频P2,取出Z′<sub>k</sub>的第<img file="FDA0000498677330000024.GIF" wi="220" he="83" />个符号到第<img file="FDA0000498677330000025.GIF" wi="170" he="85" />个符号的数据作为Z′的后段导频P3;(4e)利用前段导频P2和后段导频P3,分别得到前段导频去调制信号P2<sup>d</sup>=P2×(P2')<sup>*</sup>和后段导频去调制信号P3<sup>d</sup>=P3×(P3')<sup>*</sup>,其中P2'和P3'为收发双方事先约定的发送信号的前段导频和后段导频,(P2')<sup>*</sup>为P2'的共轭,(P3')<sup>*</sup>为P′<sub>3</sub>的共轭;(4f)利用前段导频去调制信号P2<sup>d</sup>和后段导频去调制信号P3<sup>d</sup>,通过时域相关算法,得到频偏估计值f<sub>C</sub>;(4g)将频偏估计值f<sub>R</sub>与频偏估计值f<sub>C</sub>相加,得到频偏粗估值f<sub>co</sub>=f<sub>R</sub>+f<sub>C</sub>;(4h)利用前段导频去调制信号P2<sup>d</sup>,通过最大似然算法,得到相位粗估值φ<sub>co</sub>;(5)细估计步骤(5a)利用频偏粗估值f<sub>co</sub>和相偏粗估值φ<sub>co</sub>,通过复共轭相乘的方法,对接收信号Z进行校正,得到频偏校正信号Z″;(5b)利用频偏校正信号Z″,通过成形偏移正交相移键控SOQPSK解调,得到解调后的信号Z″′;(5c)取出解调后的信号Z″′的前<img file="FDA0000498677330000026.GIF" wi="54" he="82" />个数据,作为Z″′的前段导频P4,取出Z″′的第<img file="FDA0000498677330000027.GIF" wi="112" he="79" />个符号到第<img file="FDA0000498677330000028.GIF" wi="148" he="80" />个符号的数据作为Z″′的第一段数据Q<sub>1</sub>,取出Z″′的第<img file="FDA0000498677330000029.GIF" wi="216" he="84" />个符号到第<img file="FDA00004986773300000210.GIF" wi="172" he="76" />个符号的数据作为Z″′的后段导频P5,取出Z″′的第<img file="FDA00004986773300000211.GIF" wi="235" he="80" />个符号到第<img file="FDA00004986773300000212.GIF" wi="260" he="77" />个符号的数据作为Z″′的第二段数据Q<sub>2</sub>;(5d)利用第一段数据Q<sub>1</sub>和第二段数据Q<sub>2</sub>,通过TPC译码,得到译码后的二进制比特信息A;(5e)利用二进制比特信息A,通过Turbo乘积码TPC编码,得到编码后的N个TPC块;(5f)将编码后的N个Turbo乘积码TPC块、前段导频P4、后段导频P5,通过复用模块,得到复用信号S1,即将前B<sub>co</sub>个Turbo乘积码TPC块插入到前段导频P4和后段导频P5之间,将剩下N‑B<sub>co</sub>个Turbo乘积码TPC块插入到后段导频P5后面;(5g)利用复用信号S1,通过成形偏移正交相移键控SOQPSK调制,得到调制后的信号S2;(5h)利用调制后的信号S2及频偏校正信号Z″,得到频偏校正信号Z″的去调制信号S3=Z″×(S2)<sup>*</sup>,其中(S2)<sup>*</sup>为S2的共轭;(5i)利用去调制信号S3,通过时域相关算法,得到频偏细估值f<sub>in</sub>;(5j)利用去调制信号S3及频偏细估值f<sub>in</sub>,通过最大似然算法,得到相偏细估值φ<sub>in</sub>;(6)将频偏粗估值f<sub>co</sub>与频偏细估值f<sub>in</sub>相加,得到最终的频偏估计值:f=f<sub>co</sub>+f<sub>in</sub>;(7)将相偏粗估值φ<sub>co</sub>与相偏细估值φ<sub>in</sub>相加,得到最终的相偏估计值:φ=φ<sub>co</sub>+φ<sub>in</sub>。
地址 710071 陕西省西安市太白南路2号