CN113777599A - 一种基于fda-mimo雷达的速度解模糊方法 - Google Patents

一种基于fda-mimo雷达的速度解模糊方法 Download PDF

Info

Publication number
CN113777599A
CN113777599A CN202111067790.1A CN202111067790A CN113777599A CN 113777599 A CN113777599 A CN 113777599A CN 202111067790 A CN202111067790 A CN 202111067790A CN 113777599 A CN113777599 A CN 113777599A
Authority
CN
China
Prior art keywords
speed
array
signal
fda
data
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
CN202111067790.1A
Other languages
English (en)
Other versions
CN113777599B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202111067790.1A priority Critical patent/CN113777599B/zh
Publication of CN113777599A publication Critical patent/CN113777599A/zh
Application granted granted Critical
Publication of CN113777599B publication Critical patent/CN113777599B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S13/00—Systems 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • 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/28—Details of pulse systems
    • G01S7/285—Receivers
    • G01S7/292—Extracting wanted echo-signals
    • G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods

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

本发明公开了一种基于FDA‑MIMO雷达的速度解模糊方法，通过构建FDA‑MIMO雷达的发射端阵列与接收端阵列；通过预设频率增量计算FDA‑MIMO雷达发射端阵列的阵元发射信号，并通过阵元发射信号计算接收端阵列的目标回波信号，并得到目标所在距离门的信号数据；根据信号数据分别计算存在多普勒模糊的各载频对应的估计速度与无多普勒模糊的估计速度，并利用一维集法得到目标速度估计值；本发明在不发射多个不同脉冲重复周期的信号的前提下，实现高速运动目标的无模糊速度估计，具有计算量小，获取速度精度高的优点，简化FDA‑MIMO雷达的系统复杂度，提高测量速度范围，有效解决了低重频下运动目标速度出现多普勒模糊的问题。

Description

