首页> 中国专利> 一种多个GPS卫星微弱回波信号检测的方法

一种多个GPS卫星微弱回波信号检测的方法

摘要

本发明公开了一种多个GPS卫星微弱回波信号检测的方法,所述方法包括以下步骤:利用GPS公开的调制参数和C/A码对参考通道接收到的GPS信号进行分离重构得到参考信号;对回波信号进行直达波及其多径的抑制;对回波通道信号和不同的参考信号分别进行循环互模糊函数处理得到多个时延‑多普勒谱;对多个时延‑多普勒谱经过TF‑LV坐标变换得到多个检测量并将这些检测量进行叠加;根据检测量的概率分布得到最佳的检测门限,并与检验量的峰值进行比较判决,从而实现对多个GPS卫星微弱回波信号检测。本发明可以对低信噪比环境下的多个GPS回波信号具有良好的检测性能。

著录项

  • 公开/公告号CN105866750A

    专利类型发明专利

  • 公开/公告日2016-08-17

    原文格式PDF

  • 申请/专利权人 西安电子科技大学;

    申请/专利号CN201610343896.2

  • 发明设计人 刘明骞;李雨亭;李兵兵;

    申请日2016-05-23

  • 分类号G01S7/36(20060101);G01S7/41(20060101);G01S7/292(20060101);

  • 代理机构北京科亿知识产权代理事务所(普通合伙);

  • 代理人汤东凤

  • 地址 710071 陕西省西安市太白南路2号西安电子科技大学

  • 入库时间 2023-06-19 00:17:55

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2022-06-17

    专利权的转移 IPC(主分类):G01S 7/36 专利号:ZL2016103438962 登记生效日:20220606 变更事项:专利权人 变更前权利人:西安电子科技大学 变更后权利人:桂林长海发展有限责任公司 变更事项:地址 变更前权利人:710071 陕西省西安市太白南路2号西安电子科技大学 变更后权利人:541000 广西壮族自治区桂林市桂林长海路3号

    专利申请权、专利权的转移

  • 2022-03-01

    著录事项变更 IPC(主分类):G01S 7/36 专利申请号:2016103438962 变更事项:发明人 变更前:刘明骞李雨亭李兵兵 变更后:李明刘明骞李雨亭李兵兵

    著录事项变更

  • 2018-10-12

    授权

    授权

  • 2016-09-14

    实质审查的生效 IPC(主分类):G01S7/36 申请日:20160523

    实质审查的生效

  • 2016-08-17

    公开

    公开

说明书

技术领域

本发明属于通信技术和卫星信号处理技术领域,尤其涉及一种多个GPS卫星微弱回波信号检测的方法。

背景技术

随着空间技术的不断发展,全球导航卫星系统GNSS(Global NavigationSatellite System)的信号形式已经越来越多样和成熟。利用分布广泛的GNSS导航卫星作为外辐射源来对目标进行探测具有很重大的意义,以最成熟的GPS卫星系统为例,国内外的一些文献都表明GPS作为外辐射源进行目标探测是可行的而且探测性能优良。然而,传统GPS外辐射源回波检测方法都是基于单个GPS卫星信号进行分析处理,实际上,GPS卫星信号都工作在一个频段,接收天线不可避免的接收到3-4个以上GPS信号,这导致传统单个GPS卫星外辐射源回波检测方法在该场景下性能较差。因此研究多GPS卫星的微弱回波信号检测方法具有重要意义。

Wang lei等人介绍了地面多个CDMA辐射源的定位系统,但是该系统直达波是直接用电缆从基站引来的,因此不适用于GPS卫星直达波提取(WangLei,Wang Jun,Xiao Long.Passive location and precision analysis based on multipleCDMA base stations[C],Radar Conference,2009IET International,2009:1-4)。Michael Edrich等人介绍了基于地面FM/DAB/DVB-T辐射源的回波检测系统,但是该系统通过信号频段的不同来分离信号,所以不适合频段重叠的GPS信号(Michael Edrich,Alexander Schroeder,Fabienne Meyer.Design and performanceevaluation of a matureFM/DAB/DVB-T multi-illuminator passive radar system[J],IET Radar Sonar Navig,2014,8(2):114–122)。黄磊等人依据GPS信号公开的PN码来重构出纯净的参考信号,从而消除了参考通道噪声对检测的影响,但是 该论文只重构了接收到的多个GPS信号中的一个来进行检测,本质上还是单个GPS回波检测(黄磊,李立萍.基于GPS的外辐射源雷达关键技术研究[D],电子科技大学,2013)。

