CN101635936B - 对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法 - Google Patents
对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法 Download PDFInfo
- Publication number
- CN101635936B CN101635936B CN200910034158XA CN200910034158A CN101635936B CN 101635936 B CN101635936 B CN 101635936B CN 200910034158X A CN200910034158X A CN 200910034158XA CN 200910034158 A CN200910034158 A CN 200910034158A CN 101635936 B CN101635936 B CN 101635936B
- Authority
- CN
- China
- Prior art keywords
- sequence
- tau
- overbar
- prime
- kasami
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了一种对子集为Gold-like序列的Kasami大集合序列的FHT(Fast HadamardTransform)快速相关检测方法,属于适用于极低谱密度条件下使用扩频技术的无线通信系统的多目标相位检测方法。该方法在维度k=0,1,...,N′-2时,计算接收信号与Gold-like-Kasami大集合序列码的互相关函数:先固定两个维度,使用一个m序列构成Walsh-Hadamard矩阵,做FHT得到一维相关值,然后分别对两个维度取所有可能值得到所有二维相关值,重排二维相关值得到三维相关值,最后确定峰值所在的三个维度的取值;当k=N′-1时采用Kasami小集合序列的FHT快速相关检测方法。本方法只包含加法运算,没有乘法运算,极大地减少了运算量。
Description
技术领域
本发明涉及一种对子集为Gold-like序列的Kasami大集合序列的FHT(Fast HadamardTransform)快速相关检测方法,是一种适用于极低谱密度条件下使用扩频技术的无线通信系统的多目标相位检测方法,属于移动通信中的同步技术领域。
背景技术
在极低谱密度条件下的扩频无线通信中,例如深空通信、全球导航卫星系统(GNSS)、军事通信和无线传感器网络等领域,他们的主要特点是利用长PN码的良好的相关特性来实现扩频通信。对长PN码的快速捕获是极低谱密度扩频无线通信系统的一项关键技术,目前比较广为人知的PN码有最大移位寄存器序列(简称为m序列)、戈尔德(Gold)序列、Kasami序列等。
常规实现极长PN码快速捕获的方法有时域并行相关法和FFT(Fast Fourier Transform)频域相关法。
时域并行相关法的运算量极大,复杂度和成本极高:对于一个长度为p的PN码,运算量是p2量级的乘法和加法;对于K个长度为p的PN码,运算量是Kp2量级的乘法和加法。
FFT频域相关法具有方便快捷、易于实现PN码相位的并行搜索和载波频偏的并行搜索等优点,使K个长度为p的PN码快速相关的运算量降低为Kplog2p量级的复数乘法和加法,成为国内外研究长PN码快速捕获的热点。不过,FFT频域相关法存在以下问题:1、FFT进行复数乘法运算多,对于长PN码频域相关,FFT运算中存在很小数值的复权值,导致复数乘法的比特位数多,运算复杂度很高,其实现需要采用昂贵的高精度高速AD器件;2、FFT相关算法必须使用IFFT,使PN码捕获最终回到时域以确定PN码相位,这也增加了运算量;3、受当前DSP器件限制,FFT核长度已成为进一步提高搜索速度的瓶颈。
发明内容
本发明针对背景技术中传统的相关检测方法计算复杂度高的缺陷而提出一种不包含乘法运算的对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法。
本发明对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法,包括如下步骤:
当维度k=0,1,...,N′-2时,计算{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ),步骤如下:
步骤1:先固定两个维度,使用组成Gold-like-Kasami大集合序列的一个m序列构成Walsh-Hadamard矩阵,做FHT得到一维相关值;
步骤2:将步骤1中固定的两个维度中的一个维度取所有可能值,重复步骤(1),得到二维相关值;
步骤3:将步骤1中固定的两个维度中的另一个维度取所有可能值,重复步骤(1),得到所有的二维相关值,将所有二维相关值重排得到三维相关值,确定三维相关值峰值所在的三个维度的取值,峰值所在的三个维度的取值确定了在伽罗瓦域GF(q)中的Gold-like-Kasami大集合序列的码字组成及其码相位;
当k=N′-1时,{ak,k′(n)}即变成Kasami小序列集合,计算{r(n)}和{ak,k′(n)}的互相关函数C(k,k′;τ)采用Kasami小集合序列的FHT快速相关检测方法。
其中:Gold-like-Kasami大集合序列即是在伽罗瓦域GF(q)中,q=pn,p=2,当n≡0mod 4时的Kasami大集合序列,q-1为Gold-like-Kasami大集合序列码的码长,n为码长的阶数;m序列即最大移位寄存器序列;k、k′和τ是Gold-like-Kasami大集合序列码字组成及其码相位的三个维度;{r(n)}为接收信号;{ak,k′(n)}为Gold-like-Kasami大集合序列码;N′为v(c)的周期长度,v(c)是对序列Tcu每2n/2+1比特抽取1比特组成的周期为N′=(2n-1)/3的周期序列,c=0,1,2,T表示向量向左循环移位,其上标表示移位数。
本发明的有益效果如下:
本发明方法可以用FHT快速相关检测伽罗瓦域GF(q)中Gold-like-Kasami大集合序列的码字组成及其初始码相位。由于FHT中不包含乘法运算,只有加法运算,因此相对于用FFT进行相关检测的方法能大大降低计算复杂度,提高捕获速度。
根据一次q=2n点的FFT需要(q/2)log2q次复数乘法和qlog2q次复数加法,一次复数乘法需要三次实数乘法和三次实数加法,一次复数加法需要两次实数加法。将FFT频域相关法中的乘法运算折算成加法运算,折算方法为:q点FFT中使用的乘法折算成log2q次加法,具体效果如下:
当Gold-like-Kasami大集合序列是为类I时,用FFT频域相关法检测出序列的码字组成及其码相位,所需的实数加法运算量为:
用本发明提出的算法1.1和1.2完成同样的相关检测运算所需的实数加法运算量为:
用本发明提出的算法1.3完成同样的相关检测运算所需的实数加法运算量为:
当Gold-like-Kasami大集合序列是为类II时,用FFT频域相关法检测出序列的码字组成及其码相位,所需的实数加法运算量为:
用本发明提出的算法2完成同样的相关检测运算所需的实数加法运算量为:
综上,可见本发明方法相比FFT频域相关法:类I的算法1.3能节省90%以上的运算量,并且节约的运算量随着码长的增加而增大;类I的算法1.1、1.2和类II的算法2在码字较短的情况下能节约大量的运算量。
附图说明
图1是算法1.1的二维计算过程示意图。
图2是算法1.1的相关值三维排列示意图。
图3是算法1.2的二维计算过程示意图。
图4是算法1.2的相关值三维排列示意图。
图5是算法1.3.1的二维计算过程示意图。
图6是算法1.3.1的相关值三维排列示意图。
图7是算法1.3.2的二维计算过程示意图。
图8是算法2的二维计算过程示意图。
图9是本发明的方法流程示意图。
在图2、图4、图6中:k0为k轴上的任意值;k′0为k′轴上的任意值。
具体实施方式
文中的字母及符号定义:q-1为Gold-like-Kasami大集合序列码的码长;n为码长的阶数;m序列即最大移位寄存器序列;u是周期为2n-1的m序列;{r(n)}为接收信号;r*(·)表示r(·)的共轭;{ak,k′(n)}为Gold-like-Kasami大集合序列码;N′为v(c)的周期长度,v(c)(c=0,1,2)是对序列Tcu(c=0,1,2)每t(n)=2n/2+1比特抽取1比特组成的周期为N′=(2n-1)/3的周期序列;N=q-1;w是对u每2n/2+1比特抽取1比特组成的周期为2n/2-1的m序列;k、k′和τ是Gold-like-Kasami大集合序列码字组成及其码相位的三个维度;T表示向量向左循环移位,其上标表示移位数,下标为长度;A为维度k与维度τ的差值;B为维度k′与维度τ的差值;j为将向量分割成块的编号, q=2n;·表示取模,其下标为模值;为模2加法;Cx为任意值符号。
Gold-like-Kasami大集合序列即是在伽罗瓦域GF(q)中,q=pn,p=2,当n≡0mod4时的Kasami大集合序列,可以表述为:
Ht(n)(u)是由u、v(c)(c=0,1,2)两个序列组成的Gold-like序列,Ht(n)(u)序列集合可以表述为:
本发明的方法流程示意图如图9所示。
根据Gold-like-Kasami大集合序列的定义,{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ)可分为两类:(q=2n,N′=N/3,c=0,1,2)
类I:当Gold-like-Kasami大集合序列是由2个m序列和1个周期序列组成的周期为2n-1的序列时,{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ)可以表示为:
观察C(k,k′;τ)有三个变量:k、k′和τ,用FHT(Fast Hadamard Transform)计算三维相关值时,需要先固定两个变量,再求得二维相关值,最后得到三维相关值。
由式(1)可以看出,当k=N′- 1,{ak,k′(n)}即变成Kasami小序列集合,计算{r(n)}和{ak,k′(n)}的互相关函数c(k,k′;τ)采用Kasami小集合序列的FHT快速相关检测方法。
当k=0,1,...,N′-2时,计算c(k,k′;τ)的算法如下:
算法1.1:
先固定k=0,1,...,N′-2和τ=0,1,...,q-2两个维度,使用w生成Walsh-Hadamard矩阵,先得到k′和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
令n=n-τ q-1,式(1)可得到:
k=0,1,...,N′-2,
使用w生成Walsh-Hadamard矩阵可得:
生成方法参见文献[1][2][3],下同:
[1]A.Lempel,‘Hadamard and M-sequence transforms are permutationally similar’,Applied optics,vol.18,no.25,Dec.79
[2]Srdjan Z.Budisin,‘Fast PN sequence correlation by using FWT’,MediterraneanElectrotechnical Conference Proc.,1989,pp.513-515.
[3]Martin Cohn and Abraham Lempel,‘On Fast M-sequence Transforms’,IEEE Transactionson Information Theory,1977
(4)
其中
然后把个变换后的向量累加起来即可得到相关阵C(k,k′;τ)的第τ列,即:
使用这种方法,计算的复杂度为
算法1.1的二维计算过程如图1所示,算法1.1的三维相关值排列过程如图2所示。
当k和k′均与τ的差值固定时,即:
τ-k N′-1=A(A=0,1,2,...,N′-2)或 由式(1)可得到:
另外两种计算相关值的方法:
算法1.2:
先固定A=0,1,...,N′-2和τ=0,1,...,q-2的值,使用w生成Walsh-Hadamard矩阵,先得到k′和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
使用w生成Walsh-Hadamard矩阵可得:
(8)
其中
使用这种方法,计算的复杂度为
算法1.2的二维计算过程如图3所示,算法1.2的三维相关值排列过程如图4所示。
算法1.3:
先固定A=0,1,...,N′-2和 的值,使用u生成Walsh-Hadamard矩阵,先得到到k和τ或k′和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
使用w生成Walsh-Hadamard矩阵可得:
(10)
接下来分两种情况:
算法1.3.1:先计算k和τ两个维度的相关值。算法1.3.1的二维计算过程如图5所示(方框内是A值),算法1.3.1的三维相关值排列过程如图6所示。
算法1.3.2:先计算k′和τ两个维度的相关值。算法1.3.2的二维计算过程如图7所示,(方框内是B值),算法1.3.2的三维相关值排列过程与图6类似,只需将k与k′的维度互换。
使用这种方法,计算的复杂度为
类II:当Gold-like-Kasami大集合序列是由1个m序列和1个周期序列组成的周期组成的周期为2n-1的序列时,{r(n)}与{ak(n)}的互相关函数C(k,k′;τ)可以表示为:
算法2:
令n=n-τ q-1,式(11)可得到:
使用w生成Walsh-Hadamard矩阵可得:
其中
然后把个变换后的向量累加起来即可得到相关阵C(k;τ)的第τ列,即:
使用这种方法,计算的复杂度为
算法2的计算过程如图8所示。
下面分别为类I和类II算法与FFT频域相关法检测Gold-like-Kasami大集合序列时的运算比较:
表1:类I算法与FFT频域相关法检测Gold-like-Kasami大集合序列时的运算量比较
表2:类II算法与FFT频域相关法检测Gold-like-Kasami大集合序列时的运算量比较
Claims (3)
1.一种对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法,其特征在于:
当维度k=0,1,...,N′-2时,计算{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ),步骤如下:
步骤1:先固定两个维度,使用组成Gold-like-Kasami大集合序列的一个m序列构成Walsh-Hadamard矩阵,做FHT得到一维相关值;
步骤2:将步骤1中固定的两个维度中的一个维度取所有可能值,重复步骤(1),得到二维相关值;
步骤3:将步骤1中固定的两个维度中的另一个维度取所有可能值,重复步骤(1),得到所有的二维相关值,将所有二维相关值重排得到三维相关值,确定三维相关值峰值所在的三个维度的取值,峰值所在的三个维度的取值确定了在伽罗瓦域GF(q)中的Gold-like-Kasami大集合序列的码字组成及其码相位;
当k=N′-1时,{ak,k′(n)}即变成Kasami小序列集合,计算{r(n)}和{ak,k′(n)}的互相关函数C(k,k′;τ)采用Kasami小集合序列的FHT快速相关检测方法;
在上述内容中:Gold-like-Kasami大集合序列即是在伽罗瓦域GF(q)中,q=2n,当n≡0 mod 4时的Kasami大集合序列,q-1为Gold-like-Kasami大集合序列码的码长,n为码长的阶数;m序列即最大移位寄存器序列;k、k′和τ是Gold-like-Kasami大集合序列码字组成及其码相位的三个维度;{r(n)}为接收信号;r*(n)表示r(n)的共轭;{ak,k′(n)}为Gold-like-Kasami大集合序列码;N′为v(c)的周期长度,v(c)是对序列Tcu每2n/2+1比特抽取1比特组成的周期为N′=(2n-1)/3的周期序列,T表示向量向左循环移位,T的上标表示移位数,c=0,1,2;N=q-1;u为周期为2n-1的m序列;w是对u每2n/2+1比特抽取1比特组成的周期为2n/2-1的m序列。
2.根据权利要求1所述的对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法,其特征在于:所述{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ)分为两类:
q=2n,N′=N/3,c=0,1,2
类I:当Gold-like-Kasami大集合序列是由2个m序列和1个周期序列组成的周期为2n-1的序列时,{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ)表示为:
类II:当Gold-like-Kasami大集合序列是由1个m序列和1个周期序列组成的周期为2n-1的序列时,{r(n)}与{ak(n)}的互相关函数C(k;τ)表示为:
将C(k;τ)看作是C(k,k′;τ)的特殊情况,所以只用FHT计算{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ)。
3.根据权利要求1所述的对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法,其特征在于:所述Gold-like-Kasami大集合序列的构成有两种情况:
情况1:Gold-like-Kasami大集合序列是由2个m序列和1个周期序列组成,记做{ak,k′(n)},此时序列的码字组成及其码相位的不确定度为三维不确定度,由于FHT算法需要一个m序列生成Walsh-Hadamard矩阵,因此存在3种快速相关检测的方法;
情况2:Gold-like-Kasami大集合序列是由1个m序列和1个周期序列组成,记做{ak(n)},此时序列的码字组成及其码相位的不确定度为二维不确定度,由于FHT算法需要一个m序列生成Walsh-Hadamard矩阵,因此存在1种快速相关检测的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910034158XA CN101635936B (zh) | 2009-09-01 | 2009-09-01 | 对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910034158XA CN101635936B (zh) | 2009-09-01 | 2009-09-01 | 对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101635936A CN101635936A (zh) | 2010-01-27 |
CN101635936B true CN101635936B (zh) | 2012-02-01 |
Family
ID=41594947
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910034158XA Expired - Fee Related CN101635936B (zh) | 2009-09-01 | 2009-09-01 | 对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101635936B (zh) |
-
2009
- 2009-09-01 CN CN200910034158XA patent/CN101635936B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101635936A (zh) | 2010-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100553168C (zh) | 一种长周期扩频码的频率折叠时频并行搜索方法 | |
Kong | A deterministic compressed GNSS acquisition technique | |
CN104753561B (zh) | 一种抑制水声通信中多途干扰的直接序列扩频调制方法 | |
CN104168233B (zh) | 基于特征分解和梅西算法的dsss/uqpsk信号的伪码序列估计方法 | |
CN101561484B (zh) | 一种gnss信号伪码捕获方法 | |
CN104977504A (zh) | 一种电缆故障在线检测及定位装置 | |
CN100553169C (zh) | 基于频率折叠的短周期扩频码时频并行搜索方法 | |
CN100501441C (zh) | 定位全球定位系统相关峰值信号的方法和系统 | |
CN101635581B (zh) | 对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法 | |
CN101063717B (zh) | 一种用fft进行gps信号捕获的方法 | |
CN101320084A (zh) | 卫星定位信号的实时处理方法 | |
CN103926604A (zh) | 基于重叠差分循环相干积分的弱信号捕获方法 | |
Sagiraju et al. | Fast acquisition implementation for high sensitivity global positioning systems receivers based on joint and reduced space search | |
Choi et al. | Design of binary sequences with optimal cross-correlation values | |
CN101635936B (zh) | 对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法 | |
CN103969664B (zh) | 一种基于差分循环相干积分的gps弱信号捕获方法 | |
CN103441813B (zh) | 一种用于cdma系统的低相关二元序列集生成方法 | |
Fatima et al. | New chaotic binary sequences with good correlation property using logistic maps | |
CN105842714A (zh) | 卫星信号捕获方法及装置 | |
CN101762803B (zh) | 一种应用于全球导航卫星系统接收机的差分相关器 | |
CN101969353A (zh) | 卫星导航系统中扩频码的生成方法 | |
Akansu et al. | Improved correlation of Generalized DFT with nonlinear phase for OFDM and CDMA Communications | |
Akopian et al. | Fast Time-Recursive Block Correlators for Pseudorandom Sequences | |
CN101609139A (zh) | 卫星导航信号捕获中移除电文位影响的方法 | |
Zhu et al. | Performance evaluation of three BOC signal acquisition techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120201 Termination date: 20150901 |
|
EXPY | Termination of patent right or utility model |