一种基于FDA-MIMO雷达的速度解模糊方法
技术领域
本发明涉及雷达信号处理技术领域，具体是一种基于FDA-MIMO雷达的速度解模糊方法。
背景技术
频率分集阵列(英文全称：Frequency Diverse Array，英文缩写：FDA)的概念最早由Antonik和Wicks提出，这种阵列主要通过调节各阵元之间的载频差来实现新的系统功能。频率分集阵列发射的多载频信号给目标的检测与参数估计带来了新的优势。
多输入多输出(Multiple-Input Multiple-Output，MIMO)技术可以充分利用频率分集阵列发射多载频信号的优势，使其在接收端可以很好地分离出发射导向矢量。所以，常常将两种技术结合，形成FDA-MIMO体制雷达。
对于传统的脉冲多普勒雷达来说，由于奈崔斯特采样频率的限制，高速运动的目标往往会引起速度模糊。为了解决这一问题，可以利用多个不同脉冲重复周期的信号脉冲来估计运动目标的速度，提高测量速度的范围。但是这样会增加雷达系统的复杂性。对于FDA-MIMO体制雷达来说，这样无疑会大大增加该体制雷达的系统复杂性。
发明内容
针对现有技术中的上述：高速运动目标速度模糊，本发明提供了一种基于FDA-MIMO雷达的速度解模糊方法。
为了达到上述发明目的，本发明采用的技术方案为：
一种基于FDA-MIMO雷达的速度解模糊方法，包括以下步骤：
S1、构建FDA-MIMO雷达的发射端阵列与接收端阵列；
S2、根据预设频率增量计算步骤S1中发射端阵列的阵元发射信号；
S3、根据步骤S2中阵元发射信号得到接收端阵列的目标回波信号，并得到目标所在距离门的信号数据；
S4、根据步骤S3中信号数据计算存在多普勒模糊的各载频对应的估计速度；
S5、根据步骤S3中信号数据计算无多普勒模糊的估计速度；
S6、利用一维集法根据步骤S4中存在多普勒模糊的各载频对应的估计速度与步骤S5中无多普勒模糊的估计速度得到目标速度估计值。
本发明具有以下有益效果：
通过构建FDA-MIMO雷达的发射端阵列与接收端阵列，并通过预设频率增量计算发射端阵列的阵元发射信号，根据得到的发射端阵列的阵元发射信号计算接收端阵列的目标回波信号，并得到目标所在距离门的信号数据，利用信号数据分别计算存在多普勒模糊的各载频对应的估计速度与无多普勒模糊的估计速度，并结合一维集法得到目标速度估计值，在不发射多个不同脉冲重复周期的信号的前提下，实现高速运动目标的无模糊速度估计，具有计算量小，获取速度精度高等优点，有效解决了低重频下运动目标速度出现多普勒模糊的问题。
进一步地，步骤S1具体为：
预设FDA-MIMO雷达的发射端与接收端中阵元数，并根据发射信号的中心频率设置发射阵元与接收阵元的阵元间距，阵元间距的计算式表示为：
Figure BDA0003259050600000021
其中，dT为发射阵元间距，dR为接收阵元间距，
Figure BDA0003259050600000022
f0为发射信号的中心频率，c为光速。
进一步地，步骤S2具体包括以下分步骤：
S21、根据预设频率增量计算步骤S1中发射端阵列的阵元载频，计算式表示为：
fm＝f0+mΔf
其中，fm为发射端阵列中第m个阵元载频，Δf是频偏；
S22、根据步骤S21中阵元载频计算发射端阵元发射信号，计算式表示为：
Figure BDA0003259050600000031
其中，sm(t)为发射端第m个阵元的发射信号，φm(t)为发射端第m个阵元发射的基带信号，j为虚数，e、π均为常数，t为时间变量。
进一步地，S31、根据步骤S2中阵元发射信号计算接收端阵列接收的目标回波信号，表示为：
Figure BDA0003259050600000032
其中，yk(·)为接收端阵列接收的第k个脉冲的目标回波信号，θ为目标方位角，r为目标与接收端间距离，v为目标速度，(·)T为转置，ξ为目标的复反射特性参数，Tpr为发射信号中脉冲串的脉冲重复周期，N(t)为各接收端阵列中接收阵元的高斯白噪声，aT(·)为发射端阵列的导向矢量，aR(·)为接收端阵列的导向矢量，ΩD(·)为FDA-MIMO雷达的多普勒频移矩阵，e(·)为FDA-MIMO雷达的载波矢量，sm(·)为发射端第m个阵元的发射信号，MT为发射阵元总数；
S32、对步骤S31中目标回波信号进行多通道混频与匹配滤波，得到匹配滤波后的回波信号；
S33、采用FDA-MIMO雷达信号处理方法对步骤S32中匹配滤波后的回波信号进行处理，得到目标所在距离门的信号数据，表示为：
Figure BDA0003259050600000041
其中，Yk为目标所在距离门的信号数据，aT(·)为距离-角度联合导向矢量，MR为接收阵元总数。
进一步地，步骤S4具体包括以下分步骤：
步骤S4具体包括以下分步骤：
S41、提取步骤S3中信号数据的各列数据，构建数据矩阵；
S42、根据步骤S41中数据矩阵计算协方差矩阵，表示为：
Figure BDA0003259050600000042
其中，Rm为协方差矩阵，K为发射信号的脉冲总数，Zm为第m个阵元载频对应的数据矩阵，(·)H为共轭转置；
S43、根据步骤S41中数据矩阵与步骤S42中协方差矩阵计算各载频对应的多普勒频率，表示为：
Figure BDA0003259050600000043
其中，
Figure BDA0003259050600000044
为第m个阵元载频对应的多普勒频率，
Figure BDA0003259050600000045
为取该函数最大值所对应的多普勒频率fd，
Figure BDA0003259050600000046
ω(fd)为所有可能的多普勒频率矢量，(·)*为共轭；
S44、根据步骤S43中多普勒频率计算存在多普勒模糊的各载频对应的估计速度，表示为：
Figure BDA0003259050600000051
其中，vm为第m个阵元载频对应的估计速度。
进一步地，步骤S41具体为：
提取步骤S3中信号数据各列数据，得到各载频对应的数据矢量，并在得到全部脉冲数据后，构建各载频对应的数据矩阵。
进一步地，步骤S5具体包括以下分步骤：
步骤S5具体包括以下分步骤：
S51、根据步骤S3中信号数据计算无多普勒模糊的初始估计速度，表示为：
Figure BDA0003259050600000052
其中，xk为无多普勒模糊的初始估计速度，即为：各个载频回波数据间的相位差，
Figure BDA0003259050600000053
为步骤S32中目标所在距离门的信号数据Yk的矩阵元素yj,i,k共轭，yj,i+1,k为步骤S32中目标所在距离门的信号数据Yk的矩阵元素；
S52、对步骤S51中初始估计速度进行快速傅里叶变换，并选取变换后频谱最大值所对应的频率计算无多普勒模糊的估计速度，表示为：
Figure BDA0003259050600000054
其中，vx为无多普勒模糊的估计速度，FFT(xk)为对初始估计速度xk进行快速傅里叶变换。
进一步地，步骤S6具体包括以下分步骤：
S61、根据步骤S4中存在多普勒模糊的各载频对应的估计速度与步骤S5中无多普勒模糊的估计速度计算模糊参数，并确定搜索范围；
S62、根据步骤S61确定的搜索范围，将搜索范围内目标速度值进行排序，并分组计算均值与方差；
S63、利用一维集法根选取步骤S62中各组方差内最小值所对应的均值作为目标速度估计值。
进一步地，步骤S61中模糊参数计算式表示为：
Figure BDA0003259050600000061
其中，
Figure BDA0003259050600000062
为模糊参数，[·]为取整函数，
Figure BDA0003259050600000063
为第m个阵元载频所对应的无模糊估计的最大速度。
进一步地，步骤S62中方差的计算式表示为：
Figure BDA0003259050600000064
其中，Cv(j)为第j数据组的方差，
Figure BDA0003259050600000065
为第j数据组的第m*个速度值，Vj为搜索范围内MT个连续速度值构建的第j数据组，
Figure BDA0003259050600000066
为第j数据组的平均值。
该进一步方案具有以下有益效果：
1、通过对目标回波信号进行多通道混频与匹配滤波提高了信噪比，便于后续处理；
2、通过一维集法进一步提高速度的估计精度。
附图说明
图1为本发明提供的一种基于FDA-MIMO雷达的速度解模糊方法的步骤流程图；
图2为本发明实施例提供的FDA-MIMO雷达结构示意图；
图3为本发明中步骤S2的分步骤流程图；
图4为本发明中步骤S3中分步骤流程图；
图5为本发明实施例提供的FDA-MIMO雷达接收机的结构示意图；
图6为本发明中步骤S4中分步骤流程图；
图7为本发明中步骤S5中分步骤流程图；
图8为本发明中步骤S6中分步骤流程图；
图9为本发明实施例中仿真实验1提供的方法对不同速度目标的速度估计RMSE随信噪比变化曲线；
图10为本发明实施例中仿真实验1提供的方法对不同速度目标的速度估计Psuccess随信噪比变化曲线；
图11为本发明实施例中仿真实验2提供的方法对不同速度目标的速度估计RMSE随脉冲数变化曲线；
图12为本发明实施例中仿真实验2提供的方法对不同速度目标的速度估计Psuccess随脉冲数变化曲线。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，本发明实施例提供了一种基于FDA-MIMO雷达的速度解模糊方法，包括以下步骤S1至步骤S6：
S1、构建FDA-MIMO雷达的发射端阵列与接收端阵列；
本实施例中，步骤S1具体为：
预设FDA-MIMO雷达的发射端与接收端中阵元数，并根据发射信号的中心频率设置发射阵元与接收阵元的阵元间距，阵元间距的计算式表示为：
Figure BDA0003259050600000081
其中，dT为发射阵元间距，dR为接收阵元间距，
Figure BDA0003259050600000082
f0为发射信号的中心频率，c为光速。
实际中，通过构建FDA-MIMO雷达的发射端阵列与接收端阵列得到的FDA-MIMO雷达结构，如图2所示。
S2、根据预设频率增量计算步骤S1中发射端阵列的阵元发射信号；
如图3所示，本发明实施例中，步骤S2具体包括以下分步骤：
S21、根据预设频率增量计算步骤S1中发射端阵列的阵元载频，计算式表示为：
fm＝f0+mΔf
其中，fm为发射端阵列中第m个阵元载频，Δf是频偏；
实际中，亦可表示为：fm＝(m-1)f0。
S22、根据步骤S21中阵元载频计算发射端阵元发射信号，计算式表示为：
Figure BDA0003259050600000091
其中，sm(t)为发射端第m个阵元的发射信号，φm(t)为发射端第m个阵元发射的基带信号，j为虚数，e、π均为常数，t为时间变量。
实际中，各发射端阵元发射信号为脉冲串信号，脉冲的重复周期为Tpr。
S3、根据步骤S2中阵元发射信号得到接收端阵列的目标回波信号，并得到目标所在距离门的信号数据；
如图4所示，本发明实施例中，步骤S3具体包括以下分步骤：
步骤S3具体包括以下分步骤：
S31、根据步骤S2中阵元发射信号计算接收端阵列接收的目标回波信号，表示为：
Figure BDA0003259050600000092
其中，yk(·)为接收端阵列接收的第k个脉冲的目标回波信号，θ为目标方位角，r为目标与接收端间距离，v为目标速度，(·)T为转置，ξ为目标的复反射特性参数，Tpr为发射信号中脉冲串的脉冲重复周期，N(t)为各接收端阵列中接收阵元的高斯白噪声，aT(·)为发射端阵列的导向矢量，aR(·)为接收端阵列的导向矢量，ΩD(·)为FDA-MIMO雷达的多普勒频移矩阵，e(·)为FDA-MIMO雷达的载波矢量，sm(·)为发射端第m个阵元的发射信号，MT为发射阵元总数；
实际中，对于一个方位角为θ，距离为r，以速度v径向匀速靠近雷达的目标计算得到接收端阵列接收的目标回波信号
Figure BDA0003259050600000093
其中，
Figure BDA0003259050600000101
为接收端第MR-1个阵元接收的第k个脉冲的目标回波信号；
其中，发射端阵列的导向矢量
Figure BDA0003259050600000102
MT为发射阵元总数，计算式表示为：
Figure BDA0003259050600000103
接收端阵列的导向矢量
Figure BDA0003259050600000104
MR为接收阵元总数，计算式表示为:
Figure BDA0003259050600000105
FDA-MIMO雷达的多普勒频移矩阵
Figure BDA0003259050600000106
计算式表示为：
Figure BDA0003259050600000107
diag{·}为对角矩阵；
FDA-MIMO雷达的载波矢量
Figure BDA0003259050600000108
计算式表示为：
Figure BDA0003259050600000109
S32、对步骤S31中目标回波信号进行多通道混频与匹配滤波，得到匹配滤波后的回波信号；
实际中，目标回波信号进行多通道混频与匹配滤波，得到匹配滤波后的回波信号的过程，可表示为：
Figure BDA00032590506000001010
其中，y'n,m,k(t)为接收端第n个接收阵元接收到的第k个目标回波信号经过与第m个发射阵元的载频混频以及匹配滤波后的回波信号，yn,k(t)为接收端第n个接收阵元接收到的第k个目标回波信号，φm'(t)为发射端的基带信号φm(t)所匹配冲激响应函数，
Figure BDA00032590506000001011
表示卷积运算。
S33、采用FDA-MIMO雷达信号处理方法对步骤S32中匹配滤波后的回波信号进行处理，得到目标所在距离门的信号数据，表示为：
Figure BDA0003259050600000111
其中，Yk为目标所在距离门的信号数据，aT(·)为距离-角度联合导向矢量，MR为接收阵元总数。
实际中，将接收端接收的回波信号进行多通道混频与匹配滤波，其中接收机的结构，如图5所示，回波信号通过上述滤波处理后，对每个脉冲进行采样，便可得到目标所在距离门的信号数据Yk，其中，距离-角度联合导向矢量aT(r,θ)，定义表示为：
Figure BDA0003259050600000112
⊙为哈达玛积。
S4、根据步骤S3中信号数据计算存在多普勒模糊的各载频对应的估计速度；
如图6所示，本发明实施例中，步骤S4具体包括以下分步骤：
S41、提取步骤S3中信号数据的各列数据，构建数据矩阵；
本实施例中，步骤S41具体为：
提取步骤S3中信号数据各列数据，得到各载频对应的数据矢量，并在得到全部脉冲数据后，构建各载频对应的数据矩阵。
实际中，根据全部K个脉冲的目标所在距离门的信号数据Yk得到不同载频所估计的速度vm，但此时估计得到的速度是存在多普勒模糊，因此可计算多普勒模糊干扰下各载频对应的估计速度，需将目标所在距离门的信号数据Yk中每一列提取出来作为第m个载频对应的数据质量，表示为：
Figure BDA0003259050600000113
在经过总数为K个脉冲数据后，可得到第m个载频对应的数据矩阵Zm＝[zm,0 zm,1… zm,K-1]。
S42、根据步骤S41中数据矩阵计算协方差矩阵，表示为：
Figure BDA0003259050600000121
其中，Rm为协方差矩阵，K为发射信号的脉冲总数，Zm为第m个阵元载频对应的数据矩阵，(·)H为共轭转置；
S43、根据步骤S41中数据矩阵与步骤S42中协方差矩阵计算各载频对应的多普勒频率，表示为：
Figure BDA0003259050600000122
其中，
Figure BDA0003259050600000123
为第m个阵元载频对应的多普勒频率，
Figure BDA0003259050600000124
为取该函数最大值所对应的多普勒频率fd，
Figure BDA0003259050600000125
ω(fd)为所有可能的多普勒频率矢量，可表示为：
Figure BDA0003259050600000126
(·)*为共轭；
实际中，通过搜索多普勒频率可以得到该载频对应的多普勒频率。
S44、根据步骤S43中多普勒频率计算存在多普勒模糊的各载频对应的估计速度，表示为：
Figure BDA0003259050600000127
其中，vm为第m个阵元载频对应的估计速度。
S5、根据步骤S3中信号数据计算无多普勒模糊的估计速度；
实际中，根据全部K个脉冲的目标所在距离门的信号数据Yk得到无多普勒模糊的速度vx，此时的vx虽然没有多普勒模糊但是由于噪声的影响不够准确，与真实值之间的误差较大，因此需要进行下一步处理。
如图7所示，本实施例中，步骤S5具体包括以下分步骤：
S51、根据步骤S3中信号数据计算无多普勒模糊的初始估计速度，表示为：
Figure BDA0003259050600000131
其中，xk为无多普勒模糊的初始估计速度，即为：各个载频回波数据间的相位差，
Figure BDA0003259050600000132
为步骤S32中目标所在距离门的信号数据Yk的矩阵元素yj,i,k共轭，yj,i+1,k为步骤S32中目标所在距离门的信号数据Yk的矩阵元素；
S52、对步骤S51中初始估计速度进行快速傅里叶变换，并选取变换后频谱最大值所对应的频率计算无多普勒模糊的估计速度，表示为：
Figure BDA0003259050600000133
其中，vx为无多普勒模糊的估计速度，FFT(xk)为对初始估计速度xk进行快速傅里叶变换。
实际中，对初始估计速度xk进行FFT(快速傅里叶变换)，接着找到该频谱最大值对应的频率fdx，计算式表示为：
Figure BDA0003259050600000134
结合得到的该频谱最大值对应的频率fdx计算无多普勒模糊下各载频对应的估计速度，可表示为：
Figure BDA0003259050600000135
S6、利用一维集法根据步骤S4中存在多普勒模糊的各载频对应的估计速度与步骤S5中无多普勒模糊的估计速度得到目标速度估计值。
如图8所示，本实施例中，步骤S6具体包括以下分步骤：
S61、根据步骤S4中存在多普勒模糊的各载频对应的估计速度与步骤S5中无多普勒模糊的估计速度计算模糊参数，并确定搜索范围；
本实施例中，步骤S61中模糊估计参数计算式表示为：
Figure BDA0003259050600000141
其中，
Figure BDA0003259050600000149
为模糊估计参数，[·]为取整函数，
Figure BDA0003259050600000142
为第m个阵元载频所对应的无模糊估计的最大速度，计算式表示为：
Figure BDA0003259050600000143
S62、根据步骤S61确定的搜索范围，将搜索范围内目标速度值进行排序，并分组计算均值与方差；
本实施例中，步骤S62中方差的计算式表示为：
Figure BDA0003259050600000144
其中，Cv(j)为第j数据组的方差，
Figure BDA0003259050600000145
为第j数据组的第m*个速度值，Vj为搜索范围内MT个连续速度值构建的第j数据组，
Figure BDA0003259050600000146
为第j数据组的平均值。
实际中，由于噪声的存在会导致模糊数估计的错误，所以需要再模糊数附近进行搜素得到准确的无多普勒模糊速度，对于第m个阵元载频估计目标的可能速度可以表示为：
Figure BDA0003259050600000147
其中，L为预计的搜索范围，在
Figure BDA0003259050600000148
范围内进行搜索，l为单次搜索预计的模糊数，将搜索得到的全部目标可能的速度排成一组，可得到MT组数据，并且将MT组数据中数据进行从小到大进行排列，接着依次以MT个连续的速度值为一组，共得到2LMT+1组，分别按照上式进行均值与方差计算。
S63、利用一维集法根选取步骤S62中各组方差内最小值所对应的均值作为目标速度估计值。
实际中，依据一维集法，通过对比Cv(j)数值大小，并选取方差Cv(j)数最小的那组数据作为存在目标速度的最佳估计值的数据组，并将该组数据的均值作为目标速度估计值
Figure BDA0003259050600000154
本发明实施例中，通过以下仿真实验进行进一步说明，预设如表1所示的仿真参数，将发射、接收阵列的阵元间的间距预设为二分之一个波长，添加加性高斯白噪声，并假定发射端发射信号为线性调频信号，其中，信号模型表示为：
Figure BDA0003259050600000151
其中，
Figure BDA0003259050600000152
为脉宽为Tp的矩形脉冲，Bp为该基带信号的带宽。
表1仿真参数
Figure BDA0003259050600000153
采用速度的平均均方根误差(RMSE)来测量最优目标无多普勒模糊速度获取方法的精度。此外，为了测量该最优目标无多普勒模糊获取方法在解决多普勒模糊时的成功率，本发明实施例中假定当估计速度和真实速度之间的误差小于真实速度的1％时，该估计是正确的；接着计算多次试验中成功试验的百分比Psuccess。
仿真实验1中，限定目标运动速度分别为28.3m/s、275.7m/s和5675.3m/s，分别代表汽车、飞机和高速运动飞机的运动速度，脉冲数K为64，其余仿真参数如表1所示，实验结果如图9，图10所示，可知，由于目标的运动速度较大，在该仿真参数的条件下均存在速度模糊，若使用传统方法难以准确估计出目标的速度。但是通过本专利中的方法可以正确地估计目标的速度，并且从图中可以看出对于不同速度的目标该方法均可以得到较好的估计结果，此外，随着信噪比的增加，目标速度的估计精度会趋于一致。当信噪比超过-26dB时，平均均方根误差RMSE会小于0.01m/s，并且Psuccess会达到100％。
仿真实验2中，将信噪比固定为-26dB，脉冲数从8到128变化，其余仿真参数如表1所示，仿真实验结果如图11、图12所示，其中，图11展示了平均均方根误差RMSE与脉冲数的关系，图12展示了Psuccess与脉冲数量的关系，可知，当脉冲数量超过64时，速度的平均均方根误差RMSE低于0.01m/s，并且Psuccess为1。这表明在脉冲数量有限的情况下，本专利所提供的获取方法具有更好的估计精度。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (10)

