升余弦成型mpsk信号的盲符号速率估计与定时方法
【专利摘要】本发明公开了一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。该方法包括：计算输入基带信号的功率谱密度；根据所述的功率谱密度确定符号速率范围；根据定时同步环路能同步的最大定时频率偏差和所述符号速率范围确定不同的初始符号速率估计值；依次利用所述初始符号速率估计值进行定时同步，并将能使定时同步环路同步的初始符号速率估计值加上此时的定时频率偏差作为符号速率值。本发明的特点在于：（1）具有较好实时性；（2）符号速率估计与定时同步都不受载波频率偏差的影响；（3）在进行盲符号速率估计的同时完成了符号定时同步。本方法也适用于根升余弦成型的MPSK信号。
【专利说明】升余弦成型MPSK信号的盲符号速率估计与定时方法
【技术领域】
[0001]本发明涉及数字通信中一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。
【背景技术】
[0002]多进制相位键控(M-ary Phase Shift Keying，MPSK)调制信号具有良好的抗干扰性能和工程上易于实现等优点，在航天航空领域的飞行器遥测、遥控、遥感、数传、通信等领域得到广泛应用。而在非协作的通信信号分析中，符号速率作为一个解调的重要参数，需要准确估计。符号速率的准确估计，对于定时同步、调制方式的识别、特定信号搜索以及解调等均有重要意义。
[0003]目前，常见的实现PSK信号盲符号速率估计的几种代表性算法主要有:(I)基于循环谱的载频和码速率估计算法，该算法需要较大的采样数据量，不能满足工程上实时性的要求。(2)基于小波变换的突变点和码速率估计算法，该方法存在最优尺度和母函数的选取问题。(3)基于最大似然估计的PSK信号参数估计算法，该算法需要多维搜索，实时性仍较差。
【发明内容】
[0004]鉴于上述现有技术所存在的问题，本发明的目的是提供一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。特点在于:(1)具有较好实时性；(2)符号速率估计与定时同步都不受载波频率偏差的影响；(3)在进行符号速率估计的同时完成了符号定时同步。本方法也适用于根升余弦成型的MPSK信号。
[0005]本发明的目的通过以下技术方案实现的:
[0006]一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，包括:
[0007]A、根据输入的基带信号求信号功率谱密度；
[0008]B、根据所述的功率谱密度求出符号速率范围；
[0009]C、根据定时同步环路的最大定时频率偏差,在所述符号速率范围内，等间隔的确定η个符号速率值作为初始符号速率值；
[0010]D、依次选取所述η个初始符号速率值，计算定时同步环路之前的抽取系数，然后对抽取后的信号进行定时同步。当定时同步后，输出定时同步信号和符号速率。
[0011 ] 所述步骤A进一步包括:
[0012]Al、设置采样点数Nsaniple，采样频率fs ；
[0013]A2、设置 FFT 点数 Nfft ;
[0014]A3、设置窗函数类型与窗函数长度Nwindw ；
[0015]A4、设置重叠数据长度Ntjverlap ；
[0016]A5、设置抽取模块I的抽取系数为S1=I ；[0017]A6、计算基带信号的功率谱密度Sy⑴。
[0018]所述步骤B进一步包括:
[0019]B1、搜索Sy (f)找到其最大值，其值记为V ；
[0020]B2、在第一奈奎斯特频率范围(O~0.5fs)内搜索满足Sy (f)≥0.5v的最大频率点，记为fL;
[0021]B3、在第二奈奎斯特频率范围(0.5fs~fs)内搜索满足Sy (f)≥0.5v的最小频率点，记为fH ;
[0022]BMf (fs-fH+f1)作为基带信号带宽的估计值Bg
[0023]B5、根据基带信号带宽的估计值Bg，如果满足Bg <fx/8，则设置S1=S1*2,并返
回步骤A6 ;如果满足Bg>fg/2，则设置S1=S1 /2,并返回步骤A6 ;若都不满足，则执行步骤B6 ；
[0024]B6、根据信号带宽估计值&，求得信号速率估计范围为(0.8Bg ~1.2Bg)。
[0025]所述步骤C进一步包括:
[0026]Cl、确定定时同步环路支持的最大定时频率偏差ΔfTMAX;
[0027]C2、根据所述符号速率估计范围(0.8Bg~1.2Bg)和最大定时频率偏差Δfimx，确定
η个长度为2 Δ fTMAX的频率区间，每个区间中心频率值为fn(i-1,2,...,n)，将其作为定时同步
环路的初始符号速率。
[0028]所述步骤D进一步包括:
[0029]D1、设置i=l，设置每个区间的定时检测时间为tp ；
[0030]D2、选取初始符号速率估计值/n_，计算抽取模块2的抽取因子S2 ；
[0031]D3、等待tp时间后，检测环路滤波器的积分支路，若其值不收敛，令i=i+l，返回步骤D2 ;若其值收敛，则表示定时已同步，输出符号定时同步信号和相对定时频率偏差AfT，进入步骤D4 ；
[0032]D4、将此时的初始符号速率估计值定时频率偏差Λ fT相加的结果作为基带信号的符号速率估计值，并输出该值。
[0033]由上述本发明提供的技术方案可以看出，本发明原理简单，计算量小，具有较好实时性，利于工程实现，在进行符号速率估计的同时实现了符号定时同步，并且符号速率估计与定时同步都不受载波频率偏差干扰；通过对实际采集的信号分析表明，由该方法对符号速率进行估计时，具有优越的性能。
【专利附图】
【附图说明】
[0034]图1为本发明的系统框图
[0035]图2为本发明的符号速率范围分割图
[0036]图3为本发明的处理流程图【具体实施方式】
[0037]本发明提供了一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法。本发明的核心为:对接收到的基带信号进行功率谱密度估计，通过其功率谱密度估计出符号速率范围，根据定时同步环路能同步的最大定时频率偏差在符号速率范围中确定不同的符号速率估计值，将能使定时同步环路同步的符号速率估计值加上相对定时频率偏差作为符号速率输出，同时输出定时同步信号。
[0038]本发明所述针对升余弦成型MPSK信号的盲符号速率估计与定时装置的结构示意图如图1所示,其中包括:
[0039]1、抽取模块1:用于对基带信号进行抽取，其抽取系数受到后面信号带宽估计模块的控制。
[0040]2、功率谱估计模块:用于估计接收到基带信号的功率谱密度，并且将估计出来的功率谱密度送入后面的信号带宽估计模块；本发明中设置FFT的点数为256点，设置的窗函数类型为海明窗，窗函数长度为255，参与计算的采样点数为20000，重叠数据的长度128。
[0041]3、信号带宽估计模块:用于估计信号的带宽，并根据信号的带宽控制抽取模块I的抽取系数。
[0042]设输入的功率谱密度为Sy (f)，则搜索其最大值，记为V，然后在第一奈奎斯特频率范围(O~0.5fs)内搜索满足Sy(f)≥0.5v的最大频率点，记为4 ;在第二奈奎斯特频率范围(0.5fs~fs)内搜索满足Sy(f)≥0.5v的最小频率点，记为fH，则信号的带宽估计
值巧=(/；-Zb +Λ>当输入信号带宽估计值A満足Β,<.?,β，则设置抽取模块I的抽取系数S1=S1XZ,当输入信号带宽估计值巧满足％ >又/2，则设置抽取模块I的抽取系数S1=S1 + ],反之则输出信号带宽估计值Ig。
[0043]4、符号速率范围估计模块:用于估计信号的速率范围，根据输入的信号带宽估计值，符号的速率范围为(0.8瓦~1.2瓦
S'--S5 /
[0044]5、符号速率估计值选取模块:用于从符号速率范围中确定符号速率估计值，同时根据选取的符号速率估计值计算抽取模块2的抽取系数。
[0045]确定符号速率估计值的标准如图2所示，为从最小符号速率开始，下一个符号速率估计值距离前一个符号速率估计值的频率间隔为定时同步环路的最大定时频率偏差的两倍。
[0046]6、抽取模块2:用于对基带信号进行抽取，使采样率满足所选符号速率估计值的16倍过采样，其抽取系数受到符号速率估计值选取模块的控制。
[0047]7、定时同步模块:用于对信号进行定时同步，同时输出定时同步信号和矫正后的符号速率估计值作为符号速率输出。
[0048]本发明所述方法的处理流程图如图3所示，包括如下步骤:
[0049]步骤31:设置采样点数Nsample,采样频率fs。
[0050]步骤32:设置FFT点数Nfft。
[0051]步骤33:设置窗函数类型为海明窗与窗函数长度Nwind?。
[0052]步骤34:设置重叠数据长度NQverlap。[0053]步骤35:设置抽取模块I的抽取系数为S1=I。
[0054]步骤36:计算基带信号的功率谱密度Sy (f)。
[0055]步骤37:搜索Sy (f)找到其最大值，其值记为V。
[0056]步骤38:在第一奈奎斯特频率范围(O~0.5fs)内搜索满足Sy (f)≥0.5v的最大频率点，记为4。
[0057]步骤39:在第二奈奎斯特频率范围(0.5fs~fs)内搜索满足Sy (f)≥0.5v的最小频率点，记为fH。
[0058]步骤310:将(fs-fH+fL)作为基带信号带宽的估计值&
[0059]步骤311:根据基带信号带宽的估计值&，如果满足lg</s/8，则执行步骤312 ；
如果满足A 2/、/2，则执行步骤313 ;若都不满足，则执行步骤314。
[0060]步骤312:设置抽取模块I抽取系数S1=S1 X 2，执行步骤36。
[0061]步骤313:设置抽取模块I抽取系数SpSi + 2，执行步骤36。
[0062]步骤314: 根据信号带宽估计值&，求得信号速率估计范围为(0.8&~1.2；^)。
[0063]步骤315:选择定时同步环路的最大定时频率偏差Λ fTMAX。
[0064]步骤316:根据所述符号速率估计范围(0.8&~1.2:^)和最大定时频率偏差Λ fTMAX,确定η个长度为2Λ fTMAX的频率区间，此时η的值为n = round[{l.2-0.8?/(2x AfTMAX)],每个区间中心频率值为/ -1,2,...,?)，将其作为定时同步环路的初始符号速率。
[0065]步骤317:设置i=l，设置每个区间的定时检测时间为tp。
[0066]步骤318:选取初始符号速率估计值f\计算抽取模块2的抽取因子S2。
[0067]步骤319:等待tp时间后，检测环路滤波器的积分支路，若其值不收敛，令i=i+l，执行步骤320 ;若其值收敛，则表示定时已同步，输出符号定时同步信号和相对定时频率偏差Λ fT，进入步骤321。
[0068]步骤320:令i=i+l，选取下一个初始符号速率估计值。
[0069]步骤321:输出定时同步信号和相对定时频率偏差AfT。
[0070]步骤322:输出基带信号的速率估计值/r =/fn? +Afr。
【权利要求】
1.一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，包括: A、根据输入的基带信号求信号功率谱密度； B、根据所述的功率谱密度求出符号速率范围； C、根据定时同步环路的最大定时频率偏差，在所述符号速率范围内，等间隔的确定η个符号速率值作为初始符号速率值； D、依次选取所述η个初始符号速率值，计算定时同步环路之前的抽取系数，然后对抽取后的信号进行定时同步。当定时同步后，输出定时同步信号和符号速率。
2.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤A进一步包括: Al、设置采样点数N—采样频率fs ； A2、设置FFT点数Nfft ； A3、设置窗函数类型与窗函数长度Nwindmt ； A4、设置重叠数据长度NOTOTlap ； A5、设置抽取模块I的抽取系数为S1=I ； A6、计算基带信号的功率谱密度Sy (f)。
3.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤B进一步包括: B1、搜索Sy(f)找到其最大值，其值记为V ； B2、在第一奈奎斯特频率范围(O~0.5fs)内搜索满足Sy(f)≥0.5v的最大频率点，记为fL ； B3、在第二奈奎斯特频率范围(0.5fs~fs)内搜索满足Sy (f)≥0.5v的最小频率点，记为fH ； B4、将(fs-fH+f^)作为基带信号带宽的估计值& B5、根据基带信号带宽的估计值豆，如果满足巧^/:/8，则设置S1=S1XZ,并返回步骤A6 ;如果满足& >./:/2，则设置S1=S1S,并返回步骤A6 ;若都不满足，则执行步骤B6 ； B6、根据信号带宽估计值瓦，求得信号速率估计范围为(0.8&~\巩)。
4.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤C进一步包括: Cl、确定定时同步环路支持的最大定时频率偏差Afmx ；C2、根据所述符号速率估计范围(60.8?~1.2?)和最大定时频率偏差Afimx，确定η个长度为2Afmx的频率区间，每个区间中心频率值为/),(/ = 1,2,.η)将其作为定时同步环路的初始符号速率。
5.根据权利要求1所述的一种针对升余弦成型MPSK信号的盲符号速率估计与定时方法，其特征在于，所述步骤D进一步包括: D1、设置i=l，设置每个区间的定时检测时间为tp ；D2、选取初始符号速率估计值/π，计算抽取模块2的抽取因子S2 ； D3、等待tp时间后，检测环路滤波器的积分支路，若其值不收敛，令i=i+l，返回步骤D2 ;若其值收敛，则表示定时已同步，输出符号定时同步信号和相对定时频率偏差AfT，进入步骤D4 ； D4、将此时的初始符号速率估计值1和定时频率偏差Λ fT相加的结果作为基带信号的符号速率估计值，并输出该值。
【文档编号】H04L25/02GK103916341SQ201410152732
【公开日】2014年7月9日 申请日期:2014年4月16日 优先权日:2014年4月16日
【发明者】欧静兰, 吴皓威, 刘晓明, 王牌, 邹玉涛, 童玄, 杨力生, 王韬 申请人:重庆大学