[发明专利]一种长周期扩频码的频率折叠时频并行搜索方法

摘要

基于频率折叠的长周期扩频码时频并行搜索方法属于扩频通信技术领域，其特征在于，通过将同一段接收信号映射到不同的频率单元格下进行折叠，实现了对频率单元的并行搜索，再利用多个码相位并行搜索方法实现对多个码相位单元的并行搜索，在实现时，根据待搜索的频率范围以及频率搜索步长，确定待搜索的频率单元数，再将映射到不同频率单元格上的同一段接收信号进行折叠，然后把折叠后的信号用FFT技术实现该信号与本地伪码间的并行相关，实现多个码相位的并行搜索，本方法在不显著增加算法复杂度和硬件实现规模的条件下，同时实现了载波频率单元和多个码相位单元的并行搜索，提高了搜索速度，减少了伪码捕获时间。

主权项

1、一种长周期扩频码的频率折叠时频并行搜索方法，其特征在于，所述方法是在接收端数字信号处理器中按以下步骤实现的：步骤(1)设定待搜索的频率范围fu和频率搜索步长fstep，得到待搜索的频率单元数K＝fu/fstep；步骤(2)接收N/2点经前端A/D采样后的数字信号r(n)，n＝0，1，2，...，N/2-1，再把这N/2点数字信号分别映射到不同频率单元的频率格上，记作rk(n)，k＝1，2，...，K，k为不同的频率格；然后把这K组信号进行折叠得到新的信号b(n)：$b\left(n\right)=\underset{k=1}{\overset{K}{\Sigma }}{r}^{\left(k\right)}\left(n\right),n=\mathrm{0,1,2},...,N/2-1;$ 步骤(3)对步骤(2)得到的折叠后的信号b(n)补零至N点得到新的补零信号c(n)：$c\left(n\right)=\left\{\begin{array}{cc}b\left(n\right),& 0\le n\le (N/2-1)\\ 0,& N/2\le n\le (N-1)\end{array}\right.,$ 再对c(n)进行FFT运算，得到c(m)：$C\left(m\right)=\underset{n=0}{\overset{N-1}{\Sigma }}c\left(n\right)e^{-j\frac{2\pi }{N}\mathrm{nm}},$ m＝0，1，2，...，N-1，$j=\sqrt{-1};$ 步骤(4)对本地伪码l(n)作FFT运算，得到L(m)：$L\left(m\right)=\underset{n=0}{\overset{N-1}{\Sigma }}l\left(n\right)e^{-j\frac{2\pi }{N}\mathrm{nm}},$ m＝0，1，2，...，N-1，$j=\sqrt{-1},$ 对C(m)取复共轭，得到C*(m)，把C*(m)与L(m)相乘，得到Y(m)：Y(m)＝C*(m)L(m)，m＝0，1，2，...，N-1，对Y(m)进行IFFT变换后，得到相关结果d(n)：$d\left(n\right)=\frac{1}{N}\underset{m=0}{\overset{N-1}{\Sigma }}Y\left(m\right)e^{j\frac{2\pi }{N}\mathrm{mn}},$ n＝0，1，2，...，N-1；步骤(5)存储前N/2+1点结果‖d(n)‖2的值，n＝0，1，2，...，N/2，“‖‖”表示取模；步骤(6)重复步骤(2)至步骤(5)，共F次，将步骤(5)所得结果记为‖di(n)‖2，i表示非相干累加脚标，然后进行非相干累加，得到结果e(n)：$e\left(n\right)=\underset{i=0}{\overset{F-1}{\Sigma }}{\left|\right|d_i\left(n\right)\left|\right|}^2,$ n＝0，1，2，...，N/2；步骤(7)从步骤(6)得到的e(n)中选出最大值emax：$e_{\max }=\underset{0\le n\le N/2}{\max }\left\{e\left(n\right)\right\},$ 将emax与设定的门限Po进行比较，$P_o=\sqrt{-\mathrm{\lambda N}({\delta }^2+A^2)\ln P_{\mathrm{fa}}},$ 其中Pfa表示单个待搜索单元的虚惊概率、δ2表示噪声方差、A表示接收信号幅度、N表示FFT长度，λ为一系数与频率折叠次数有关；步骤(8)若最大值emax大于门限，则根据该最大值所在位置来确定伪码相位，若最大值emax 小于门限则移动N/2+1点本地伪码相位，对另外一段伪码相位进行搜索，重复步骤(2)至步骤(7)，直到找到大于门限的emax或搜索完整个码相位空间为止；步骤(9)根据步骤(8)所得到的伪码相位，再依次对每个频率单元格进行搜索，得到载波频率。