1.一种基于FDA-MIMO雷达的速度解模糊方法，其特征在于，包括以下步骤：
S1、构建FDA-MIMO雷达的发射端阵列与接收端阵列；
S2、根据预设频率增量计算步骤S1中发射端阵列的阵元发射信号；
S3、根据步骤S2中阵元发射信号得到接收端阵列的目标回波信号，并得到目标所在距离门的信号数据；
S4、根据步骤S3中信号数据计算存在多普勒模糊的各载频对应的估计速度；
S5、根据步骤S3中信号数据计算无多普勒模糊的估计速度；
S6、利用一维集法根据步骤S4中存在多普勒模糊的各载频对应的估计速度与步骤S5中无多普勒模糊的估计速度得到目标速度估计值。
2.根据权利要求1所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S1具体为：
预设FDA-MIMO雷达的发射端与接收端中阵元数，并根据发射信号的中心频率设置发射阵元与接收阵元的阵元间距，阵元间距的计算式表示为：
Figure FDA0003259050590000011
其中，dT为发射阵元间距，dR为接收阵元间距，
Figure FDA0003259050590000012
f0为发射信号的中心频率，c为光速。
3.根据权利要求2所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S2具体包括以下分步骤：
S21、根据预设频率增量计算步骤S1中发射端阵列的阵元载频，计算式表示为：
fm＝f0+mΔf
其中，fm为发射端阵列中第m个阵元载频，Δf是频偏；
S22、根据步骤S21中阵元载频计算发射端阵元发射信号，计算式表示为：
Figure FDA0003259050590000021
其中，sm(t)为发射端第m个阵元的发射信号，φm(t)为发射端第m个阵元发射的基带信号，j为虚数，e、π均为常数，t为时间变量。
4.根据权利要求3所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S3具体包括以下分步骤：
S31、根据步骤S2中阵元发射信号计算接收端阵列接收的目标回波信号，表示为：
Figure FDA0003259050590000022
其中，yk(·)为接收端阵列接收的第k个脉冲的目标回波信号，θ为目标方位角，r为目标与接收端间距离，v为目标速度，(·)T为转置，ξ为目标的复反射特性参数，Tpr为发射信号中脉冲串的脉冲重复周期，N(t)为各接收端阵列中接收阵元的高斯白噪声，aT(·)为发射端阵列的导向矢量，aR(·)为接收端阵列的导向矢量，ΩD(·)为FDA-MIMO雷达的多普勒频移矩阵，e(·)为FDA-MIMO雷达的载波矢量，sm(·)为发射端第m个阵元的发射信号，MT为发射阵元总数；
S32、对步骤S31中目标回波信号进行多通道混频与匹配滤波，得到匹配滤波后的回波信号；
S33、采用FDA-MIMO雷达信号处理方法对步骤S32中匹配滤波后的回波信号进行处理，得到目标所在距离门的信号数据，表示为：
Figure FDA0003259050590000031
其中，Yk为目标所在距离门的信号数据，aT(·)为距离-角度联合导向矢量，MR为接收阵元总数。
5.根据权利要求4所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S4具体包括以下分步骤：
S41、提取步骤S3中信号数据的各列数据，构建数据矩阵；
S42、根据步骤S41中数据矩阵计算协方差矩阵，表示为：
Figure FDA0003259050590000032
其中，Rm为协方差矩阵，K为发射信号的脉冲总数，Zm为第m个阵元载频对应的数据矩阵，(·)H为共轭转置；
S43、根据步骤S41中数据矩阵与步骤S42中协方差矩阵计算各载频对应的多普勒频率，表示为：
Figure FDA0003259050590000033
其中，
Figure FDA0003259050590000034
为第m个阵元载频对应的多普勒频率，
Figure FDA0003259050590000035
为取该函数最大值所对应的多普勒频率fd，
Figure FDA0003259050590000036
ω(fd)为所有可能的多普勒频率矢量，(·)*为共轭；
S44、根据步骤S43中多普勒频率计算存在多普勒模糊的各载频对应的估计速度，表示为：
Figure FDA0003259050590000041
其中，vm为第m个阵元载频对应的估计速度。
6.根据权利要求5所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S41具体为：
提取步骤S3中信号数据各列数据，得到各载频对应的数据矢量，并在得到全部脉冲数据后，构建各载频对应的数据矩阵。
7.根据权利要求5所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S5具体包括以下分步骤：
S51、根据步骤S3中信号数据计算无多普勒模糊的初始估计速度，表示为：
Figure FDA0003259050590000042
其中，xk为无多普勒模糊的初始估计速度，即为：各个载频回波数据间的相位差，
Figure FDA0003259050590000043
为步骤S32中目标所在距离门的信号数据Yk的矩阵元素yj,i,k共轭，yj,i+1,k为步骤S32中目标所在距离门的信号数据Yk的矩阵元素；
S52、对步骤S51中初始估计速度进行快速傅里叶变换，并选取变换后频谱最大值所对应的频率计算无多普勒模糊的估计速度，表示为：
Figure FDA0003259050590000044
其中，vx为无多普勒模糊的估计速度，FFT(xk)为对初始估计速度xk进行快速傅里叶变换。
8.根据权利要求7所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S6具体包括以下分步骤：
S61、根据步骤S4中存在多普勒模糊的各载频对应的估计速度与步骤S5中无多普勒模糊的估计速度计算模糊参数，并确定搜索范围；
S62、根据步骤S61确定的搜索范围，将搜索范围内目标速度值进行排序，并分组计算均值与方差；
S63、利用一维集法根选取步骤S62中各组方差内最小值所对应的均值作为目标速度估计值。
9.根据权利要求8所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S61中模糊参数计算式表示为：
Figure FDA0003259050590000051
其中，
Figure FDA0003259050590000052
为模糊参数，[·]为取整函数，
Figure FDA0003259050590000053
为第m个阵元载频所对应的无模糊估计的最大速度。
10.根据权利要求8所述的基于FDA-MIMO雷达的速度解模糊方法，其特征在于，步骤S62中方差的计算式表示为：
Figure FDA0003259050590000054
其中，Cv(j)为第j数据组的方差，
Figure FDA0003259050590000055
为第j数据组的第m*个速度值，Vj为搜索范围内MT个连续速度值构建的第j数据组，
Figure FDA0003259050590000061
为第j数据组的平均值。
CN202111067790.1A 2021-09-13 2021-09-13 一种基于fda-mimo雷达的速度解模糊方法 Active CN113777599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111067790.1A CN113777599B (zh) 2021-09-13 2021-09-13 一种基于fda-mimo雷达的速度解模糊方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111067790.1A CN113777599B (zh) 2021-09-13 2021-09-13 一种基于fda-mimo雷达的速度解模糊方法

