CN106597407A - 多星协同下微弱回波信号时延和多普勒频移联合估计方法 - Google Patents
多星协同下微弱回波信号时延和多普勒频移联合估计方法 Download PDFInfo
- Publication number
- CN106597407A CN106597407A CN201611111559.7A CN201611111559A CN106597407A CN 106597407 A CN106597407 A CN 106597407A CN 201611111559 A CN201611111559 A CN 201611111559A CN 106597407 A CN106597407 A CN 106597407A
- Authority
- CN
- China
- Prior art keywords
- tau
- alpha
- omega
- infin
- prime
- 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
Links
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种多星协同下微弱回波信号时延和多普勒频移联合估计方法,对参考通道中的多个直达波信号进行分离;对回波通道中的直达波信号和多径信号进行自适应抑制;对回波通道中的信号和不同的直达波信号分别进行基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数处理得到多个特征量;对模糊函数峰值进行谱峰提取后对应的坐标即为时延和多普勒频移的估计值;对多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,并通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移。当信噪比大于‑10dB时,信噪比估计的归一化均方误差小于1。
Description
技术领域
本发明属于通信技术与信号处理技术领域,尤其涉及一种多星协同下微弱回波信号时延和多普勒频移联合估计方法。
背景技术
随着电磁环境的复杂,传统的目标探测技术无法满足探测精度和反应能力的需求。因此,如何利用目标与电波传播环境、杂波环境等的相互作用,通过检测目标对电波环境的扰动来对目标进行探测具有重要的意义。在目标探测后,为了对目标速度等信息的获取和对目标的定位,需要对微弱回波信号进行参数估计。时延和多普勒频移是表征目标回波信号的重要特征参数,准确迅速的估计时延及多普勒信息可以进一步确定移动目标的距离、方位、速度和移动方向等信息。因此,需要对微弱回波的时延和多普勒频移进行精确估计,可实现对空中目标的跟踪定位,进而进行目标的识别。卫星通信作为现代通信的重要组成部分,在通信、导航、定位等领域占有着重要的地位。随着卫星通信的发展,基于卫星辐射源下目标探测技术越来越受到重视,一些学者展开了将卫星作为辐射源对目标的微弱回波信号进行检测的方法研究。但是,这些研究往往以单个卫星作为辐射源进行信号的分析和处理,在这样场景下的目标探测和微弱回波参数估计方法存在覆盖范围小、可靠性低、应用场景受限,易受地理环境和设备条件限制等缺点,并且接收机不可避免的同时接收到多个卫星的辐射信号,从而导致单个卫星辐射源下的目标探测和参数估计在该场景下的容易受到其它卫星辐射信号的干扰。为了进一步提高目标探测与微弱回波信号参数估计的准确性和可靠性,研究基于多个卫星辐射源下的目标探测和微弱回波信号参数估计方法具有研究意义和使用价值。Bilal Hammoud等人提出了一种基于张量滤波的GNSS回波信号的时延估计方法,但是过程中会涉及到对矩阵的处理,运算量较大(Bilal Hammoud,FelixAntreich,Josef A.Nossek,et al.Tensor-Based Approach for Time-Delay Estimation[C].Smart Antennas(WSA 2016);Proceedings ofthe 20th International ITGWorkshop on,2016,:1-7.)。Nicola Linty等人提出迫零和双快速傅里叶变换相结合的多普勒频移估计算法。但是由于噪声和干扰的影响,在信噪比较低时估计性能较低(NicolaLinty,Letizia Lo Presti.Doppler Frequency Estimation in GNSS Receivers Basedon Double FFT[J].IEEE Transactions on Vehicular Technology,2015,65(2):509-524)。Xu Yunxiang等人提出利用GPS卫星星历和光学锁相环估计多普勒频移的方法,但是该方法随着卫星轨道和仰角的变化比较敏感,具有较差的适应性和普遍性(Yunxiang Xu,Bin Wu.Doppler shift estimation using broadcast ephemeris in satelliteoptical communication[C],201625th Wireless and Optical CommunicationConference(WOCC),2016:1-4.)。Baharak Soltanian等人提出了一种基于快速傅里叶变换的GNSS回波信号多普勒频移估计方法,但是该方法没有考虑因为噪声的不确定性对估计结果产生的影响。(Baharak Soltanian1,AliMuratDemirtas,Ali Shahed hagh ghadam,Markku Renfors.Reduced-complexity FFT-based method for Doppler estimation inGNSS receivers[J].Eurasip Journal on Applied Signal Processing,2014,2014(1):1-15.)。Christoph Enneking,Manuel Stein,Mario等人针对多GNSS辐射源下信号之间的相互干扰问题,提出了一种最大似然估计的时延估计方法并进行了相应的数学推导,为了减弱信号之间的相互干扰和噪声对估计结果的影响,依据期望最大原理构造了一个迭代函数,随着迭代次数的增加,信号之间的相互干扰减弱,但是计算复杂度也会增加(Christoph Enneking,Manuel Stein,Marioet al.Multi-satellite time-delay estimation for reliable high-resolution GNSS receivers[C].PositionLocation and Navigation Symposium(PLANS),2012IEEE/ION,2012:488-494.)。
综上所述,现有基于多个卫星辐射源下的目标探测和微弱回波信号参数估计方法存在低信噪比环境下多个卫星的微弱回波信号的时延和多普勒频移的估计性能较差的问题。
发明内容
本发明的目的在于提供一种多星协同下微弱回波信号时延和多普勒频移联合估计方法,旨在解决现有基于多个卫星辐射源下的目标探测和微弱回波信号参数估计方法存在低信噪比环境下多个卫星的微弱回波信号的时延和多普勒频移的估计性能较差的问题。
本发明是这样实现的,一种多星协同下微弱回波信号时延和多普勒频移联合估计方法,所述多星协同下微弱回波信号时延和多普勒频移联合估计方法包括以下步骤:
步骤一,利用多个不同的带通滤波器对参考通道中的多个直达波信号进行分离;根据正交投影矩阵的机理对回波通道中的直达波信号和多径信号进行自适应抑制;
步骤二,通过对回波通道中的信号和不同的直达波信号分别进行基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数处理得到多个特征量;通过对模糊函数峰值进行谱峰提取后对应的坐标即为时延和多普勒频移的估计值;
步骤三,对多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,并通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移。
进一步,所述带通滤波器的***函数为:
其中,ak和br为带通滤波器的系数,ωp为带通滤波器的通带上边界频率,ωpu为通带下边界频率,ωs为阻带上边界频率,ωsu为阻带下边界频率,αp为通带最大衰减,αs为阻带最小衰减。
进一步,经过分离之后,参考通道中的多个直达波信号xi(t)表示为:
xi(t)=risi(t)+n(t);
监测通道中的回波信号z(t)模型描述为:
其中,M为不同体制卫星下的回波信号的个数,n(t),n′(t)为均值为零的平稳高斯白噪声,si(t)是不同的直达波信号,ri为参考通道中的不同直达波信号的幅度,rη′为回波信号的幅度,Dη为回波信号相对应直达波的时延,为回波信号相对应直达波的多普勒频移,Ωη为监测通道中的直达波信号的幅度,H为监测通道中多径信道的径数,ωηj为第η个卫星的直达波信号经过第j条径之后的幅度,τηj为第η个卫星的直达波信号经过第j条多径之后的时延。
进一步,对监测通道中的直达波和多径信号的抑制采用并行序列抵消算法,具体的步骤如下:
首先对个个直达波信号和监测通道信号进行采样,采样后的信号表示如下:
直达波信号:
xi=[xi[-R+1],xi[-R+2],......,xi[N-1]]T;
其中i=1,2,3,......,M,M为接收机同时可以接收到的直达波信号的个数。xi代表不同的直达波信号,N为采样点数,R为参考信号相对于监测通道信号多出来的采样点;xi[·]表示第i个直达波的不同采样点;
监测通道信号:
z=[z[0],z[1],......,z[N-1]]T;
令NB=N/b,其中N表示原始信号的总长度,b表示信号分段后的段数,NB表示分段后每段信号的长度,分段后的信号可以表示如下:
进行直达波和杂波抑制之后的检测通道信号表示为:
y(t)为监测通道中抑制直达波和多径干扰后的信号,并且y(t)只包含回波信号和噪声;y(t)表示为:
进一步,基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数的处理得到不同的特征量具体包括:
首先,对经过分离的直达波信号xi(t)进行四阶自循环累积量的计算,计算表达式为:
其中αi为循环频率,是信号si(t)的四阶自循环矩和二阶自循环矩,其分别表示为:
其中T为周期,是有限的时间长度;
然后对经过分离后的直达波信号xi(t)和经过直达波和多径抑制的回波信号y(t)求四阶互循环累积量表示如下:
其中B=E[xi(t)y(t)],αi-f为xi(t)与y(t)的互四阶循环频率,为第i个si(t)信号的循环自相关,其表示为:
为si(t)与参考通道噪声n(t)的互四阶循环矩,为si(t)与噪声n′(t)的互四阶循环矩,其表示为:
将和带入下式四阶循环累积量互模糊函数:
得到四项加权的分数傅里叶变换如下:
其中β为分数傅里叶变换的角度,g0(x)~g3(x)分别为函数g(x)的0~3次傅里叶变换,加权系数ω0(β)~ω3(β)定义如下:
原模糊函数的自变量按照式下式变换即可:
AFβ(t,w)=AF(t cosβ-ωsinβ,t sinβ+ωcosβ);
当模糊函数为四阶循环累积量的互模糊函数时,令按照上就可以得到模糊函数0~3次的傅里叶变换。
进一步,
①当β=0时:
②当β=π/2时:
③当β=π时:
④当β=3π/2时:
当变换角度β确定后,分数傅里叶变换系数ω0(β)~ω3(β)也确定,则得到基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数表示如下:
进一步,通过对模糊函数峰值进行谱峰提取后得到多组时延和多普勒频移的初始估计值具体包括:
第一步,把模糊函数截面的数据存入数组W,其中第i个元素记为wi:
第二步,对数组W中的元素进行后向差分处理,得到新的数组元素,并将新的数组元素存入数组W′,其中,第j个元素记为wj′,其表达式为:
wj′=wi-wi-1;
第三步,在得到包含了谱峰特征的数组后,对大于零的数组元素进行如下变换:
其中n为数组中所有大于零的元素的个数;
第四步,对由第三步得到的数组中所有大于零的数组元素求和以及求平均值,把小于平均值的数组元素置零;
第五步,对第四步进行循环操作,直到得到的数组元素的和等于其中一个数组元素w″k,循环结束,则数组元素w″k的下标k+1就是二维截面中谱峰处对的横坐标;
第五步,将多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,再通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移。
进一步,所述第五步具体包括:将不同的时延-多普勒谱经过TF-LV变换使检测峰值坐标相同,使多个不同卫星信号下目标的的距离-速度统一,具体包括:
δ为双基地夹角平分线与飞机速度v之间的夹角,β为双基地角;τ与L以及fd与v的关系如下:
求解上式可得:
根据下式:
可得:
式中λ=1/f,β可由下式求得:
经过TF-LV的变换后,对应着多对L和v的值,为了得到更加准确和可靠的估计结果,需要对多对L、v进行加权融合。
进一步,所述加权融合为:按照一定的准则,为每一个估计值分配对应的加权系数;由所有的局部估计和加权系数得到全局估计;对于变换之后得到的L1、L2、L3和v1、v2、v3分别使用如下的加权数据融合算法进行融合:
首先,对于L1、L2、L3,对应的加权系数分别为w1、w2、w3,则融合后的应该满足如下条件:
根据上式得到不同估计值总的均方误差和为:
因为当i≠j时,Li与Lj相互独立,所以有以下等式成立:
因此可得:
其中Pi=E[(L-Li)2],为每个估计值的方差;
总均方误差是以各个加权系数为变量的多元二次函数,存在极小值:
通过对上式进行求解,可得当取得最小值时的加权系数为:
可得此时均方误差和的最小值为:
得到L1、L2、L3的加权系数后,求出经过融合后的值;速度的融合过程与距离的融合过程相同;
经过多对距离L和速度v加权融合,得到一组和通过下式可得时延的估计值:
频偏的估计值:
通过将每个卫星对应的θi和βi带入上面两式,即可估计出对应的时延τi′和多普勒
本发明的另一目的在于提供一种利用所述多星协同下微弱回波信号时延和多普勒频移联合估计方法的卫星通信***。
本发明提供的多星协同下微弱回波信号时延和多普勒频移联合估计方法,可以对低信噪比环境下,多个卫星协同下微弱回波信号的时延和多普勒频移具有良好的估计性能。当信噪比大于-10dB时,信噪比估计的归一化均方误差小于1。由此可见,本发明在低信噪比条件下,对多星协同下微弱回波信号的时延和多普勒频移具有良好的估计性能;利用本发明技术可以利用多个卫星辐射源信号对运动目标进行跟踪定位。
附图说明
图1是本发明实施例提供的多星协同下微弱回波信号时延和多普勒频移联合估计方法流程图。
图2是本发明实施例提供的实施例1的流程图。
图3是本发明实施例提供的对多个卫星辐射源信号在不同信噪比下的多普勒频移和时延估计性能示意图。
图4是本发明实施例提供的接收***的几何结构示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示,本发明实施例提供的多星协同下微弱回波信号时延和多普勒频移联合估计方法包括以下步骤:
S101:利用多个不同的带通滤波器对参考通道中的多个直达波信号进行分离;根据正交投影矩阵的机理对回波通道中的直达波信号和多径信号进行自适应抑制;
S102:通过对回波通道中的信号和不同的直达波信号分别进行基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数处理得到多个特征量;通过对模糊函数峰值进行谱峰提取后对应的坐标即为时延和多普勒频移的估计值;
S103:对多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,并通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移。
下面结合具体实施例对本发明的应用原理作进一步的描述。
如图2所示,本发明为多星协同下微弱回波信号的时延和多普勒频移的联合估计方法,所述方法包括以下步骤:
S1利用多个不同的数字带通滤波器对参考通道中的多个直达波信号进行分离。
需要说明的是,步骤S1中利用多个不同的数字带通滤波器对参考通道中的多个直达波信号进行分离原理为:
假设参考通道中接收机同时接收到多个不同的直达波信号,由于不同卫星的直达波信号的频率不同,因此可以设计多个不同的带通滤波器将多个直达波信号同时并行分离,其中,带通滤波器的***函数为:
其中,ak和br为带通滤波器的系数,ωp为带通滤波器的通带上边界频率,ωpu为通带下边界频率,ωs为阻带上边界频率,ωsu为阻带下边界频率,αp为通带最大衰减,αs为阻带最小衰减。根据所占用不同的频谱范围设置不同滤波器的参数,并确定不同带通滤波器***函数中Q,S,ak和br的值,然后将多个带通滤波器并联以同时对多个不同直达波信号分离,得到多个不同的直达波信号。经直达波信号分离后,在参考通道中的信号为相互独立的多个直达波信号。
S2利用基于正交投影矩阵的自适应滤波方法对回波通道中的直达波和多径进行自适应抑制;
需要说明的是,步骤S2中利用基于正交投影矩阵的自适应滤波方法对回波通道中的直达波和多径进行自适应抑制原理为:
经过分离之后,参考通道中的多个直达波信号xi(t)表示为:
xi(t)=risi(t)+n(t);
监测通道中的回波信号z(t)模型描述为:
其中,M为不同体制卫星下的回波信号的个数,n(t),n′(t)为均值为零的平稳高斯白噪声,si(t)是不同的直达波信号,ri为参考通道中的不同直达波信号的幅度,rη′为回波信号的幅度,Dη为回波信号相对应直达波的时延,为回波信号相对应直达波的多普勒频移,Ωη为监测通道中的直达波信号的幅度,H为监测通道中多径信道的径数,ωηj为第η个卫星的直达波信号经过第j条径之后的幅度,τηj为第η个卫星的直达波信号经过第j条多径之后的时延。
对监测通道中的直达波和多径信号的抑制可以采用并行序列抵消算法,该算法的步骤如下:
首先对个个直达波信号和监测通道信号进行采样,采样后的信号表示如下:
直达波信号:
xi=[xi[-R+1],xi[-R+2],......,xi[N-1]]T;
其中i=1,2,3,......,M,M为接收机同时可以接收到的直达波信号的个数。xi代表不同的直达波信号,N为采样点数,R为参考信号相对于监测通道信号多出来的采样点。xi[·]表示第i个直达波的不同采样点。
监测通道信号:
z=[z[0],z[1],......,z[N-1]]T;
由于原始信号长度较长,为了解决对较长信号进行矩阵运算时计算量大的问题,对原始直达波信号和监测通道信号进行分段处理,令NB=N/b,其中N表示原始信号的总长度,b表示信号分段后的段数,NB表示分段后每段信号的长度,分段后的信号可以表示如下:。
其中
zk=[zk[kNB],......,zk[(k+1)NB-1]]T;
k=0,1,......,b-1。
则进行直达波和杂波抑制之后的检测通道信号(即回波信号)可以表示为:
其中由下式计算
其中:
PM=IN;
Pj=Qj+1Qj+2……QM;
其中xj-1为矩阵Xk的第j列:
Tmax为最大时延,其中Tmax可由最大探测距离Lmax得到,即Tmax=Lmax/c。通过以上算法就可以得到对监测通道中直达波和多径信号抑制之后的回波信号zMSCA-B。
设y(t)为监测通道中抑制直达波和多径干扰后的信号,并且y(t)只包含回波信号和噪声。y(t)可以表示为:
S3基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数的处理得到不同的特征量;
需要说明的是,步骤S3中基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数的处理得到不同的特征量原理为:
首先,对经过分离的直达波信号xi(t)进行四阶自循环累积量的计算,计算表达式为:
其中αi为循环频率,是信号si(t)的四阶自循环矩和二阶自循环矩,其分别表示为:
其中T为周期,是有限的时间长度。
然后对经过分离后的直达波信号xi(t)和经过直达波和多径抑制的回波信号y(t)求四阶互循环累积量表示如下:
其中B=E[xi(t)y(t)],αi-f为xi(t)与y(t)的互四阶循环频率,为第i个si(t)信号的循环自相关,其表示为:
为si(t)与参考通道噪声n(t)的互四阶循环矩,为si(t)与噪声n′(t)的互四阶循环矩,其表示为:
将和带入下式四阶循环累积量互模糊函数:
如果对函数进行0次、1次、2次、3次傅里叶变换,再进一步乘以加权系数求和,可以得到分数傅里叶变换的结果。从傅里叶变换的基本性质出发,结合傅里叶变换以4为周期的特点,可以得到四项加权的分数傅里叶变换如下:
其中β为分数傅里叶变换的角度,g0(x)~g3(x)分别为函数g(x)的0~3次傅里叶变换,加权系数ω0(β)~ω3(β)定义如下:
模糊函数的分数傅里叶变换相当于对模糊函数在原坐标系进行旋转变换到分数域,在表达式上的体现为原模糊函数的自变量按照式下式变换即可:
AFβ(t,w)=AF(t cosβ-ωsinβ,t sinβ+ωcosβ);
当模糊函数为四阶循环累积量的互模糊函数时,令按照上就可以得到模糊函数0~3次的傅里叶变换。
①当β=0时:
②当β=π/2时:
③当β=π时:
④当β=3π/2时:
当变换角度β确定后,分数傅里叶变换系数ω0(β)~ω3(β)也可以确定,则可以得到基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数表示如下:
S4通过对模糊函数峰值进行谱峰提取后得到多组时延和多普勒频移的初始估计值;
需要说明的是,步骤S4中通过对模糊函数峰值进行谱峰提取后得到多组时延和多普勒频移的初始估计值原理为:
Step1:把模糊函数截面的数据存入数组W,其中第i个元素记为wi:
Step2:对数组W中的元素进行后向差分处理,得到新的数组元素,并将新的数组元素存入数组W′,其中,第j个元素记为wj′,其表达式为:
wj′=wi-wi-1;
当wj′<0时,取wj′=0。经过后向差分运算后,谱峰处的突变特性以一个突变的最大值的形式存储在数组中,所以在这个突变的数组元素处能够体现出谱峰的特征。
Step3:在得到包含了谱峰特征的数组后,对大于零的数组元素进行如下变换:
其中n为数组中所有大于零的元素的个数。上式通过一个减法操作使所有的数组元素都不同程度的减小,但是较大的数组元素减小的少,较小的数组元素减小的多。因为包含了谱峰特征的数组元素为数组元素的最大值,所以通过上式换后,这个数组元素减小的最少,几乎不变,这使数组中更加明显的体现出谱峰的信息。式中减去含有均值平方项的作用为使数组中各个较小的元素更快的收敛到零,从而使谱峰的特征更加突出。
Step4:对由Step3得到的数组中所有大于零的数组元素求和以及求平均值,把小于平均值的数组元素置零。
Step5:对Step4进行循环操作,直到得到的数组元素的和等于其中一个数组元素w″k,循环结束,则数组元素w″k的下标k+1就是二维截面中谱峰处对的横坐标。
S5将多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,再通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移。
步骤S5中将多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,再通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移原理为:
不同卫星分布的位置不同,因此通过互模糊函数的峰值所对应的时延τ和多普勒频移fd不同,但是不同的τ和fd所对应目标的L和速度v是相同的。对于这个问题,将不同的时延-多普勒谱经过TF-LV变换使检测峰值坐标相同,从而使多个不同卫星信号下目标的的距离-速度统一。
图4是接收***的几何结构,其中,δ为双基地夹角平分线与飞机速度v之间的夹角,β为双基地角(即卫星与目标连线和接收天线与目标连线之间的夹角)。τ与L以及fd与v的关系如下:
求解上式可得:
根据下式:
可得:
式中λ=1/f,β可由下式求得:
经过TF-LV的变换后,对应着多对L和v的值,为了得到更加准确和可靠的估计结果,需要对多对L、v进行加权融合。加权融合的思想为:按照一定的准则,为每一个估计值分配对应的加权系数。最后,由所有的局部估计和加权系数得到全局估计。对于变换之后得到的L1、L2、L3和v1、v2、v3分别使用如下的加权数据融合算法进行融合:
首先,对于L1、L2、L3,假设其对应的加权系数分别为w1、w2、w3,则融合后的应该满足如下条件:
根据上式得到不同估计值总的均方误差和为:
因为当i≠j时,Li与Lj相互独立,所以有以下等式成立:
因此可得:
其中Pi=E[(L-Li)2],为每个估计值的方差。
从上式可以看出,总均方误差是以各个加权系数为变量的多元二次函数,所以一定存在极小值:
通过对上式进行求解,可得当取得最小值时的加权系数为:
可得此时均方误差和的最小值为:
得到L1、L2、L3的加权系数后,可求出经过融合后的值。速度的融合过程与距离的融合过程相同。
经过多对距离L和速度v加权融合,得到一组和包含了更加准确的运动目标的信息,通过下式可得时延的估计值:
频偏的估计值:
通过将每个卫星对应的θi和βi带入上面两式,即可估计出对应的时延τi′和多普勒
下面结合仿真实验对本发明的应用效果作详细的描述。
为了评估方法的性能,下面的仿真实验采用信号的类型为C/A码调制的BPSK信号、BPSK信号和QPSK信号,并进行2000次Monte Carlo实验。检测的评估标准为归一化均方误差(NMSE)。
为了测试本方法的检验统计量的性能,下面的仿真实验采用信号的类型为三个不同体制卫星(地球同步卫星,全球导航卫星,地球同步卫星)下的微弱回波信号,参数设置如下:采样频率设为 采样持续时间为100ms,GPS卫星直达波信号中载频1575.42MHZ,DVB_S信号的载频为12.38GHZ,卫星移动信号的载频为14.5GHZ,相对于直达波的回波信号的时延为τ1=10ms、τ2=22ms、τ3=34ms,相对于直达波的回波信号的频偏为分数傅里叶变换的角度β=45°,直达波信号与回波信号的功率比为40dB。仿真结果如图3所示,当信噪比SNR≥-10dB时,多普勒频移估计归一化最小均方误差达到10-4;当信噪比SNR≥-15dB时,时延估计归一化最小均方误差达到10-4,说明本发明的多星协同下微弱回波信号的时延和多普勒频移的联合估计方法是有效可行的。由此说明本发明方法在低信噪比条件下,对多星协同下微弱回波信号的时延和多普勒频移的联合估计具有较好的估计性能。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,所述多星协同下微弱回波信号时延和多普勒频移联合估计方法包括以下步骤:
步骤一,利用多个不同的带通滤波器对参考通道中的多个直达波信号进行分离;根据正交投影矩阵的机理对回波通道中的直达波信号和多径信号进行自适应抑制;
步骤二,通过对回波通道中的信号和不同的直达波信号分别进行基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数处理得到多个特征量;通过对模糊函数峰值进行谱峰提取后对应的坐标即为时延和多普勒频移的估计值;
步骤三,对多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,并通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移。
2.如权利要求1所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,所述带通滤波器的***函数为:
其中,ak和br为带通滤波器的系数,ωp为带通滤波器的通带上边界频率,ωpu为通带下边界频率,ωs为阻带上边界频率,ωsu为阻带下边界频率,αp为通带最大衰减,αs为阻带最小衰减。
3.如权利要求1所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,经过分离之后,参考通道中的多个直达波信号xi(t)表示为:
xi(t)=risi(t)+n(t);
监测通道中的回波信号z(t)模型描述为:
其中,M为不同体制卫星下的回波信号的个数,n(t),n′(t)为均值为零的平稳高斯白噪声,si(t)是不同的直达波信号,ri为参考通道中的不同直达波信号的幅度,rη′为回波信号的幅度,Dη为回波信号相对应直达波的时延,为回波信号相对应直达波的多普勒频移,Ωη为监测通道中的直达波信号的幅度,H为监测通道中多径信道的径数,ωηj为第η个卫星的直达波信号经过第j条径之后的幅度,τηj为第η个卫星的直达波信号经过第j条多径之后的时延。
4.如权利要求1所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,对监测通道中的直达波和多径信号的抑制采用并行序列抵消算法,具体的步骤如下:
首先对个个直达波信号和监测通道信号进行采样,采样后的信号表示如下:
直达波信号:
xi=[xi[-R+1],xi[-R+2],......,xi[N-1]]T;
其中i=1,2,3,......,M,M为接收机同时接收到的直达波信号的个数;xi代表不同的直达波信号,N为采样点数,R为参考信号相对于监测通道信号多出来的采样点;xi[·]表示第i个直达波的不同采样点;
监测通道信号:
z=[z[0],z[1],......,z[N-1]]T;
令NB=N/b,其中N表示原始信号的总长度,b表示信号分段后的段数,NB表示分段后每段信号的长度,分段后的信号表示如下:
进行直达波和杂波抑制之后的检测通道信号表示为:
y(t)为监测通道中抑制直达波和多径干扰后的信号,并且y(t)只包含回波信号和噪声;y(t)表示为:
5.如权利要求1所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数的处理得到不同的特征量具体包括:
首先,对经过分离的直达波信号xi(t)进行四阶自循环累积量的计算,计算表达式为:
其中αi为循环频率,是信号si(t)的四阶自循环矩和二阶自循环矩,其分别表示为:
其中T为周期,是有限的时间长度;
然后对经过分离后的直达波信号xi(t)和经过直达波和多径抑制的回波信号y(t)求四阶互循环累积量表示如下:
其中B=E[xi(t)y(t)],αi-f为xi(t)与y(t)的互四阶循环频率,为第i个si(t)信号的循环自相关,其表示为:
为si(t)与参考通道噪声n(t)的互四阶循环矩,为si(t)与噪声n′(t)的互四阶循环矩,其表示为:
将和带入下式四阶循环累积量互模糊函数:
得到四项加权的分数傅里叶变换如下:
其中β为分数傅里叶变换的角度,g0(x)~g3(x)分别为函数g(x)的0~3次傅里叶变换,加权系数ω0(β)~ω3(β)定义如下:
原模糊函数的自变量按照式下式变换即可:
AFβ(t,w)=AF(t cosβ-ωsinβ,t sinβ+ωcosβ);
当模糊函数为四阶循环累积量的互模糊函数时,令得到模糊函数0~3次的傅里叶变换。
6.如权利要求5所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,
①当β=0时:
②当β=π/2时:
③当β=π时:
④当β=3π/2时:
当变换角度β确定后,分数傅里叶变换系数ω0(β)~ω3(β)也确定,则得到基于四项加权分数傅里叶变换的四阶循环累积量互模糊函数表示如下:
7.如权利要求1所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,通过对模糊函数峰值进行谱峰提取后得到多组时延和多普勒频移的初始估计值具体包括:
第一步,把模糊函数截面的数据存入数组W,其中第i个元素记为wi:
第二步,对数组W中的元素进行后向差分处理,得到新的数组元素,并将新的数组元素存入数组W′,其中,第j个元素记为wj′,其表达式为:
wj′=wi-wi-1;
第三步,在得到包含了谱峰特征的数组后,对大于零的数组元素进行如下变换:
其中n为数组中所有大于零的元素的个数;
第四步,对由第三步得到的数组中所有大于零的数组元素求和以及求平均值,把小于平均值的数组元素置零;
第五步,对第四步进行循环操作,直到得到的数组元素的和等于其中一个数组元素w″k,循环结束,则数组元素w″k的下标k+1就是二维截面中谱峰处对的横坐标;
第五步,将多组时延和多普勒频移的估计值进行转化为目标到接收机的距离和速度进行数据融合,再通过数据融合后的距离和速度值再转化为多组时延和多普勒的估计值,从而估计出多星协同下微弱回波信号的时延和多普勒频移。
8.如权利要求7所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,所述第五步具体包括:将不同的时延-多普勒谱经过TF-LV变换使检测峰值坐标相同,使多个不同卫星信号下目标的的距离-速度统一,具体包括:
δ为双基地夹角平分线与飞机速度v之间的夹角,β为双基地角;τ与L以及fd与v的关系如下:
求解上式可得:
根据下式:
可得:
式中λ=1/f,β可由下式求得:
经过TF-LV的变换后,对应着多对L和v的值,为了得到更加准确和可靠的估计结果,需要对多对L、v进行加权融合。
9.如权利要求8所述的多星协同下微弱回波信号时延和多普勒频移联合估计方法,其特征在于,所述加权融合为:按照一定的准则,为每一个估计值分配对应的加权系数;由所有的局部估计和加权系数得到全局估计;对于变换之后得到的L1、L2、L3和v1、v2、v3分别使用如下的加权数据融合算法进行融合:
首先,对于L1、L2、L3,对应的加权系数分别为w1、w2、w3,则融合后的应该满足如下条件:
根据上式得到不同估计值总的均方误差和为:
因为当i≠j时,Li与Lj相互独立,所以有以下等式成立:
因此可得:
其中为每个估计值的方差;
总均方误差是以各个加权系数为变量的多元二次函数,存在极小值:
通过对上式进行求解,可得当取得最小值时的加权系数为:
可得此时均方误差和的最小值为:
得到L1、L2、L3的加权系数后,求出经过融合后的值;速度的融合过程与距离的融合过程相同;
经过多对距离L和速度v加权融合,得到一组和通过下式可得时延的估计值:
频偏的估计值:
通过将每个卫星对应的θi和βi带入上面两式,即可估计出对应的时延τi′和多普勒
10.一种利用权利要求1~9任意一项所述多星协同下微弱回波信号时延和多普勒频移联合估计方法的卫星通信***。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611111559.7A CN106597407B (zh) | 2016-12-06 | 2016-12-06 | 多星协同下微弱回波信号时延和多普勒频移联合估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611111559.7A CN106597407B (zh) | 2016-12-06 | 2016-12-06 | 多星协同下微弱回波信号时延和多普勒频移联合估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106597407A true CN106597407A (zh) | 2017-04-26 |
CN106597407B CN106597407B (zh) | 2019-05-07 |
Family
ID=58596073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611111559.7A Active CN106597407B (zh) | 2016-12-06 | 2016-12-06 | 多星协同下微弱回波信号时延和多普勒频移联合估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106597407B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107064903A (zh) * | 2017-05-18 | 2017-08-18 | 西安电子科技大学 | 一种多个异构卫星下基于glrt的目标探测方法 |
CN109031311A (zh) * | 2018-05-28 | 2018-12-18 | 西北工业大学 | 一种基于时间滑动窗处理的双基地大分置角主动探测方法 |
CN109188344A (zh) * | 2018-08-23 | 2019-01-11 | 北京邮电大学 | 脉冲噪声环境下基于互循环相关music算法信源个数与来波方向角估计方法 |
CN109633586A (zh) * | 2018-12-21 | 2019-04-16 | 电子科技大学 | 一种消除相位模糊度的时延估计方法 |
CN109991628A (zh) * | 2019-03-07 | 2019-07-09 | 深圳开阳电子股份有限公司 | 一种弱导航信号载噪比实时估计方法、装置及计算机存储介质 |
CN110365611A (zh) * | 2019-05-29 | 2019-10-22 | 北京邮电大学 | 一种多普勒频移估计方法及装置 |
CN110554378A (zh) * | 2019-09-10 | 2019-12-10 | 深圳锐越微技术有限公司 | 单通道多普勒雷达径向运动方向识别方法及装置 |
CN110831147A (zh) * | 2019-11-08 | 2020-02-21 | 北京邮电大学 | 一种适用于毫米波星间链路的载波频率同步方法 |
CN111198374A (zh) * | 2020-01-13 | 2020-05-26 | 哈尔滨工程大学 | 基于时空频联合干扰抑制的多普勒敏感信号动目标水声探测方法 |
CN111948619A (zh) * | 2020-06-30 | 2020-11-17 | 西安电子科技大学 | 一种多类型外辐射源照射下目标被动协同探测方法及*** |
CN111948620A (zh) * | 2020-06-30 | 2020-11-17 | 西安电子科技大学 | 一种基于多类型外辐射源的目标被动协同探测方法及*** |
WO2021251902A1 (en) * | 2020-06-12 | 2021-12-16 | Nanyang Technological University | Radar device and method for performing radar detection |
CN114338297A (zh) * | 2021-11-26 | 2022-04-12 | 河南工程学院 | 一种非相干LoRa***下的联合定时同步与频偏估计方法 |
CN117289256A (zh) * | 2023-11-24 | 2023-12-26 | 成都本原星通科技有限公司 | 一种基于低轨通信卫星的目标远距离高精度追踪方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006220532A (ja) * | 2005-02-10 | 2006-08-24 | Mazeran Systems Japan Kk | 衛星測位方法とそのシステム及び周波数誤差検出方法 |
CN105866750A (zh) * | 2016-05-23 | 2016-08-17 | 西安电子科技大学 | 一种多个gps卫星微弱回波信号检测的方法 |
CN106100769A (zh) * | 2016-05-23 | 2016-11-09 | 西安电子科技大学 | 一种多个不同体制卫星下微弱回波信号联合检测方法 |
-
2016
- 2016-12-06 CN CN201611111559.7A patent/CN106597407B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006220532A (ja) * | 2005-02-10 | 2006-08-24 | Mazeran Systems Japan Kk | 衛星測位方法とそのシステム及び周波数誤差検出方法 |
CN105866750A (zh) * | 2016-05-23 | 2016-08-17 | 西安电子科技大学 | 一种多个gps卫星微弱回波信号检测的方法 |
CN106100769A (zh) * | 2016-05-23 | 2016-11-09 | 西安电子科技大学 | 一种多个不同体制卫星下微弱回波信号联合检测方法 |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107064903B (zh) * | 2017-05-18 | 2020-06-09 | 西安电子科技大学 | 一种多个异构卫星下基于glrt的目标探测方法 |
CN107064903A (zh) * | 2017-05-18 | 2017-08-18 | 西安电子科技大学 | 一种多个异构卫星下基于glrt的目标探测方法 |
CN109031311A (zh) * | 2018-05-28 | 2018-12-18 | 西北工业大学 | 一种基于时间滑动窗处理的双基地大分置角主动探测方法 |
CN109188344A (zh) * | 2018-08-23 | 2019-01-11 | 北京邮电大学 | 脉冲噪声环境下基于互循环相关music算法信源个数与来波方向角估计方法 |
CN109188344B (zh) * | 2018-08-23 | 2020-11-27 | 北京邮电大学 | 脉冲噪声环境下基于互循环相关music算法信源个数与来波方向角估计方法 |
CN109633586A (zh) * | 2018-12-21 | 2019-04-16 | 电子科技大学 | 一种消除相位模糊度的时延估计方法 |
CN109633586B (zh) * | 2018-12-21 | 2022-07-12 | 电子科技大学 | 一种消除相位模糊度的时延估计方法 |
CN109991628A (zh) * | 2019-03-07 | 2019-07-09 | 深圳开阳电子股份有限公司 | 一种弱导航信号载噪比实时估计方法、装置及计算机存储介质 |
CN109991628B (zh) * | 2019-03-07 | 2022-10-11 | 深圳开阳电子股份有限公司 | 一种弱导航信号载噪比实时估计方法、装置及计算机存储介质 |
CN110365611A (zh) * | 2019-05-29 | 2019-10-22 | 北京邮电大学 | 一种多普勒频移估计方法及装置 |
CN110365611B (zh) * | 2019-05-29 | 2020-12-04 | 北京邮电大学 | 一种多普勒频移估计方法及装置 |
CN110554378A (zh) * | 2019-09-10 | 2019-12-10 | 深圳锐越微技术有限公司 | 单通道多普勒雷达径向运动方向识别方法及装置 |
CN110831147A (zh) * | 2019-11-08 | 2020-02-21 | 北京邮电大学 | 一种适用于毫米波星间链路的载波频率同步方法 |
CN111198374A (zh) * | 2020-01-13 | 2020-05-26 | 哈尔滨工程大学 | 基于时空频联合干扰抑制的多普勒敏感信号动目标水声探测方法 |
WO2021251902A1 (en) * | 2020-06-12 | 2021-12-16 | Nanyang Technological University | Radar device and method for performing radar detection |
CN111948619A (zh) * | 2020-06-30 | 2020-11-17 | 西安电子科技大学 | 一种多类型外辐射源照射下目标被动协同探测方法及*** |
CN111948620A (zh) * | 2020-06-30 | 2020-11-17 | 西安电子科技大学 | 一种基于多类型外辐射源的目标被动协同探测方法及*** |
CN111948620B (zh) * | 2020-06-30 | 2024-01-16 | 西安电子科技大学 | 一种基于多类型外辐射源的目标被动协同探测方法及*** |
CN111948619B (zh) * | 2020-06-30 | 2024-01-30 | 西安电子科技大学 | 一种多类型外辐射源照射下目标被动协同探测方法及*** |
CN114338297A (zh) * | 2021-11-26 | 2022-04-12 | 河南工程学院 | 一种非相干LoRa***下的联合定时同步与频偏估计方法 |
CN117289256A (zh) * | 2023-11-24 | 2023-12-26 | 成都本原星通科技有限公司 | 一种基于低轨通信卫星的目标远距离高精度追踪方法 |
CN117289256B (zh) * | 2023-11-24 | 2024-01-30 | 成都本原星通科技有限公司 | 一种基于低轨通信卫星的目标远距离高精度追踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106597407B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106597407A (zh) | 多星协同下微弱回波信号时延和多普勒频移联合估计方法 | |
CN103901395B (zh) | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 | |
CN103364772B (zh) | 基于实数域广义多重信号分类算法的目标低仰角估计方法 | |
CN105929378A (zh) | 基于外辐射源联合时延与多普勒频率的直接跟踪方法 | |
CN108449123B (zh) | 对地扩频通信***多目标检测、识别与二维角度估计方法 | |
CN104155662B (zh) | 基于gnss相关峰值探测器的自适应互干扰抑制方法 | |
CN104049262A (zh) | 一种基于矢量跟踪环路的波束形成抗干扰方法 | |
CN107272030B (zh) | 一种基于多个gps卫星信息融合的目标探测方法 | |
CN107255822A (zh) | 多径环境下gnss接收机信号参数估计方法 | |
CN106291543A (zh) | 一种运动平台扫描雷达超分辨成像方法 | |
Liang et al. | Multi-Bernoulli filter for target tracking with multi-static Doppler only measurement | |
RU2661357C1 (ru) | Способ обзорной пассивной однопозиционной моноимпульсной трёхкоординатной угломерно-разностно-доплеровской локации перемещающихся в пространстве радиоизлучающих объектов | |
CN106249216A (zh) | 基于高频地波雷达的静止目标双路径回波信息实现电离层高度估计的方法 | |
Cheng et al. | An EM-based multipath interference mitigation in GNSS receivers | |
CN106100769A (zh) | 一种多个不同体制卫星下微弱回波信号联合检测方法 | |
Lee et al. | Array-based GNSS signal tracking with a reduced state signal model | |
Yang et al. | RF emitter geolocation using amplitude comparison with auto-calibrated relative antenna gains | |
Qiao et al. | A survey of GNSS interference monitoring technologies | |
CN101788675B (zh) | 基于单通道多延迟互相关处理的gps干扰抑制方法 | |
CN107918126B (zh) | 基于多特征自动分割的多通道近海岸模糊杂波抑制方法 | |
CN115616628B (zh) | 基于角跟踪环路的gnss天线阵接收机盲波束形成方法 | |
Huang et al. | Navigation domain multipath characterization using GNSS direct position estimation in urban canyon environment | |
Xue et al. | Overview of multipath mitigation technology in global navigation satellite system | |
Giremus et al. | A fixed-lag particle filter for the joint detection/compensation of interference effects in GPS navigation | |
Zhang et al. | GNSS multipath mitigation algorithm with antenna arrays based on matrix reconstruction |
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 |