CN109901149B - 一种基于fda-mimo雷达的目标参数估计方法 - Google Patents
一种基于fda-mimo雷达的目标参数估计方法 Download PDFInfo
- Publication number
- CN109901149B CN109901149B CN201910227674.8A CN201910227674A CN109901149B CN 109901149 B CN109901149 B CN 109901149B CN 201910227674 A CN201910227674 A CN 201910227674A CN 109901149 B CN109901149 B CN 109901149B
- Authority
- CN
- China
- Prior art keywords
- transmitting
- target
- array element
- mth
- signal
- 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.)
- Active
Links
Images
Abstract
本发明属于雷达信号处理技术领域,公开了一种基于FDA‑MIMO雷达的目标参数估计方法。该方法包括如下步骤:通过对导致波束耦合的频率增量和阵元间距进行设计,打破它们之间线性递增的同步性,达到解耦和的效果,为目标距离角度的唯一估计提供了条件。通过设计使得各阵元频率增量满足sinc函数递增的关系,再进行角度、距离的估计,得到精确的目标角度和距离信息。本发明提供的方法能够解决发射波束在距离‑角度上的耦合性,并且能够在多目标的距离角度都不同时,较好的估计目标的距离角度,避免出现假目标的情况,以提升雷达的估计性能。
Description
技术领域
本发明涉及雷达信号处理技术领域,尤其涉及一种基于 FDA-MIMO雷达的目标参数估计方法,可以用于FDA-MIMO体制雷达的参数估计。
背景技术
多输入多输出(英文全称:Multiple-Input Multiple-Output,英文缩写:MIMO)雷达能够根据不同的工作模式提供灵活的发射波束,然而与传统的相控阵类似,MIMO雷达的发射波束只与目标的角度有关,无法区分角度相同但距离不同的目标。此外,若存在角度相同而距离不同的干扰时,MIMO雷达的定位性能会受到很大影响。
频率分集阵列(英文全称:Frequency Diverse Array,英文缩写: FDA)的概念最早由Antonik和Wicks提出,这种阵列主要通过调节各阵元之间的载频差来实现新的系统功能。与传统相控阵不同,FDA的波束不仅具有方位角的依赖性,而且具有距离依赖性,因而能够填补传统相控阵波束没有距离指向分辨能力的不足,有效抑制距离依赖干扰。由于FDA的波束是与距离、角度和时间都相关的,及其发射波束是时变的,单纯的利用FDA并没有太大意义。
现有许多学者开始关注FDA-MIMO体制雷达,该体制雷达既能利用MIMO雷达的发射自由度,又能利用FDA雷达的距离维自由度,为解决新的问题提供了可能。但FDA-MIMO体制雷达的发射波束在距离 -角度上是耦和的,因此解决其耦和性是该体制雷达面临的主要问题。现有一些学者利用非线性频偏法来解决该问题,比如平方法、立方法、 log法等,在一定程度上这些方法达到了解决FDA-MIMO体制雷达发射波束的耦和,但是不同频偏法的目标发射波束会影响波束的使用性能,即雷达仅能较好的估计单目标参数或者同一角度下的多目标参数,当多目标的距离角度都不同时,就会出现假目标的情况,导致雷达的估计性能大大降低。
发明内容
本发明的实施例提供的一种基于FDA-MIMO雷达的目标参数估计方法,能够解决发射波束在距离-角度上的耦合性,并且能够在多目标的距离角度都不同时,较好的估计目标的距离角度,避免出现假目标的情况,以提升雷达的估计性能。
为达到上述目的,本发明的实施例采用如下技术方案:
一种基于FDA-MIMO雷达的目标参数估计方法,包括以下步骤:
其中,M=N,λ为雷达发射信号的波长;
步骤2,设置参考频率增量△f,计算第m个发射阵元的发射信号的的频率增量和第m个发射阵元的发射信号的载频,m=1,2,…,M;
步骤3,根据第m个发射阵元的载频,计算第m个发射阵元的发射信号;
步骤4,根据全部M个发射阵元的发射信号,计算每个接收阵元的接收到的目标回波信号,进一步得到全部N个接收阵元的接收到的目标回波信号矩阵;
步骤5,根据第m个发射阵元的发射信号,计算第m路匹配滤波函数;
步骤6,利用全部M路匹配滤波函数对全部N个接收阵元的接收到的目标回波信号矩阵进行匹配滤波,得到滤波后的回波数据矩阵Z;
步骤7,计算滤波后的回波数据矩阵Z的协方差矩阵,根据协方差矩阵,计算目标的角度和距离。
本发明提供的一种基于FDA-MIMO雷达的目标参数估计方法,通过对导致波束耦合的频率增量和阵元间距进行设计,打破它们之间线性递增的同步性,达到解耦和的效果,为目标距离角度的唯一估计提供了条件。通过设计使得各阵元频率增量满足sinc函数递增的关系,再利用Capon法进行角度、距离的估计,得到精确的目标角度和距离信息。本发明提供的方法能够解决发射波束在距离-角度上的耦合性,并且能够在多目标的距离角度都不同时,较好的估计目标的距离角度,避免出现假目标的情况,以提升雷达的估计性能。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种基于FDA-MIMO雷达的目标参数估计方法的流程示意图;
图2为本发明实施例提供的方法的发射端频率增量设置图;
图3为本发明实施例提供的方法的发射波束形成图;
图4为现有方法的非线性频偏法的发射波束形成图;
图5为单目标情况下,本发明实施例提供的方法计算得到的距离和角度;图5(a)为距离估计结果图;图5(b)为角度估计结果;
图6为同一角度、多目标情况下,本发明实施例提供的方法计算得到的距离和角度;图6(a)为距离估计结果,图6(b)为多目标的角度估计结果;
图7为不同角度、不同距离且多目标情况下,本发明实施例提供的方法计算得到的距离和角度;图7(a)为距离估计结果图,图7(b) 为多目标的角度估计结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例提供的一种基于FDA-MIMO雷达的目标参数估计方法的流程示意图,参照图1,本发明实施例提供的方法包括以下步骤:
其中,M=N,λ为雷达发射信号的波长;
步骤2,设置参考频率增量△f,计算第m个发射阵元的发射信号的的频率增量和第m个发射阵元的发射信号的载频,m=1,2,…,M;
进一步的,步骤2具体包括:
步骤2.1,计算第m个发射阵元的频率增量△fm=sinc((m-1)π/M)·△f;
步骤2.2,根据所述第m个发射阵元的频率增量△fm,计算第m个发射阵元的发射信号的载频fm=f0+△fm=f0+sinc((m-1)π/M)·△f;
其中,sinc为辛克函数,f0为参考频率,△f<<f0。
步骤3,根据第m个发射阵元的载频,计算第m个发射阵元的发射信号;
进一步的,在步骤3中,利用预设公式计算第m个发射阵元的发射信号sm(t),预设公式为:
步骤4,根据全部M个发射阵元的发射信号,计算每个接收阵元的接收到的目标回波信号,进一步得到全部N个接收阵元的接收到的目标回波信号矩阵;
进一步的,步骤4具体包括:
步骤4.3,计算全部N个接收阵元的接收到的目标回波信号矩阵 y=[y1,y2,…,yn,…,yN]T。
步骤5,根据第m个发射阵元的发射信号,计算第m路匹配滤波函数;
优选的,在步骤5中,计算第m路匹配滤波函数具体包括:
其中,hm(t)表示第m路匹配滤波函数,fm为第m个发射阵元的发射信号的载频,c为光速,r为目标距离发射端的参考距离,为第m 个发射阵元的发射信号对应的基带调制信号,(·)*表示复共轭操作符, 0≤t≤T,T为发射信号脉冲持续时间。
步骤6,利用全部M路匹配滤波函数对全部N个接收阵元的接收到的目标回波信号矩阵进行匹配滤波,得到滤波后的回波数据矩阵Z;
优选的,步骤6具体包括:
根据公式
计算得到滤波后的回波数据矢量Z;
其中,y=[y1,y2,…,yn,…,yN]T,h(t)=[h1(t),h2(t),…,hm(t),…,hM(t)]T,[·]T为转置操作符,
为接收导向矢量,
步骤7,计算滤波后的回波数据矩阵Z的协方差矩阵,根据协方差矩阵,计算目标的角度和距离。
进一步的,步骤7具体包括:
步骤7.1,计算滤波后的回波数据矢量Z的协方差矩阵R=E{ZZH};
其中,(·)H表示共轭转置,E{·}表示求期望值;
其中,a(θ,r)为发射导向矢量,且a(θ,r)∈CM×1, b(θ)为接收导向矢量,且b(θ)∈CN×1,(·)-1表示求矩阵的逆,表示求使得i最大时θ的值,表示求使得·最大时r的值,arg{·}表示满足·时的变量值。
本发明提供的一种基于FDA-MIMO雷达的目标参数估计方法,通过对导致波束耦合的频率增量和阵元间距进行设计,打破它们之间线性递增的同步性,达到解耦和的效果,为目标距离角度的唯一估计提供了条件。通过设计使得各阵元频率增量满足sinc函数递增的关系,再利用Capon法进行角度、距离的估计,得到精确的目标角度和距离信息。本发明提供的方法能够解决发射波束在距离-角度上的耦合性,并且能够在多目标的距离角度都不同时,较好的估计目标的距离角度,避免出现假目标的情况,以提升雷达的估计性能。
进一步的,以下通过仿真实验来对本发明上述有益效果进行验证。
1.仿真参数:
本发明实施例提供的方法的仿真参数设置如表1所示。
表1
参数 | 参数值 | 参数 | 参数值 |
阵元个数 | 16 | 参考频率增量△f(KHz) | 70 |
脉冲重复频率 | 50000 | 阵元间距d(m) | 0.015 |
参考载频(GHz) | 10 |
2.仿真内容:
图3为本发明实施例提供的方法中发射端个阵元对应的频率增量设置图,如图所示个阵元的频率增量变化趋势符合sinc函数形状。
图4是t=0us时本发明实施例提供的方法的发射天线方向图,图5 是t=0us时基于log频偏法的发射天线方向图。图4中的发射波束宽度类似一个圆点,波束能量集中在目标附近,图5中的发射波束宽度呈椭圆状,由图4和图5对比可以看出,相对而言图4的发射波束要窄的多,在进行位置比较近的多目标识别的时候更具优势。
图6为3个同一角度的情况下,本发明实施例提供的方法的目标信息估计结果,其中,目标1的距离为28km,角度为0°;目标2的距离为30km,角度为0°;目标3的距离为32km,角度为0°:图6(a) 为距离估计结果,图中有三个峰值,对应的距离分别为(28 30 32)km,与目标的真实距离分布一致;图6(b)为多目标的角度估计结果,图中只有一个峰值,对应的角度为0°,由于三个目标的真实角度相同,所以角度维无法区分出来。由此可以说明在解决传统FDA-MIMO雷达距离模糊的前提下,基于sinc频偏的FDA-MIMO体制雷达能够准确实现同一角度下的多目标的无模糊定位。
图7为3个不同角度、不同距离的情况下,本发明实施例提供的方法的目标信息估计结果,其中,目标1的距离为28km,角度为0°;目标2的距离为32km,角度为30°;目标3的距离为37km,角度为45°:图7(a)为距离估计结果图,图中对应的三个距离峰值分别为(283237.1)km;与目标的真实距离分布一致;图7(b)为多目标的角度估计结果,图中对应的三个角度峰值分别为(0°,30°,45°),与目标的真实角度分布一致。由此可以看出,本发明提供的方法具有较好的距离、角度分辨率,能够在目标分布情况未知的情况下实现多目标参数的无模糊估计。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于一计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种基于FDA-MIMO雷达的目标参数估计方法,其特征在于,包括以下步骤:
其中,M=N,λ为雷达发射信号的波长;
步骤2,设置参考频率增量△f,计算第m个发射阵元的发射信号的频率增量和第m个发射阵元的发射信号的载频,m=1,2,...,M;
所述步骤2具体包括:
步骤2.1,计算第m个发射阵元的频率增量△fm=sinc((m-1)π/M)·△f;
步骤2.2,根据所述第m个发射阵元的频率增量△fm,计算第m个发射阵元的发射信号的载频fm=f0+△fm=f0+sinc((m-1)π/M)·△f;
其中,sin c为辛克函数,f0为参考频率,△f<<f0;
步骤3,根据第m个发射阵元的载频,计算第m个发射阵元的发射信号;
步骤4,根据全部M个发射阵元的发射信号,计算每个接收阵元的接收到的目标回波信号,进一步得到全部N个接收阵元的接收到的目标回波信号矩阵;
步骤5,根据所述第m个发射阵元的发射信号,计算第m路匹配滤波函数;
步骤6,利用所述全部M路匹配滤波函数对所述全部N个接收阵元的接收到的目标回波信号矩阵进行匹配滤波,得到滤波后的回波数据矩阵Z;
步骤7,计算所述滤波后的回波数据矩阵Z的协方差矩阵,根据所述协方差矩阵,计算目标的角度和距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910227674.8A CN109901149B (zh) | 2019-03-25 | 2019-03-25 | 一种基于fda-mimo雷达的目标参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910227674.8A CN109901149B (zh) | 2019-03-25 | 2019-03-25 | 一种基于fda-mimo雷达的目标参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109901149A CN109901149A (zh) | 2019-06-18 |
CN109901149B true CN109901149B (zh) | 2022-10-28 |
Family
ID=66953452
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910227674.8A Active CN109901149B (zh) | 2019-03-25 | 2019-03-25 | 一种基于fda-mimo雷达的目标参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109901149B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110346764B (zh) * | 2019-08-21 | 2022-03-29 | 上海无线电设备研究所 | 一种mimo雷达目标距离角度解耦方法 |
CN111352078B (zh) * | 2019-12-20 | 2020-11-10 | 湖北工业大学 | 杂波下基于admm的低截获频控阵mimo雷达系统的设计方法 |
CN111355517B (zh) * | 2020-03-17 | 2021-04-06 | 电子科技大学 | 一种面向高速移动用户的频控阵基站协作传输方法 |
CN112534299B (zh) * | 2020-08-05 | 2022-03-29 | 华为技术有限公司 | 一种基于雷达信号的发射方法和装置 |
CN112363151B (zh) * | 2020-11-05 | 2022-09-23 | 西安电子科技大学 | 频率分集阵多输入多输出雷达的自适应目标检测方法 |
CN113777599B (zh) * | 2021-09-13 | 2023-06-16 | 电子科技大学 | 一种基于fda-mimo雷达的速度解模糊方法 |
CN114047500B (zh) * | 2021-09-28 | 2023-03-07 | 电子科技大学 | 一种大规模频控阵非线性频偏产生电路 |
CN114035180B (zh) * | 2021-11-05 | 2023-03-07 | 电子科技大学 | 一种大规模频控阵非线性频偏产生电路 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105699953A (zh) * | 2016-01-28 | 2016-06-22 | 西安电子科技大学 | 频率分集mimo雷达距离-角度解耦合波束形成方法 |
CN108196239A (zh) * | 2018-01-18 | 2018-06-22 | 西安电子科技大学 | 一种频率分集mimo雷达的无模糊参数估计方法 |
CN109375213A (zh) * | 2018-08-08 | 2019-02-22 | 西安电子科技大学 | 基于子阵划分的频率分集阵列信号处理方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106559371B (zh) * | 2015-09-24 | 2021-01-12 | 索尼公司 | 用于无线通信的电子设备以及无线通信方法 |
CN107356921B (zh) * | 2017-08-11 | 2020-09-11 | 桂林电子科技大学 | 一种频率分集阵列雷达基于一次频偏目标定位的方法 |
CN108627806A (zh) * | 2018-05-17 | 2018-10-09 | 西安电子科技大学 | 基于对数频偏频率分集阵列自适应抗空间分布干扰的方法 |
-
2019
- 2019-03-25 CN CN201910227674.8A patent/CN109901149B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105699953A (zh) * | 2016-01-28 | 2016-06-22 | 西安电子科技大学 | 频率分集mimo雷达距离-角度解耦合波束形成方法 |
CN108196239A (zh) * | 2018-01-18 | 2018-06-22 | 西安电子科技大学 | 一种频率分集mimo雷达的无模糊参数估计方法 |
CN109375213A (zh) * | 2018-08-08 | 2019-02-22 | 西安电子科技大学 | 基于子阵划分的频率分集阵列信号处理方法 |
Non-Patent Citations (5)
Title |
---|
Frequency Diverse Array Radar With Logarithmically Increasing Frequency Offset;Waseem Khan;《IEEE Antennas Wireless Propagation Letters》;20141110;全文 * |
MIMO双基地机载雷达距离模糊杂波抑制方法;王宇卓等;《雷达学报》;20180428(第02期);全文 * |
前视阵FDA-MIMO雷达距离模糊杂波抑制方法;徐义正等;《系统工程与电子技术》;20181116(第12期);全文 * |
频控阵雷达技术及其应用研究进展;王文钦等;《雷达学报》;20180428(第02期);全文 * |
频率分集阵雷达技术探讨;许京伟等;《雷达学报》;20180420(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109901149A (zh) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109901149B (zh) | 一种基于fda-mimo雷达的目标参数估计方法 | |
CN102156279B (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
CN110456334B (zh) | 基于优化稀疏阵的tdm-mimo雷达系统及其信号处理方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN102213761B (zh) | 双基地共址多输入多输出雷达多目标定位方法 | |
CN109407055B (zh) | 基于多径利用的波束形成方法 | |
CN111880171B (zh) | 一种消除雷达目标盲速的脉冲分段编码方法 | |
CN103076596B (zh) | 基于先验信息的mimo雷达发射方向图设计方法 | |
CN110146871B (zh) | 基于双频偏fda-mimo雷达的目标参数估计方法 | |
CN105044684B (zh) | 基于射频隐身的mimo跟踪雷达发射波束的形成方法 | |
CN108196239B (zh) | 一种频率分集mimo雷达的无模糊参数估计方法 | |
CN110412570B (zh) | 基于空间脉冲相位编码的hrws-sar成像方法 | |
CN112034444B (zh) | 基于循环编码阵列的多波束雷达通信一体化方法 | |
CN103983952A (zh) | 一种非圆信号双基地mimo雷达低复杂度收发角度联合估计方法 | |
CN110471033A (zh) | 基于相干fda的二维空时编码干扰抑制方法 | |
Wang et al. | Range-azimuth decouple beamforming for frequency diverse array with Costas-sequence modulated frequency offsets | |
CN103777197A (zh) | 单基地mimo雷达中降维传播算子的方位估计方法 | |
CN109946655B (zh) | Mimo雷达lfm信号置零波形设计方法 | |
CN104678362B (zh) | Mimo天波超视距雷达波形优化方法 | |
CN114114188A (zh) | 一种低旁瓣的fda雷达通信一体化波形设计方法 | |
CN110988816A (zh) | 基于空时相位编码的相干mimo雷达发射波形设计方法 | |
CN108896974B (zh) | 一种改进的mimo阵列高分辨空间谱估计方法 | |
CN113238212B (zh) | 基于空时编码的频率分集阵列雷达距离分辨率增强方法 | |
CN111352077B (zh) | 基于多比例分式规划的低截获频控阵mimo雷达系统的设计方法 | |
CN111352079B (zh) | 一种基于频控阵mimo雷达的低截获系统的设计方法 |
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 |