CN101635581B - 对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法 - Google Patents
对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法 Download PDFInfo
- Publication number
- CN101635581B CN101635581B CN2009100341575A CN200910034157A CN101635581B CN 101635581 B CN101635581 B CN 101635581B CN 2009100341575 A CN2009100341575 A CN 2009100341575A CN 200910034157 A CN200910034157 A CN 200910034157A CN 101635581 B CN101635581 B CN 101635581B
- Authority
- CN
- China
- Prior art keywords
- sequence
- tau
- overbar
- kasami
- fht
- 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
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公布了一种对子集为Gold序列的Kasami大集合序列的FHT(Fast HadamardTransform)快速相关检测方法,属于适用于极低谱密度条件下使用扩频技术的无线通信***的多目标相位检测方法。该方法在维度k=0,1,...,q-2时,计算接收信号与Gold-Kasami大集合序列码的互相关函数:先固定两个维度,使用一个m序列构成Walsh-Hadamard矩阵,做FHT得到一维相关值,然后分别对两个维度取所有可能值得到所有二维相关值,重排二维相关值得到三维相关值,最后确定峰值所在的三个维度的取值;当k=q-1时采用Kasami小集合序列的FHT快速相关检测方法。本方法中使用的FHT算法,只包含加法运算,没有乘法运算,因而极大地减少了运算量。
Description
技术领域
本发明涉及一种对子集为Gold序列的Kasami大集合序列的FHT(Fast HadamardTransform)快速相关检测方法,是一种适用于极低谱密度条件下使用扩频技术的无线通信***的多目标相位检测方法,属于移动通信中的同步技术领域。
背景技术
在极低谱密度条件下的扩频无线通信中,例如深空通信、全球导航卫星***(GNSS)、军事通信和无线传感器网络等领域,他们的主要特点是利用长PN码的良好的相关特性来实现扩频通信。对长PN码的快速捕获是极低谱密度扩频无线通信***的一项关键技术,目前比较广为人知的PN码有最大移位寄存器序列(或简称为m序列)、戈尔德(Gold)序列、Kasami序列等等。
常规实现极长PN码快速捕获的方法有时域并行相关法和FFT(Fast FourierTransform)频域相关法。
时域并行相关法的运算量极大,复杂度和成本极高:对于一个长度为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序列的Kasami大集合序列的FHT快速相关检测方法。
技术方案:本发明为实现上述目的,采用如下技术方案:
本发明所提出的对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法。当维度k=0,1,...,q-2时,计算{r(n)}与{ak,k’(n)}的互相关函数C(k,k′;τ)采用的相关步骤如下:
步骤(1):先固定两个维度,使用组成Gold-Kasami大集合序列的一个m序列构成Walsh-Hadamard矩阵,FHT得到一维相关值;
步骤(2):将步骤1中固定的两个维度中的一个维度取所有可能值,重复步骤(1),得到二维相关值;
步骤(3):将步骤1中固定的两个维度中的另一个维度取所有可能值,重复步骤(1),得到所有的二维相关值,将所有二维相关值重排得到三维相关值,确定三维相关值峰值所在的三个维度的取值,峰值所在的三个维度的取值确定了在伽罗瓦域GF(q)中的Gold-Kasami大集合序列的码字组成及其码相位。
当维度k=q-1时,{ak,k’(n)}即变成Gold-Kasami小序列集合,这时计算{r(n)}和{ak, k’(n)}的互相关函数C(k,k′;τ)采用Kasami小集合序列的FHT快速相关检测方法。
其中:Gold-Kasami大集合序列即是在伽罗瓦域GF(q)中,q=pn,p=2,当n≡2mod4时的Kasami大集合序列,q-1为Gold-Kasami大集合序列码的码长,n为码长的阶数;m序列即最大移位寄存器序列;k、k′和τ是Gold-Kasami大集合序列码字组成及其码相位的三个维度;{r(n)}为接收信号;r*(n)表示r(n)的共轭;{ak,k′(n)}为Gold-Kasami大集合序列码;u是周期为2n-1的m序列,v是对序列u每2n/2+1+1比特抽取1比特组成的周期为2n-1的m序列,w是对序列u每2n/2+1比特抽取1比特组成的周期为2n/2-1的m序列。
有益效果:通过本发明提出的算法,可以用FHT(Fast Hadamard Transform)快速相关检测出伽罗瓦域GF(q)中Gold-Kasami大集合序列的码字组成及其初始码相位。由于FHT中不包含乘法运算,只有加法运算,因此相对于用FFT进行相关检测的算法能大大降低时间复杂度,提高捕获速度。
根据一次q=2n点的FFT需要(q/2)log2q次复数乘法和qlog2q复数加法,一次复数乘法需要三次实数乘法和三次实数加法,一次复数加法需要两次实数加法;Gold-Kasami大集合序列由三个不同的m序列组成,用FFT频域相关法检测出码字组成及其码相位,将乘法运算折算成加法运算,折算方法为:q点FFT中使用的乘法折算成log2q次加法,具体效果如下:
用FFT频域相关法检测出码字组成及其码相位所需的实数加法运算量为:
用本专利提出的V-FHT类和U-FHT类算法完成同样的相关检测运算所需的实数加法运算量为:
用本专利提出的W-FHT类算法完成同样的相关检测运算所需的实数加法运算量为:
综上,可见本发明提出的用FHT快速相关检测出Gold-Kasami大集合序列的码字组成及其初始码相位的算法,相比用FFT频域相关法进行相关检测的算法,V-FHT类和U-FHT类能节省90%以上的运算量,并且节约的运算量随着码长的增加而增大;W-FHT类仅在码字较短的情况下能节约相当大的运算量。
附图说明
图1是算法V-FHT(a)的二维计算过程。
图2是算法V-FHT(a)的相关值三维排列示意图。
图3是算法V-FHT(b)的二维计算过程。
图4是算法V-FHT(b)的相关值三维排列示意图。
图5是算法U-FHT情况I的二维计算过程。
图6是算法U-FHT情况I从二维到三维的相关值排列示意图。
图7是算法U-FHT情况II的二维计算过程。
图8是算法W-FHT(a)的二维计算过程。
图9是算法W-FHT(a)的相关值三维排列示意图。
图10是算法W-FHT(b)的二维计算过程。
图11是算法W-FHT(b)的相关值三维排列示意图。
图12是算法W-FHT(c)从二维到三维的相关值排列示意图。
图13是本发明的方法流程示意图。
图2、4、6、9、11、12中:k0为k轴上的任意值;k′0为k′轴上的任意值。
具体实施方式
文中的字母及符号定义:u、v和w为m序列,u是周期为2n-1的m序列,v是对序列u每2n/2+1+1比特抽取1比特组成的周期为2n-1的m序列,w是对序列u每2n/2+1比特抽取1比特组成的周期为2n/2-1的m序列;N=q-1为Gold-Kasami大集合序列码的码长;n为码长的阶数;T表示向量向左循环移位,其上标表示移位数,下标为长度;G(u,v)是Gold序列;{r(n)}为接收信号;r*(·)表示r(·)的共轭;{ak,k’(n)}是Gold-Kasami大集合序列码;C(k,k′;τ)是互相关函数;k、k′和τ是Gold-Kasami大集合序列码字组成及其码相位组成的三个维度,A为维度k与维度τ的差值,B为维度k′与维度τ的差值;j为将向量分割成块的编号, q=2n;·表示取模,其下标为模值;为模2加法;Cx为任意值符号;A0为A轴上的任意值;B0为B轴上的任意值。
子集为Gold序列的Kasami大集合序列,即为在伽罗瓦域GF(q)中,q=pn,p=2,当n≡2mod 4时由3个m序列组成的周期为N=2n-1的Kasami大集合序列,也称为Gold-Kasami大集合序列,可以表述为
其中,G(u,v)是由u、v两个m序列组成的Gold码,可以表述为:
根据Gold-Kasami大集合序列的定义,{r(n)}与{ak,k’(n)}的互相关函数C(k,k′;τ)可以表示为(q=2n):
观察相关函数C(k,k′;τ)有三个变量:k、k′和τ,用FHT(Fast Hadamard Transform)计算三维相关值时,需要先固定两个变量,使用组成Gold-Kasami大集合序列的一个m序列构成Walsh-Hadamard矩阵,FHT得到一维相关值,再求得二维相关值,最后重排得到三维相关值。确定三维相关值峰值所在的三个维度的取值,峰值所在的三个维度的取值确定了在伽罗瓦域GF(q)中Gold-Kasami大集合序列的码字组成及其码相位。
观察式(1),当k=q-1时,{ak,k’(n)}即变成Kasami小序列集合,计算{r(n)}和{ak,k’(n)}的互相关函数C(k,k′;τ)采用Kasami小集合序列的FHT快速相关检测方法。
当k=0,1,...,q-2时,计算{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ)的3类共7种算法如下:
V-FHT类——使用m序列v生成Walsh-Hadamard矩阵的FHT快速相关检测方法类
a)先固定维度 和维度τ=0,1,...,q-2,得到一维相关值,再得到k和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
令n=n-τ q-1,式(1)可得到
τ=0,1,...,q-2,k=0,1,...,q-2,
使用m序列v生成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’,IEEETransactions on Information Theory,1977
使用这种相关方法,计算复杂度为
算法V-FHT(a)的二维的计算过程如图1所示。
算法V-FHT(a)的三维的相关值排列过程如图2所示。
b):先固定维度 和维度τ=0,1,...,q-2,得到一维相关值,再得到k和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
令τ-k q-1=A(A=0,1,2,...,q-2),即k=τ-A q-1,式(1)可得到
使用m序列v生成Walsh-Hadamard矩阵可得:
使用这种相关方法,计算复杂度为
算法V-FHT(b)二维的计算过程如图3所示。
当 时,用一次Gold码的FHT得到接收信号{r(n)}与{ak,k’(n)}的互相关函数C(k,k′;τ)。
算法V-FHT(b)三维的相关值排列如图4所示。
c):当维度k和维度k′均与维度τ的差值固定时,即τ-k q-1=A(A=0,1,2,...,q-2)或 由式(1)可得到
先固定 和维度τ=0,1,...,q-2的值,得到一维相关值,再得到k和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
使用m序列v生成Walsh-Hadamard矩阵可得:
使用这种相关方法,计算复杂度为
当 时,用一次Gold码的FHT得到接收信号{r(n)}与{ak,k′(n)}的互相关函数C(k,k′;τ)。
观察发现(7)与(5)类似,计算k和τ两个维度的相关值的方法是一样的,只是在第三维相关值的排列上不同。
算法V-FHT(c)二维的计算过程如图3。
算法V-FHT(c)三维的相关值排列过程与方法五类似,如图4,只需要将维度k与维度k′互换即可。
注:以上算法中用序列v生成Walsh-Hadamard矩阵时,序列v不是对序列u每2n/2+1+1比特抽取任1比特组成的周期为2n-1的m序列,因为用这种方法生成的m序列组成的M矩阵不一定能等效为Walsh-Hadamard矩阵。
U-FHT类——使用m序列u生成Walsh-Hadamard矩阵的FHT快速相关检测方法类
当维度k和维度k′均与维度τ的差值固定时,即τ-k q-1=A(A=0,1,2,...,q-2)或 由式(1)可得到
先固定A=0,1,...,q-2和 的值,得到一维相关值,再得到k和τ或k′和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
使用m序列u生成Walsh-Hadamard矩阵可得:
接下来分两种情况:
情况I:如果先计算维度k和维度τ的二维相关值。
算法U-FHT情况I的二维的计算过程如图10所示(方框内是A的值)。
算法U-FHT情况I的三维的相关值排列如图11所示。
情况II:如果先计算维度k′和维度τ的二维相关值。
算法U-FHT情况II的二维的计算过程如图12所示(方框内是B的值)。
算法U-FHT情况II的三维的相关值排列过程与图11类似,之需要将维度k与维度k′的互换即可。
使用U-FHT相关方法,计算复杂度为
当维度 时,用一次Gold码的FHT得到接收信号{r(n)}与{ak,k’(n)}的互相关函数C(k,k′;τ)。
W-FHT类——使用m序列w生成Walsh-Hadamard矩阵的FHT快速相关检测方法类
a):先固定维度k=0,1,...,q-2和维度τ=0,1,...,q-2,得到一维相关值,再得到k′和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
令n=n-τ q-1,式(1)可得到
τ=0,1,...,q-2,k=0,1,...,q-2,
使用m序列w生成Walsh-Hadamard矩阵可得:
(12)
其中
使用这种相关方法,计算复杂度为
算法W-FHT(a)二维的计算过程如图8所示。
算法W-FHT(a)三维的相关值排列如图9所示。
当 时,用一次Gold码的FHT得到接收信号{r(n)}与{ak,k’(n)}的互相关函数C(k,k′;τ)。
b):先固定维度k=0,1,...,q-2和维度τ=0,1,...,q-2,得到一维相关值,再得到k′和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
令 即 式(1)可得到
使用m序列w生成FHT矩阵可得:
(15)
(16)
其中
使用这种相关方法,计算复杂度为
算法W-FHT(b)二维的计算过程如图10所示。
算法W-FHT(b)三维的相关值排列如图11所示。
当 时,用一次Gold码的FHT得到接收信号{r(n)}与{ak,k’(n)}的互相关函数C(k,k′;τ)。
c):当维度k和维度k′均与维度τ的差值固定时,即τ-k q-1=A(A=0,1,2,...,q-2)或 由式(1)可得到
先固定A=0,1,...,q-2和维度τ=0,1,...,q-2的值,得到一维相关值,再得到k′和τ两个维度的相关值,最后得到k、k′和τ三个维度的相关值:
使用m序列w生成Walsh-Hadamard矩阵可得:
使用这种相关方法,计算复杂度为
观察发现式(19)与(15)类似,计算维度k′和维度τ二维相关值的方法是一样的,只是在第三维相关值的排列上不同。
算法W-FHT(c)二维的计算过程如图10。
算法W-FHT(c)三维的相关值排列如图12所示。
本发明的方法流程示意图如图13所示。
当 时,用一次Gold码的FHT得到接收信号{r(n)}与{ak,k’(n)}的互相关函数C(k,k′;τ)。
表1是算法类V-FHT、U-FHT与FFT频域相关法检测Gold-Kasami大集合序列时的运算量比较。
表2是算法类W-FHT与FFT频域相关法检测Gold-Kasami大集合序列时的运算量比较。
表1
表2
Claims (3)
1.一种对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法,其特征在于:当维度k=0,1,...,q-2时,计算{r(n)}与{ak,k(n)}——Gold-Kasami大集合序列码的互相关函数C(k,k′;τ)采用的相关步骤如下:
步骤(1):先固定两个维度,使用组成Gold-Kasami大集合序列的一个m序列构成Walsh-Hadamard矩阵,FHT得到一维相关值;
步骤(2):将步骤1中固定的两个维度中的一个维度取所有可能值,重复步骤(1),得到二维相关值;
步骤(3):将步骤1中固定的两个维度中的另一个维度取所有可能值,重复步骤(1),得到所有的二维相关值,将所有二维相关值重排得到三维相关值,确定三维相关值峰值所在的三个维度的取值,峰值所在的三个维度的取值确定了在伽罗瓦域GF(q)中的Gold-Kasami大集合序列的码字组成及其码相位;
当维度k=q-1时,{ak,k,(n)}即变成Gold-Kasami小序列集合,这时计算{r(n)}和{ak,k,(n)}的互相关函数C(k,k′;τ)采用Kasami小集合序列的FHT快速相关检测方法,其中{r(n)}是接收信号,{ak,k’(n)}是Gold-Kasami大集合序列码;C(k,k′;τ)是互相关函数,Gold-Kasami大集合序列即是在伽罗瓦域GF(q)中,q=pn,p=2,当n≡2mod4时的Kasami大集合序列,q-1为Gold-Kasami大集合序列码的码长,n为码长的阶数;m序列即最大移位寄存器序列;k、k′和τ是Gold-Kasami大集合序列码字组成及其码相位的三个维度。
2.根据权利要求1所述的关于对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法,其特征在于:{r(n)}和{ak,k’(n)}——Gold-Kasami大集合序列码的互相关函数C(k,k′;τ)表示为(q=2n):
r*(n)表示r(n)的共轭;u是周期为2n-1的m序列;v是对序列u每2n/2+1+1比特抽取1比特组成的周期为2n-1的m序列;w是对序列u每2n/2+1比特抽取1比特组成的周期为2n/2-1的m序列。
3.根据权利要求1所述的关于对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法,其特征在于:Gold-Kasami大集合序列由三个不同的m序列组成,记做{ak,k(n)},此时序列的码字组成及其码相位的不确定度为三维不确定度,由于FHT算法需要一个m序列生成Walsh-Hadamard矩阵,因此所有可能的快速相关检测的方法有3类,共7种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100341575A CN101635581B (zh) | 2009-09-01 | 2009-09-01 | 对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100341575A CN101635581B (zh) | 2009-09-01 | 2009-09-01 | 对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101635581A CN101635581A (zh) | 2010-01-27 |
CN101635581B true CN101635581B (zh) | 2013-07-03 |
Family
ID=41594657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100341575A Expired - Fee Related CN101635581B (zh) | 2009-09-01 | 2009-09-01 | 对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101635581B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103488612B (zh) * | 2013-09-17 | 2016-08-24 | 浙江工业大学 | 一种应用于数字滤波的沃什-新梅森数快速变换方法 |
CN103532892A (zh) * | 2013-09-27 | 2014-01-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种海洋石油勘探声学定位应答器通信编码方法 |
CN108507596B (zh) * | 2018-03-05 | 2020-06-30 | 中国科学院上海光学精密机械研究所 | 二维gold矩阵绝对位置编码方法及其解码方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1533045A (zh) * | 2000-06-12 | 2004-09-29 | 三星电子株式会社 | 用于编码和解码传输格式组合指示符的设备和方法 |
-
2009
- 2009-09-01 CN CN2009100341575A patent/CN101635581B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1533045A (zh) * | 2000-06-12 | 2004-09-29 | 三星电子株式会社 | 用于编码和解码传输格式组合指示符的设备和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101635581A (zh) | 2010-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100553168C (zh) | 一种长周期扩频码的频率折叠时频并行搜索方法 | |
Kong | A deterministic compressed GNSS acquisition technique | |
CN102223226B (zh) | Sc-fdm***时频同步装置及同步方法 | |
CN104977504A (zh) | 一种电缆故障在线检测及定位装置 | |
CN100553169C (zh) | 基于频率折叠的短周期扩频码时频并行搜索方法 | |
CN101635581B (zh) | 对子集为Gold序列的Kasami大集合序列的FHT快速相关检测方法 | |
CN101702827A (zh) | Toa定位的估计方法及基于该估计方法的精确定位方法 | |
CN101056140A (zh) | 用于同步的前置码 | |
CN101063717B (zh) | 一种用fft进行gps信号捕获的方法 | |
CN106817326B (zh) | 多用户周期长短码直扩信号的伪码盲估计方法 | |
CN103926604A (zh) | 基于重叠差分循环相干积分的弱信号捕获方法 | |
CN101330318B (zh) | 一种下行同步***中辅同步信道序列的加扰和解扰方法 | |
CN102801459A (zh) | 多模卫星信号接收装置及其方法 | |
CN103885073B (zh) | 基于数字水印和压缩感知的导航信号通信方法 | |
CN101635936B (zh) | 对子集为Gold-like序列的Kasami大集合序列的FHT快速相关检测方法 | |
CN103969664A (zh) | 一种基于差分循环相干积分的gps弱信号捕获方法 | |
CN104280748B (zh) | 先累加后相关的批处理与差分相干的导航信号捕获方法 | |
Alvarez et al. | Real-time implementation of an efficient correlator for complementary sets of four sequences applied to ultrasonic pulse compression systems | |
CN101762803B (zh) | 一种应用于全球导航卫星***接收机的差分相关器 | |
CN114114335A (zh) | 一种弱信号gnss接收机快速捕获方法 | |
CN102324996B (zh) | 多值伪随机序列生成方法及装置 | |
CN103124438B (zh) | 一种随机接入前导的生成方法及装置 | |
CN102035567B (zh) | 可变长度数字扩频信号快速相关处理方法 | |
KR101812325B1 (ko) | 메모리 어레이에 기반한 p 코드 획득 방법 및 메모리 어레이에 기반하여 p 코드를 획득하는 수신기 | |
Zhang et al. | Fast correlation for gold large sets of Kasami sequences |
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: 20130703 Termination date: 20160901 |
|
CF01 | Termination of patent right due to non-payment of annual fee |