CN109239683A - 宽带被动mimo雷达的克拉美罗界分析方法 - Google Patents
宽带被动mimo雷达的克拉美罗界分析方法 Download PDFInfo
- 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
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/40—Means 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雷达的克拉美罗界分析。
背景技术
近年来,无源被动雷达在抗“四大威胁”方面的优势正受到各国雷达研究者的广泛关注,其基本理论也取得了突破性进展。由于被动雷达系统只存在接收机,自身不发射电磁信号,而是利用外部辐射源信号进行目标探测,因此具有非常强的隐蔽性、“四抗”能力,无电磁污染、无探测盲区、便携和低成本等优势。
多输入多输出(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雷达的克拉美罗界。
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)
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)
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雷达克拉美罗界计算方法 |
-
2018
- 2018-08-28 CN CN201810990597.7A patent/CN109239683B/zh active Active
Patent Citations (4)
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)
Title |
---|
QING WANG 等: "A Low-Complexity Method for Two-Dimensional Direction-of-Arrival Estimation Using an L-Shaped Array", 《SENSORS》 * |
何茜: "MIMO雷达检测与估计理论研究", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》 * |
陈新亮 等: "高分辨雷达扩展目标检测算法研究", 《中国科学: 信息科学》 * |
Cited By (7)
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 |