发明内容

本发明针对现有技术的不足,本发明旨在提供一种有效的多个GPS卫星的微弱回波信号检测方法,以提高在低信噪比环境下多个GPS卫星的微弱回波信号的检测概率。

本发明是这样实现的,一种多个GPS卫星微弱回波信号检测的方法包括以下步骤:

S1利用GPS的调制参数和C/A码对参考通道接收到的GPS信号进行分离重构得到各个独立的GPS参考信号;

S2利用步骤S1得到的GPS参考信号对回波信号进行直达波及其多径的抑制,得到无直达波和多径影响的回波信号;

S3对回波通道信号和不同的参考信号分别进行循环互模糊函数处理得到多个时延-多普勒谱;

S4对S3得到的多个时延多普勒谱经过TF-LV坐标变换得到多个检测量并将这些检测量进行叠加;

S5根据叠加后的距离速度检测量的概率分布得到最佳的检测门限,通过将该检测门限与步骤S4中得到的检验量的峰值进行比较判决,从而实现对多个GPS卫星微弱回波信号检测。

需要说明的是,步骤S1中,对参考通道接收到的GPS信号进行分离重构得到各个独立的GPS参考信号按以下进行:

参考通道接收到n个GPS信号的直达波信号模型表示为:

其中,nr(t)为直达波通道噪声,xr(t)为参考通道总信号,为第i个卫星的直达波信号,Pi是第i个卫星的发射信号功率,Ci(t)是带宽为1.023MHz的第i个GPS卫星的粗码(C/A码),D(t)i是带宽为50Hz的第i个GPS卫星的导航数据,f0为1575.42MHz的载频,φi是第i个GPS卫星信号的初始相位,nr(t)为直达波通道噪声。

利用本地产生的1575.42MHz载频对接收到的xr(t)进行下变频处理后为D*(t),其表示为:

D*(t)=Σi=1nPi·Ci(t)·Di(t)+nr(t);

其中n'r(t)为直达波通道噪声的下变频,可表示为:

nr(t)=nr(t)·e-j2πf0t;

得到基带混合信号D*(t)之后,选取信号l的本地C/A码与D*(t)相乘,解调出卫星l的信息Dl(t)为:

Dl(t)=[Σi=1nPi·Ci(t)·Di(t)+nr(t)]·Cl(t)=Pl·Cl(t)·Cl(t)·Dl(t)+[Σi=1ilnPi·Ci(t)·Di(t)]·Cl(t)+nr(t)·Cl(t)=Pl·Dl(t)+nr(t);

其中为常数,n″r(t)为噪声项,可表示为:

n″r(t)=n'r(t)·Cl(t);

由于n'r(t)与Cl(t)不相关,所以n″r(t)很小,不会影响电平判断。在获得信息Dl(t)之后,将该信息调制到与直达波相位同步的本地C/A码上,再将其乘上本地载波,就可得到了纯净的参考信号l,其表示为:

需要说明的是,步骤S2中,利用参考信号对回波通道信号进行直达波及其多径抑制按以下进行:

回波通道中信号的模型为:

xs(t)=Σi=0Nωixr1(t-τi)+Σi=0Nωixr1(t-τi)+.......+Σi=0Nωixrn(t-τi)+Sother(t);

其中,ωi是第i条径的增益,τi是第i条径带来的时延,Sother为回波和噪声信号,表示为:

Sother(t)=xr1(t-τ1)exp(j2πfd1)+xr2(t-τ2)exp(j2πfd2)+...+xrn(t-τn)exp(j2πfdn)+n(t);

其中,为多普勒偏移,Sother的求解步骤如下:

首先构造n个信号的多径空间为:

其中,Q为采样点数,K为最大时延,可以由最大探测距离Rmax/c得到,其中xreci(t)为参考信号,Rmax为最大探测距离,c为光速。然后利用最小二乘法原理抑制直达波及其多径,其实施步骤为:将求min||Ssur-Xref·α||2转化为求从而得出:

代入αestim,解得:

Sother=ssur(t)-Xrefαestim=Ssur-Xref(XrefHXref)-1XrefSsur;