Publications (2)

Publication Number Publication Date
CN113777599A true CN113777599A (zh) 2021-12-10
CN113777599B CN113777599B (zh) 2023-06-16

Family

ID=78842965

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111067790.1A Active CN113777599B (zh) 2021-09-13 2021-09-13 一种基于fda-mimo雷达的速度解模糊方法

Country Status (1)

Country Link
CN (1) CN113777599B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113341418A (zh) * 2021-05-21 2021-09-03 南京航空航天大学 基于dbf机载气象雷达前视快速扫描下的解模糊方法
CN114814817A (zh) * 2022-04-25 2022-07-29 南京隼眼电子科技有限公司 运动目标速度解模糊方法、装置、电子设备及存储介质
CN116299303A (zh) * 2023-05-19 2023-06-23 南京隼眼电子科技有限公司 速度解模糊方法、装置、雷达设备及存储介质

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106353744A (zh) * 2016-10-28 2017-01-25 中国人民解放军信息工程大学 基于双基地fda‑mimo雷达的多参数联合估计方法
WO2017050246A1 (zh) * 2015-09-24 2017-03-30 索尼公司 用于无线通信的电子设备以及无线通信方法
CN109901149A (zh) * 2019-03-25 2019-06-18 西安电子科技大学 一种基于fda-mimo雷达的目标参数估计方法
EP3514569A1 (en) * 2018-01-22 2019-07-24 Infineon Technologies AG System and method for human behavior modelling and power control using a millimeter-wave radar sensor
CN110109078A (zh) * 2019-05-05 2019-08-09 西安电子科技大学 基于相参频率捷变雷达的速度解模糊方法
CN110133631A (zh) * 2019-06-11 2019-08-16 电子科技大学 一种基于模糊函数的频控阵mimo雷达目标定位方法
CN110146871A (zh) * 2019-05-21 2019-08-20 西安电子科技大学 基于双频偏fda-mimo雷达的目标参数估计方法
CN110221283A (zh) * 2019-05-14 2019-09-10 电子科技大学 一种基于fda-mimo的运动目标检测方法
CN113093137A (zh) * 2021-04-02 2021-07-09 电子科技大学 一种基于fda-mimo雷达最优频偏估计的杂波抑制方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017050246A1 (zh) * 2015-09-24 2017-03-30 索尼公司 用于无线通信的电子设备以及无线通信方法
CN106353744A (zh) * 2016-10-28 2017-01-25 中国人民解放军信息工程大学 基于双基地fda‑mimo雷达的多参数联合估计方法
EP3514569A1 (en) * 2018-01-22 2019-07-24 Infineon Technologies AG System and method for human behavior modelling and power control using a millimeter-wave radar sensor
CN109901149A (zh) * 2019-03-25 2019-06-18 西安电子科技大学 一种基于fda-mimo雷达的目标参数估计方法
CN110109078A (zh) * 2019-05-05 2019-08-09 西安电子科技大学 基于相参频率捷变雷达的速度解模糊方法
CN110221283A (zh) * 2019-05-14 2019-09-10 电子科技大学 一种基于fda-mimo的运动目标检测方法
CN110146871A (zh) * 2019-05-21 2019-08-20 西安电子科技大学 基于双频偏fda-mimo雷达的目标参数估计方法
CN110133631A (zh) * 2019-06-11 2019-08-16 电子科技大学 一种基于模糊函数的频控阵mimo雷达目标定位方法
CN113093137A (zh) * 2021-04-02 2021-07-09 电子科技大学 一种基于fda-mimo雷达最优频偏估计的杂波抑制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XU, ZT 等: "Performance Prediction of FDA-MIMO Radar Detector", IEEE RADAR CONFERENCE (RADARCONF) *
王文钦 等: "频控阵雷达技术及其应用研究进展", 雷达学报, vol. 7, no. 2, pages 153 - 166 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113341418A (zh) * 2021-05-21 2021-09-03 南京航空航天大学 基于dbf机载气象雷达前视快速扫描下的解模糊方法
CN114814817A (zh) * 2022-04-25 2022-07-29 南京隼眼电子科技有限公司 运动目标速度解模糊方法、装置、电子设备及存储介质
CN114814817B (zh) * 2022-04-25 2023-10-20 南京隼眼电子科技有限公司 运动目标速度解模糊方法、装置、电子设备及存储介质
CN116299303A (zh) * 2023-05-19 2023-06-23 南京隼眼电子科技有限公司 速度解模糊方法、装置、雷达设备及存储介质
CN116299303B (zh) * 2023-05-19 2023-08-04 南京隼眼电子科技有限公司 速度解模糊方法、装置、雷达设备及存储介质

