CN113126076A - 基于双模式雷达的协作式距离和角度参数估计方法及装置 - Google Patents

基于双模式雷达的协作式距离和角度参数估计方法及装置 Download PDF

Info

Publication number
CN113126076A
CN113126076A CN202110262416.0A CN202110262416A CN113126076A CN 113126076 A CN113126076 A CN 113126076A CN 202110262416 A CN202110262416 A CN 202110262416A CN 113126076 A CN113126076 A CN 113126076A
Authority
CN
China
Prior art keywords
radar
distance
estimation
angle
estimation value
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
CN202110262416.0A
Other languages
English (en)
Other versions
CN113126076B (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN202110262416.0A priority Critical patent/CN113126076B/zh
Publication of CN113126076A publication Critical patent/CN113126076A/zh
Application granted granted Critical
Publication of CN113126076B publication Critical patent/CN113126076B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/42Simultaneous measurement of distance and other co-ordinates
    • 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/41Details 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

Landscapes

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

Abstract

本发明公开了一种基于双模式雷达的协作式距离和角度参数估计方法及装置,包括:构建发射‑接收信号模型,所述发射‑接收信号模型包括PA雷达、FDA雷达;根据所述发射‑接收信号模型获取接收信号,所述接收信号包括PA雷达接收数据、FDA雷达接收数据;根据所述PA雷达接收数据进行角度估计得到第一角度估计值;根据所述FDA雷达接收数据、所述第一角度估计值进行联合距离估计得到无模糊距离估计值;根据所述FDA雷达接收数据、所述第一角度估计值进行联合角度估计得到最终角度估计值。本发明方法,将PA雷达与FDA雷达的优点相结合,不同于分时模式,本实施例占用时间资源更少,可以同时实现角度、距离的高精度估计。

Description

基于双模式雷达的协作式距离和角度参数估计方法及装置
技术领域
本发明属于雷达信号处理技术领域,具体涉及一种基于双模式雷达的协作式距离和角度参数估计方法及装置。
背景技术
脉冲雷达是雷达的一种。能够辐射较短的高频脉冲,然后天线转接到接收机接收信号,因此发射和接收信号在时间上是分开的。脉冲雷达用于测距,尤其适于同时测量多个目标的距离。当前常用的雷达大多数是脉冲雷达。
脉冲雷达接收信号中包含了感兴趣信号(Signal Of Interest,简称SOI)的距离、角度信息,通过波束形成输出可以估计目标距离和角度参数。但由于阵元稀疏分布或高脉冲重复频率,估计得到的参数信息可能会出现模糊。传统相控阵雷达通过波束形成实现高天线增益因此具有高角度分辨率,但由于其不具有距离维自由度,无法直接估计距离模糊度,进而导致距离估计精度降低。同时,频率分集阵列雷达具有距离维自由度,具有距离模糊分辨能力,但是其距离和角度信息在发射域内耦合。如何解决模糊并同时实现高精度的距离、角度参数估计是当前的一个难点。“W.Wang and H.Shao,“Range-anglelocalization of targets by a double-pulse frequency diverse array radar,”IEEEJournal of Selected Topics in Signal Processing,vol.8,no.1,pp.106-114,February 2014.”中提出了一种利用双脉冲FDA雷达进行参数估计的方法(以下简称:双脉冲法),其算法为:结合PA和FDA两种雷达,使用分时模式,在第一个脉冲使用零频率增量估计SOI的角度,在第二个脉冲使用非零频率增量估计距离参数;“J.Xu,G.Liao,S.Zhu,L.Huang and H.C.So,“Joint range and angle estimation using MIMO radar withfrequency diverse array,”IEEE Transactions on Signal Processing,vol.63,no.13,pp.3396-3410,July 2015.”中提出了一种基于FDA-MIMO体制下联合距离角度估计,其算法为:利用MIMO体制优势,在发射域进行波束形成,实现无模糊估计;“C.Cui,J.Xu,R.Gui,W.Wang and W.Wu,“Search-free DOD,DOA and range estimation for bistatic FDA-MIMO radar,”IEEE Access,vol.6,pp.15431-15445,March 2018.”中提出了一种基于旋转不变技术的双基地FDA-MIMO雷达无需搜索的参数估计方法,其算法为:在双基地模式下,利用非线性频率增量,使用旋转不变技术,得到参数估计结果。
但是,双脉冲法由于PA和FDA雷达的分时工作,需要较长的时间资源;FDA-MIMO体制下的参数估计由于发射域内距离自由度与角度耦合,仅能直接在接收域估计角度,导致角度估计精度低;若FDA雷达采用非线性的频率增量,会导致发射域内波束能量损失,使得估计性能下降;最后,传统PA雷达由于不具有距离维自由度,无法直接由接收端波束形成输出估计距离模糊数,使得距离估计精度低。
发明内容
为了解决现有技术中存在的上述问题,本发明提供了一种基于双模式雷达的协作式距离和角度参数估计方法及装置。
本发明的一个实施例提供了一种基于双模式雷达的协作式距离和角度参数估计方法,包括:
构建发射-接收信号模型,所述发射-接收信号模型包括PA雷达、FDA雷达;
根据所述发射-接收信号模型获取接收信号,所述接收信号包括PA雷达接收数据、FDA雷达接收数据;
根据所述PA雷达接收数据进行角度估计得到第一角度估计值;
根据所述FDA雷达接收数据、所述第一角度估计值进行联合距离估计得到无模糊距离估计值;
根据所述FDA雷达接收数据、所述第一角度估计值进行联合角度估计得到最终角度估计值。
在本发明的一个实施例中,构建的发射-接收信号模型包括发射端、接收端,其中:
所述发射端包括PA雷达、FDA雷达,所述PA雷达位于左侧且包含L个阵元,所述FDA雷达位于右侧且包含M个阵元,相邻阵元间距均为d,所有阵元均发射正交信号;
所述接收端接收PA雷达、FDA雷达发射的雷达信号。
在本发明的一个实施例中,根据所述发射-接收信号模型获取的接收信号表示为:
Figure BDA0002970567400000031
其中,y表示接收信号,yPA表示PA雷达接收数据,yFDA表示FDA雷达接收数据,ξ表示目标散射系数,
Figure BDA0002970567400000032
表示Kronecker积,
Figure BDA0002970567400000033
表示发射后波束形成输出,λ0表示发射信号波长,θ0表示PA发射端波束指向,θ表示空间内目标角度,d表示阵元间距,
Figure BDA0002970567400000041
表示相位差值,a(R,θ)表示FDA雷达发射导向矢量,
Figure BDA0002970567400000042
Δf表示FDA雷达频率增量,c表示光速,R表示目标距离,[·]T表示转置操作,b(θ)表示接收导向矢量,
Figure BDA0002970567400000043
N=L+M。
在本发明的一个实施例中,根据所述PA雷达接收数据进行角度估计得到第一角度估计值包括:
根据所述接收导向矢量计算第一权矢量;
根据所述PA雷达接收数据、所述第一权矢量构建第一角度估计函数;
对所述第一角度估计函数进行角度估计得到所述第一角度估计值。
在本发明的一个实施例中,根据所述FDA雷达接收数据、所述第一角度估计值进行联合距离估计得到无模糊距离估计值包括:
对所述FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据;
根据所述补偿后的FDA雷达接收数据、所述第一角度估计值进行距离模糊数估计得到距离模糊数估计值;
根据所述补偿后的FDA雷达接收数据、所述第一角度估计值、所述距离模糊数估计值进行主值距离差估计得到主值距离差估计值;
根据所述距离模糊数估计值、所述主值距离差估计值计算得到所述无模糊距离估计值。
在本发明的一个实施例中,对所述FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据包括:
构建补偿矢量;
根据所述补偿矢量对所述FDA雷达接收数据进行补偿得到所述补偿后的FDA雷达接收数据。
在本发明的一个实施例中,根据所述补偿后的FDA雷达接收数据、所述第一角度估计值进行距离模糊数估计得到距离模糊数估计值包括:
计算补偿距离,根据所述补偿距离得到补偿后的FDA雷达发射导向矢量;
根据所述补偿后的FDA雷达发射导向矢量、所述接收导向矢量、所述第一角度估计值计算第二权矢量;
根据所述补偿后的FDA雷达接收数据、所述第二权矢量构建距离模糊数估计函数;
对所述距离模糊数估计函数进行模糊距离估计得到所述距离模糊数估计值。
在本发明的一个实施例中,根据所述补偿后的FDA雷达接收数据、所述第一角度估计值、所述距离模糊数估计值进行主值距离差估计得到主值距离差估计值包括:
构建补偿后的FDA雷达接收数据的协方差矩阵;
根据所述补偿后的FDA雷达发射导向矢量、所述接收导向矢量、所述第一角度估计值、所述距离模糊数估计值、所述补偿后的FDA雷达接收数据的协方差矩阵计算第三权矢量;
根据所述补偿后的FDA雷达接收数据、所述第三权矢量构建主值距离差估计函数;
对所述主值距离差估计函数进行主值距离差估计得到所述主值距离差估计值。
在本发明的一个实施例中,根据所述FDA雷达接收数据、所述第一角度估计值进行联合角度估计得到最终角度估计值包括:
根据所述补偿后的FDA雷达发射导向矢量、所述接收导向矢量、所述主值距离差估计值、所述距离模糊数估计值计算第四权矢量;
根据所述补偿后的FDA雷达接收数据、所述第四权矢量构建第二角度估计函数;
对所述第二角度估计函数进行角度估计得到第二角度估计值;
根据所述第一角度估计值、所述第二角度估计值计算所述最终角度估计值。
本发明的另一个实施例提供了一种基于双模式雷达的协作式距离和角度参数估计装置,包括:
数据模型构建模块,用于构建发射-接收信号模型,所述发射-接收信号模型包括PA雷达、FDA雷达;
数据获取模块,用于根据所述发射-接收信号模型获取接收信号,所述接收信号包括PA雷达接收数据、FDA雷达接收数据;
第一数据估计模块,用于根据所述PA雷达接收数据进行角度估计得到第一角度估计值;
第二数据估计模块,用于根据所述FDA雷达接收数据、所述第一角度估计值进行联合距离估计得到无模糊距离估计值;
第三数据估计模块,用于根据所述FDA雷达接收数据、所述第一角度估计值进行联合角度估计得到最终角度估计值。
与现有技术相比,本发明的有益效果:
本发明提供的基于双模式雷达的协作式距离和角度参数估计方法,将PA雷达与FDA雷达的优点相结合,不同于分时模式,本实施例占用时间资源更少,可以同时实现角度、距离的高精度估计。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的一种基于双模式雷达的协作式距离和角度参数估计方法的流程示意图;
图2是本发明实施例提供的一种基于双模式雷达的协作式距离和角度参数估计方法中发射-接收信号模型示意图;
图3是本发明实施例提供的在不同信噪比下本发明协作式雷达估计方法距离估计的均方根误差示意图;
图4是本发明实施例提供的在不同信噪比下本发明协作式雷达估计方法角度估计的均方根误差对比示意图;
图5是本发明实施例提供的一种基于双模式雷达的协作式距离和角度参数估计装置结构示意图。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述,但本发明的实施方式不限于此。
实施例一
为解决由于距离模糊和角度估计而导致的目标定位问题,请参见图1,图1是本发明实施例提供的一种基于双模式雷达的协作式距离和角度参数估计方法的流程示意图。本实施例提出了一种基于双模式雷达的协作式距离和角度参数估计方法,该基于双模式雷达的协作式距离和角度参数估计方法包括以下步骤:
步骤1、构建发射-接收信号模型,发射-接收信号模型包括PA雷达、FDA雷达。
具体而言,请参见图2,图2是本发明实施例提供的一种基于双模式雷达的协作式距离和角度参数估计方法中发射-接收信号模型示意图,本实施例构建的发射-接收信号模型包括发射端、接收端,发射端为一维线阵,具体该一维线阵被分为PA雷达、FDA雷达两个部分:PA雷达位于左侧且包含L个阵元,FDA雷达位于右侧且包含M个阵元,相邻阵元间距均为d,所有阵元均发射正交信号,以使接收端可以良好分离接收PA雷达、FDA雷达发射的雷达信号。
步骤2、根据发射-接收信号模型获取接收信号,接收信号包括PA雷达接收数据、FDA雷达接收数据。
具体而言,本实施例根据步骤1构建的发射-接收信号模型,通过回波快拍数据,重建获取双模雷达的接收信号,获取的双模雷达的接收信号表示为:
Figure BDA0002970567400000081
其中,y表示接收信号,yPA表示PA雷达接收数据,yFDA表示FDA雷达接收数据,ξ表示目标散射系数,
Figure BDA0002970567400000082
表示Kronecker积,
Figure BDA0002970567400000083
表示发射后波束形成输出,λ0表示发射信号波长,θ0表示PA发射端波束指向,θ表示空间内目标角度,d表示阵元间距,
Figure BDA0002970567400000091
表示相位差值,a(R,θ)表示FDA雷达发射导向矢量,
Figure BDA0002970567400000092
Δf表示FDA雷达频率增量,c表示光速,R表示目标距离,[·]T表示转置操作,b(θ)表示接收导向矢量,
Figure BDA0002970567400000093
N=L+M。
步骤3、根据PA雷达接收数据进行角度估计得到第一角度估计值。
具体而言,本实施例步骤3根据PA雷达接收数据进行角度估计得到第一角度估计值包括步骤3.1、步骤3.2:
步骤3.1、根据接收导向矢量计算第一权矢量。
具体而言,由公式(1)得知接收导向矢量
Figure BDA0002970567400000094
本实施例根据接收导向矢量b(θ)计算第一权矢量,具体计算的第一权矢量w(θ)表示为:
w(θ)=b(θ) (2)
步骤3.2、根据PA雷达接收数据、第一权矢量构建第一角度估计函数。
具体而言,根据公式(1)可以获取PA雷达接收数据,根据PA雷达接收数据和公式(2)计算的第一权矢量w(θ)构建第一角度估计函数,构建的第一角度估计函数表示为:
Figure BDA0002970567400000095
其中,k表示第k个快拍,K表示快拍总数,(·)H表示共轭转置操作,
Figure BDA0002970567400000096
表示PA雷达的第k次快拍PA雷达接收数据。
步骤3.3、对第一角度估计函数进行角度估计得到第一角度估计值。
具体而言,本实施例步骤3.3采用最大似然(Maximum likelihood,简称ML)方法对第一角度估计函数进行角度估计得到第一角度估计值。
步骤4、根据FDA雷达接收数据、第一角度估计值进行联合距离估计得到无模糊距离估计值。
具体而言,本实施例步骤4根据FDA雷达接收数据、第一角度估计值进行联合距离估计得到无模糊距离估计值包括步骤4.1、步骤4.2、步骤4.3、步骤4.4:
步骤4.1、对FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据。
具体而言,由于FDA雷达存在距离依赖性,因此本实施例步骤4.1对FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据,具体包括:构建补偿矢量;根据补偿矢量对FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据。本实施例构造的补偿矢量gl表示为:
Figure BDA0002970567400000101
其中,rl表示主值距离,由距离单元个数以及距离单元大小决定。
然后,根据补偿矢量gl对每一个快拍得到的FDA雷达接收数据进行补偿,得到补偿后的FDA雷达接收数据,补偿后的FDA雷达接收数据表示为:
Figure BDA0002970567400000102
其中,
Figure BDA0002970567400000103
表示FDA雷达的补偿后的第k次快拍数据,diag(·)表示对角化操作,
Figure BDA0002970567400000104
表示FDA雷达的第k个快拍FDA雷达接收数据,
Figure BDA0002970567400000105
表示Kronecker积,1M表示一M维单位矩阵。
步骤4.2、根据补偿后的FDA雷达接收数据、第一角度估计值进行距离模糊数估计得到距离模糊数估计值。
具体而言,本实施例采用ML方法根据补偿后的FDA雷达接收数据、第一角度估计值进行距离模糊数估计,具体步骤4.2包括步骤4.2.1、步骤4.2.2、步骤4.2.3、步骤4.2.4:
步骤4.2.1、计算补偿距离,根据补偿距离得到补偿后的FDA雷达发射导向矢量。
具体而言,考虑到当前待检测距离门,本实施例首先构建了补偿距离方程:
R-rl=rΔ+(p-1)Ru (6)
其中,R表示目标距离,rl表示主值距离,rΔ表示主值距离差,p表示距离模糊数,Ru表示FDA雷达最大无模糊距离。则根据补偿距离得到的补偿后的FDA雷达发射导向矢量表示为:
Figure BDA0002970567400000111
步骤4.2.2、根据补偿后的FDA雷达发射导向矢量、接收导向矢量、第一角度估计值计算第二权矢量。
具体而言,由公式(7)得到补偿后的FDA雷达发射导向矢量,此时
Figure BDA0002970567400000112
中rΔ取值为0、θ取值为第一角度估计值
Figure BDA0002970567400000113
即此时的补偿后的FDA雷达发射导向矢量表示为
Figure BDA0002970567400000114
同样接收导向矢量表示为
Figure BDA0002970567400000115
则本实施例计算得到的第二权矢量表示为:
Figure BDA0002970567400000121
步骤4.2.3、根据补偿后的FDA雷达接收数据、第二权矢量构建距离模糊数估计函数。
具体而言,本实施例根据公式(5)得到了补偿后的FDA雷达接收数据,以及根据公式(8)计算的第二权矢量构建距离模糊数估计函数,构建的距离模糊数估计函数表示为:
Figure BDA0002970567400000122
其中,Na表示最大距离模糊数。
步骤4.2.4、对距离模糊数估计函数进行模糊距离估计得到距离模糊数估计值。
具体而言,本实施例步骤4.2.4采用ML方法对距离模糊数估计函数进行模糊距离估计得到距离模糊数估计值。
步骤4.3、根据补偿后的FDA雷达接收数据、第一角度估计值、距离模糊数估计值进行主值距离差估计得到主值距离差估计值。
具体而言,本实施例采用ML方法根据补偿后的FDA雷达接收数据、第一角度估计值、距离模糊数估计值进行主值距离差估计,具体步骤4.3包括步骤4.3.1、步骤4.3.2、步骤4.3.3:
步骤4.3.1、构建补偿后的FDA雷达接收数据的协方差矩阵。
具体而言,本实施例首先构建了补偿后的FDA雷达接收数据的协方差矩阵方程:
Figure BDA0002970567400000123
步骤4.3.2、根据补偿后的FDA雷达发射导向矢量、接收导向矢量、第一角度估计值、距离模糊数估计值、补偿后的FDA雷达接收数据的协方差矩阵计算第三权矢量。
具体而言,本实施例在自适应波束形成(Minimum Variance DistortionlessResponse,简称MVDR)准则下计算第三权矢量,由公式(10)得到了补偿后的FDA雷达接收数据的协方差矩阵,以及此时的补偿后的FDA雷达发射导向矢量表示为
Figure BDA0002970567400000131
接收导向矢量表示为
Figure BDA0002970567400000132
则本实施例计算得到的第三权矢量表示为:
Figure BDA0002970567400000133
其中,
Figure BDA0002970567400000134
(·)-1表示矩阵取逆操作。
步骤4.3.3、根据补偿后的FDA雷达接收数据、第三权矢量构建主值距离差估计函数。
具体而言,本实施例根据公式(5)得到了补偿后的FDA雷达接收数据,以及根据公式(11)计算的第三权矢量构建主值距离差估计函数,构建的主值距离差估计函数表示为:
Figure BDA0002970567400000135
步骤4.3.4、对主值距离差估计函数进行主值距离差估计得到主值距离差估计值。
具体而言,本实施例步骤4.3.4采用ML方法对主值距离差估计函数进行主值距离差估计得到主值距离差估计值。
步骤4.4、根据距离模糊数估计值、主值距离差估计值计算得到无模糊距离估计值。
具体而言,本实施例根据公式(9)得到了距离模糊数估计值,以及根据公式(12)得到了主值距离差估计值,将其均代入公式(6)计算得到无模糊距离估计值,即表示为:
Figure BDA0002970567400000141
步骤5、根据FDA雷达接收数据、第一角度估计值进行联合角度估计得到最终角度估计值。
具体而言,本实施例步骤5根据FDA雷达接收数据、第一角度估计值进行联合角度估计得到最终角度估计值包括步骤5.1、步骤5.2、步骤5.3、步骤5.4:
步骤5.1、根据补偿后的FDA雷达发射导向矢量、接收导向矢量、主值距离差估计值、距离模糊数估计值计算第四权矢量。
具体而言,在计算第四权矢量时,补偿后的FDA雷达发射导向矢量表示为
Figure BDA0002970567400000142
接收导向矢量表示为b(θ),则本实施例计算得到的第四权矢量表示为:
Figure BDA0002970567400000143
步骤5.2、根据补偿后的FDA雷达接收数据、第四权矢量构建第二角度估计函数。
具体而言,本实施例根据公式(5)得到了补偿后的FDA雷达接收数据,以及根据公式(14)计算的第四权矢量构建第二角度估计函数,构建的第二角度估计函数表示为:
Figure BDA0002970567400000144
步骤5.3、对第二角度估计函数进行角度估计得到第二角度估计值。
具体而言,本实施例步骤5.3采用ML方法对第二角度估计函数进行角度估计得到第二角度估计值。
步骤5.4、根据第一角度估计值、第二角度估计值计算最终角度估计值。
具体而言,本实施例根据公式(3)得到了第一角度估计值
Figure BDA0002970567400000151
以及根据公式(15)得到了第二角度估计值
Figure BDA0002970567400000152
根据第一角度估计值
Figure BDA0002970567400000153
第二角度估计值
Figure BDA0002970567400000154
计算最终角度估计值,最终估计估计值表示为:
Figure BDA0002970567400000155
为了验证本实施例提出的基于双模式雷达的协作式距离和角度参数估计方法的有效性,通过以下仿真实验进一步说明。
仿真条件:
假设发射端一维线阵中分为的两部分中:PA雷达发射阵元数为10,FDA发射阵元数为10,相邻阵元间距为0.15m,脉冲重复频率PRF为5KHz,快拍数为200,最大无模糊距离为30km,距离模糊数为5。
仿真1:
请参见图3,图3是本发明实施例提供的在不同信噪比下本发明协作式雷达估计方法距离估计的均方根误差示意图。由于PA雷达不具有距离自由度,无法实现无模糊距离估计,而由图3可以看出,在发射频率增量Δf=451KHz的双模式雷达协作下,距离估计具有较小的克拉美罗界,并且在较低信噪比时,均方根误差(Root Mean Squared Error,简称RMSE)也较小。在信噪比大于-17dB情况下,双模式雷达协作的距离估计RMSE小于1m,说明本发明可以实现高精度的估计距离。同时还可看出随着信噪比的增加,RMSE逐渐逼近克拉美罗界,精度提升。
仿真2:
请参见图4,图4是本发明实施例提供的在不同信噪比下本发明协作式雷达估计方法角度估计的均方根误差示意图,此时信噪比SNR由-25dB至5dB,并且进行了200次蒙特卡洛实验。由图4可以看出,随着信噪比的增加,角度估计的RMSE逐渐逼近克拉美罗界,说明信噪比可以提升估计精度。在信噪比为5dB情况下,FDA雷达模式下对应克拉美罗界最大,PA雷达次之,本发明提出协作式雷达估计方法的克拉美罗界最小,为0.01°,并且实际情况下,三种情况RMSE也符合这一规律,说明本发明提出方法角度可以有效减小估计误差,提升角度估计精度。
综上所述,本发明提出的基于双模式雷达的协作式距离和角度参数估计方法,可解决距离模糊情况下的距离估计问题,并且给出高精度的角度估计。
综上所述,本实施例提出的基于双模式雷达的协作式距离和角度参数估计方法,将PA雷达与FDA雷达的优点相结合,不同于分时模式,本实施例占用时间资源更少,可以同时实现角度、距离的估计,由于可以有效获取目标距离模糊数,实现了无模糊距离估计,提升了距离估计精度,并且解决了现在角度估计精度不足的问题,提升了一定的角度估计精度。
实施例二
在实施例一的基础上,请参见图5,图5是本发明实施例提供的一种基于双模式雷达的协作式距离和角度参数估计装置结构示意图,本实施例提出了一种基于双模式雷达的协作式距离和角度参数估计装置,包括:
数据模型构建模块,用于构建发射-接收信号模型,发射-接收信号模型包括PA雷达、FDA雷达。
具体而言,本实施例数据模型构建模块中构建的发射-接收信号模型包括发射端、接收端,其中:
发射端包括PA雷达、FDA雷达,PA雷达位于左侧且包含L个阵元,FDA雷达位于右侧且包含M个阵元,相邻阵元间距均为d,所有阵元均发射正交信号;
接收端接收PA雷达、FDA雷达发射的雷达信号。
数据获取模块,用于根据发射-接收信号模型获取接收信号,接收信号包括PA雷达接收数据、FDA雷达接收数据。
具体而言,本实施例数据获取模块中根据发射-接收信号模型获取的接收信号表示为:
Figure BDA0002970567400000171
其中,y表示接收信号,yPA表示PA雷达接收数据,yFDA表示FDA雷达接收数据,ξ表示目标散射系数,
Figure BDA0002970567400000172
表示Kronecker积,
Figure BDA0002970567400000173
表示发射后波束形成输出,λ0表示发射信号波长,θ0表示PA雷达发射端指向,θ表示空间内目标角度,d表示阵元间距,
Figure BDA0002970567400000174
表示相位差值,a(R,θ)表示FDA雷达发射导向矢量,
Figure BDA0002970567400000175
Δf表示FDA雷达频率增量,c表示光速,R表示目标距离,[·]T表示转置操作,b(θ)表示接收导向矢量,
Figure BDA0002970567400000176
N=L+M。
第一数据估计模块,用于根据PA雷达接收数据进行角度估计得到第一角度估计值。
具体而言,本实施例第一数据估计模块中根据PA雷达接收数据进行角度估计得到第一角度估计值包括:
根据接收导向矢量计算第一权矢量;
根据PA雷达接收数据、第一权矢量构建第一角度估计函数;
对第一角度估计函数进行角度估计得到第一角度估计值。
第二数据估计模块,用于根据FDA雷达接收数据、第一角度估计值进行联合距离估计得到无模糊距离估计值。
具体而言,本实施例第二数据估计模块中根据FDA雷达接收数据、第一角度估计值进行联合距离估计得到无模糊距离估计值包括:
对FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据;
根据补偿后的FDA雷达接收数据、第一角度估计值进行距离模糊数估计得到距离模糊数估计值;
根据补偿后的FDA雷达接收数据、第一角度估计值、距离模糊数估计值进行主值距离差估计得到主值距离差估计值;
根据距离模糊数估计值、主值距离差估计值计算得到无模糊距离估计值。
进一步地,本实施例对FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据包括:
构建补偿矢量;
根据补偿矢量对FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据。
进一步地,本实施例根据补偿后的FDA雷达接收数据、第一角度估计值进行距离模糊数估计得到距离模糊数估计值包括:
计算补偿距离,根据补偿距离得到补偿后的FDA雷达发射导向矢量;
根据补偿后的FDA雷达发射导向矢量、接收导向矢量、第一角度估计值计算第二权矢量;
根据补偿后的FDA雷达接收数据、第二权矢量构建距离模糊数估计函数;
对距离模糊数估计函数进行模糊距离估计得到距离模糊数估计值。
进一步地,本实施例根据补偿后的FDA雷达接收数据、第一角度估计值、距离模糊数估计值进行主值距离差估计得到主值距离差估计值包括:
构建补偿后的FDA雷达接收数据的协方差矩阵;
根据补偿后的FDA雷达发射导向矢量、接收导向矢量、第一角度估计值、距离模糊数估计值、补偿后的FDA雷达接收数据的协方差矩阵计算第三权矢量;
根据补偿后的FDA雷达接收数据、第三权矢量构建主值距离差估计函数;
对主值距离差估计函数进行主值距离差估计得到主值距离差估计值。
第三数据估计模块,用于根据FDA雷达接收数据、第一角度估计值进行联合角度估计得到最终角度估计值。
具体而言,本实施例第三数据估计模块中根据FDA雷达接收数据、第一角度估计值进行联合角度估计得到最终角度估计值包括:
根据补偿后的FDA雷达发射导向矢量、接收导向矢量、主值距离差估计值、距离模糊数估计值计算第四权矢量;
根据补偿后的FDA雷达接收数据、第四权矢量构建第二角度估计函数;
对第二角度估计函数进行角度估计得到第二角度估计值;
根据第一角度估计值、第二角度估计值计算最终角度估计值。
本实施例提供的一种基于双模式雷达的协作式距离和角度参数估计装置,可以执行上述基于双模式雷达的协作式距离和角度参数估计方法实施例,其实现原理和技术效果类似,在此不再赘述。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (10)

1.一种基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,包括:
构建发射-接收信号模型,所述发射-接收信号模型包括PA雷达、FDA雷达;
根据所述发射-接收信号模型获取接收信号,所述接收信号包括PA雷达接收数据、FDA雷达接收数据;
根据所述PA雷达接收数据进行角度估计得到第一角度估计值;
根据所述FDA雷达接收数据、所述第一角度估计值进行联合距离估计得到无模糊距离估计值;
根据所述FDA雷达接收数据、所述第一角度估计值进行联合角度估计得到最终角度估计值。
2.根据权利要求1所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,构建的发射-接收信号模型包括发射端、接收端,其中:
所述发射端包括PA雷达、FDA雷达,所述PA雷达位于左侧且包含L个阵元,所述FDA雷达位于右侧且包含M个阵元,相邻阵元间距均为d,所有阵元均发射正交信号;
所述接收端接收PA雷达、FDA雷达发射的雷达信号。
3.根据权利要求2所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,根据所述发射-接收信号模型获取的接收信号表示为:
Figure FDA0002970567390000011
其中,y表示接收信号,yPA表示PA雷达接收数据,yFDA表示FDA雷达接收数据,ξ表示目标散射系数,
Figure FDA0002970567390000021
表示Kronecker积,
Figure FDA0002970567390000022
表示发射后波束形成输出,λ0表示发射信号波长,θ0表示PA雷达发射端波束指向,θ表示空间内目标角度,d表示阵元间距,
Figure FDA0002970567390000023
表示相位差值,a(R,θ)表示FDA雷达发射导向矢量,
Figure FDA0002970567390000024
Δf表示FDA雷达频率增量,c表示光速,R表示目标距离,[·]T表示转置操作,b(θ)表示接收导向矢量,
Figure FDA0002970567390000025
N=L+M。
4.根据权利要求3所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,根据所述PA雷达接收数据进行角度估计得到第一角度估计值包括:
根据所述接收导向矢量计算第一权矢量;
根据所述PA雷达接收数据、所述第一权矢量构建第一角度估计函数;
对所述第一角度估计函数进行角度估计得到所述第一角度估计值。
5.根据权利要求3所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,根据所述FDA雷达接收数据、所述第一角度估计值进行联合距离估计得到无模糊距离估计值包括:
对所述FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据;
根据所述补偿后的FDA雷达接收数据、所述第一角度估计值进行距离模糊数估计得到距离模糊数估计值;
根据所述补偿后的FDA雷达接收数据、所述第一角度估计值、所述距离模糊数估计值进行主值距离差估计得到主值距离差估计值;
根据所述距离模糊数估计值、所述主值距离差估计值计算得到所述无模糊距离估计值。
6.根据权利要求5所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,对所述FDA雷达接收数据进行补偿得到补偿后的FDA雷达接收数据包括:
构建补偿矢量;
根据所述补偿矢量对所述FDA雷达接收数据进行补偿得到所述补偿后的FDA雷达接收数据。
7.根据权利要求5所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,根据所述补偿后的FDA雷达接收数据、所述第一角度估计值进行距离模糊数估计得到距离模糊数估计值包括:
计算补偿距离,根据所述补偿距离得到补偿后的FDA雷达发射导向矢量;
根据所述补偿后的FDA雷达发射导向矢量、所述接收导向矢量、所述第一角度估计值计算第二权矢量;
根据所述补偿后的FDA雷达接收数据、所述第二权矢量构建距离模糊数估计函数;
对所述距离模糊数估计函数进行模糊距离估计得到所述距离模糊数估计值。
8.根据权利要求7所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,根据所述补偿后的FDA雷达接收数据、所述第一角度估计值、所述距离模糊数估计值进行主值距离差估计得到主值距离差估计值包括:
构建补偿后的FDA雷达接收数据的协方差矩阵;
根据所述补偿后的FDA雷达发射导向矢量、所述接收导向矢量、所述第一角度估计值、所述距离模糊数估计值、所述补偿后的FDA雷达接收数据的协方差矩阵计算第三权矢量;
根据所述补偿后的FDA雷达接收数据、所述第三权矢量构建主值距离差估计函数;
对所述主值距离差估计函数进行主值距离差估计得到所述主值距离差估计值。
9.根据权利要求8所述的基于双模式雷达的协作式距离和角度参数估计方法,其特征在于,根据所述FDA雷达接收数据、所述第一角度估计值进行联合角度估计得到最终角度估计值包括:
根据所述补偿后的FDA雷达发射导向矢量、所述接收导向矢量、所述主值距离差估计值、所述距离模糊数估计值计算第四权矢量;
根据所述补偿后的FDA雷达接收数据、所述第四权矢量构建第二角度估计函数;
对所述第二角度估计函数进行角度估计得到第二角度估计值;
根据所述第一角度估计值、所述第二角度估计值计算所述最终角度估计值。
10.一种基于双模式雷达的协作式距离和角度参数估计装置,其特征在于,包括:
数据模型构建模块,用于构建发射-接收信号模型,所述发射-接收信号模型包括PA雷达、FDA雷达;
数据获取模块,用于根据所述发射-接收信号模型获取接收信号,所述接收信号包括PA雷达接收数据、FDA雷达接收数据;
第一数据估计模块,用于根据所述PA雷达接收数据进行角度估计得到第一角度估计值;
第二数据估计模块,用于根据所述FDA雷达接收数据、所述第一角度估计值进行联合距离估计得到无模糊距离估计值;
第三数据估计模块,用于根据所述FDA雷达接收数据、所述第一角度估计值进行联合角度估计得到最终角度估计值。
CN202110262416.0A 2021-03-10 2021-03-10 基于双模式雷达的协作式距离和角度参数估计方法及装置 Active CN113126076B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110262416.0A CN113126076B (zh) 2021-03-10 2021-03-10 基于双模式雷达的协作式距离和角度参数估计方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110262416.0A CN113126076B (zh) 2021-03-10 2021-03-10 基于双模式雷达的协作式距离和角度参数估计方法及装置

Publications (2)

Publication Number Publication Date
CN113126076A true CN113126076A (zh) 2021-07-16
CN113126076B CN113126076B (zh) 2023-12-26

Family

ID=76773014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110262416.0A Active CN113126076B (zh) 2021-03-10 2021-03-10 基于双模式雷达的协作式距离和角度参数估计方法及装置

Country Status (1)

Country Link
CN (1) CN113126076B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103823217A (zh) * 2014-02-26 2014-05-28 南京信息工程大学 基于双频发射的双基地mimo雷达高速运动目标参数估计方法
WO2015109870A1 (zh) * 2014-01-24 2015-07-30 深圳大学 Mimo雷达系统及其目标端相位同步方法
CN109143228A (zh) * 2018-08-28 2019-01-04 海南大学 互耦条件下双基地mimo雷达非圆目标的角度估计方法
CN110146871A (zh) * 2019-05-21 2019-08-20 西安电子科技大学 基于双频偏fda-mimo雷达的目标参数估计方法
WO2019158254A1 (de) * 2018-02-15 2019-08-22 Robert Bosch Gmbh WINKELSCHÄTZUNG UND MEHRDEUTIGKEITSAUFLÖSUNG VON RADARSENSOREN FÜR KRAFTFAHRZEUGE MIT GROßEM ANTENNENARRAY

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015109870A1 (zh) * 2014-01-24 2015-07-30 深圳大学 Mimo雷达系统及其目标端相位同步方法
CN103823217A (zh) * 2014-02-26 2014-05-28 南京信息工程大学 基于双频发射的双基地mimo雷达高速运动目标参数估计方法
WO2019158254A1 (de) * 2018-02-15 2019-08-22 Robert Bosch Gmbh WINKELSCHÄTZUNG UND MEHRDEUTIGKEITSAUFLÖSUNG VON RADARSENSOREN FÜR KRAFTFAHRZEUGE MIT GROßEM ANTENNENARRAY
CN109143228A (zh) * 2018-08-28 2019-01-04 海南大学 互耦条件下双基地mimo雷达非圆目标的角度估计方法
CN110146871A (zh) * 2019-05-21 2019-08-20 西安电子科技大学 基于双频偏fda-mimo雷达的目标参数估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
田语柔;王晶琦;谭晶;薛伟;吴文;: "MIMO雷达角度超分辨估计算法研究", 微波学报, no. 1 *
赵智昊;王志民;孙扬;: "双基地FDA-MIMO雷达角度、距离及速度无模糊估计方法", 四川大学学报(自然科学版), no. 06 *
陈宝欣;关键;董云龙;黄勇;陈小龙;: "多频连续波雷达与角度-距离联合估计方法", 电子学报, no. 02 *
陈金立;李家强;顾红;: "双基地MIMO雷达高速运动目标的DOD和DOA联合估计", 电子与信息学报, no. 04 *

Also Published As

Publication number Publication date
CN113126076B (zh) 2023-12-26

Similar Documents

Publication Publication Date Title
CN111308437B (zh) 一种毫米波mimo交通雷达求熵解速度模糊的方法
CN108459307B (zh) 基于杂波的mimo雷达收发阵列幅相误差校正方法
CN109975807B (zh) 一种适用于毫米波车载雷达的降维子空间测角方法
CN101251597B (zh) 一种多输入多输出雷达系统阵列误差自校正的方法
CN109782243B (zh) 基于块Hankel矩阵填充的阵元故障MIMO雷达角度估计方法
US9562968B2 (en) Sensor system and method for determining target location using sparsity-based processing
CN109581352B (zh) 一种基于毫米波雷达的超分辨测角系统
CN110488255B (zh) 一种相控阵雷达单脉冲高分辨测角系统及方法
CN108693511B (zh) 时分复用mimo雷达的运动目标角度计算方法
US11454702B2 (en) Synthetic aperture radar method and synthetic aperture radar device
CN110412559A (zh) 分布式无人机mimo雷达的非相参融合目标检测方法
CA3083033A1 (en) Synthetic aperture radar apparatus and methods
CN109188386B (zh) 基于改进二维esprit算法的mimo雷达高分辨率参数估计方法
US10247815B1 (en) Phased array radar system with specular multipath mitigation
CN113189592B (zh) 考虑幅相互耦误差的车载毫米波mimo雷达测角方法
CN108828504B (zh) 基于部分相关波形的mimo雷达目标方向快速估计方法
CN114895260A (zh) 一种俯仰向空时编码星载sar系统回波分离方法
CN110196417A (zh) 基于发射能量集中的双基地mimo雷达角度估计方法
CN113960583A (zh) 基于发射波束域的机载mimo雷达发射接收稳健联合优化方法
CN108828586A (zh) 一种基于波束域的双基地mimo雷达测角优化方法
CN109597034B (zh) 一种基于欧几里得距离的空时自适应处理方法
CN113126076B (zh) 基于双模式雷达的协作式距离和角度参数估计方法及装置
CN114779199A (zh) 端射阵机载雷达自适应互耦补偿与杂波抑制方法
Cassidy et al. Analysis of MIMO-DBS performance for wide FOV automotive imaging
CN111142097A (zh) 一种调幅广播外辐射源雷达中的目标方向估计方法

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