其中,Ssur为回波通道信号,α为自适应权值,αestim为α的估计值,为Xref的转置,Sother为回波通道中最终所剩的回波和噪声。

需要说明的是,步骤S3中,所述对参考信号和回波进行循环互模糊函数处理按以下进行:

首先对参考信号xreci(t)做循环自相关运算:

Rxrecixreciρ(u)=lim1TT/2T/2xreci(t+u/2)xreci(t-u/2)*e-j2πρtdt=PiRrrρ(u);

其中,ρ为循环频率,为两个直达波信号的循环自相关,为两个直达波信号的幅值归一化的循环自相关,可表述为:

然后,对参考信号和回波通道信号做循环互相关运算,得到:

Rxrecixsρ-f(u)=lim1TT/2T/2xreci(t+u/2)xs(t-u/2)*e-j2π(ρ-f)tdt=riPie-jπfdiτie-jπ(ρ-f+fdi)τiRrrρ-f+fdi(u-τi)+Rrnρ(u);

其中,n(t)为回波通道的噪声,为直达波与噪声的循环互相关,其表达式为:

Rrnρ(u)=lim1TT/2T/2xreci(t+u/2)n(t-u/2)*e-j2π(ρ-f)tdt;

最后将和进行互模糊函数处理,得到时延-多普勒谱,其表示为:

S(τ,f)=Rxrecixsρ-f(u)Rxrecixreciρ(u-τ)*ejπfτdu=riPi23e-jπατiRrrρ-f+fdi(u-τi)Rrrρ(u-τ)*ejπ(f-fdi)(u+τi)du+PiRrnρ(u)Rrrρ(u-τ)ejπfudu=SRR(τ,f)+SNR(τ,f);

其中,SRR(τ′,f)为参考信号的循环自相关和参考信号与回波信号的循环互相关的互模糊函数,其可表示为:

SRR(τ,f)=riPi23e-jπατiRrrρ-f+fdi(u-τi)Rrrρ(u-τ)*ejπ(f-fdi)(u+τi)du;

SNR(τ′,f)为参考信号与噪声的循环互相关和直达波信号的循环自相关的互模糊函数,其可表示为:

SNR(τ,f)=PiRrnρ(u)Rrrρ(u-τ)ejπfudu.

需要说明的是,步骤S4中,所述对不同卫星得到的检测量的TF-LV变换按以下进行:

在发明中基线距离L和目标视角θ是固定的,根据几何关系可以列出方程组为:

Rr+Rt=L+Rt2=Rr2+L2-2RrL>cosθ;

其中,Rr为回波通道中目标到接收端的距离,,Rt为卫星到目标的距离,τ′为时延。通过解方程组得到距离和时延的关系为:

Rr=c2τ2+2Lcτ2(L+-L>cosθ)=f(τ)

多普勒频移和速度的关系为:

f2vλcosβ2=g(v)

其中,v为目标速度,λ为信号的波长,未知量β可以由以下公式计算得出:

sin(β)=2Rrsinθ(Rr+-Rrcosθ)Rr2-2Rr(Rr+)cos>Rr+(L+)2;

通过直达波与回波的循环互模糊函数得到的时延—多普勒谱为Si(τ′,f),i=1,2,3,…,该谱经过变换得到距离—速度谱Si(R,v),然后将多个GPS卫星的Si(R,v)谱进行叠加得到检测量为

需要说明的是,步骤5中所述根据检测量的概率分布得到最优判决门限,并与检测量的峰值进行比较判决,按以下进行:

检测量S(R,v)在不同的假设下分布分别为:

H0假设:

H1假设:

其中,H0假设为回波通道无回波信号,H1假设为回波通道存在回波信号, N(·,·)为高斯分布,Q为采样点数,为回波通道的噪声功率,σreci为参考信号功率。

根据上述假设设计检测器为:

max(S(R,v))H1><H0T;

其中,检测器的最佳检测门限T为:

T=-Σi=1nσn2σreci3Qln(PFA);

式中,PFA为虚警概率。通过检测量S(R,v)的峰值与该门限T进行比较从而得到检测结果。

本发明提供的多个GPS卫星微弱回波信号检测的方法,具有检测多个GPS回波信号的效果,在信噪比-25dB以上,联合检测概率在90%以上,可见本发明在低信噪比环境下具有良好的检测性能。

附图说明

图1是本发明实施例提供的多个GPS卫星微弱回波信号检测的方法流程图。

图2是本发明实施例提供的在不同信噪比下的检测性能示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。

下面结合附图对本发明的应用原理作详细的描述。

如图1所示,本发明为一种多个GPS卫星微弱回波信号检测的方法,其特征在于:所述方法包括以下步骤:

S1利用GPS的调制参数和C/A码对参考通道接收到的GPS信号进行分离 重构得到各个独立的GPS参考信号;

需要说明的是,步骤S1中,对参考通道接收到的GPS信号进行分离重构得到纯净的参考信号按以下进行:

参考通道接收到n个GPS信号的直达波信号模型表示为:

其中,nr(t)为直达波通道噪声,xr(t)为参考通道总信号,为第i个卫星的直达波信号,Pi是第i个卫星的发射信号功率,Ci(t)是带宽为1.023MHz的第i个GPS卫星的粗码(C/A码),D(t)i是带宽为50Hz的第i个GPS卫星的导航数据,f0为1575.42MHz的载频,φi是第i个GPS卫星信号的初始相位,nr(t)为直达波通道噪声。

利用本地产生的1575.42MHz载频对接收到的xr(t)进行下变频处理后为D*(t),其表示为:

D*(t)=Σi=1nPi·Ci(t)Di(t)+nr(t);

其中n'r(t)为直达波通道噪声的下变频,可表示为:

nr(t)=nr(t)·e-j2πf0t;

得到基带混合信号D*(t)之后,选取信号l的本地C/A码与D*(t)相乘,解调出卫星l的信息Dl(t)为:

Dl(t)=[Σi=1nPi·Ci(t)·Di(t)+nr(t)]·Cl(t)=Pl·Cl(t)·Cl(t)·Dl(t)+[Σi=1ilnPi·Ci(t)·Di(t)]·Cl(t)+nr(t)·Cl(t)=Pl·Dl(t)+nr(t);

其中为常数,n″r(t)为噪声项,可表示为:

n″r(t)=n'r(t)·Cl(t);

由于n'r(t)与Cl(t)不相关,所以n″r(t)很小,不会影响电平判断。在获得信息Dl(t)之后,将该信息调制到与直达波相位同步的本地C/A码上,再将其乘上本地载波,就可得到了纯净的参考信号l,其表示为:

S2利用步骤S1得到的GPS参考信号对回波信号进行直达波及其多径的抑制,得到无直达波和多径影响的回波信号;

需要说明的是,步骤S2中,利用参考信号对回波通道信号进行直达波及其多径抑制按以下进行:

回波通道中信号的模型为:

xs(t)=Σi=0Nωixr1(t-τi)+Σi=0Nωixr1(t-τi)+.......+Σi=0Nωixrn(t-τi)+Sother(t);

其中,ωi是第i条径的增益,τi是第i条径带来的时延,Sother为回波和噪声信号,表示为:

Sother(t)=xr1(t-τ1)exp(j2πfd1)+xr2(t-τ2)exp(j2πfd2)+...+xrn(t-τn)exp(j2πfdn)+n(t);

其中,为多普勒偏移,Sother的求解步骤如下:

首先构造n个信号的多径空间为:

其中,Q为采样点数,K为最大时延,可以由最大探测距离Rmax/c得到,其中xreci(t)为参考信号,Rmax为最大探测距离,c为光速。然后利用最小二乘法原理抑制直达波及其多径,其实施步骤为:将求min||Ssur-Xref·α||2转化为求从而得出:

代入αestim,解得:

Sother=ssur(t)-Xrefαestim=Ssur-Xref(XrefHXref)-1XrefSsur;

其中,Ssur为回波通道信号,α为自适应权值,αestim为α的估计值,为Xref的转置,Sother为回波通道中最终所剩的回波和噪声。

S3对回波通道信号和不同的参考信号分别进行循环互模糊函数处理得到多个时延-多普勒谱;

需要说明的是,步骤S3中,所述对参考信号和回波进行循环互模糊函数处理按以下进行:

首先对参考信号xreci(t)做循环自相关运算:

Rxrecixreciρ(u)=lim1TT/2T/2xreci(t+u/2)xreci(t-u/2)*e-j2πρtdt=PiRrrρ(u);

其中,ρ为循环频率,为两个直达波信号的循环自相关,为两个直达波信号的幅值归一化的循环自相关,可表述为:

然后,对参考信号和回波通道信号做循环互相关运算,得到:

