一种扩频信号速率自适应捕获方法与流程

1.本发明属于通信技术领域，具体涉及一种扩频信号速率自适应捕获方法。


背景技术：

2.直接序列扩频信号因其具备良好的抗干扰能力和隐蔽性，在现有移动通信系统、卫星通信系统中得到了广泛应用。在卫星通信系统中，测控链路、用户信令信道经常使用扩频通信体制。由于地面终端大小、能力、通信信道状态差异较大，因此星载扩频接收机通常需要同时处理多种速率的同时接入。传统的方式是将每路星载扩频接收机配置为固定速率的信道，每路扩频接收机只处理固定速率的信号，这种僵化的分配方式会导致处理资源的浪费和效率的降低，从而浪费了宝贵的星载处理资源。


技术实现要素：

3.为了克服现有技术的不足，本发明提供了一种扩频信号速率自适应捕获方法，首先在体制设计时保证不同速率档所使用的扩频码之间具有良好的互相关特性；在此基础上通过频域多路并行积分矩阵实现多速率档pn码与输入信号的时域相关；提出一种基于多路相关的扩频信号信噪比估计方法，采用基于多路相关的信噪比估计算法，计算出各支路信噪比，并进行门限判决，超出门限的信噪比最大值对应的支路即为正确的速率支路，最大相关峰值对应的时域位置即为码相位起始位置，当前搜索的频点即为对应的正确频偏，从而完成扩频信号的自适应捕获。该方法可在未知信号速率的条件下，自适应完成扩频信号的捕获，给出正确的速率、频偏与码相位，极大增加了扩频接收机的灵活性和鲁棒性。
4.本发明解决其技术问题所采用的技术方案包括如下步骤：
5.步骤1：定义输入信号为a/d采样后的复数型中频信号s
if
(n)，表示为：
[0006][0007]
其中，x(n)为输入扩频序列；fi为未知的、包含多普勒频移在内的中频频率；a为中频信号的幅度；n为样点序列计数值；
[0008]
步骤2：采用下变频模块对输入信号进行变频，具体如下：
[0009][0010][0011]
其中，s
ddc
(n)为经过下变频后的零频信号，f
l
为本地搜索频偏，δf为输入频偏与本地搜索频偏的差值；
[0012]
步骤3：采用fft变换模块对下变频模块的输出数据继续fft运算，具体如下：
[0013][0014]
其中，n为fft样点个数，s
ddc
(k)为s
ddc
(n)在频域上的表达；
[0015]
步骤4：将本地pn序列经过fft运算后的结果存储在本地pn码存储器中，输出m路信号描述为：
[0016][0017]
其中，ym(n)为第m支路本地pn序列，ym(k)为第m支路本地pn序列经过fft运算后的频域信号表达式；
[0018]
步骤5：采用多路并行积分矩阵将步骤3输出结果与本地pn序列存储器输出的多路信号进行并行积分运算，计算过程如下所示：
[0019][0020]
其中，c(m,k)表示第m支路并行积分运算结果；
[0021]
步骤6：采用m路并行ifft模块将多路并行积分矩阵输出的m路数据进行并行ifft运算，计算过程如下：
[0022][0023]
式中，cm(n)表示第m支路并行ifft运算的结果；
[0024]
对式(6)进行离散相位变换，从而得到：
[0025][0026]
考虑实际传输系统中的噪声，式(7)进一步描述为：
[0027][0028]
其中，r(τ)为本地pn序列与输入扩频序列的码相位差，nm为第m支路的噪声；
[0029]
步骤7：采用基于多路相关的信噪比估计算法，计算正确速率档对应的支路，并估计出信噪比，具体过程如下：
[0030]
步骤7-1：分别求取m个支路上的功率；
[0031]
步骤7-2：分别找到m个支路上的功率最大值点，对应该支路的自相关峰值，判定此时相位差τ≈0；
[0032]
步骤7-3：在m个自相关峰值中选择最大值，将该支路作为真正的信号支路，该支路对应的本地扩频码与输入扩频序列一致，此时信号峰值功率描述为：
[0033]
p＝(r(0)e
j2πδfn
)2，ym(n)＝x(n)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0034]
步骤7-4：由于不同支路的扩频码是不相关的，因此判定其余支路的功率为噪声，取其余支路的功率平均值作为噪声功率：
[0035][0036]
式中，为除最大峰支路之外的其余支路的功率；
[0037]
步骤7-5：计算出最大峰值对应的信噪比snr＝[p]-[σ2]；
[0038]
步骤8：对最大峰值对应的信噪比snr进行判决，当信噪比snr超过预设门限snr_thre时，判定捕获成功，输出相应标志信息；此时，最大峰值对应的支路即代表正确的速率支路，最大峰值在当前支路中的位置即为码相位起始位置，当前搜索的频点即为对应的正确频偏，从而完成扩频信号的自适应捕获；
[0039]
若信噪比snr小于等于门限snr_thre，判定本次捕获未成功，此时修改下变频模块中的下变频频点，返回步骤2，进行下一个频点的搜索捕获；
[0040]
优选地，所述预设门限值snr_thre取值范围为[0:10]。
[0041]
本发明的有益效果如下：
[0042]
1)本发明的基于多路相关的信噪比估计算法只需要用到简单的比较器与累加器，消耗资源少，性能优良，非常适合处理资源有限的星载扩频接收机。
[0043]
2)本发明的扩频信号速率自适应捕获方法可在未知信号速率的条件下，自适应完成扩频信号的捕获，给出正确的速率、频偏与码相位，大大增加了星载扩频接收机的灵活性和鲁棒性。
附图说明
[0044]
图1为本发明方法的原理图。
具体实施方式
[0045]
下面结合附图和实施例对本发明进一步说明。
[0046]
本发明解决的技术问题是：克服现有技术的不足，提出一种扩频信号速率自适应捕获方法，在扩频信号捕获时自适应识别通信速率，解决了星载扩频接收机灵活自适应捕获可变通信速率扩频信号的问题。
[0047]
本发明方法首先将ad采样后的信号进行下变频，下变频频点为搜索频点，以抵消残留频偏；然后将下变频后的数据送往fft模块，转换为频域信号；将多档速率对应的不同pn码进行fft运算，将运算结果存储在本地存储资源中，从而避免了多个fft运算带来的资源消耗；采用多路并行积分矩阵，将频域数据信号与多路频域pn码信号进行多路并行积分，将积分后的结果送往m路并行ifft模块，进行ifft运算；将m个支路的运算结果送往信噪比估计模块，采用基于多路相关的信噪比估计算法，计算出各支路信噪比，并进行门限判决，超出门限的信噪比最大值对应的支路即为正确的速率支路，ifft对应的峰值即为码相位起始位置，当前搜索的频点即为对应的正确频偏。从而完成扩频信号的自适应捕获。
[0048]
步骤1：定义输入信号为a/d采样后的复数型中频信号s
if
(n)，表示为：
[0049][0050]
其中，x(n)为输入扩频序列；fi为未知的、包含多普勒频移在内的中频频率；a为中频信号的幅度；n为样点序列计数值；
[0051]
步骤2：采用下变频模块对输入信号进行变频，具体如下：
[0052][0053]
其中，s
ddc
(n)为经过下变频后的零频信号，f
l
为本地搜索频偏，δf为输入频偏与本地搜索频偏的差值；
[0054]
步骤3：采用fft变换模块对下变频模块的输出数据继续fft运算，具体如下：
[0055][0056]
其中，n为fft样点个数，s
ddc
(k)为s
ddc
(n)在频域上的表达；
[0057]
步骤4：将本地pn序列经过fft运算后的结果存储在本地pn码存储器中，输出m路信号描述为：
[0058][0059]
其中，ym(n)为第m支路本地pn序列，ym(k)为第m支路本地pn序列经过fft运算后的频域信号表达式；
[0060]
步骤5：采用多路并行积分矩阵将步骤3输出结果与本地pn序列存储器输出的多路信号进行并行积分运算，计算过程如下所示：
[0061][0062]
其中，c(m,k)表示第m支路并行积分运算结果；
[0063]
步骤6：采用m路并行ifft模块将多路并行积分矩阵输出的m路数据进行并行ifft运算，计算过程如下：
[0064][0065]
式中，cm(n)表示第m支路并行ifft运算的结果；
[0066]
对式(6)进行离散相位变换，从而得到：
[0067][0068]
考虑实际传输系统中的噪声，式(7)进一步描述为：
[0069][0070]
其中，r(τ)为本地pn序列与输入扩频序列的码相位差，nm为第m支路的噪声；
[0071]
步骤7：采用基于多路相关的信噪比估计算法，计算正确速率档对应的支路，并估计出信噪比，具体过程如下：
[0072]
步骤7-1：分别求取m个支路上的功率；
[0073]
步骤7-2：分别找到m个支路上的功率最大值点，对应该支路的自相关峰值，判定此时相位差τ≈0；
[0074]
步骤7-3：在m个自相关峰值中选择最大值，将该支路作为真正的信号支路，该支路
对应的本地扩频码与输入扩频序列一致，此时信号峰值功率描述为：
[0075]
p＝(r(0)e
j2πδfn
)2，ym(n)＝x(n)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0076]
步骤7-4：由于不同支路的扩频码是不相关的，因此判定其余支路的功率为噪声，取其余支路的功率平均值作为噪声功率：
[0077][0078]
式中，为除最大峰支路之外的其余支路的功率；
[0079]
步骤7-5：计算出最大峰值对应的信噪比snr＝[p]-[σ2]；
[0080]
步骤8：对最大峰值对应的信噪比snr进行判决，当信噪比snr超过预设门限snr_thre时，判定捕获成功，输出相应标志信息；此时，最大峰值对应的支路即代表正确的速率支路，最大峰值在当前支路中的位置即为码相位起始位置，当前搜索的频点即为对应的正确频偏，从而完成扩频信号的自适应捕获；
[0081]
若信噪比snr小于等于门限snr_thre，判定本次捕获未成功，此时修改下变频模块中的下变频频点，返回步骤2，进行下一个频点的搜索捕获；
[0082]
预设门限值根据经验值选取，通常snr_thre取值范围为[0:10]。
[0083]
具体实施例：
[0084]
本发明所提方法的输入条件如下：输入信号为ad采样后的复数型中频信号s
if
(n)，可表示为
[0085][0086]
其中，x(n)为输入重复扩频序列，不同速率档对应的扩频序列不同，不同速率档之间的扩频序列从系统设计上具有良好的互相关特性；fd为未知的、包含多普勒频移在内的中频频率；a为中频信号的幅度；n为样点序列计数值。
[0087]
本发明实现框图如图1所示。包括如下模块：下变频模块1、fft变换模块2、本地pn码(fft后)存储器3、多路并行积分矩阵4、m路并行ifft模块5、非相干累加模块6、基于多路相关的信噪比估计模块7、门限判决模块8。
[0088]
下变频模块1：对a/d采样后的带有残留频偏的中频信号进行变频，将变频后的零频信号送往后续模块进行处理。变频频点为保持一定间隔的搜索频点，当一次捕获运算完成后，若判决本次捕获失败，则需要挪动变频频点，进行下一次捕获运算。
[0089]
fft变换模块2：对下变频模块1的输出数据继续fft运算，将运算结果输出。
[0090]
本地pn码(fft后)存储器3：为节省处理资源，将本地pn码经过fft运算后的结果预先存储在本地存储器中，待需要时按照时序要求读出。
[0091]
多路并行积分矩阵4：将fft变换模块2输出的数据与本地pn码(fft后)存储器3输出的多路信号进行并行积分运算，将运算结果送往下一级模块处理。
[0092]
m路并行ifft模块5：将多路并行积分矩阵4输出的m路数据进行并行ifft运算，然后输出。
[0093]
基于多路相关的信噪比估计模块6：根据m路并行ifft模块5输出的m路数据，采用基于多路相关的信噪比估计算法，计算出正确速率档对应的支路，并估计出信噪比，输出给门限判决模块7模块。
[0094]
门限判决模块7模块：对输入的信噪比估计结果进行判决，当超出门限时，认为捕
获成功，输出相应标志信息。
[0095]
具体信号处理流程如下：
[0096]
1、下变频模块1对ad采样后的带有残留频偏的中频信号进行变频。其处理流程描述如下：
[0097][0098]
2、fft变换模块2对输入的s
ddc
(n)信号进行fft(傅里叶变换)操作。其处理流程描述如下：
[0099][0100]
3、本地pn码(fft后)存储器3：将本地pn码经过fft运算后的结果预先存储在本地存储器中，按时序要求输出。其输出信号可描述为：
[0101][0102]
4、多路并行积分矩阵4：将fft变换模块2输出的数据与本地pn码(fft后)存储器3输出的多路信号进行并行积分运算，其计算过程如下所示：
[0103][0104]
5、m路并行ifft模块5：将多路并行积分矩阵4输出的m路数据进行并行ifft运算，其计算过程如下所述：
[0105][0106]
从而得到
[0107][0108]
当输入扩频序列与本地扩频序列一致时，二者积分可以用自相关函数r(τ)来描述；当输入扩频序列与本地扩频序列无关时，根据扩频序列的特性，互相关结果约等于0。同时，考虑实际传输系统中的噪声，上式可以进一步描述为：
[0109][0110]
6、基于多路相关的信噪比估计模块6：根据m路并行ifft模块5输出的m路计算结
果，采用基于多路相关的信噪比估计算法，计算出正确速率档对应的支路。具体过程如下：
[0111]
1)分别求取m个支路上的功率；
[0112]
2)分别找到m个支路上的功率最大值点，对应该支路的自相关峰值，认为此时相位差τ≈0；
[0113]
3)在m个自相关峰值中选择最大值，认为该支路是真正的信号支路，该支路对应的本地扩频码与输入扩频序列一致。此时信号峰值功率可以描述为：
[0114]
p＝(r(0)e
j2πδfn
)2，ym(n)＝x(n)
[0115]
4)由于不同支路的扩频码是不相关的，因此认为其余支路的功率为噪声，取其余支路的功率平均值作为噪声功率：
[0116][0117]
5)计算出最大峰对应的信噪比snr＝[p]-[σ2]，送往后级模块进行判决。
[0118]
7、门限判决模块7模块
[0119]
对输入的信噪比估计结果进行判决，当超出门限snr_thre时，认为捕获成功，输出相应标志信息。此时，最大峰对应的支路即代表正确的速率支路，最大峰值在当前支路中的位置即为码相位起始位置，当前搜索的频点即为对应的正确频偏。从而完成扩频信号的自适应捕获。
[0120]
若信噪比估计结果小于门限snr_thre，认为本次捕获未成功。此时改变下变频模块1中的下变频频点，进行下一个频点的搜索捕获。
[0121]
门限值snr_thre根据经验值选取，通常取值范围为[0:10]。