CN112731486A - 信号角度和信号频率的估计方法、装置、设备及存储介质 - Google Patents

信号角度和信号频率的估计方法、装置、设备及存储介质 Download PDF

Info

Publication number
CN112731486A
CN112731486A CN202011556223.8A CN202011556223A CN112731486A CN 112731486 A CN112731486 A CN 112731486A CN 202011556223 A CN202011556223 A CN 202011556223A CN 112731486 A CN112731486 A CN 112731486A
Authority
CN
China
Prior art keywords
signal
satellite
frequency
angle
estimating
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.)
Granted
Application number
CN202011556223.8A
Other languages
English (en)
Other versions
CN112731486B (zh
Inventor
李强
黄磊
赵博
黄敏
孙维泽
张沛昌
刘仕奇
赵源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen University
Original Assignee
Shenzhen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen University filed Critical Shenzhen University
Priority to CN202011556223.8A priority Critical patent/CN112731486B/zh
Publication of CN112731486A publication Critical patent/CN112731486A/zh
Priority to PCT/CN2021/123920 priority patent/WO2022134762A1/zh
Priority to US18/041,734 priority patent/US20230314622A1/en
Application granted granted Critical
Publication of CN112731486B publication Critical patent/CN112731486B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction
    • 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/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/37—Hardware or software details of the signal processing chain
    • 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/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/36—Constructional details or hardware or software details of the signal processing chain relating to the receiver frond end
    • 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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开一种信号角度和信号频率的估计方法、装置、设备及存储介质，其中，所述信号角度和信号频率的估计方法包括步骤：通过天线阵元采集卫星信号的初始数据；将初始数据代入信号函数中，求解得到卫星信号；根据所述卫星信号，获取卫星信号的角度和频率，考虑在非高斯噪声环境下，同时考虑卫星信号空域稀疏特性和频率域稀疏特性，进行卫星信号角度和频率联合估计，有利于卫星导航接收机后续捕获跟踪卫星信号。

Description

信号角度和信号频率的估计方法、装置、设备及存储介质
技术领域
本发明涉及卫星信号处理领域，尤其涉及一种信号角度和信号频率的估计方法、装置、设备及存储介质。
背景技术
研究卫星信号角度估计或卫星信号频率估计通常是在高斯噪声背景下进行研究的。然而，在卫星导航接收机工作中，不可避免的出现非高斯噪声情况，如脉冲噪声，目前还缺少在非高斯噪声背景下估计卫星信号的角度和卫星信号的频率。
发明内容
本发明的主要目的在于提出一种信号角度和信号频率的估计方法、装置、设备及计算机可读存储介质，旨在提供一种在非高斯噪声的背景下估计卫星信号的角度和频率。信号角度和信号频率的估计方法包括以下步骤：
通过天线阵元采集卫星信号的初始数据；
将所述初始数据代入信号函数中，求解得到卫星信号；
根据所述卫星信号，获取所述卫星信号的角度和频率。
在一种实施方式中，将所述初始数据代入信号函数中，求解得到卫星信号的步骤包括：
对所述信号函数进行重加权处理，得到目标函数；
将所述初始数据代入所述目标函数，求解得到卫星信号。
在一种实施方式中，所述将所述初始数据代入所述目标函数，求解得到卫星信号的步骤包括：
将所述初始数据代入所述目标函数，采用梯度求导的方法，求解所述目标函数，得到卫星信号。
在一种实施方式中，所述根据所述卫星信号，获取所述卫星信号的角度和频率的步骤包括：
对所述卫星信号进行傅里叶变换，得到所述卫在一种实施方式中，所述信号函数的表达式为：
Figure BDA0002855967630000021
其中，所述X表示所述卫星信号，所述Y表示初始数据，所述A表示导向矢量矩阵，所述λ1所述λ2表示正则化因子，所述F表示傅里叶变换。
在一种实施方式中，所述目标函数的表达式为：
Figure BDA0002855967630000022
其中，ρi和κi分别为加权向量ρ和κ中的第i个元素，i＝1，2，…，J，所述J为所述A的列数。
此外，为实现上述目的，本发明还提供一种信号角度和信号频率的估计装置，所述信号角度和信号频率的估计装置包括：
采集模块，用于通过天线阵元采集卫星信号的初始数据；
代入求解模块，用于将所述反馈数据发送至所述共享平台；
获取模块，用于根据所述卫星信号，获取所述卫星信号的角度和频率。
此外，为实现上述目的，本发明还提供一种信号角度和信号频率的估计设备，所述信号角度和信号频率的估计设备包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的信号角度和信号频率的估计程序，所述信号角度和信号频率的估计程序被所述处理器执行时实现如上所述的信号角度和信号频率的估计方法的步骤。
此外，为实现上述目的，本发明还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有所述信号角度和信号频率的估计程序，所述信号角度和信号频率的估计程序被处理器执行时实现如上所述的信号角度和信号频率的估计方法的步骤。
本发明通过天线阵元采集卫星信号的初始数据，将初始数据代入信号函数中求解得到卫星信号，根据卫星信号获取卫星信号的角度和频率，考虑在非高斯噪声环境下，同时考虑卫星信号空域稀疏特性和频率域稀疏特性，进行卫星信号角度和频率联合估计，有利于卫星导航接收机后续捕获跟踪卫星信号。
附图说明
图1为实现本发明各个实施例一种设备的硬件结构示意图；
图2为本发明信号角度和信号频率的估计方法第一实施例的流程示意图；
图3为本发明空域和时域示意图；
图4为本发明空域与频率域示意图；
图5为卫星信号估计值的MSE随着迭代次数仿真图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做说明。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明提供了一种信号角度和信号频率的估计设备，参照图1，图1是本发明实施例方案涉及的硬件运行环境的结构示意图。
需要说明的是，图1即可为信号角度和信号频率的估计设备的硬件运行环境的结构示意图。本发明实施例信号角度和信号频率的估计设备可以是PC(Personal Computer，个人电脑)，便携计算机，服务器等设备。
如图1所示，该信号角度和信号频率的估计设备可以包括：处理器1001，例如CPU，存储器1005，用户接口1003，网络接口1004，通信总线1002。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(Display)、输入单元比如键盘(Keyboard)，可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如WI-FI接口)。存储器1005可以是高速RAM存储器，也可以是稳定的存储器(non-volatile memory)，例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
可选地，信号角度和信号频率的估计设备还可以包括RF(Radio Frequency，射频)电路，传感器、WiFi模块等等。
本领域技术人员可以理解，图1中示出的信号角度和信号频率的估计设备结构并不构成信号角度和信号频率的估计设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
如图1所示，作为一种计算机存储可读存储介质的存储器1005中可以包括操作设备、网络通信模块、用户接口模块以及信号角度和信号频率的估计程序。其中，操作设备是管理和控制信号角度和信号频率的估计设备硬件和软件资源的程序，支持信号角度和信号频率的估计程序以及其它软件或程序的运行。
图1所示的信号角度和信号频率的估计设备，用于提供一种在非高斯噪声的背景下估计卫星信号的角度和频率，用户接口1003主要用于侦测或者输出各种信息，如输入初始数据和输出卫星信号等；网络接口1004主要用于与后台服务器交互，进行通信；处理器1001可以用于调用存储器1005中存储的信号角度和信号频率的估计程序，并执行以下操作：
通过天线阵元采集卫星信号的初始数据；
将所述初始数据代入信号函数中，求解得到卫星信号；
根据所述卫星信号，获取所述卫星信号的角度和频率。
本发明通过天线阵元采集卫星信号的初始数据，将初始数据代入信号函数中求解得到卫星信号，根据卫星信号获取卫星信号的角度和频率，考虑在非高斯噪声环境下，同时考虑卫星信号空域稀疏特性和频率域稀疏特性，进行卫星信号角度和频率联合估计，有利于卫星导航接收机后续捕获跟踪卫星信号。
本发明移动终端具体实施方式与下述信号角度和信号频率的估计方法各实施例基本相同，在此不再赘述。
基于上述结构，提出本发明信号角度和信号频率的估计方法的各个实施例。
本发明提供一种信号角度和信号频率的估计方法。
参照图2，图2为本发明信号角度和信号频率的估计方法第一实施例的流程示意图。
在本实施例中，提供了信号角度和信号频率的估计方法的实施例，需要说明的是，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
在本实施例中，信号角度和信号频率的估计方法包括：
步骤S10，通过天线阵元采集卫星信号的初始数据；
本实施例以阵列天线卫星导航接收机为研究对象，考虑非高斯噪声环境下，卫星信号波入射角度与信号频率联合估计方法。卫星导航接收机在工作时，一般能观察到6-10颗可视卫星，可视卫星入射信号角度在天线可接收信号空间范围是稀疏的。此外，由于载体与不同卫星之间的相对运动不同，因此，不同卫星信号到达导航接收机处的多普勒频率值是不一样的，位于-10kHz到10kHz之间。由于卫星信号是单频信号，其在频率域上只具有单个明显的峰值，即卫星信号在频率域也可以认为是具有稀疏特性的。因此，本实施例利用卫星信号空域和频率域稀疏特性，设计优化目标函数，进行卫星信号角度与频率联合估计。
卫星导航接收机采用均匀线性阵列天线，天线阵元数目为M，相邻阵元间距为d，每个天线阵元接收到L个信号，对每个信号进行K个采样点的时域采样，得到卫星信号的初始数据，可以理解的是初始数据是包括非高斯噪声的卫星信号。初始数据中的角度区间定义为Θ，被等间隔分成J个角度。
步骤S20，将所述初始数据代入信号函数中，求解得到卫星信号；
本实施例中用Y表示初始数据一般的，阵列天线接收到的初始数据Y模型可以表示为Y＝AX+N，X表示的是不含非高斯噪声的卫星信号，即需要求的卫星信号，是具有行稀疏结构的J×K维矩阵，其中包含L个实际卫星信号，其他部分无数据。阵列天线接收到的初始数据Y＝[y(1) y(2)…y(K)]为M×K维矩阵y(k)，i＝1，2，…，K，表示第k个采样时刻天线阵列采集到的包含非高斯噪声的数据，即初始数据，N为非高斯噪声，为M×N维矩阵。A为M×J维导向矢量矩阵，具体可以表示为A＝[a(θ1) a(θ2)…a(θJ)]，其中，a(θj)是角度为θj的导向矢量，表示为：
Figure BDA0002855967630000051
(1)式中，λ表示卫星信号波长。
本实施例的目的是根据阵列天线所采集到的初始数据，准确恢复卫星信号，从而估计出卫星信号角度和频率。本发明同时考虑非高斯噪声背景条件和卫星信号的空域、频率域稀疏特性，建立如下函数：
Figure BDA0002855967630000061
(1)式中，(·)T表示转置操作，λ1和λ2表示正则化因子，F为K×K维傅里叶变换矩阵，||·||2，1表示l2，1范数，定义为：
Figure BDA0002855967630000062
式中，因为F为傅里叶变换矩阵，存在F＝FT特性。所以，可以进一步简化上式得最终的信号函数
Figure BDA0002855967630000063
将初始数据代入(4)式中，求解得到卫星信号。
在一些实施例中，步骤S20还包括：
步骤a，对所述信号函数进行重加权处理，得到目标函数；
步骤b，将所述初始数据代入所述目标函数，求解得到卫星信号。
在上述信号函数中，空域约束项||X||2，1和频率域约束项||X F||2，1均采用l2，1范数表示。为了使得这两项具有更加稀疏的特性，引入重加权方法，将信号函数转化为目标函数，得到：
Figure BDA0002855967630000064
上式中，Xi表示矩阵X的第i行数据，(XF)i表示矩阵(X F)的第i行数据，||·||2表示l2范数，ρi和κi分别为加权向量ρ和κ中的第i个元素，i＝1，2，…，J，其中，ρ和κ均为J×1维加权向量。
将初始函数代入目标函数中，得到卫星信号。
在一些实施例中，步骤b还包括：
步骤b1，将所述初始数据代入所述目标函数，采用梯度求导的方法，求解所述目标函数，得到卫星信号。
令上述目标函数对矩阵X求导并令倒数等于0，有
Figure BDA0002855967630000065
(6)式中，(·)H表示转置操作，矩阵C，D和G均为对角矩阵，它们对角线上的元素分别可以表示为
Figure BDA0002855967630000071
Figure BDA0002855967630000072
Figure BDA0002855967630000073
(9)式中，(AX-Y)i表示矩阵(AX-Y)的第i行数据，ρi和κi分别为：
Figure BDA0002855967630000074
Figure BDA0002855967630000075
(11)式中，η1和η2均为大于0的较小常数值，为了防止上式分母为0。
在式(6)中，因为F为傅里叶变换矩阵，则有：
Figure BDA0002855967630000076
(12)式中，K为傅里叶变换矩阵F的维度。因此，式(6)可以简写为
Figure BDA0002855967630000077
(13)式可以进一步调整为：
Figure BDA0002855967630000078
根据(14)式，很容易得到：
Figure BDA0002855967630000079
步骤S30，根据所述卫星信号，获取所述卫星信号的角度和频率。
可以理解的是卫星信号X的每一行刚好对应阵列天线接收信号Θ的一个角度，还需要根据卫星信号得到卫星信号的频率。
在一些实施例中，步骤S30还包括：
步骤c，对所述卫星信号进行傅里叶变换，得到所述卫星信号的频率。
根据图3可知，所恢复的卫星信号X的每一行刚好对应阵列天线接收信号Θ的一个角度。根据X的稀疏结构，存在数据的行所对应的角度，即为卫星信号的入射角度，并且该行的数据就是卫星信号的采样值。再根据附图4，对每一行数据进行傅里叶变换，便可以找到卫星信号的频率。所以，根据所恢复的卫星信号X，便可以同时得到卫星信号的入射角度和频率。
图3为卫星信号角度与频率联合估计方法原理示意图，其中，图3为空域与时域示意图。图3中，纵向区间代表阵列天线接收信号角度范围，横向区间代表对接收到的卫星信号进行时域采样。彩色区域代表有接收到的卫星信号，目标所在的角度具有空间稀疏特性。图4为空域与频率域示意图，图4中纵向区间仍表示阵列天线接收信号角度范围，横向区间代表不同卫星信号经过傅里叶变换后的频率域数据。由于不同卫星与导航接收机间相对运动产生的多普勒频率不同，因此，不同卫星信号经过傅里叶变换后的频率也是不同的。如附图4所示，每行均有一个数值大的数据，该数据所对应的频率即为该颗卫星信号的频率。每行中数值大的数据所处位置不同，即对应的信号频率不同。
本实施例通过天线阵元采集卫星信号的初始数据，将初始数据代入信号函数中求解得到卫星信号，根据卫星信号获取卫星信号的角度和频率，考虑在非高斯噪声环境下，同时考虑卫星信号空域稀疏特性和频率域稀疏特性，进行卫星信号角度和频率联合估计，有利于卫星导航接收机后续捕获跟踪卫星信号。
本实施例以阵列天线导航接收机为研究对象，考虑在非高斯噪声环境下，充分利用卫星信号空域和频率域稀疏特性，设计优化目标函数，进行卫星信号角度与频率联合估计，使得卫星导航接收机天线能够在卫星信号方向形成波束，增强信号接收增益；同时，从恢复信号中提取卫星信号的频率，为导航接收机后续基带信号处理中的频率锁定环路提供先验信息。
对本实施例提出的方法进行测试。均匀线性天线阵列具有80个天线阵元，相邻阵元间距d为卫星信号半波长，卫星接收信号角度区间为Θ＝[-20°，20°]，以1°等间隔划分，即J＝41，假定有L＝4个卫星信号，入射角度在Θ区间随机分布，信号采样点K＝50，正则化因子λ1＝λ2＝0.1，噪声采用混合高斯模型，信噪比为20dB，公式(11)和(12)中参数η1和η2均为0.2。由于本实施例需要采用迭代算法，初始化对角矩阵C，D和G为单位矩阵。图5给出了卫星信号估计值X的均方误差(Mean Square Error,MSE)随着迭代次数仿真图。从图5中可以看出，随着迭代次数的增加，卫星信号估计值X的MSE逐渐降低。当迭代次数到第7次时，MSE值收敛至3×10-5，具有很好估计性能。
此外，本发明实施例还提出一种信号角度和信号频率的估计装置，所述信号角度和信号频率的估计装置包括：
采集模块，用于通过天线阵元采集卫星信号的初始数据；
代入求解模块，用于将所述反馈数据发送至所述共享平台；
获取模块，用于根据所述卫星信号，获取所述卫星信号的角度和频率本发明所述信号角度和信号频率的估计装置实施方式与上述信号角度和信号频率的估计各实施例基本相同，在此不再赘述。
此外，本发明实施例还提出一种计算机可读存储介质，所述计算机可读存储介质上存储有信号角度和信号频率的估计程序，所述信号角度和信号频率的估计程序被处理器执行时实现如上所述的信号角度和信号频率的估计方法的各个步骤。
需要说明的是，计算机可读存储介质可设置在信号角度和信号频率的估计设备中。
本发明计算机可读存储介质具体实施方式与上述信号角度和信号频率的估计方法各实施例基本相同，在此不再赘述。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其它变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其它要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者设备中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其它相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (9)

