CN107290761B - 一种并行码相位搜索装置及实现并行码相位搜索的方法 - Google Patents
一种并行码相位搜索装置及实现并行码相位搜索的方法 Download PDFInfo
- Publication number
- CN107290761B CN107290761B CN201610228762.6A CN201610228762A CN107290761B CN 107290761 B CN107290761 B CN 107290761B CN 201610228762 A CN201610228762 A CN 201610228762A CN 107290761 B CN107290761 B CN 107290761B
- Authority
- CN
- China
- Prior art keywords
- assembly
- result
- data
- preset
- judgment
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
Abstract
一种并行码相位搜索装置及实现并行码相位搜索的方法，包括：根据采样率将混频复数数据和本地码分别组装为预设质数的幂次方尺寸的第一组装数据和第二组装数据；傅里叶变换第一组装数据获得第一组装变换结果；对第二组装数据进行傅里叶变换和复数共轭处理获得组装共轭结果；将第一组装变换结果和组装共轭结果相乘获得组装乘积，傅里叶反变换组装乘积获得组装反变换结果；取模处理组装反变换结果获得取模结果后，按照预设策略进行相位判决获得进行相位搜索的判决输出。本发明实施例技术方案，按照预设质数的幂次方尺寸对混频复数数据和本地码进行组装，按照预设策略进行相位判决输出，降低了并行码相位搜索的复杂度，简化了傅里叶变换单元的设置。
Description
技术领域
本发明涉及信号处理技术，尤指一种并行码相位搜索装置及实现并行码相位搜索的方法。
背景技术
全球卫星导航系统(GNSS，Global Navigation Satellite System)在人们的日常生活中发挥着越来越不可替代的重要作用，尤其在导航、定时、测绘等领域得到越来越多的应用。目前，全球卫星导航系统主要包括美国的全球定位系统(GPS)、中国的北斗(BD)系统、俄罗斯的全球导航卫星定位系统(GLONASS，GLObal NAvigation Satellite System)，以及欧洲的伽利略(Galileo)系统。在中国和亚太地区，GPS和北斗系统应用较为广泛；而在俄罗斯，以GPS和GLONASS应用较多。由于伽利略系统远未成熟，尚不可提供正式服务。利用全球卫星导航系统进行定位、定时等业务时，首先需要捕获到至少四颗可见卫星的无线信号，通过捕获的无线信号实现三维搜索算法，三维搜索算法包括：卫星伪码、码相位和多普勒频移。
常见的线性搜索方法按照一维、二维、三维的顺序逐步搜索，耗时太长。为了较少耗时，并行搜索方法得到了发展，例如、并行码相位搜索算法，可明显减小用时。图1为现有的并行码相位搜索电路原理图，如图1所示，当数字中频输入信号分别与同相(I)支路和正交(Q)支路上第一频带的复制正弦和复制余弦载波信号混频后，以同相和正交混频结果的复数形式的混频复数数据通过第一傅里叶变换单元进行傅里叶变换获得傅里叶变换结果；将傅里叶变换结果与本地码共轭结果(复制粗捕获(C/A)码发生器产生的本地码，将复制的本地码经第二傅里叶变换单元和复数共轭单元处理获得本地码共轭结果)通过乘法器相乘，将通过乘法器相乘获得的乘积经傅里叶反变换单元进行处理获得在时域内的相关结果，将获得的时域内的相关结果通过取模单元进行取模后检测判断并行码相位信号是否存在。在完成了对当前频带的搜索与检测后，接收机接着让载波数控振荡器(NCO，NumericalControl Oscillator)进行第二频带正弦载波和余弦载波复制，然后类似地完成对其他频带的搜索与检测，这里，第一频带、第二频带、及其他频带的数值为并行码相位搜索过程中使用的遍历频带，为本领域技术人员的公知常识。在对同一个卫星信号不同频带内的搜索过程中，复制C/A码的相位可保持不变，相应地其傅里叶变换及其共轭值也保持不变。当搜索另一个卫星信号时，接收机可让C/A码发生器复制相应的另一个C/A码，然后重复上述在各个频带中的信号搜索过程。
并行码相位搜索算法实际上利用傅里叶变换这种数字信号处理技术来替代数字相关器的相关运算，下面证明两者的等价性。记两个长度均为N点的周期性序列为x(n)和y(n)，其相关值序列为z(n)，n＝0,1,…,N-1，可表示为公式(1)：
对上述相关值序列z(n)进行离散傅里叶变换，得z(n)的离散傅里叶变换Z(k)如公式(2)所示：
将公式(2)带入公式(1)得到公式(3)
对公式(3)做进一步变换得：
最终得：
其中，X(k)与Y(k)分别为x(n)与y(n)的离散傅里叶变换，
代表复数Y(k)的共轭。上式表明：两个序列x(n)与y(n)在时域内做相关运算，相当于它们的离散傅里叶变换X(k)与Y(k)(确切地讲是Y(k)的共轭
)在频域内做乘积运算。于是倒过来，乘积
的离散傅里叶反变换正好是接收机需要进行检测的在各个码相位处的相关值z(n)。一旦接收机通过傅里叶反变换计算得到相关值z(n)，那么接下来的信号检测就同线性搜索捕获法一样，即找出在所有搜索单元中自相关幅值|z(n)|的峰值，并将该峰值与捕获门限值相比较。若峰值超过捕获门限值，则接收机捕获到了信号，并且也从中获得了该信号的频率和码相位两个参数值。
对于GPS导航系统来说，上述方案中的y(n)序列一般是C/A码序列，而C/A码序列的序列长度是N＝1023个码片，时间长度为1ms；即傅里叶变换和傅里叶反变换的尺寸通常为1023或其倍数；而进行傅里叶变换时需要通过设置傅里叶变换单元中的蝶形装置进行；以1023个码片为例，最少设置出尺寸为3、11和31三个蝶形装置才可以进行；大尺寸的蝶形装置也造成傅里叶变换处理过程更为复杂(尺寸越小的蝶形装置，其复杂度越低)；而设置多种尺寸的蝶形装置使傅里叶变换单元设置复杂。
发明内容
为了解决上述技术问题，本发明提供一种并行码相位搜索的装置及实现并行码相位搜索的方法，能够降低傅里叶变换的复杂程度。
为了达到本发明目的，本发明提供了一种并行码相位搜索装置，包括：第一组装单元、第一傅里叶变换单元、第二组装单元、第二傅里叶变换单元、复数共轭单元、乘法器、傅里叶反变换单元、取模单元和判决单元；其中，
第一组装单元，用于将发往并行码相位搜索电路第一傅里叶变换单元的混频复数数据根据采样率组装为预设质数的幂次方尺寸的第一组装数据；
第一傅里叶变换单元用于，接收第一组装数据，对第一组装数据进行傅里叶变换后输出第一组装变换结果，将第一组装变换结果发往乘法器；
第二组装单元用于，将并行码相位搜索电路中粗捕获CA码发生器产生的本地码根据采样率组装为预设质数的幂次方尺寸的第二组装数据；
第二傅里叶变换单元用于，接收第二组装数据，将第二组装数据进行傅里叶变换处理后获得第二组装变换结果，并发送第二组装变换结果至复数共轭单元；
复数共轭单元用于，接收第二组装变换结果，对第二组装变换结果进行复数共轭处理获得组装共轭结果，将组装共轭结果发往乘法器；
乘法器用于，接收第一组装变换结果和组装共轭结果，将第一组装变换结果和组装共轭结果相乘后获得组装乘积；
傅里叶反变换单元用于，对组装乘积进行傅里叶反变换处理，获得组装反变换结果；
取模单元用于，对获得的组装反变换结果进行取模处理，获得取模结果；
判决单元用于，对获得的取模结果按照预设策略进行相位判决后，获得进行相位搜索的判决输出。
可选的，所述预设质数为2。
可选的，所述第一组装单元具体用于，
确定所述预设质数的幂次方尺寸与所述混频复数数据的尺寸的第一组装差值；
在所述混频复数数据的末端添加确定的第一组装差值个随机数以组装完成所述第一组装数据。
可选的，所述第二组装单元具体用于，
确定所述预设质数的幂次方尺寸与所述本地码的尺寸的第二组装差值；
在所述本地码的末端添加确定的第二组装差值个随机数以组装完成所述第二组装数据。
可选的，所述判决单元具体用于，
获取所述取模结果中的序列最大值；
获取所述取模结果中大于预设的第一门限的索引序列；
所述索引序列的长度大于1时，选择索引序列中的最小值作为判决结果；
所述索引序列的长度等于1时，
如果所述序列最大值大于预设的第二门限，且索引序列大于预设的索引门限，将所述索引序列减去1后获得的索引差值作为判决结果；如果所述序列最大值小于或等于第二门限、或所述索引序列小于或等于预设的索引门限时，将所述索引序列作为判决结果；
将所述混频复数数据的尺寸减去获得的所述判决结果获得判决差值后，将获得的判决差值对混频复数数据的尺寸进行求余，将求余结果作为所述判决输出；
所述第二门限大于第一门限。
另一方面，本发明实施例方法还提供一种实现并行码相位搜索的方法，包括：
根据采样率将发往并行码相位搜索电路第一傅里叶变换单元的混频复数数据组装为预设质数的幂次方尺寸的第一组装数据；
对第一组装数据进行傅里叶变换获得第一组装变换结果；
根据采样率将并行码相位搜索电路中粗捕获CA码发生器产生的本地码组装为预设质数的幂次方尺寸的第二组装数据；
对第二组装数据进行傅里叶变换处理获得第二组装变换结果；
对第二组装变换结果进行复数共轭处理获得组装共轭结果；
将第一组装变换结果和组装共轭结果相乘后获得组装乘积；
对组装乘积进行傅里叶反变换处理，获得组装反变换结果；
对获得的组装反变换结果进行取模处理，获得取模结果；
对获得的取模结果按照预设策略进行相位判决后，获得进行相位搜索的判决输出。
可选的，组装为预设质数的幂次方尺寸的第一组装数据包括：
确定所述预设质数的幂次方尺寸与所述混频复数数据的尺寸的第一组装差值；
在所述混频复数数据的末端添加确定的第一组装差值个随机数以组装完成所述第一组装数据。
可选的，组装为预设质数的幂次方尺寸的第二组装数据包括：
确定所述预设质数的幂次方尺寸与所述本地码的尺寸的第二组装差值；
在所述本地码的末端添加确定的第二组装差值个随机数以组装完成所述第二组装数据。
可选的，按照预设策略进行相位判决包括：
获取所述取模结果中的序列最大值；
获取所述取模结果中大于预设的第一门限的索引序列；
所述索引序列的长度大于1时，选择索引序列中的最小值作为判决结果；
所述索引序列的长度等于1时，
如果所述序列最大值大于预设的第二门限，且索引序列大于预设的索引门限，将所述索引序列减去1后获得的索引差值作为判决结果；如果所述序列最大值小于或等于第二门限、或所述索引序列小于或等于预设的索引门限时，将所述索引序列作为判决结果；
将所述混频复数数据的尺寸减去获得的所述判决结果获得判决差值后，将获得的判决差值对混频复数数据的尺寸进行求余，将求余结果作为所述判决输出；
所述第二门限大于第一门限。
可选的，预设质数为2。
可选的，组装完成所述第一组装数据的随机数和组装完成所述第二组装数据的随机数相同。
可选的，随机数包括全为0的数。
与相关技术相比，本申请技术方案包括：根据采样率将发往并行码相位搜索电路第一傅里叶变换单元的混频复数数据组装为预设质数的幂次方尺寸的第一组装数据；对第一组装数据进行傅里叶变换获得第一组装变换结果；根据采样率将并行码相位搜索电路中粗捕获(CA)码发生器产生的本地码组装为预设质数的幂次方尺寸的第二组装数据；对第二组装数据进行傅里叶变换处理获得第二组装变换结果；对第二组装变换结果进行复数共轭处理获得组装共轭结果；将第一组装变换结果和组装共轭结果相乘后获得组装乘积；对组装乘积进行傅里叶反变换处理，获得组装反变换结果；对获得的组装反变换结果进行取模处理，获得取模结果；对获得的取模结果按照预设策略进行相位判决后，获得进行相位搜索的判决输出。本发明实施例技术方案，按照预设质数的幂次方尺寸对混频复数数据和本地码进行组装，按照预设策略进行相位判决输出，降低了并行码相位搜索的复杂度，简化了傅里叶变换单元的设置。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为现有的并行码相位搜索电路原理图；
图2为本发明实施例并行码相位搜索装置的结构框图；
图3为本发明实施例实现并行码相位搜索的方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下文中将结合附图对本发明的实施例进行详细说明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
图2为本发明实施例并行码相位搜索装置的结构框图，如图2所示，包括：第一组装单元、第一傅里叶变换单元、第二组装单元、第二傅里叶变换单元、复数共轭单元、乘法器、傅里叶反变换单元、取模单元和判决单元；其中，
第一组装单元，用于将发往并行码相位搜索电路第一傅里叶变换单元的混频复数数据根据采样率组装为预设质数的幂次方尺寸的第一组装数据；
需要说明的是，这里的采样率包括并行码相位搜索电路中的采样率。
可选的，本发明实施例预设质数可以为2。
需要说明的是，预设质数为2时，第一傅里叶变换单元中的蝶形装置尺寸为2，此时第一傅里叶变换单元进行傅里叶变变换处理时，实质是快速傅里叶变换。
可选的，第一组装单元具体用于，
确定预设质数的幂次方尺寸与混频复数数据的尺寸的第一组装差值；
在混频复数数据的末端添加确定的第一组装差值个随机数以组装完成第一组装数据。
需要说明的是，随机数可以是任意的整数、例如、0、1或其他数。另外、预设质数的幂次方尺寸与混频复数数据的尺寸是数值相近，例如、尺寸长度为1023时，组装为尺寸长度为1024，而不会将1023封装为尺寸长度为2048的数据。
第一傅里叶变换单元用于，接收第一组装数据，对第一组装数据进行傅里叶变换后输出第一组装变换结果，将第一组装变换结果发往乘法器；
第二组装单元用于，将并行码相位搜索电路中粗捕获CA码发生器产生的本地码根据采样率组装为预设质数的幂次方尺寸的第二组装数据；
需要说明的是，这里的采样率包括并行码相位搜索电路中的采样率。
可选的，第二组装单元具体用于，
确定预设质数的幂次方尺寸与本地码的尺寸的第二组装差值；
在本地码的末端添加确定的第二组装差值个随机数以组装完成第二组装数据。
可选的，本发明实施例搜索装置中，组装完成所述第一组装数据的随机数和组装完成所述第二组装数据的随机数可以相同。
可选的，随机数包括全为0的数。
第二傅里叶变换单元用于，接收第二组装数据，将第二组装数据进行傅里叶变换处理后获得第二组装变换结果，并发送第二组装变换结果至复数共轭单元；
需要说明的是，预设质数为2时，第二傅里叶变换单元中的蝶形装置尺寸为2，此时第二傅里叶变换单元进行傅里叶变变换处理时，实质是快速傅里叶变换。
复数共轭单元用于，接收第二组装变换结果，对第二组装变换结果进行复数共轭处理获得组装共轭结果，将组装共轭结果发往乘法器；
乘法器用于，接收第一组装变换结果和组装共轭结果，将第一组装变换结果和组装共轭结果相乘后获得组装乘积；
傅里叶反变换单元用于，对组装乘积进行傅里叶反变换处理，获得组装反变换结果；
取模单元用于，对获得的组装反变换结果进行取模处理，获得取模结果；
判决单元用于，对获得的取模结果按照预设策略进行相位判决后，获得进行相位搜索的判决输出。
可选的，判决单元具体用于，
获取取模结果中的序列最大值；
获取取模结果中大于预设的第一门限的索引序列；
索引序列的长度大于1时，选择索引序列中的最小值作为判决结果；
索引序列的长度等于1时，
如果序列最大值大于预设的第二门限，且索引序列大于预设的索引门限，将索引序列减去1后获得的索引差值作为判决结果；如果序列最大值小于或等于第二门限、或索引序列小于或等于预设的索引门限时，将索引序列作为判决结果；
将混频复数数据的尺寸减去获得的判决结果获得判决差值后，将获得的判决差值对混频复数数据的尺寸进行求余，将求余结果作为判决输出；
第二门限大于第一门限。
需要说明的是，通过预设策略对取模结果的多个峰值序列进行判断获得进行相位搜索的信号。另外、第一门限、第二门限、索引门限主要根据系统信噪比等信号属性通过仿真获得的数值；一般的，混频复数数据为1023个码片时，第一门限可以取值295～305之间的数值，第二门限可以取值495～505之间的数值、索引门限可以取值395～405之间的数值；仿真方式为本领域技术人员的惯用技术手段，对不同种类的信号，获得的第一门限、第二门限、索引门限可以不同。
本发明实施例装置，按照预设质数的幂次方尺寸对混频复数数据和本地码进行组装，按照预设策略进行相位判决输出，降低了并行码相位搜索的复杂度，简化了傅里叶变换单元的设置。
图3为本发明实施例实现并行码相位搜索的方法的流程图，如图3所示，包括：
步骤300、根据采样率将发往并行码相位搜索电路第一傅里叶变换单元的混频复数数据组装为预设质数的幂次方尺寸的第一组装数据；
可选的，本发明实施例中，预设质数可以为2。
可选的，组装为预设质数的幂次方尺寸的第一组装数据包括：
确定预设质数的幂次方尺寸与混频复数数据的尺寸的第一组装差值；
在混频复数数据的末端添加确定的第一组装差值个随机数以组装完成第一组装数据。
步骤301、对第一组装数据进行傅里叶变换获得第一组装变换结果；
步骤302、根据采样率将并行码相位搜索电路中粗捕获(CA)码发生器产生的本地码组装为预设质数的幂次方尺寸的第二组装数据；
可选的，组装为预设质数的幂次方尺寸的第二组装数据包括：
确定预设质数的幂次方尺寸与本地码的尺寸的第二组装差值；
在本地码的末端添加确定的第二组装差值个随机数以组装完成第二组装数据。
可选的，本发明实施例搜索方法中，组装完成第一组装数据的随机数和组装完成第二组装数据的随机数可以相同。
可选的，随机数包括全为0的数。
步骤303、对第二组装数据进行傅里叶变换处理获得第二组装变换结果；
步骤304、对第二组装变换结果进行复数共轭处理获得组装共轭结果；
需要说明的是，步骤300、301获得第一组装变换结果与步骤302～步骤304获得组装共轭结果并不存在时间先后顺序。
步骤305、将第一组装变换结果和组装共轭结果相乘后获得组装乘积；
步骤306、对组装乘积进行傅里叶反变换处理，获得组装反变换结果；
步骤307、对获得的组装反变换结果进行取模处理，获得取模结果；
步骤308、对获得的取模结果按照预设策略进行相位判决后，获得进行相位搜索的判决输出。
可选的，按照预设策略进行相位判决包括：
获取取模结果中的序列最大值；
获取取模结果中大于预设的第一门限的索引序列；
索引序列的长度大于1时，选择索引序列中的最小值作为判决结果；
索引序列的长度等于1时，
如果序列最大值大于预设的第二门限，且索引序列大于预设的索引门限，将索引序列减去1后获得的索引差值作为判决结果；如果序列最大值小于或等于第二门限、或索引序列小于或等于预设的索引门限时，将索引序列作为判决结果；
将混频复数数据的尺寸减去获得的判决结果获得判决差值后，将获得的判决差值对混频复数数据的尺寸进行求余，将求余结果作为判决输出；
第二门限大于第一门限。
本发明实施例方法，按照预设质数的幂次方尺寸对混频复数数据和本地码进行组装，按照预设策略进行相位判决输出，降低了并行码相位搜索的复杂度，简化了傅里叶变换单元的设置。
以下通过应用示例对本发明方法进行清楚详细的说明，应用示例仅用于陈述本发明，并不用于限定本发明方法的保护范围。
应用示例
本应用示例根据混频复数数据的采样率将1毫秒(ms)的混频复数数据组装为2的幂次方尺寸的第一组装数据；将本地码组装为2的幂次方尺寸的第二组装数据；
假设采样率为1.023兆赫兹(MHz)，那么1ms的数据包含1023个点数；组装处理后，第一组装数据的尺寸为1024长度。组装是可以采用在数据末端添加随机数的方式实现，例如在末端添加数据0。
第一傅里叶变换单元、第二傅里叶变换单元、复数共轭单元、乘法器、取模单元和并行码相位搜索电路中的单元工作原理相同。
假设混频复数数据的采样率为1.023MHz，取1ms的数据段共1023个数据，记为x(0)、x(1)、…、x(1022)，组装为第一组装数据的处理过程包括：
在x(0)、x(1)、…、x(1022)后面添加一个值为0的数据，即新序列x2(n)为：x2(n)＝x(0)、x(1)、…、x(1022)、0，共1024点。
对x2(n)做傅里叶变换，记为X2(k)。
C/A码发生器产生1ms的数据，共1023个，记为y(0)、y(1)、…、y(1022)，同样经组装处理，获得的新序列y2(n)为
y2(n)＝y(0)、y(1)、…、y(1022)、0，共1024点。
同样，对y2(n)序列做傅里叶变换和复数共轭处理，记为
傅里叶反变换单元的输入为：
输出为z2(n)序列，n＝0、1、2、…、1023。
对z2(n)序列经过取模单元处理，即进行|z2(n)|计算获得取模结果；
获取取模结果中的序列最大值；
获取取模结果中大于预设的第一门限的索引序列；
索引序列的长度大于1时，选择索引序列中的最小值作为判决结果；
索引序列的长度等于1时，
如果序列最大值大于预设的第二门限，且索引序列大于预设的索引门限，将索引序列减去1后获得的索引差值作为判决结果；如果序列最大值小于或等于第二门限、或索引序列小于或等于预设的索引门限时，将索引序列作为判决结果；这里，第二门限大于第一门限；
将混频复数数据的尺寸减去获得的判决结果获得判决差值后，将获得的判决差值对混频复数数据的尺寸进行求余，将求余结果作为判决输出；
本应用示例中，通过仿真获得的第一门限等于300、第二门限等于500、索引门限等于400。
虽然本发明所揭露的实施方式如上，但所述的内容仅为便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式及细节上进行任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。
Claims (12)
1.一种并行码相位搜索装置，其特征在于，包括：第一组装单元、第一傅里叶变换单元、第二组装单元、第二傅里叶变换单元、复数共轭单元、乘法器、傅里叶反变换单元、取模单元和判决单元；其中，
第一组装单元，用于将发往并行码相位搜索电路第一傅里叶变换单元的混频复数数据根据采样率组装为预设质数的幂次方尺寸的第一组装数据；
第一傅里叶变换单元用于，接收第一组装数据，对第一组装数据进行傅里叶变换后输出第一组装变换结果，将第一组装变换结果发往乘法器；
第二组装单元用于，将并行码相位搜索电路中粗捕获CA码发生器产生的本地码根据采样率组装为预设质数的幂次方尺寸的第二组装数据；
第二傅里叶变换单元用于，接收第二组装数据，将第二组装数据进行傅里叶变换处理后获得第二组装变换结果，并发送第二组装变换结果至复数共轭单元；
复数共轭单元用于，接收第二组装变换结果，对第二组装变换结果进行复数共轭处理获得组装共轭结果，将组装共轭结果发往乘法器；
乘法器用于，接收第一组装变换结果和组装共轭结果，将第一组装变换结果和组装共轭结果相乘后获得组装乘积；
傅里叶反变换单元用于，对组装乘积进行傅里叶反变换处理，获得组装反变换结果；
取模单元用于，对获得的组装反变换结果进行取模处理，获得取模结果；
判决单元用于，对获得的取模结果按照预设策略进行相位判决后，获得进行相位搜索的判决输出。
2.根据权利要求1所述的并行码相位搜索装置，其特征在于，所述预设质数为2。
3.根据权利要求1或2所述的并行码相位搜索装置，其特征在于，所述第一组装单元具体用于，
确定所述预设质数的幂次方尺寸与所述混频复数数据的尺寸的第一组装差值；
在所述混频复数数据的末端添加确定的第一组装差值个随机数以组装完成所述第一组装数据。
4.根据权利要求1或2所述的并行码相位搜索装置，其特征在于，所述第二组装单元具体用于，
确定所述预设质数的幂次方尺寸与所述本地码的尺寸的第二组装差值；
在所述本地码的末端添加确定的第二组装差值个随机数以组装完成所述第二组装数据。
5.根据权利要求1或2所述的并行码相位搜索装置，其特征在于，所述判决单元具体用于，
获取所述取模结果中的序列最大值；
获取所述取模结果中大于预设的第一门限的索引序列；
所述索引序列的长度大于1时，选择索引序列中的最小值作为判决结果；
所述索引序列的长度等于1时，
如果所述序列最大值大于预设的第二门限，且索引序列大于预设的索引门限，将所述索引序列减去1后获得的索引差值作为判决结果；如果所述序列最大值小于或等于第二门限、或所述索引序列小于或等于预设的索引门限时，将所述索引序列作为判决结果；
将所述混频复数数据的尺寸减去获得的所述判决结果获得判决差值后，将获得的判决差值对混频复数数据的尺寸进行求余，将求余结果作为所述判决输出；
所述第二门限大于第一门限。
6.一种实现并行码相位搜索的方法，其特征在于，包括：
根据采样率将发往并行码相位搜索电路第一傅里叶变换单元的混频复数数据组装为预设质数的幂次方尺寸的第一组装数据；
对第一组装数据进行傅里叶变换获得第一组装变换结果；
根据采样率将并行码相位搜索电路中粗捕获CA码发生器产生的本地码组装为预设质数的幂次方尺寸的第二组装数据；
对第二组装数据进行傅里叶变换处理获得第二组装变换结果；
对第二组装变换结果进行复数共轭处理获得组装共轭结果；
将第一组装变换结果和组装共轭结果相乘后获得组装乘积；
对组装乘积进行傅里叶反变换处理，获得组装反变换结果；
对获得的组装反变换结果进行取模处理，获得取模结果；
对获得的取模结果按照预设策略进行相位判决后，获得进行相位搜索的判决输出。
7.根据权利要求6所述的方法，其特征在于，所述组装为预设质数的幂次方尺寸的第一组装数据包括：
确定所述预设质数的幂次方尺寸与所述混频复数数据的尺寸的第一组装差值；
在所述混频复数数据的末端添加确定的第一组装差值个随机数以组装完成所述第一组装数据。
8.根据权利要求6所述的方法，其特征在于，所述组装为预设质数的幂次方尺寸的第二组装数据包括：
确定所述预设质数的幂次方尺寸与所述本地码的尺寸的第二组装差值；
在所述本地码的末端添加确定的第二组装差值个随机数以组装完成所述第二组装数据。
9.根据权利要求6所述的方法，其特征在于，所述按照预设策略进行相位判决包括：
获取所述取模结果中的序列最大值；
获取所述取模结果中大于预设的第一门限的索引序列；
所述索引序列的长度大于1时，选择索引序列中的最小值作为判决结果；
所述索引序列的长度等于1时，
如果所述序列最大值大于预设的第二门限，且索引序列大于预设的索引门限，将所述索引序列减去1后获得的索引差值作为判决结果；如果所述序列最大值小于或等于第二门限、或所述索引序列小于或等于预设的索引门限时，将所述索引序列作为判决结果；
将所述混频复数数据的尺寸减去获得的所述判决结果获得判决差值后，将获得的判决差值对混频复数数据的尺寸进行求余，将求余结果作为所述判决输出；
所述第二门限大于第一门限。
10.根据权利要求6～9中任一项所述的方法，其特征在于，所述预设质数为2。
11.根据权利要求7或8所述的方法，其特征在于，所述组装完成所述第一组装数据的随机数和组装完成所述第二组装数据的随机数相同。
12.根据权利要求7或8所述的方法，其特征在于，所述随机数包括全为0的数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610228762.6A CN107290761B (zh) | 2016-04-13 | 2016-04-13 | 一种并行码相位搜索装置及实现并行码相位搜索的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610228762.6A CN107290761B (zh) | 2016-04-13 | 2016-04-13 | 一种并行码相位搜索装置及实现并行码相位搜索的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107290761A CN107290761A (zh) | 2017-10-24 |
| CN107290761B true CN107290761B (zh) | 2020-07-28 |
Family
ID=60095921
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610228762.6A Active CN107290761B (zh) | 2016-04-13 | 2016-04-13 | 一种并行码相位搜索装置及实现并行码相位搜索的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107290761B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101072045A (zh) * | 2007-06-15 | 2007-11-14 | 清华大学 | 基于频率折叠的短周期扩频码时频并行搜索方法 |
| CN101464508A (zh) * | 2008-12-19 | 2009-06-24 | 苏州莱迪斯特电子有限公司 | 一种gps的c/a码信号的捕获方法 |
| CN102495416A (zh) * | 2011-10-31 | 2012-06-13 | 北京中微星通电子有限公司 | 卫星导航系统信号快速捕获的方法及装置 |
| CN103499825A (zh) * | 2013-09-24 | 2014-01-08 | 航天恒星科技有限公司 | 一种高轨航天器bds弱信号快速捕获方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7702002B2 (en) * | 2004-01-28 | 2010-04-20 | Qualcomm Incorporated | Rapid acquisition methods and apparatus for GPS signals |
| US8094702B2 (en) * | 2008-04-28 | 2012-01-10 | Qualcomm Incorporated | System and/or method for detecting multi-tone jamming |
-
2016
- 2016-04-13 CN CN201610228762.6A patent/CN107290761B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101072045A (zh) * | 2007-06-15 | 2007-11-14 | 清华大学 | 基于频率折叠的短周期扩频码时频并行搜索方法 |
| CN101464508A (zh) * | 2008-12-19 | 2009-06-24 | 苏州莱迪斯特电子有限公司 | 一种gps的c/a码信号的捕获方法 |
| CN102495416A (zh) * | 2011-10-31 | 2012-06-13 | 北京中微星通电子有限公司 | 卫星导航系统信号快速捕获的方法及装置 |
| CN103499825A (zh) * | 2013-09-24 | 2014-01-08 | 航天恒星科技有限公司 | 一种高轨航天器bds弱信号快速捕获方法 |
Non-Patent Citations (1)
| Title |
|---|
| "北斗B1I信号的捕获算法";王丽黎 等;《计算机系统应用》;20160331;第25卷(第3期);第194-198页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107290761A (zh) | 2017-10-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4757425B2 (ja) | コード位相トラッキング方法及び受信機 | |
| RU2510134C2 (ru) | Способ и система для обнаружения сигнала с расширенным спектром | |
| US8462616B2 (en) | Apparatus and method for estimating a frequency shift and a time shift | |
| EP1712010A1 (en) | Rapid acquisition methods and apparatus for gps signals background | |
| CN104111465B (zh) | 连续波干扰信号的估计器、估计方法、消除器及消除方法 | |
| Borio | M-sequence and secondary code constraints for GNSS signal acquisition | |
| Tamazin et al. | Robust fine acquisition algorithm for GPS receiver with limited resources | |
| CN103760578B (zh) | 一种gnss卫星导航信号的无模糊跟踪方法 | |
| Zheng | A software-based frequency domain parallel acquisition algorithm for GPS signal | |
| Ahamed et al. | Fast acquisition of GPS signal using Radix-2 and Radix-4 FFT algorithms | |
| CN102798871A (zh) | 基于伪码重构的伪码捕获方法及装置 | |
| CN107290761B (zh) | 一种并行码相位搜索装置及实现并行码相位搜索的方法 | |
| CN107132552B (zh) | 一种并行码相位搜索装置及实现并行码相位搜索的方法 | |
| KR100778096B1 (ko) | 확산 대역 신호를 이용하는 다중 접속 수신기의 신호 동기획득 방법 및 그 장치 | |
| CN112764063A (zh) | 一种实现捕获处理的方法及接收机 | |
| CN107346028B (zh) | 一种实现码相位搜索的方法及装置 | |
| JP2006217601A (ja) | 地理的局所化システムの位置決め信号を捕捉するための方法、地理的局所化システム用受信機及び上記の方法を実行するためのプログラム命令を含んでいるコンピュータデータ媒体 | |
| US20120274512A1 (en) | Signal Processing Method, Device and System | |
| CN108226968B (zh) | 一种导航信号快速捕获方法 | |
| CN107132555A (zh) | 一种并行码相位搜索装置及实现并行码相位搜索的方法 | |
| Lin et al. | Acquisition of GPS software receiver using split-radix FFT | |
| CN107132554A (zh) | 一种并行码相位搜索装置及实现并行码相位搜索的方法 | |
| Sivani et al. | Development of Acquisition Module for a Software Based IRNSS Receiver | |
| GB2481575A (en) | Improvements to reception of spread spectrum signals | |
| CN106896384A (zh) | 一种频率搜索的装置及实现频率搜索的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |