CN109239683A - 宽带被动mimo雷达的克拉美罗界分析方法 - Google Patents

宽带被动mimo雷达的克拉美罗界分析方法 Download PDF

Info

Publication number
CN109239683A
CN109239683A CN201810990597.7A CN201810990597A CN109239683A CN 109239683 A CN109239683 A CN 109239683A CN 201810990597 A CN201810990597 A CN 201810990597A CN 109239683 A CN109239683 A CN 109239683A
Authority
CN
China
Prior art keywords
target
broadband
signal
cramér
lower bound
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
CN201810990597.7A
Other languages
English (en)
Other versions
CN109239683B (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.)
Tianjin University
Original Assignee
Tianjin 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 Tianjin University filed Critical Tianjin University
Priority to CN201810990597.7A priority Critical patent/CN109239683B/zh
Publication of CN109239683A publication Critical patent/CN109239683A/zh
Application granted granted Critical
Publication of CN109239683B publication Critical patent/CN109239683B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及雷达领域与通信领域,为提出一种宽带的CRB求解算法,既保证距离分辨率,也能保证检测性能的提升。本发明,宽带被动MIMO雷达的克拉美罗界分析方法,步骤如下:步骤一:确定信号模型;步骤二:得到目标回波信;步骤三:求待估计参量的似然比,得到接收信号的似然函数;步骤四:推导费歇尔信息矩阵FIM公式,通过对似然函数求得二阶导数;步骤五:根据CRB界的定义,求得宽带被动MIMO雷达的克拉美罗界,并利用MATLAB软件仿真出其结果。本发明主要应用于雷达设计制造场合。

Description

宽带被动MIMO雷达的克拉美罗界分析方法
技术领域
本发明涉及雷达领域与通信领域,具体讲,宽带被动MIMO雷达的克拉美罗界分析。
背景技术
近年来,无源被动雷达在抗“四大威胁”方面的优势正受到各国雷达研究者的广泛关注,其基本理论也取得了突破性进展。由于被动雷达系统只存在接收机,自身不发射电磁信号,而是利用外部辐射源信号进行目标探测,因此具有非常强的隐蔽性、“四抗”能力,无电磁污染、无探测盲区、便携和低成本等优势。
多输入多输出(Multiple Input Multiple Output,MIMO)雷达是利用系统发射端和接收端的多副天线获得空间分集增益[1],对被动雷达性能的评价可以从克拉美罗界和模糊函数两方面介绍,文献[2]在相关条件下推导了正交信号的速度位移联合估计的克拉美罗界(Cramer-Rao Bounds,CRB),文献[3]推导了对高斯最小频移键控(GaussianFiltered Minimum Shift Keying,GMSK)更加普遍的CRB,文献[4]分别计算了在相关和非相关条件下修正后的CRB。当前的参数估计CRB主要是针对窄带信号以及点目标,具有局限性,因此,可以扩展到距离扩展目标,提高距离分辨率,而当前对MIMO雷达往往要求发射波形相互正交,使得其在接收端能够被区分,除此之外,5G非正交技术也发展迅猛,可以将两者结合起来,用以提高频谱利用率。综上所述,当前的得动雷达研究呈现出从窄带信号到宽带信号,从正交到非正交的发展趋势。
[1]陈浩文,黎湘,庄钊文.一种新兴的雷达体制—MIMO雷达[J].电子学报,2012,40(6):1190-1198.
[2]P.L.Sun,J.Tang and S.Wan,“Cramer-Rao bound of joint estimation oftarget location and velocity for coherent MIMO radar,”Journal of SystemsEngineering and Electronics,vol.25,no.4,pp.566-572,2014.
[3]S.Gogineni,M.Rangaswamy,B.D.Rigling and A.Nehorai,“Cramr-RaoBounds for UMTS-Based Passive Multistatic Radar,”IEEE Transactions on SignalProcessing,vol.62,no.1,pp.95-106,2014.
[4]Q.He,J.Hu,R.S.Blum and Y.Wu,“Generalized Cramer-Rao Bound forJoint Estimation of Target Position and Velocity for Active and Passive RadarNetworks,”IEEE Transactions on Signal Processing,vol.64,no.8,pp.2078-2089,2016.
发明内容
为克服现有技术的不足,本发明旨在提出提出一种宽带的CRB求解算法,既保证距离分辨率,也能保证检测性能的提升。为此,本发明采用的技术方案是,宽带被动MIMO雷达的克拉美罗界分析方法,步骤如下:
步骤一:确定信号模型,由于距离分辨率和带宽之间是反比的关系,使得宽带信号的距离分辨率提高,此时检测目标不能被认为是点目标,在这个前提下,检测目标是距离扩展目标;
步骤二:得到目标回波信号,r(t)为接收到的目标回波信号,数学模型为
服从零均值,其中,其中E是发送端的传输能量,f(t)是发送信号的复包络,其被归一化为λ表示目标在分辨率下不同不同部位的延时,b(λ)是目标反射函数,w(t)为噪声,服从零均值,协方差为σω 2的复高斯分布,并且与目标反射系数是相互独立的;
步骤三:求待估计参量的似然比,得到接收信号的似然函数为
SR(λ)是距离散射函数,f*(t)是f(t)的共轭函数,λ,fD分别表示信号在经过目标反射后所受到的时延和多普勒。在可分辨多径条件下表示为
其中N表示可分辨多径的数目,bi是每条径的目标反射系数,服从零均值复高斯分布;
步骤四:推导费歇尔信息矩阵FIM公式,通过对似然函数求得二阶导数,得到FIM为
其中为未知参量,E()表示求期望值;
步骤五:根据CRB界的定义,求得宽带被动MIMO雷达的克拉美罗界,并利用MATLAB软件仿真出其结果。
本发明的一个实例中:
其中的步骤一:检测目标是距离扩展目标,待估计参量为目标的速度位移,通过检测目标的时延和多普勒来估计速度位移θ=[x,y,vx,vy]T
其中的步骤三:求待估计参量的似然比,当信号所有特征值ei均小于噪声时,特别当满足
n表示信号特征值的个数,称之为低能量相关条件LEC(Low-Energy-Coherence),当LEC条件满足时,距离扩展目标的接收信号的似然函数为
SR(λ)是目标散射函数,在可分辨多径条件下表示为:
其中N表示可分辨多径的数目,bi是每条径的目标反射系数,服从零均值复高斯分布;当LEC条件不满足时,为了分析检测性能,依然把它当做接收模型,这时是次最优的接收模型;
步骤四:推导FIM公式,通过对似然函数求得二阶导数期望值的相反数,得到FIM为:
其中为未知参量。以下,E()或E[]表示求期望值,Re[]表示取实部的值。
因此,
其中
得到
得到:
似然函数对多普勒的求导过程如下
因此,得到
所以FIM表示为
步骤五:求解得到CRB=J-1,根据链式法则:
θ=[x,y,vx,vy]T
因此,可得
根据克拉美罗界的定义,求得宽带被动MIMO雷达的克拉美罗界。
本发明的特点及有益效果是:
本发明基于宽带被动MIMO雷达的CRB界分析,不同于以往窄带点目标视距条件下的速度位移估计,在宽带信号下,提高了距离分辨率,可以扩展到距离扩展目标,分析距离扩展目标模型并推导了CRB,因此具有更好的距离检测性能。除此之外,5G非正交技术也发展迅猛,可以将两者结合起来,用以提高频谱利用率,能够为被动MIMO雷达检测估计提供新的指导。
附图说明:
图1距离扩展目标模型。
图2距离扩展目标的最佳接收机。
图3距离扩展目标下30M带宽OFDM信号CRB。
图4距离扩展目标下30M带宽线性调频信号CRB。
具体实施方式
当前的速度位移联合参数估计CRB主要是针对窄带信号以及点目标,具有局限性,因此,可以扩展到距离扩展目标,提高距离分辨率,而当前对MIMO雷达往往要求发射波形相互正交,使得其在接收端能够被区分,除此之外,5G非正交技术也发展迅猛,可以将两者结合起来,用以提高频谱利用率。本发明的目的在于提出一种宽带的CRB求解算法,既保证距离分辨率,也能保证检测性能的提升。
为克服现有技术的不足,本发明旨在分析宽带被动MIMO雷达的克拉美罗界。首先假设点目标和视距范围内的条件下,对于多载波宽带信号,我们同样可以使用窄带条件下最大似然估计的方法分析CRB。如果假设在距离扩展目标下,最佳检测的模型需要重新设计,而费歇尔信息矩阵(Fisher Information Matrix,FIM)和CRB将重新推导。进而探索出针对求解一般宽带信号CRB的分析方法。本发明采用的技术方案是,宽带被动MIMO雷达的克拉美罗界分析,步骤如下:
步骤一:确定信号模型,由于距离分辨率和带宽之间是反比的关系,使得宽带信号的距离分辨率提高,此时检测目标不能被认为是点目标。在这个前提下,假设检测目标是距离扩展目标。
步骤二:得到目标回波信号,r(t)为接收到的目标回波信号,数学模型为
其中E是发送端的传输能量,f(t)是发送信号的复包络,假设其被归一化为w(t)为噪声,服从零均值,协方差为σω 2的复高斯分布,并且与目标反射系数是相互独立的,λ,fD分别表示信号在经过目标反射后所受到的时延和多普勒。
步骤三:求待估计参量的似然比,可以得到接收信号的似然函数为
SR(λ)是距离散射函数,在可分辨多径条件下可以表示为
其中N表示可分辨多径的数目,bi是每条径的目标反射系数,服从零均值复高斯分布。
步骤四:推导FIM公式,通过对似然函数求得二阶导数的期望值的相反数,可以得到FIM为
其中为未知参量。
步骤五:根据CRB界的定义,求得宽带被动MIMO雷达的克拉美罗界,并利用MATLAB软件仿真出其结果。
本发明的一个实例如下:
本发明属雷达领域与通信领域,用宽带信号来分析被动MIMO雷达的克拉美罗界,以提高雷达系统的整体检测性能。详细步骤如下:
步骤一:确定信号模型,由于距离分辨率和带宽之间是反比的关系,使得宽带信号的距离分辨率提高,此时检测目标不能被认为是点目标。在这个前提下,假设检测目标是距离扩展目标,如图1。待估计参量为目标的速度位移,通过检测目标的时延和多普勒来估计速度位移θ=[x,y,vx,vy]T
步骤二:得到目标回波信号,r(t)为接收到的目标回波信号,数学模型为
b(λ)是目标反射函数,λ表示目标在分辨率下不同不同部位的延时,w(t)为噪声,服从零均值,协方差为σω 2的复高斯分布,并且与目标反射系数是相互独立的。
步骤三:求待估计参量的似然比。在许多复杂情况下,当信号所有特征值ei均小于噪声时,特别当满足
n表示信号特征值的个数,我们称之为低能量相关条件(Low-Energy-Coherence,LEC),尽管每个信号特征值很小,但结合大量的输出接收,最终在两种不同假设检验下,有不同的概率密度结果。当LEC条件满足时,距离扩展目标的最佳接收机可以表示如图2,距离扩展目标的接收信号的似然函数为
SR(λ)是目标散射函数,在可分辨多径条件下可以表示为
其中N表示可分辨多径的数目,bi是每条径的目标反射系数,服从零均值复高斯分布。当LEC条件不满足时,为了分析检测性能,我们依然把它当做接收模型,这时是次最优的接收模型。
步骤四:推导FIM公式,通过对似然函数求得二阶导数,可以得到FIM为
其中为未知参量,E()表示求期望,Re[]表示取实部的值。
因此,
其中
可以得到
类似的,可得
似然函数对多普勒的求导过程如下
因此,可以得到
所以FIM可以表示为
步骤五:求解得到CRB=J-1,根据链式法则:
θ=[x,y,vx,vy]T
因此,可得
根据克拉美罗界的定义,求得宽带被动MIMO雷达的克拉美罗界,并利用MATLAB软件分别仿真出OFDM信号与线性调频信号的结果,如图3和4。

Claims (2)

1.一种宽带被动MIMO雷达的克拉美罗界分析方法,其特征是,步骤如下:
步骤一:确定信号模型,由于距离分辨率和带宽之间是反比的关系,使得宽带信号的距离分辨率提高,此时检测目标不能被认为是点目标,在这个前提下,检测目标是距离扩展目标;
步骤二:得到目标回波信号,r(t)为接收到的目标回波信号,数学模型为
服从零均值,其中,其中E是发送端的传输能量,f(t)是发送信号的复包络,其被归一化为λ表示目标在分辨率下不同不同部位的延时,b(λ)是目标反射函数,w(t)为噪声,服从零均值,协方差为σω 2的复高斯分布,并且与目标反射系数是相互独立的;
步骤三:求待估计参量的似然比,得到接收信号的似然函数为
SR(λ)是距离散射函数,f*(t)是f(t)的共轭函数,λ,fD分别表示信号在经过目标反射后所受到的时延和多普勒。在可分辨多径条件下表示为
其中N表示可分辨多径的数目,bi是每条径的目标反射系数,服从零均值复高斯分布;
步骤四:推导费歇尔信息矩阵FIM公式,通过对似然函数求得二阶导数,得到FIM为
其中为未知参量,E()表示求期望值;
步骤五:根据CRB界的定义,求得宽带被动MIMO雷达的克拉美罗界,并利用MATLAB软件仿真出其结果。
2.如权利要求1所述的宽带被动MIMO雷达的克拉美罗界分析方法,其特征是,
其中的步骤一进一步地:检测目标是距离扩展目标,待估计参量为目标的速度位移,通过检测目标的时延和多普勒来估计速度位移θ=[x,y,vx,vy]T
其中的步骤三进一步地:求待估计参量的似然比,当信号所有特征值ei均小于噪声时,特别当满足
n表示信号特征值的个数,称之为低能量相关条件LEC(Low-Energy-Coherence),当LEC条件满足时,距离扩展目标的接收信号的似然函数为
SR(λ)是目标散射函数,在可分辨多径条件下表示为:
其中N表示可分辨多径的数目,bi是每条径的目标反射系数,服从零均值复高斯分布;当LEC条件不满足时,为了分析检测性能,依然把它当做接收模型,这时是次最优的接收模型;
步骤四进一步地:推导FIM公式,通过对似然函数求得二阶导数期望值的相反数,得到FIM为:
其中为未知参量。以下,E()或E[]表示求期望值,Re[]表示取实部的值。
因此,
其中
得到
得到:
似然函数对多普勒的求导过程如下
因此,得到
所以FIM表示为
步骤五进一步地:求解得到CRB=J-1,根据链式法则:
θ=[x,y,vx,vy]T
因此,可得
根据克拉美罗界的定义,求得宽带被动MIMO雷达的克拉美罗界。
CN201810990597.7A 2018-08-28 2018-08-28 宽带被动mimo雷达的克拉美罗界分析方法 Active CN109239683B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810990597.7A CN109239683B (zh) 2018-08-28 2018-08-28 宽带被动mimo雷达的克拉美罗界分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810990597.7A CN109239683B (zh) 2018-08-28 2018-08-28 宽带被动mimo雷达的克拉美罗界分析方法

Publications (2)

Publication Number Publication Date
CN109239683A true CN109239683A (zh) 2019-01-18
CN109239683B CN109239683B (zh) 2022-08-02

Family

ID=65067884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810990597.7A Active CN109239683B (zh) 2018-08-28 2018-08-28 宽带被动mimo雷达的克拉美罗界分析方法

Country Status (1)

Country Link
CN (1) CN109239683B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109782249A (zh) * 2019-02-22 2019-05-21 南京航空航天大学 一种两目标相关时延估计算法
CN113075621A (zh) * 2021-03-30 2021-07-06 电子科技大学 针对分布式网络化雷达的信号级定位算法精度界计算方法
CN113189574A (zh) * 2021-04-02 2021-07-30 电子科技大学 基于量化时延的云mimo雷达目标定位克拉美罗界计算方法
CN113359095A (zh) * 2021-04-27 2021-09-07 电子科技大学 一种相干被动mimo雷达克拉美罗界的计算方法
CN113701759A (zh) * 2021-08-27 2021-11-26 杭州腓腓科技有限公司 基于可重构超材料的室内同步定位与地图构建方法与系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140347211A1 (en) * 2013-05-24 2014-11-27 Robert Bosch Gmbh Method for operating a mimo radar
US20150295628A1 (en) * 2014-04-10 2015-10-15 Robert Bosch Gmbh Method for determining a time multiplexing sequence for a mimo radar
CN105068049A (zh) * 2015-07-27 2015-11-18 电子科技大学 一种分置天线mimo雷达的克拉美罗界计算方法
CN106909779A (zh) * 2017-02-17 2017-06-30 电子科技大学 基于分布式处理的mimo雷达克拉美罗界计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140347211A1 (en) * 2013-05-24 2014-11-27 Robert Bosch Gmbh Method for operating a mimo radar
US20150295628A1 (en) * 2014-04-10 2015-10-15 Robert Bosch Gmbh Method for determining a time multiplexing sequence for a mimo radar
CN105068049A (zh) * 2015-07-27 2015-11-18 电子科技大学 一种分置天线mimo雷达的克拉美罗界计算方法
CN106909779A (zh) * 2017-02-17 2017-06-30 电子科技大学 基于分布式处理的mimo雷达克拉美罗界计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
QING WANG 等: "A Low-Complexity Method for Two-Dimensional Direction-of-Arrival Estimation Using an L-Shaped Array", 《SENSORS》 *
何茜: "MIMO雷达检测与估计理论研究", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》 *
陈新亮 等: "高分辨雷达扩展目标检测算法研究", 《中国科学: 信息科学》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109782249A (zh) * 2019-02-22 2019-05-21 南京航空航天大学 一种两目标相关时延估计算法
CN113075621A (zh) * 2021-03-30 2021-07-06 电子科技大学 针对分布式网络化雷达的信号级定位算法精度界计算方法
CN113189574A (zh) * 2021-04-02 2021-07-30 电子科技大学 基于量化时延的云mimo雷达目标定位克拉美罗界计算方法
CN113359095A (zh) * 2021-04-27 2021-09-07 电子科技大学 一种相干被动mimo雷达克拉美罗界的计算方法
CN113359095B (zh) * 2021-04-27 2022-10-14 电子科技大学 一种相干被动mimo雷达克拉美罗界的计算方法
CN113701759A (zh) * 2021-08-27 2021-11-26 杭州腓腓科技有限公司 基于可重构超材料的室内同步定位与地图构建方法与系统
CN113701759B (zh) * 2021-08-27 2024-05-03 杭州腓腓科技有限公司 基于可重构超材料的室内同步定位与地图构建方法与系统

Also Published As

Publication number Publication date
CN109239683B (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
CN109239683A (zh) 宽带被动mimo雷达的克拉美罗界分析方法
CN105303137B (zh) 一种环境反向散射系统的读写器的门限值的确定方法
CN106707257B (zh) 基于嵌套阵列的mimo雷达波达方向估计方法
Liu et al. Persymmetric adaptive target detection with distributed MIMO radar
CN103353591B (zh) 基于mimo的双基地雷达局域化降维杂波抑制方法
CN105807267B (zh) 一种mimo雷达扩展目标的检测方法
CN106909779B (zh) 基于分布式处理的mimo雷达克拉美罗界计算方法
CN102156279A (zh) 基于mimo的双基地雷达地面动目标检测方法
CN103353592A (zh) 基于mimo的双基地雷达多通道联合降维杂波抑制方法
CN103852759A (zh) 扫描雷达超分辨成像方法
CN104991232A (zh) 目标信号相关下信号级融合的组网雷达抗欺骗式干扰方法
US20150318881A1 (en) Communications system
CN103728607A (zh) 机载mimo雷达空时码三维自适应杂波对消方法
Li et al. MIMO radar moving target detection against compound-Gaussian clutter
CN108828504A (zh) 基于部分相关波形的mimo雷达目标方向快速估计方法
CN106154242A (zh) 基于分数低阶类相关熵的目标参数联合估计新算法
Colone et al. Reference-free amplitude-based WiFi passive sensing
Li et al. Adaptive filter bank multi-carrier waveform design for joint communication-radar system
Jameson et al. Reconnaissance using adaptive multi-carrier radar
Jin et al. Time reversal adaptive waveform in MIMO radar
Sturm et al. Deterministic propagation modeling for joint radar and communication systems
Shome et al. Moving target detection and Doppler extraction using digital spread spectrum radar
Santra et al. Towards Adaptive MIMO Radar–Receiver Processing for Orthogonally Coded FMCW Waveforms
Dong et al. Integrated communication waveform design for bistatic radar
Giza Fusion of Sensing and Backscatter Communications via OFDM

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