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

权利要求书

1、一种长周期扩频码的频率折叠时频并行搜索方法，其特征在于，所述方法是在接收端数字信号处理器中按以下步骤实现的：

步骤(1)设定待搜索的频率范围f_u和频率搜索步长f_step，得到待搜索的频率单元数K＝f_u/f_step；

步骤(2)接收N/2点经前端A/D采样后的数字信号r(n)，n＝0，1，2，...，N/2-1，再把这N/2点数字信号分别映射到不同频率单元的频率格上，记作r^k(n)，k＝1，2，..，K，k为不同的频率格；然后把这K组信号进行折叠得到新的信号b(n)：

b(n)=Σk=1Kr(k)(n),n=0,1,2,...,N/2-1;]]>

步骤(3)对步骤(2)得到的折叠后的信号b(n)补零至N点得到新的补零信号c(n)：

c(n)=b(n),0≤n≤(N/2-1)0,N/2≤n≤(N-1),]]>

再对c(n)进行FFT运算，得到C(m)：

C(m)=Σn=0N-1c(n)e-j2πNnm,m=0,1,2,...,N-1,j=-1;]]>

步骤(4)对本地伪码l(n)作FFT运算，得到L(m)：

L(m)=Σn=0N-1l(n)e-j2πNnm,m=0,1,2,...,N-1,j=-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(n)=1NΣm=0N-1Y(m)ej2πNmn,n=0,1,2,...,N-1;]]>

步骤(5)存储前N/2+1点结果‖d(n)‖²的值，n＝0，1，2，...，N/2，“‖‖”表示取模；

步骤(6)重复步骤(2)至步骤(5)，共F次，将步骤(5)所得结果记为‖d_i(n)‖²，i表示非相干累加脚标，然后进行非相干累加，得到结果e(n)：

e(n)=Σi=0F-1||di(n)||2,n=0,1,2,...,N/2;]]>

步骤(7)从步骤(6)得到的e(n)中选出最大值e_max：

emax=max0≤n≤N/2{e(n)},]]>

将e_max与设定的门限P_o进行比较，

Po=-λN(δ2+A2)lnPfa,]]>其中P_fa表示单个待搜索单元的虚警概率、δ²表示噪声方差、A表示接收信号幅度、N表示FFT长度，λ为一系数与频率折叠次数有关；

步骤(8)若最大值e_max大于门限，则根据该最大值所在位置来确定伪码相位，若最大值e_max小于门限则移动N/2+1点本地伪码相位，对另外一段伪码相位进行搜索，重复步骤(2)至步骤(7)，直到找到大于门限的e_max或搜索完整个码相位空间为止；

步骤(9)根据步骤(8)所得到的伪码相位，再依次对每个频率单元的频率格进行搜索，得到载波频率。