Rxrecixsρ-f(u)=lim1TT/2T/2xreci(t+u/2)xs(t-u/2)*e-j2π(ρ-f)tdt=riPie-jπfdiτie-jπ(ρ-f+fdi)τiRrrρ-f+fdi(u-τi)+Rrnρ(u);

其中,n(t)为回波通道的噪声,为直达波与噪声的循环互相关,其表达式为:

Rrnρ(u)=lim1TT/2T/2xreci(t+u/2)n(t-u/2)*e-j2π(ρ-f)tdt;

最后将和进行互模糊函数处理,得到时延-多普勒谱,其表示为:

S(τ,f)=Rxrecixsρ-f(u)Rxrecixreciρ(u-τ)*ejπfτdu=riPi23e-jπατiRrrρ-f+fdi(u-τi)Rrrρ(u-τ)*ejπ(f-fdi)(u+τi)du+PiRrnρ(u)Rrrρ(u-τ)ejπfudu=SRR(τ,f)+SNR(τ,f);

其中,SRR(τ′,f)为参考信号的循环自相关和参考信号与回波信号的循环互相关的互模糊函数,其可表示为:

SRR(τ,f)=riPi23e-jπατiRrrρ-f+fdi(u-τi)Rrrρ(u-τ)*ejπ(f-fdi)(u+τi)du;

SNR(τ′,f)为参考信号与噪声的循环互相关和直达波信号的循环自相关的互模糊函数,其可表示为:

SNR(τ,f)=PiRrnρ(u)Rrrρ(u-τ)ejπfudu;

S4对S3得到的多个时延多普勒谱经过TF-LV坐标变换得到多个检测量并将这些检测量进行叠加;

需要说明的是,步骤S4中,所述对不同卫星得到的检测量的TF-LV变换按以下进行:

在发明中基线距离L和目标视角θ是固定的,根据几何关系可以列出方程组为:

Rr+Rt=L+Rt2=Rr2+L2-2RrL>cosθ;

其中,Rr为回波通道中目标到接收端的距离,,Rt为卫星到目标的距离,τ′为时延。通过解方程组得到距离和时延的关系为:

Rr=c2τ2+2Lcτ2(L+-L>cosθ)=f(τ);

多普勒频移和速度的关系为:

f2vλcosβ2=g(v);

其中,v为目标速度,λ为信号的波长,未知量β可以由以下公式计算得出:

sin(β)=2Rrsinθ(Rr+-Rrcosθ)Rr2-2Rr(Rr+)cos>Rr+(L+)2;

通过直达波与回波的循环互模糊函数得到的时延—多普勒谱为Si(τ′,f),i=1,2,3,…,该谱经过变换得到距离—速度谱Si(R,v),然后将多个GPS卫 星的Si(R,v)谱进行叠加得到检测量为

S5根据叠加后的距离速度检测量的概率分布得到最佳的检测门限,通过将该检测门限与步骤S4中得到的检验量的峰值进行比较判决,从而实现对多个GPS卫星微弱回波信号检测。

需要说明的是,步骤5中所述根据检测量的概率分布得到最优判决门限,并与检测量的峰值进行比较判决,按以下进行:

检测量S(R,v)在不同的假设下分布分别为:

H0假设:

H1假设:

其中,H0假设为回波通道无回波信号,H1假设为回波通道存在回波信号,N(·,·)为高斯分布,Q为采样点数,为回波通道的噪声功率,σreci为参考信号功率。

根据上述假设设计检测器为:

max(S(R,v))H1><H0T;

其中,检测器的最佳检测门限T为:

T=-Σi=1nσn2σreci3Qln(PFA);

式中,PFA为虚警概率。通过检测量S(R,v)的峰值与该门限T进行比较从而得到检测结果。

下面结合试验对本发明的应用效果作详细的描述。

为了测试本发明的检验性能,参数设置如下:采样频率设为10GHz,采样持续时间为20ms,GPS卫星直达波信号中粗码带宽为1.023MHz,导航数据码带宽为50Hz,载频为1575.42MHz,回波信号相对直达波时延分别为0.4ms、0.6ms 和0.8ms,回波信号相对于直达波多普勒偏移为500Hz、1000Hz和1500Hz,直达波相对于回波功率比为40dB。对上述参数数据进行2000次蒙特卡洛实验仿真,得到如图2所示的检测概率。

以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号