1.一种信号角度和信号频率的估计方法，其特征在于，所述信号角度和信号频率的估计方法包括以下步骤：
通过天线阵元采集卫星信号的初始数据；
将所述初始数据代入信号函数中，求解得到卫星信号；
根据所述卫星信号，获取所述卫星信号的角度和频率。
2.如权利要求1所述的信号角度和信号频率的估计方法，其特征在于，将所述初始数据代入信号函数中，求解得到卫星信号的步骤包括：
对所述信号函数进行重加权处理，得到目标函数；
将所述初始数据代入所述目标函数，求解得到卫星信号。
3.如权利要求2所述的信号角度和信号频率的估计方法，其特征在于，所述将所述初始数据代入所述目标函数，求解得到卫星信号的步骤包括：
将所述初始数据代入所述目标函数，采用梯度求导的方法，求解所述目标函数，得到卫星信号。
4.如权利要求1所述的信号角度和信号频率的估计方法，所述卫星信号的采样值表示所述卫星信号的入射角度，其特征在于，所述根据所述卫星信号，获取所述卫星信号的角度和频率的步骤包括：
对所述卫星信号进行傅里叶变换，得到所述卫星信号的频率。
5.如权利要求1所述的信号角度和信号频率的估计方法，其特征在于，所述信号函数的表达式为：
Figure FDA0002855967620000011
其中，所述X表示所述卫星信号，所述Y表示初始数据，所述A表示导向矢量矩阵，所述λ1所述λ2表示正则化因子，所述F表示傅里叶变换。
6.如权利要求2所述的信号角度和信号频率的估计方法，其特征在于，所述目标函数的表达式为：
Figure FDA0002855967620000021
其中，ρi和κi分别为加权向量ρ和κ中的第i个元素，i＝1，2，…，J，所述J为所述A的列数。
7.一种信号角度和信号频率的估计装置，其特征在于，所述信号角度和信号频率的估计装置包括：
采集模块，用于通过天线阵元采集卫星信号的初始数据；
代入求解模块，用于将所述反馈数据发送至所述共享平台；
获取模块，用于根据所述卫星信号，获取所述卫星信号的角度和频率。
8.一种信号角度和信号频率的估计设备，其特征在于，所述信号角度和信号频率的估计设备包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的信号角度和信号频率的估计程序，所述信号角度和信号频率的估计程序被所述处理器执行时实现如权利要求1至6中任一项所述的信号角度和信号频率的估计的步骤。
9.一种存储介质，其特征在于，所述存储介质为计算机可读存储介质，所述存储介质上存储有信号角度和信号频率的估计程序，所述信号角度和信号频率的估计程序被处理器执行时实现如权利要求1至6中任一项所述的信号角度和信号频率的估计方法的步骤。
CN202011556223.8A 2020-12-24 2020-12-24 信号角度和信号频率的估计方法、装置、设备及存储介质 Active CN112731486B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202011556223.8A CN112731486B (zh) 2020-12-24 2020-12-24 信号角度和信号频率的估计方法、装置、设备及存储介质
PCT/CN2021/123920 WO2022134762A1 (zh) 2020-12-24 2021-10-14 信号角度和信号频率的估计方法、装置、设备及存储介质
US18/041,734 US20230314622A1 (en) 2020-12-24 2021-10-14 Signal angle and signal frequency estimation method, apparatus and device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011556223.8A CN112731486B (zh) 2020-12-24 2020-12-24 信号角度和信号频率的估计方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN112731486A true CN112731486A (zh) 2021-04-30
CN112731486B CN112731486B (zh) 2023-07-04

