CN112444830B - 基于奇异值分解的压缩感知gnss信号捕获方法 - Google Patents
基于奇异值分解的压缩感知gnss信号捕获方法 Download PDFInfo
- Publication number
- CN112444830B CN112444830B CN202011245625.6A CN202011245625A CN112444830B CN 112444830 B CN112444830 B CN 112444830B CN 202011245625 A CN202011245625 A CN 202011245625A CN 112444830 B CN112444830 B CN 112444830B
- Authority
- CN
- China
- Prior art keywords
- matrix
- signal
- singular value
- compressed sensing
- value decomposition
- 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
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 70
- 239000011159 matrix material Substances 0.000 claims abstract description 145
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 65
- 239000013598 vector Substances 0.000 claims abstract description 24
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 238000012935 Averaging Methods 0.000 claims abstract description 9
- 230000006835 compression Effects 0.000 claims description 13
- 238000007906 compression Methods 0.000 claims description 13
- 125000004122 cyclic group Chemical group 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 4
- 238000012856 packing Methods 0.000 claims description 4
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical group [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 3
- 239000010931 gold Substances 0.000 claims description 3
- 229910052737 gold Inorganic materials 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000008447 perception Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000004806 packaging method and process Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 230000000875 corresponding effect Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Mathematical Optimization (AREA)
- Radar, Positioning & Navigation (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种基于奇异值分解的压缩感知GNSS信号捕获方法，包括：基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵，形成基于奇异值分解的压缩感知噪声信号捕获算法；所述基于奇异值分解的压缩感知噪声信号捕获算法包括：对高斯随机测量矩阵进行奇异值平均，产生新的高斯随机观测矩阵，使高斯随机观测矩阵的列向量分布比高斯随机测量矩阵分散，以使测量后的高斯随机观测矩阵信号保留更多信息；根据稀疏矩阵和新的高斯随机观测矩阵完成信号的压缩感知捕获，以提高基于奇异值分解的压缩感知噪声信号捕获算法在低信噪比条件下的捕获性能。
Description
技术领域
本发明涉及信号捕获技术领域，特别涉及一种基于奇异值分解的压缩感知GNSS信号捕获方法。
背景技术
信号的捕获是GNSS信号接收的关键步骤，其搜索过程运算量较大。压缩感知作为一种信号采集新方法，可以减少捕获过程的数据量，但对捕获性能有一定影响。
全球导航卫星定位系统(Global Navigation Satellite System，GNSS)可为用户提供全方位、高精度的导航信息和授时信息。随着GPS、北斗卫星导航系统在智能设备中的普及，人们对其接收模块的性能及功耗提出了更高的要求。在GNSS信号接收处理中，捕获模块数据运算量大、资源消耗多，是影响接收模块性能及功耗的重要因素。
GNSS一般采用扩频调制技术，在卫星通信高动态环境下，需考虑多普勒频移，因此卫星信号的捕获需要进行多个码相位的搜索及大动态多普勒频率的搜索，运算量较大。传统的码分多址系统常采用串行搜索方法，对所有卫星可能的码相位和频率进行二维搜索，需消耗着大量计算资源，且带来捕获时间的增长。因此，降低捕获模块数据运算量减少捕获时间是接收机性能优化的关键。基于GNSS信号的稀疏性，压缩感知可用于信号捕获，大大降低捕获过程的数据运算量，但其算法在低信噪比条件下的捕获概率较小。
发明内容
本发明的目的在于提供一种基于奇异值分解的压缩感知GNSS信号捕获方法，以解决现有的压缩感知算法在低信噪比条件下的捕获概率较小的问题。
为解决上述技术问题，本发明提供一种基于奇异值分解的压缩感知GNSS信号捕获方法，包括：
基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵，形成基于奇异值分解的压缩感知噪声信号捕获算法；
所述基于奇异值分解的压缩感知噪声信号捕获算法包括：对高斯随机测量矩阵进行奇异值平均，产生新的高斯随机观测矩阵，使高斯随机观测矩阵的列向量分布比高斯随机测量矩阵分散，以使测量后的高斯随机观测矩阵信号保留更多信息；
根据稀疏矩阵和新的高斯随机观测矩阵完成信号的压缩感知捕获，以提高基于奇异值分解的压缩感知噪声信号捕获算法在低信噪比条件下的捕获性能。
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，所述基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵包括：
GPS的GNSS信号由载波频率、测距码和导航信息组成，在卫星发射端，先将测距码与导航信息调制为组合码，再通过扩频通信形成最终信号；接收信号时，射频前端先接收信号并通过采样下变频后成为离散的中频信号后，再进行处理，中频信号的数学表达式为：
式中，A为信号幅度；D(n)为导航信息，C(n)为伪随机的测距码，即C/A码，wIF为中频频率，为载波相位；wd为卫星和接收相对运动产生的多普勒频移，v(n)为均值为0，方差为2σ的高斯白噪声。
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，所述基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵还包括：
C/A码是码率为1.023Mbps的Gold码，码长L＝1023；每个卫星有唯一的C/A码，不同卫星的C/A码相互正交，同一卫星的C/A码及其循环移位序列也相互正交；选择C/A码的循环移位矩阵作为信号的稀疏变换域；
设c(0)＝[c0,c1,…cL-1]∈RL×1为任意卫星的C/A码，将其循环移位m个码片后，得到序列c(m)＝[cm,cm+1,…cN-1,c0,…cm-1]∈RL×1,m＝0,1,…L-1，则构造稀疏矩阵如下：
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，所述基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵还包括：
进行多普勒频率和C/A码相位的二维搜索，根据相关峰位置获得频率和相位估计值；在对某颗GPS卫星的采样信号进行载波及多普勒剥离和数据打包等处理后，一个C/A码周期的信号表示为：
或表示为向量形式：r(N×1)，当进行半码片打包时N＝2L＝2046；根据公式(2)对剥离了载波和多普勒的信号r进行稀疏表示：
其中，θ包含C/A码的相位信息，若多普勒频率搜索成功，则θ是一个稀疏向量，其峰值位置为C/A码相位估计值。
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，还包括：
信号的观测过程完成信号由高维到低维的映射，使用和稀疏矩阵G不相关的观测矩阵对信号进行压缩采样，和稀疏矩阵G不相关的观测矩阵为B(M×N)，M<<N，压缩比a＝M/N；
得到观测序列：
Y(M×1)＝B(M×N)·r(N×1) (5)。
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，奇异值分解包括：
设M是秩为r(r>0)的m×n阶实矩阵，则存在一个分解使得：
其中U是m×m阶酉矩阵；VT即V的共轭转置，是n×n阶酉矩阵，σ1,σ2,…σr为M的全部奇异值，且σ1≥σ2≥σ3≥…σr；
奇异值σi为映射后的正交基的模的大小，即/>向量/>为/>方向上的单位向量。
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，产生新的高斯随机观测矩阵包括：
得到高斯随机测量矩阵B，B是一个M×N M<<N的矩阵；
i＝1,2,…M；j＝1,2,…N
对其进行奇异值分解，得到奇异值；
其中
通过求奇异值均值得到改进的观测矩阵使矩阵的列向量分布均匀，从而含有更多信息；
对Σ对角线上的值求平均：Mean＝(σ1+σ2+…σM)/M；令
其中σ'1＝σ'2＝…σ'M＝Mean；
得到改进的观测矩阵：
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，还包括：
进行信号的重构；包括：
由测量信号Y及感知矩阵A(M×N)＝BG求解稀疏信号θ，在信号稀疏或可压缩的前提下，求解欠定方程组的问题转化为最小l0范数问题，最小l0范数下在一定条件下和最小l1范数具有等价性，得到相同的解；
转化为最小l1范数下的最优化问题：
min||θ||1s.t.Aθ＝Y (11)
通过重构算法重构出原稀疏信号；所述重构算法包括贪婪算法、凸优化或迭代阈值。
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，基于奇异值分解的压缩感知噪声信号捕获算法的流程为：
初始化捕获参数，基于本地伪码的循环相关获得稀疏矩阵G，
根据公式(7)到公式(10)获得高斯随机观测矩阵B，并计算出传感矩阵A；
在本地载波基础上设置多普勒估计值，剥离中频信号x(n)载波和多普勒频率，得到信号r；
测量重构信号；
遍历所有多普勒估计值，找出峰值，门限判别捕获是否成功。
可选的，在所述的基于奇异值分解的压缩感知GNSS信号捕获方法中，还包括：
利用观测量Y及感知矩阵A(M×N)对得到的测量信号进行重构，重构算法为正交匹配追踪，包括：输入为测量向量Y、传感矩阵A、稀疏度K；输出为信号的稀疏估计
初始化：
循环下列步骤：
找索引λt：先计算传感矩阵A各列与残差的内积：使得：λt＝argmaxj＝1,2,…NP；
令
用最小二乘法求解获得信号的近似解：
更新残差：
t＝t+1，若t≤K，则返回至循环的第一步；否则执行下一步；
为最后一次迭代所得/>
rt表示残差，t表示迭代次数，表示空集；Λt表示第t次迭代的索引(列索引)的集合；λt表示第t次迭代找到的索引(列序号)；aj表示矩阵A的第j列；At表示按索引选出的矩阵A的列集合；∪符号表示集合并运算；<·，·>表示求向量内积；
重构得到N点稀疏数据后，这一频率点的搜索结束；
改变多普勒频率，进行下一个频率点的搜索；
将所有频率范围搜索完毕后，找出相关值最高峰，该最高峰值与该频率下的平均相关值的比值大于门限Vt时，则捕获成功，记录该恢复信号对应的多普勒频移和码相位；否则捕获失败。
在本发明提供的基于奇异值分解的压缩感知GNSS信号捕获方法中，在GNSS信号稀疏性的基础上，构造出基于奇异值分解的高斯随机观测矩阵，并利用该矩阵对压缩感知捕获算法进行仿真分析，该观测矩阵比传统高斯观测矩阵具有更好的非相关性和重构性能；与高斯压缩感知捕获算法对比，本发明算法在低信噪比情况下捕获概率有明显提升，提高压缩感知用于捕获时的性能。
本发明基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵，提出了一种基于奇异值分解的压缩感知噪声信号捕获方法，该方法通过对高斯随机测量矩阵进行奇异值平均，使其列向量分布更分散，测量后的信号保留了更多信息，进而使得压缩感知算法在低信噪比条件下的捕获性能得到提高。
附图说明
图1是本发明一实施例SVD的数学意义示意图；
图2是本发明一实施例基于奇异值分解的压缩感知GNSS信号捕获方法示意图；
图3是本发明一实施例基于奇异值分解的压缩感知GNSS信号捕获方法捕获结果示意图；
图4是本发明一实施例基于奇异值分解的压缩感知GNSS信号捕获方法不同压缩比下的捕获结果对比示意图；
图5(a)是本发明一实施例基于奇异值分解的压缩感知GNSS信号捕获方法的压缩比为0.75捕获结果示意图；
图5(b)是本发明一实施例基于奇异值分解的压缩感知GNSS信号捕获方法的压缩比为0.5捕获结果示意图；
图5(c)是本发明一实施例基于奇异值分解的压缩感知GNSS信号捕获方法的压缩比为0.25捕获结果示意图；
图5(d)是本发明一实施例基于奇异值分解的压缩感知GNSS信号捕获方法的压缩比为0.125捕获结果示意图。
具体实施方式
以下结合附图和具体实施例对本发明提出的基于奇异值分解的压缩感知GNSS信号捕获方法作进一步详细说明。根据下面说明和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
另外，除非另行说明，本发明的不同实施例中的特征可以相互组合。例如，可以用第二实施例中的某特征替换第一实施例中相对应或功能相同或相似的特征，所得到的实施例同样落入本申请的公开范围或记载范围。
本发明的核心思想在于提供一种基于奇异值分解的压缩感知GNSS信号捕获方法，以解决现有的压缩感知算法在低信噪比条件下的捕获概率较小的问题。
现有技术的一种基于圆周移位的并行码相位算法，将搜索中多次FFT/IFFT运算用计算量很小的频谱序列的搬移操作来代替，使运算量大幅下降。以及在基于FFT的并行码相位捕获算法基础上，基于稀疏傅里叶的方法，通过提高捕获过程中相关运算的效率，减少了运算量。现有技术还有一种双块零扩展截断相关的捕获算法，利用复序列频谱的非对称性对信号进行了截断预处理，结合双块零扩展与圆周移位进行分块相关以减少运算量。压缩感知(Compressed Sensing，CS)是由Donoho、Candes、Romberg、Tao等人在2004年提出的一种新的信号处理理论。该理论的提出克服了传统信号处理中奈奎斯特采样定理的限制，充分利用信号的稀疏特性，大大降低信号的采样速率，为进一步降低捕获所需的数据量提供了新思路。压缩感知用于GNSS信号捕获的关键点主要有两个：信号稀疏性和合适的测量重构算法。基于这两个关键点，现有技术以C/A码的正交基矩阵为稀疏域，利用Walsh-Hadamard矩阵的双阶段确定性压缩感知GNSS捕获技术，但需要两次Walsh-Hadamard变换，过程复杂。或基于峰值位置的压缩感知方案，但稀疏矩阵存储量巨大。以及将压缩感知理论融入交错方向乘子法(ADMM)的框架，提出一种高效的并行捕获算法。基于压缩感知改进的CSPMF－FFT算法，降低了数据量和硬件资源消耗，但没有考虑低信噪比问题。
为实现上述思想，本发明提供了一种基于奇异值分解的压缩感知GNSS信号捕获方法，包括：基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵，形成基于奇异值分解的压缩感知噪声信号捕获算法；所述基于奇异值分解的压缩感知噪声信号捕获算法包括：对高斯随机测量矩阵进行奇异值平均，产生新的高斯随机观测矩阵，使高斯随机观测矩阵的列向量分布比高斯随机测量矩阵分散，以使测量后的高斯随机观测矩阵信号保留更多信息；根据稀疏矩阵和新的高斯随机观测矩阵完成信号的压缩感知捕获，以提高基于奇异值分解的压缩感知噪声信号捕获算法在低信噪比条件下的捕获性能。
信号稀疏性以GPS为例，GNSS信号由载波频率、测距码和导航信息组成。在卫星发射端，先将测距码与导航信息调制为组合码，再通过扩频通信形成最终信号。接收信号时，射频前端先接收信号并通过采样下变频后成为离散的中频信号后，再进行处理，中频信号的数学表达式为：
式中，A为信号幅度。D(n)为导航信息，C(n)为伪随机的测距码，即C/A码，wIF为中频频率，为载波相位。wd为卫星和接收相对运动产生的多普勒频移，v(n)为均值为0，方差为2σ的高斯白噪声。
C/A码是码率为1.023Mbps的Gold码，码长L＝1023。每个卫星有唯一的C/A码，不同卫星的C/A码相互正交，同一卫星的C/A码及其循环移位序列也相互正交。卫星接收信号本身并不具有稀疏性，但由于C/A码良好的自相关性和互相关性，在码相位搜索域上，只有本地码相位与接收信号码相位对齐时相关结果才具有峰值，因此其码相关结果具有稀疏性，因此可选择C/A码的循环移位矩阵作为信号的稀疏变换域。
假设c(0)＝[c0,c1,…cL-1]∈RL×1为任意卫星的C/A码，将其循环移位m个码片后，得到序列c(m)＝[cm,cm+1,…cN-1,c0,…cm-1]∈RL×1,m＝0,1,…L-1。则构造稀疏矩阵如下。
捕获过程实质是进行多普勒频率和C/A码相位的二维搜索，根据相关峰位置获得频率和相位估计值。在对某颗GPS卫星的采样信号进行载波及多普勒剥离和数据打包等处理后，一个C/A码周期的信号可以表示为：
也可表示为向量形式：r(N×1)，当进行半码片打包时N＝2L＝2046。根据式(2)可对剥离了载波和多普勒的信号r进行稀疏表示：
其中，θ包含了C/A码的相位信息，如果多普勒频率搜索成功，即则θ是一个稀疏向量，仅有少数较大的非零值，其峰值位置为C/A码相位估计值。如果多普勒估计值存在误差，则θ中仍有多普勒分量，其稀疏性就会受到影响。
具体的，基于奇异值分解的压缩感知捕获算法包括：压缩感知的过程分为信号的稀疏表示、信号观测以及信号重构三个关键部分，其中信号的稀疏性是压缩感知应用的前提。
当长度为N的信号中只有K(K＜＜N)个非零的值，或其他值都接近于零时，我们认为信号具有稀疏性，其稀疏度为K。但通常实际处理信号不稀疏，这就需要对原始信号进行稀疏表示，即对原始信号进行某种变换，变换到某个域，使其在该域的表示是稀疏的。基于式(4)中GPS信号的稀疏表示，可利用压缩感知理论完成对码相位的搜索。
进一步的，基于奇异值分解的观测矩阵构造包括：信号的观测过程完成了信号由高维到低维的映射，使用和稀疏矩阵G不相关的观测矩阵B(M×N)，M<<N，对信号进行压缩采样，压缩比a＝M/N。得到观测序列：即：
Y(M×1)＝B(M×N)·r(N×1) (5)
观测矩阵B的构造是压缩感知的关键。要想使信号完全重构，必须保证观测矩阵不会把两个不同的稀疏信号映射到同一个采样集中，这就要求观测矩阵中任M个列向量组成的矩阵行列式不为零，即为非奇异矩阵。
具有以上特征的矩阵都可以作为压缩感知观测矩阵。目前观测矩阵主要有随机矩阵、确定性矩阵和部分随机矩阵等。随机高斯观测矩阵与大多数正交稀疏矩阵都不相关，是压缩感知中适应性最广的观测矩阵，但其自身的不确定性会造成捕获结果的不确定性，使捕获成功率下降。本发明根据奇异值分解对其进行改进。
设M是秩为r(r>0)的m×n阶实矩阵，则存在一个分解使得：
其中
U是m×m阶酉矩阵；VT即V的共轭转置，是n×n阶酉矩阵，这样的分解就称作M的奇异值分解(Singular Value Decomposition，SVD)，σ1,σ2,…σr为M的全部奇异值，且σ1≥σ2≥σ3≥…σr。
奇异值分解在代数上表现就是M将一组正交基V转化为另一组正交基U，以二维矩阵的图解为例，在图1中：奇异值σi就是映射后的正交基的模的大小，即/>向量/>为/>方向上的单位向量。几何上来说SVD就是分解，将一个复杂的空间变换分解为三个简单的基本变换，V是旋转，∑是缩放，U是投影。
从物理意义方面来说，矩阵奇异值往往对应着矩阵中隐含的重要信息，且重要性和奇异值大小正相关。奇异值衡量每个列向量对于M的权重。结合其数学意义，可以说奇异值σi代表矩阵M的列向量在子空间的分布。如果矩阵M的r个奇异值的大小差距越小，则表明其列向量分布越分散，其所含信息越多，这正好符合压缩感知对于观测矩阵的要求。
本发明基于以上分析，对高斯随机观测矩阵提出改进思路：先得到随机高斯观测矩阵B，B是一个M×N(M＜＜N)的矩阵；然后对其进行奇异值分解，得到奇异值；最后通过求奇异值均值得到改进的观测矩阵使矩阵的列向量分布均匀，从而含有更多信息，使测量后的信号更完整，便于之后更好的重构。算法如下：
1)生成高斯随机观测矩阵BM×N：
i＝1,2,…M；j＝1,2,…N
2)对B进行奇异值分解：
其中
3)对Σ对角线上的值求平均：Mean＝(σ1+σ2+…σM)/M；令
其中σ'1＝σ'2＝…σ'M＝Mean。
4)得到改进的观测矩阵：
根据矩阵奇异值的分析可知：改进的观测矩阵比原高斯随机观测矩阵B，含有更多更分散的信息，能更好的观测信号，且在相同范数意义下有更高的重构精度。得到新的观测矩阵之后，就可以利用式(5)对信号进行观测，得到观测量Y。
进一步的，重构算法包括：信号的重构是由测量信号Y及感知矩阵A(M×N)＝BG求解稀疏信号θ的过程，这是一个求解欠定方程的问题，在信号稀疏或可压缩的前提下，求解欠定方程组的问题转化为最小l0范数问题，最小l0范数下在一定条件下和最小l1范数具有等价性，可得到相同的解。则可以转化为最小l1范数下的最优化问题：
min||θ||1s.t.Aθ＝Y (11)
该方程有无穷多组解，理论分析已经证明，在满足不相关性和稀疏性的条件下，可以通过重构算法重构出原稀疏信号。常用的重构算法有贪婪算法，凸优化，迭代阈值等。
利用观测量Y及感知矩阵A(M×N)对得到的测量信号进行重构，基于信号捕获的特殊性，重构算法只需重构稀疏度很小的稀疏信号且主要关注峰值位置。重构算法选择计算量小，稳定的正交匹配追踪(Orthogonal Matching Pursuit，OMP)算法。该算法流程如下：
表1 OMP算法流程
重构得到N点稀疏数据后，这一频率点的搜索即结束。改变多普勒频率，进行下一个频率点的搜索。将所有频率范围搜索完毕后，找出相关值最高峰，该最高峰值与该频率下的平均相关值的比值大于门限Vt时，则捕获成功，记录该恢复信号对应的多普勒频移和码相位；否则捕获失败。
其中，基于以上分析，本发明提出了基于奇异值分解的压缩感知(SVD-CS)捕获算法，其算法流程如图2所示：首先初始化捕获参数，基于本地伪码的循环相关获得稀疏矩阵G，根据式(7)到式(10)获得随机观测矩阵B，并计算出传感矩阵A；然后在本地载波基础上设置多普勒估计值，剥离载波和多普勒频率，得到信号r；接着测量重构信号；最后遍历所以多普勒估计值，找出峰值，门限判别捕获是否成功。
对本发明提出的压缩感知捕获算法进行仿真。随机产生码相位移位和多普勒频移。频率搜索步长500Hz，码相位搜索步长半码片。门限值Vt设为2。
仿真一是本发明算法的捕获实现。用压缩比a为0.25，重构稀疏度K为10的SVD-CS捕获算法对信噪比为-10dB的模拟信号进行捕获，接收信号的C/A码相位设为501，多普勒频率设为1050Hz。由图3的捕获结果可以看出：在多普勒频率和码相位的二维搜索域内，出现了明显的峰值，表明捕获到了信号。同时，相关峰出现的位置对应的多普勒频率为1000Hz，码相位为501，捕获结果在精度范围要求内，实现了正确捕获。结果表明，本发明算法能成功捕获信号。
仿真二是本发明算法特性分析，如图4所示，通过比较不同压缩比下捕获概率随信噪比的变化分析本发明算法性能。仿真结果表明：本发明算法能够在低于码相位精度要求的数据运算下完成GPS信号的捕获，代价是捕获概率略有下降；当压缩比a越小时，捕获概率越小，但算法相关数目从N降到M，大大地降低了捕获的数据运算量。
仿真三是不同算法捕获性能对比。以传统串行算法为对照组，主要对比高斯随机测量矩阵压缩感知算法(GS-CS)与本发明SVD-CS算法的捕获性能，由图5(a)、(b)、(c)、(d)对比可知：不同压缩比下本发明算法的捕获概率均有明显提高，且压缩比越接近1，提高越明显，表明了本发明算法改进的有效性。基于GNSS信号的稀疏性，压缩感知可用于信号捕获，大大降低捕获过程的数据运算量，但其算法在低信噪比条件下的捕获概率较小。本发明以提高捕获概率为目标，从改进观测矩阵出发，提出了一种基于奇异值分解的压缩感知捕获算法。首先基于GNSS信号稀疏性构造稀疏矩阵；其次，对高斯随机观测矩阵进行奇异值平均，产生新的观测矩阵，最后利用稀疏矩阵和新的观测矩阵完成信号的压缩感知捕获。仿真结果验证了本发明算法的有效性，提高了低信噪比条件下压缩感知算法的捕获性能。
综上，上述实施例对基于奇异值分解的压缩感知GNSS信号捕获方法的不同构型进行了详细说明，当然，本发明包括但不局限于上述实施中所列举的构型，任何在上述实施例提供的构型基础上进行变换的内容，均属于本发明所保护的范围。本领域技术人员可以根据上述实施例的内容举一反三。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
上述描述仅是对本发明较佳实施例的描述，并非对本发明范围的任何限定，本发明领域的普通技术人员根据上述揭示内容做的任何变更、修饰，均属于权利要求书的保护范围。
Claims (10)
1.一种基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，包括：
基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵，形成基于奇异值分解的压缩感知噪声信号捕获算法；
所述基于奇异值分解的压缩感知噪声信号捕获算法包括：对高斯随机测量矩阵进行奇异值平均，产生新的高斯随机观测矩阵，使高斯随机观测矩阵的列向量分布比高斯随机测量矩阵分散，以使测量后的高斯随机观测矩阵信号保留更多信息；
根据稀疏矩阵和新的高斯随机观测矩阵完成信号的压缩感知捕获，以提高基于奇异值分解的压缩感知噪声信号捕获算法在低信噪比条件下的捕获性能。
2.如权利要求1所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，基于奇异值分解的压缩感知噪声信号捕获算法的流程为：
初始化捕获参数，基于本地伪码的循环相关获得稀疏矩阵G，
获得高斯随机观测矩阵B，并计算出传感矩阵A；
在本地载波基础上设置多普勒估计值，剥离中频信号x(n)的载波和多普勒频率，得到信号r；
测量重构信号；
遍历所有多普勒估计值，找出峰值，门限判别捕获是否成功。
3.如权利要求2所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，所述基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵包括：
GPS的GNSS信号由载波频率、测距码和导航信息组成，在卫星发射端，先将测距码与导航信息调制为组合码，再通过扩频通信形成最终信号；接收信号时，射频前端先接收信号并通过采样下变频后成为离散的中频信号后，再进行处理，中频信号的数学表达式为：
式中，A为信号幅度；D(n)为导航信息，C(n)为伪随机的测距码，即C/A码，wIF为中频频率，为载波相位；wd为卫星和接收机相对运动产生的多普勒频移，v(n)为均值为0，方差为2σ的高斯白噪声。
4.如权利要求3所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，所述基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵还包括：
C/A码是码率为1.023Mbps的Gold码，码长L＝1023；每个卫星有唯一的C/A码，不同卫星的C/A码相互正交，同一卫星的C/A码及其循环移位序列也相互正交；选择C/A码的循环移位矩阵作为信号的稀疏变换域；
设c(0)＝[c0,c1,…cL-1]∈RL×1为任意卫星的C/A码，将其循环移位m个码片后，得到序列c(m)＝[cm,cm+1,…cL-1,c0,…cm-1]∈RL×1,m＝0,1,…L-1，则构造稀疏矩阵如下：
5.如权利要求4所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，所述基于对GNSS信号稀疏性的分析，利用C/A码构造稀疏矩阵还包括：
进行多普勒频率和C/A码相位的二维搜索，根据相关峰位置获得频率和相位估计值；在对某颗GPS卫星的采样信号进行载波及多普勒剥离和数据打包处理后，一个C/A码周期的信号表示为：
或表示为向量形式：r(N×1)，当进行半码片打包时N＝2L＝2046；根据公式(2)对剥离了载波和多普勒的信号r进行稀疏表示：
其中，θ包含C/A码的相位信息，若多普勒频率搜索成功，则θ是一个稀疏向量，其峰值位置为C/A码相位估计值。
6.如权利要求5所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，还包括：
信号的观测过程完成信号由高维到低维的映射，使用和稀疏矩阵G不相关的观测矩阵对信号进行压缩采样，和稀疏矩阵G不相关的观测矩阵为B(M×N)，M<<N，压缩比a＝M/N；
得到观测序列：
Y(M×1)＝B(M×N)·r(N×1) (5)。
7.如权利要求6所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，奇异值分解包括：
设M是秩为r的m×n阶实矩阵，r>0，则存在一个分解使得：
其中U是m×m阶酉矩阵；VT即V的共轭转置，是n×n阶酉矩阵，σ1,σ2,…σr是M的全部奇异值，且σ1≥σ2≥σ3≥…σr；
奇异值σi为映射后的正交基的模的大小，即/>向量/>为/>方向上的单位向量。
8.如权利要求7所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，产生新的高斯随机观测矩阵包括：
得到高斯随机测量矩阵B，B是一个M×N M<<N的矩阵；
对其进行奇异值分解，得到奇异值；
其中
通过求奇异值均值得到改进的观测矩阵B，使矩阵的列向量分布均匀，从而含有更多信息；
对Σ对角线上的值求平均：Mean＝(σ1+σ2+…σM)/M；令
其中σ'1＝σ'2＝…σ'M＝Mean；
得到改进的观测矩阵：
9.如权利要求2所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，还包括：
进行信号的重构；包括：
由测量信号Y及感知矩阵A(M×N)＝BG求解稀疏信号θ，在信号稀疏或可压缩的前提下，求解欠定方程组的问题转化为最小l0范数问题，最小l0范数下在一定条件下和最小l1范数具有等价性，得到相同的解；
转化为最小l1范数下的最优化问题：
min||θ||1s.t.Aθ＝Y (11)
通过重构算法重构出原稀疏信号；所述重构算法包括贪婪算法、凸优化或迭代阈值。
10.如权利要求9所述的基于奇异值分解的压缩感知GNSS信号捕获方法，其特征在于，还包括：
利用观测量Y及感知矩阵A(M×N)对得到的测量信号进行重构，重构算法为正交匹配追踪，包括：输入为测量向量Y、传感矩阵A、稀疏度K；输出为信号的稀疏估计
初始化：r0＝Y；
循环下列步骤：
找索引λt：先计算传感矩阵A各列与残差的内积：P＝|<rt-1,aj>|2，使得：λt＝argmaxj＝1,2,...NP；
令Λt＝Λt-1∪{λt},
用最小二乘法求解获得信号的近似解：
更新残差：
t＝t+1，若t≤K，则返回至循环的第一步；否则执行下一步；
为最后一次迭代所得/>
rt表示残差，t表示迭代次数，表示空集；Λt表示第t次迭代的索引的集合；λt表示第t次迭代找到的索引；aj表示矩阵A的第j列；At表示按索引选出的矩阵A的列集合；∪符号表示集合并运算；<·，·>表示求向量内积；
重构得到N点稀疏数据后，这一频率点的搜索结束；
改变多普勒频率，进行下一个频率点的搜索；
将所有频率范围搜索完毕后，找出相关值最高峰，该最高峰值与该频率下的平均相关值的比值大于门限Vt时，则捕获成功，记录恢复信号对应的多普勒频移和码相位；否则捕获失败。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011245625.6A CN112444830B (zh) | 2020-11-10 | 2020-11-10 | 基于奇异值分解的压缩感知gnss信号捕获方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011245625.6A CN112444830B (zh) | 2020-11-10 | 2020-11-10 | 基于奇异值分解的压缩感知gnss信号捕获方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112444830A CN112444830A (zh) | 2021-03-05 |
| CN112444830B true CN112444830B (zh) | 2023-12-29 |
Family
ID=74736772
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011245625.6A Active CN112444830B (zh) | 2020-11-10 | 2020-11-10 | 基于奇异值分解的压缩感知gnss信号捕获方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112444830B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113159331B (zh) * | 2021-05-24 | 2023-06-30 | 同济大学 | 一种网络化机器学习系统的自适应稀疏度量化方法 |
| CN113376630A (zh) * | 2021-07-05 | 2021-09-10 | 中国人民解放军空军工程大学 | 雷达成像的方法、装置、电子设备和计算机可读存储介质 |
| CN114442032B (zh) * | 2022-04-07 | 2022-06-14 | 中国电子科技集团公司第二十九研究所 | 一种基于多极化矢量天线阵列压缩采样的测向方法 |
| CN116347329B (zh) * | 2022-12-16 | 2024-03-29 | 中交星宇科技有限公司 | 定位信号捕获方法及装置、计算设备及计算机存储介质 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7289049B1 (en) * | 2006-08-21 | 2007-10-30 | L3 Communications Integrated Systems L.P. | Method and apparatus for compressed sensing |
| CN102830409A (zh) * | 2012-08-30 | 2012-12-19 | 西安电子科技大学 | 一种基于压缩感知的导航信号采集方法 |
| KR20130087982A (ko) * | 2012-01-30 | 2013-08-07 | 한국과학기술원 | Ｇｎｓｓ 및 대역확산 신호의 신속한 신호 획득을 위한 압축 감지 방법 및 그 장치 |
| CN104199057A (zh) * | 2014-09-11 | 2014-12-10 | 重庆大学 | 一种基于压缩感知技术的boc信号无模糊捕获算法 |
| CN105763225A (zh) * | 2016-04-25 | 2016-07-13 | 山东大学 | 一种脉冲超宽带信号的稀疏字典的产生方法 |
| CN106291610A (zh) * | 2015-06-12 | 2017-01-04 | 北京信息科技大学 | 一种用于gnss信号压缩捕获处理装置的压缩并行相关模块及其实现方法 |
| CN109765586A (zh) * | 2019-01-09 | 2019-05-17 | 安徽师范大学 | 基于压缩感知的glonass信号快速捕获方法 |
| CN111624632A (zh) * | 2020-06-30 | 2020-09-04 | 青岛杰瑞自动化有限公司 | 一种卫星导航信号捕获方法及装置 |
-
2020
- 2020-11-10 CN CN202011245625.6A patent/CN112444830B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7289049B1 (en) * | 2006-08-21 | 2007-10-30 | L3 Communications Integrated Systems L.P. | Method and apparatus for compressed sensing |
| KR20130087982A (ko) * | 2012-01-30 | 2013-08-07 | 한국과학기술원 | Ｇｎｓｓ 및 대역확산 신호의 신속한 신호 획득을 위한 압축 감지 방법 및 그 장치 |
| CN102830409A (zh) * | 2012-08-30 | 2012-12-19 | 西安电子科技大学 | 一种基于压缩感知的导航信号采集方法 |
| CN104199057A (zh) * | 2014-09-11 | 2014-12-10 | 重庆大学 | 一种基于压缩感知技术的boc信号无模糊捕获算法 |
| CN106291610A (zh) * | 2015-06-12 | 2017-01-04 | 北京信息科技大学 | 一种用于gnss信号压缩捕获处理装置的压缩并行相关模块及其实现方法 |
| CN105763225A (zh) * | 2016-04-25 | 2016-07-13 | 山东大学 | 一种脉冲超宽带信号的稀疏字典的产生方法 |
| CN109765586A (zh) * | 2019-01-09 | 2019-05-17 | 安徽师范大学 | 基于压缩感知的glonass信号快速捕获方法 |
| CN111624632A (zh) * | 2020-06-30 | 2020-09-04 | 青岛杰瑞自动化有限公司 | 一种卫星导航信号捕获方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| A SINS-assisted Fast Acquisition Method for GNSS Signals Based on Compressed Sensing;Junbing Cheng et al.;《2018 Sixth International Conference on Advanced Cloud and Big Data》;第83-88页 * |
| Energy Efficient GNSS Signal Acquisition Using Singular Value Decomposition (SVD);Juan Carlos Bermúdez Ordoñez et al.;《Sensors》;第1-27页 * |
| 一种适用于GNSS信号的新颖捕获技术研究;欧松林 等;《第六届中国卫星导航学术年会论文集—S09PNT体系与导航新技术》;第102-106页 * |
| 基于压缩感知ITNSL0算法的BOC调制导航信号的稀疏性分析;张益东 等;《第九届中国卫星导航学术年会论文集——S03 卫星导航信号及抗干扰技术》;第130-134页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112444830A (zh) | 2021-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112444830B (zh) | 基于奇异值分解的压缩感知gnss信号捕获方法 | |
| Marques et al. | A review of sparse recovery algorithms | |
| CN102830409B (zh) | 一种基于压缩感知的导航信号采集方法 | |
| Wang et al. | The robust sparse Fourier transform (RSFT) and its application in radar signal processing | |
| Li et al. | GPS signal acquisition via compressive multichannel sampling | |
| Borio | M-sequence and secondary code constraints for GNSS signal acquisition | |
| He et al. | GPS signal acquisition based on compressive sensing and modified greedy acquisition algorithm | |
| Quan et al. | FPGA implementation of real-time compressive sensing with partial Fourier dictionary | |
| CN114879231A (zh) | Gnss信号压缩捕获方法及装置、电子设备、存储介质 | |
| CN114019449A (zh) | 信号源的波达方向估计方法、装置、电子设备及存储介质 | |
| He et al. | Sparse GLONASS signal acquisition based on compressive sensing and multiple measurement vectors | |
| Taghavi et al. | High resolution compressed sensing radar using difference set codes | |
| Sagiraju et al. | Fast acquisition implementation for high sensitivity global positioning systems receivers based on joint and reduced space search | |
| Park et al. | High‐resolution time delay estimation via sparse parameter estimation methods | |
| Liu et al. | General similar sensing matrix pursuit: An efficient and rigorous reconstruction algorithm to cope with deterministic sensing matrix with high coherence | |
| CN116125502B (zh) | 一种导航信号生成方法、装置、设备及存储介质 | |
| Al Irkhis et al. | Compressive chirp transform for estimation of chirp parameters | |
| CN111679123A (zh) | 适用于多模调制体制的符号边沿与频率估计方法与系统 | |
| Chang | Modified compressive sensing approach for GNSS signal reception in the presence of interference | |
| Xiong et al. | High resolution TOA estimation based on compressed sensing | |
| Xu et al. | A fast acquisition algorithm of GNSS receiver based on SFFT | |
| CN112764063A (zh) | 一种实现捕获处理的方法及接收机 | |
| CN107132552B (zh) | 一种并行码相位搜索装置及实现并行码相位搜索的方法 | |
| Cheng et al. | A SINS‐aided two‐step fast acquisition method for GNSS signal based on compressive sensing | |
| Zhang et al. | Optimizing the Gabor Bandwidth of satellite navigation signals by MCS signal expression |
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 |