Also Published As

Publication number Publication date
CN113777599B (zh) 2023-06-16

Similar Documents

Publication Publication Date Title
CN113777599B (zh) 一种基于fda-mimo雷达的速度解模糊方法
CN110412558B (zh) 基于tdm mimo的解车载fmcw雷达速度模糊方法
CN106353744B (zh) 基于双基地fda-mimo雷达的多参数联合估计方法
CN105807267B (zh) 一种mimo雷达扩展目标的检测方法
CN109188387B (zh) 基于插值补偿的分布式相参雷达目标参数估计方法
CN102520395B (zh) 基于双基地多输入多输出雷达的杂波抑制方法
CN103383449B (zh) 基于esprit算法的机载雷达近程杂波抑制方法
CN105137410B (zh) 基于ofdm的高分辨率雷达通信一体化的波形优化方法
CN109581362B (zh) 合成孔径雷达高度计在可变脉冲簇模式下的信号处理方法
DE102019112469A1 (de) Auflösen von doppler-mehrdeutigkeiten in einem multi-input-multi-output-radar unter verwendung digitaler mehrfach-impulsfolgefrequenzen
CN104730517B (zh) 双基地多输入多输出雷达多目标跟踪方法
CN112444810B (zh) 一种雷达对空多目标超分辨方法
CN108693511A (zh) 时分复用mimo雷达的运动目标角度计算方法
CN106199600A (zh) 基于多普勒估计的方位多通道合成孔径雷达成像方法
CN104297740A (zh) 基于相位分析的雷达目标多普勒谱估计方法
CN109001687A (zh) 基于广义旁瓣相消结构的机载雷达空时自适应滤波方法
CN104849711A (zh) 基于频域的i-ofdm mimo雷达信号的多普勒补偿方法
CN109001671B (zh) 一种跳频信号的目标检测和参数估计方法及装置
CN105572668B (zh) 一种基于Radon变换的运动目标径向速度估计方法
CN107356908A (zh) 一种捷变频信号相参积累方法
CN113189592A (zh) 考虑幅相互耦误差的车载毫米波mimo雷达测角方法
CN106802408B (zh) 基于稀疏恢复的机载非正侧阵近程杂波距离模糊抑制方法
CN114609595A (zh) 一种频分正交mimo雷达信号处理方法
CN114545342A (zh) 利用多通道侦察接收机的雷达脉冲信号参数测量方法
CN104202103A (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