本发明专利技术涉及一种宽范围的码辅助载波同步实现方法，将码辅助载波同步分为粗同步和细同步两部分。首先以最大化评价函数为准则，利用二维搜索窗算法将载波偏移限制在一定范围内，然后再选择具有较大参数估计范围且复杂度较低的精估计算法来准确实现载波恢复。本发明专利技术的优势在于能够大幅提高同步范围，具有较高的估计精度，可以有效实现载波同步，且获得接近理想同步的译码性能。
【技术实现步骤摘要】
【技术保护点】
一种宽范围的码辅助载波同步实现方法，其特征在于步骤如下：步骤1：将接收信号序列rk可偏移的最大相位、频率值作为相位?频率二维搜索空间的边界，在此空间内等间隔的选择n个离散的相位?频率点(θi,Δfj)，并分别利用下式对接收信号进行载波偏差的补偿：rk′=rk·e-j(θi+2πkΔfjT)式中，k=0,1,2…K?1，K是信道编码的码元长度，T是码元间隔周期；对n个补偿后的信号rk′依次进行解调、译码；步骤2：将各组译码输出的K个软信息L(ck)取绝对值后进行门限判决，选择大于门限值的|L(ck)|进行累加并求平均得然后比较n个值的大小，选择最大时对应的那个相位?频率点作为粗同步的估计值而其对应的已补偿信号rk′即为载波细同步的输入信号yk；步骤3：对载波细同步进行初始化，令迭代次数i=0，计数器设细同步的最大迭代次数为I；步骤4：利用和更新载波细同步的输入信号向量：yk(i)=yk·e-j(θ^x(i)+2πkΔf^x(i)T)对信号进行解调、译码得输出软信息Li(ck)，然后采用与步骤2同样的方法得步骤5：当i>1时，计算第i次和第i?1次迭代所得的差值，如果计数器count加1，否则计数器清零；若i=1则直接跳过此步，直接进行步骤6；步骤6：利用步骤4中大于门限值的|Li(ck)|所对应的Li(ck)通过下式求调制星座图上所有可能符号对后验概率的均值式中，M是星座点的个数，ak是接收到的调制码元序列，Am是调制星座图上第m个点的值，tanh是双曲正切运算；步骤7：利用进行频偏估计，利用简化的EM估计方法通过最大值搜索即可得到频率偏移的估计值，搜索区间利用EM算法的恒收敛特性根据前两次的迭代结果自适应地减小下一次频率搜索区间的大小，第i+2次迭代的参数的搜索区间为：Δf^x(i+2)∈(Δf^x(i+1)-fx(i+1),Δf^x(i+1)+fx(i+1))其中，fx(i+1)=|Δf^x(i+1)-Δf^x(i)|;步骤8：利用进行相偏估计，且令n=n+1；步骤9：重复上述步骤(4)～(8)，直到count=P或i=I时，得到载波偏差精估计值和退出迭代循环，执行下一步；步骤10：利用最终得到的载波偏差精估计值和对细同步输入信号yk进行补偿得zk：再经过解调、译码判决之后即可恢复出发送端的原始信息。FDA0000336735320000011.jpg,FDA0000336735320000012.jpg,FDA0000336735320000013.jpg,FDA0000336735320000014.jpg,FDA0000336735320000015.jpg,FDA0000336735320000016.jpg,FDA0000336735320000017.jpg,FDA0000336735320000019.jpg,FDA00003367353200000110.jpg,FDA00003367353200000111.jpg,FDA00003367353200000112.jpg,FDA00003367353200000113.jpg,FDA00003367353200000114.jpg,FDA0000336735320000021.jpg,FDA0000336735320000022.jpg,FDA0000336735320000025.jpg,FDA0000336735320000026.jpg,FDA0000336735320000027.jpg,FDA0000336735320000028.jpg,FDA0000336735320000029.jpg,FDA00003367353200000210.jpg,FDA00003367353200000211.jpg...
【技术特征摘要】
【专利技术属性】
技术研发人员：许华，沈海鸥，李晓东，
申请(专利权)人：中国人民解放军空军工程大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。