Family

ID=75615686

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011556223.8A Active CN112731486B (zh) 2020-12-24 2020-12-24 信号角度和信号频率的估计方法、装置、设备及存储介质

Country Status (3)

Country Link
US (1) US20230314622A1 (zh)
CN (1) CN112731486B (zh)
WO (1) WO2022134762A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113835107A (zh) * 2021-09-22 2021-12-24 深圳大学 阵列卫星导航接收机的信号处理方法、装置及智能终端
WO2022134762A1 (zh) * 2020-12-24 2022-06-30 深圳大学 信号角度和信号频率的估计方法、装置、设备及存储介质
CN115792981A (zh) * 2023-02-06 2023-03-14 深圳大学 一种基于阵列天线的可见卫星探测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090262017A1 (en) * 2008-04-18 2009-10-22 Honeywell, Spol. S.R.O. Systems and methods for fast gnss signals acquisition
CN106597480A (zh) * 2016-12-08 2017-04-26 深圳大学 用于卫星导航rtk发射电台的抗干扰定位方法及系统
CN107493106A (zh) * 2017-08-09 2017-12-19 河海大学 一种基于压缩感知的频率和角度联合估计的方法
CN109975842A (zh) * 2019-03-04 2019-07-05 广东工业大学 一种基于小波变换的北斗卫星信号高精度盲捕获方法
CN110174656A (zh) * 2019-05-21 2019-08-27 电子科技大学 一种基于频域宽带波束形成的稀疏阵列设计方法及装置
CN110780273A (zh) * 2019-11-04 2020-02-11 电子科技大学 一种混合正则化方位超分辨成像方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3383235B2 (ja) * 1999-04-14 2003-03-04 三菱電機株式会社 人工衛星の姿勢制御装置
CN103605139B (zh) * 2013-12-10 2014-12-17 武汉大学 适用于gnss接收机的载波频率和相位估计方法及系统
CN104251997B (zh) * 2014-06-27 2016-04-20 广东工业大学 一种卫星信号抗干扰稀疏处理方法
CN106130571B (zh) * 2016-06-07 2018-12-21 烟台大学文经学院 一种基于频域和空域压缩感知的信号采样接收方法及装置
CN109425875B (zh) * 2017-09-05 2021-05-04 清华大学 卫星信号分离处理装置和方法
CN108600133B (zh) * 2018-04-19 2020-12-15 深圳大学 一种基于阵列天线的欠定相位恢复方法
CN108957387B (zh) * 2018-05-21 2022-06-21 西安电子科技大学 一种卫星信号二维到达角估计方法及系统
CN109143154A (zh) * 2018-07-24 2019-01-04 南京航空航天大学 一种应用于l型阵列的信号二维doa与频率联合估计方法
CN111679296B (zh) * 2020-05-28 2023-04-25 哈尔滨工业大学 一种基于gp迭代外推的单星定位方法
CN112731486B (zh) * 2020-12-24 2023-07-04 深圳大学 信号角度和信号频率的估计方法、装置、设备及存储介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090262017A1 (en) * 2008-04-18 2009-10-22 Honeywell, Spol. S.R.O. Systems and methods for fast gnss signals acquisition
CN106597480A (zh) * 2016-12-08 2017-04-26 深圳大学 用于卫星导航rtk发射电台的抗干扰定位方法及系统
CN107493106A (zh) * 2017-08-09 2017-12-19 河海大学 一种基于压缩感知的频率和角度联合估计的方法
CN109975842A (zh) * 2019-03-04 2019-07-05 广东工业大学 一种基于小波变换的北斗卫星信号高精度盲捕获方法
CN110174656A (zh) * 2019-05-21 2019-08-27 电子科技大学 一种基于频域宽带波束形成的稀疏阵列设计方法及装置
CN110780273A (zh) * 2019-11-04 2020-02-11 电子科技大学 一种混合正则化方位超分辨成像方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
QIANG LI等: "A_New_Satellite-to-Satellite_Passive_Locating_Method_Using_Frequency-Only_Measurements_Based_on_Particle_Swarm_Optimization", 《2007 IEEE INTERNATIONAL CONFERENCE ON CONTROL AND AUTOMATION GUANGZHOU,CHINA》, pages 1459 - 1464 *
李玮等: "阵列失效单元压缩感知诊断算法研究", 《中国博士学位论文全文数据库信息科技辑》, no. 1, pages 1 - 144 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022134762A1 (zh) * 2020-12-24 2022-06-30 深圳大学 信号角度和信号频率的估计方法、装置、设备及存储介质
CN113835107A (zh) * 2021-09-22 2021-12-24 深圳大学 阵列卫星导航接收机的信号处理方法、装置及智能终端
CN113835107B (zh) * 2021-09-22 2023-09-29 深圳大学 阵列卫星导航接收机的信号处理方法、装置及智能终端
CN115792981A (zh) * 2023-02-06 2023-03-14 深圳大学 一种基于阵列天线的可见卫星探测方法

Also Published As

Publication number Publication date
WO2022134762A1 (zh) 2022-06-30
US20230314622A1 (en) 2023-10-05
CN112731486B (zh) 2023-07-04

Similar Documents

Publication Publication Date Title
CN112731486B (zh) 信号角度和信号频率的估计方法、装置、设备及存储介质
Yang et al. Kalman filtering with nonlinear state constraints
Larew et al. Adaptive beam tracking with the unscented Kalman filter for millimeter wave communication
Mohanna et al. Optimization of MUSIC algorithm for angle of arrival estimation in wireless communications
Zhang et al. Robust Capon beamforming against large DOA mismatch
US7126533B2 (en) Direction-finding for multiple cochannel sources
Liu et al. 2-D DOA robust estimation of echo signals based on multiple satellites passive radar in the presence of alpha stable distribution noise
Dai et al. Robust Bayesian learning approach for massive MIMO channel estimation
Steinwandt et al. Beamspace direction finding based on the conjugate gradient and the auxiliary vector filtering algorithms
Wu et al. Computationally efficient DOA and carrier estimation for coherent signal using single snapshot and its time-delay replications
Kitavi et al. Lower bound of the estimation error of an emitter's direction‐of‐arrival/polarisation, for a collocated triad of orthogonal dipoles/loops that fail randomly
Wan et al. Multiband delay estimation for localization using a two-stage global estimation scheme
Viberg et al. Calibration in array processing
Li et al. Adaptive beamforming based on covariance matrix reconstruction by exploiting interferences' cyclostationarity
Teng et al. Variational bayesian multiuser tracking for reconfigurable intelligent surface aided mimo-ofdm systems
Xu et al. Tensor-based angle and range estimation method in monostatic FDA-MIMO radar
Li et al. Robust cyclic MUSIC algorithm for finding directions in impulsive noise environment
Tayem Real time implemenation for DOA estimation methods on Ni-Pxi platform
CN111427006A (zh) 一种基于相控阵的单通道空间谱测向方法
CN111241470A (zh) 基于自适应零陷展宽算法的波束合成方法及装置
Parhizgar et al. Mutual coupling compensation for a practical VHF/UHF Yagi‐Uda antenna array
CN113835107B (zh) 阵列卫星导航接收机的信号处理方法、装置及智能终端
Tian et al. Modified spectral function based DOA estimation in colored noise
CN110391820B (zh) 一种基于dft的规避同频干扰的通信接收方法
Guo et al. DOA estimation of coherent wideband signals based on extended TOPS algorithm

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