CN112882007B - Realization method of single-pulse multi-target super-resolution angle measurement based on sparse array radar - Google Patents
Realization method of single-pulse multi-target super-resolution angle measurement based on sparse array radar Download PDFInfo
- Publication number
- CN112882007B CN112882007B CN202110038771.XA CN202110038771A CN112882007B CN 112882007 B CN112882007 B CN 112882007B CN 202110038771 A CN202110038771 A CN 202110038771A CN 112882007 B CN112882007 B CN 112882007B
- Authority
- CN
- China
- Prior art keywords
- target
- angle
- array radar
- matrix
- array
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 63
- 239000013598 vector Substances 0.000 claims abstract description 45
- 238000001514 detection method Methods 0.000 claims abstract description 20
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 238000003491 array Methods 0.000 claims description 16
- 238000000691 measurement method Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 21
- 238000002592 echocardiography Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000001934 delay Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
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/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
- G01S13/44—Monopulse radar, i.e. simultaneous lobing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
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
Description
技术领域technical field
本发明涉及雷达信号处理技术领域,尤其涉及一种基于稀布阵雷达的单脉冲多目标超分辨测角实现方法。The invention relates to the technical field of radar signal processing, in particular to a method for implementing single-pulse multi-target super-resolution angle measurement based on sparse array radar.
背景技术Background technique
阵列雷达在目标探测领域有着广泛的应用。在雷达探测技术中,由于目标的距离像受距离分辨率限制,角度分辨能力受雷达波束宽度约束,因此,密集多目标的高分辨技术是雷达探测领域的一项难点。Array radar has a wide range of applications in the field of target detection. In the radar detection technology, since the range image of the target is limited by the range resolution, and the angular resolution is limited by the radar beam width, the high-resolution technology of dense multi-targets is a difficult point in the field of radar detection.
目前,现有技术中的基于阵列雷达的目标检测与测角方法主要有基于一维距离像的恒虚警检测技术和基于阵列的波束扫描法等技术。上述目标检测与测角方法在密集多目标探测中都有一定的缺点。基于一维距离像的恒虚警检测技术无法实现波束内距离相近的多个目标分辨;基于阵列的波束扫描法无法对角度相近的多个目标进行准确测量,当多个目标之间的距离小于一个距离单元间隔,角度间隔小于半功率波束宽度时,上述方法都不能对空间中多个目标进行正确分辨与准确角度测量。At present, the target detection and angle measurement methods based on array radar in the prior art mainly include constant false alarm detection technology based on one-dimensional range image and beam scanning method based on array. The above target detection and angle measurement methods have certain shortcomings in dense multi-target detection. The constant false alarm detection technology based on the one-dimensional range image cannot realize the resolution of multiple targets with similar distances in the beam; the array-based beam scanning method cannot accurately measure multiple targets with similar angles, when the distance between multiple targets is less than When the interval between a range unit and the angle interval is smaller than the half-power beam width, none of the above methods can correctly distinguish and accurately measure multiple targets in space.
综上所述,如何基于阵列雷达探测体制,实现对空间密集多目标的超分辨和角度精确测量是一项亟待解决的问题。To sum up, how to achieve super-resolution and angle-accurate measurement of space-dense multi-targets based on the array radar detection system is an urgent problem to be solved.
发明内容Contents of the invention
本发明的实施例提供了一种基于稀布阵雷达的单脉冲多目标超分辨测角实现方法,以实现对雷达单脉冲回波数据进行多目标分辨与角度测量。Embodiments of the present invention provide a method for implementing single-pulse multi-target super-resolution angle measurement based on sparse array radar, so as to realize multi-target resolution and angle measurement for radar single-pulse echo data.
为了实现上述目的,本发明采取了如下技术方案。In order to achieve the above object, the present invention adopts the following technical solutions.
一种基于稀布阵雷达的单脉冲多目标超分辨测角实现方法,该方法包括:A method for implementing single-pulse multi-target super-resolution angle measurement based on sparse array radar, the method comprising:
对单脉冲稀布阵雷达多个通道的目标回波信号进行一维距离像恒虚警检测,得到各个目标所在距离单元的阵列信号,根据目标所在距离单元的阵列信号构建目标的角度观测向量;One-dimensional range image constant false alarm detection is performed on the target echo signals of multiple channels of the single pulse sparse array radar, and the array signals of the range units where each target is located are obtained, and the angular observation vector of the target is constructed according to the array signals of the range units where the target is located;
利用稀布阵雷达系统的工作参数构建角度测量字典矩阵,根据角度测量字典矩阵与稀布阵雷达的工作参数构建角度观测矩阵;The angle measurement dictionary matrix is constructed by using the working parameters of the sparse array radar system, and the angle measurement matrix is constructed according to the angle measurement dictionary matrix and the working parameters of the sparse array radar;
基于所述目标的角度观测向量和所述角度观测矩阵,利用优化算法对不同目标的相位延迟向量进行重构成像,实现所述目标回波信号的多目标超分辨与角度测量。Based on the angle observation vector of the target and the angle observation matrix, an optimization algorithm is used to reconstruct and image the phase delay vectors of different targets, so as to realize multi-target super-resolution and angle measurement of the target echo signals.
优选地,所述的对单脉冲稀布阵雷达多个通道的目标回波信号进行一维距离像恒虚警检测,得到各个目标所在距离单元的阵列信号,根据目标所在距离单元的阵列信号构建目标的角度观测向量,包括:Preferably, the one-dimensional range image constant false alarm detection is performed on the target echo signals of multiple channels of the single-pulse sparse array radar, and the array signals of the distance units where each target is located are obtained. The angle observation vector of the target, including:
对单脉冲稀布阵雷达多个通道的目标回波信号进行数字波束合成,对数字波束合成后的回波数据进行匹配滤波,得到目标一维距离像数据;Perform digital beamforming on the target echo signals of multiple channels of the single pulse sparse array radar, and perform matching filtering on the echo data after digital beamforming to obtain the target one-dimensional range image data;
对目标一维距离像数据中的每个距离单元的回波数据进行一维距离CFAR恒虚警检测,检测出不同距离的目标,得到各个目标所在距离单元的阵列信号。The one-dimensional range CFAR constant false alarm detection is performed on the echo data of each range unit in the target one-dimensional range image data to detect targets at different distances and obtain the array signals of the range units where each target is located.
利用检测到的目标距离单元各阵列信号,构建目标的角度观测向量。The angle observation vector of the target is constructed by using the detected array signals of the target range unit.
优选地,所述目标的角度观测向量表示为:Preferably, the angular observation vector of the target is expressed as:
S(k)=[s1(k) s2(k) … sP(k)] (1)S(k)=[s 1 (k) s 2 (k) … s P (k)] (1)
式(1)中,sp(k)为一维距离像CFAR检测出的第k个目标在稀布阵雷达第p个通道的距离单元复信号数据,p=1,…,P,P<N,N为阵列雷达全阵通道数。In formula (1), sp (k) is the range unit complex signal data of the kth target detected by the one-dimensional range image CFAR in the pth channel of the sparse array radar, p=1,...,P, P< N, N is the number of channels in the array radar array.
优选地,所述利用稀布阵雷达系统的工作参数构建角度测量字典矩阵,包括:Preferably, the construction of the angle measurement dictionary matrix using the working parameters of the sparse array radar system includes:
根据稀布阵雷达系统对应的全阵雷达阵列排布、工作波长、全阵子阵个数、相邻子阵间距及测角分辨间隔构建基于阵列雷达的角度测量字典矩阵。该角度测量字典矩阵表示为一N×M的矩阵:The angle measurement dictionary matrix based on the array radar is constructed according to the arrangement of the full array radar array corresponding to the sparse array radar system, the working wavelength, the number of full array sub-arrays, the distance between adjacent sub-arrays and the angular resolution interval. The angle measurement dictionary matrix is represented as a matrix of N×M:
式(2)中,为角度θm的导向矢量,θm=(m-N/2)Δθ,M为空间角度划分个数,Δθ为空间角度划分间隔,d为阵列雷达相邻两子阵间距,λ为雷达工作频率。In formula (2), is the steering vector of angle θ m , θ m = (mN/2)Δθ, M is the number of space angle divisions, Δθ is the space angle division interval, d is the distance between two adjacent sub-arrays of the array radar, and λ is the radar operating frequency.
优选地,所述的利用角度测量字典矩阵与稀布阵雷达的工作参数构建角度观测矩阵,包括:Preferably, the use of the angle measurement dictionary matrix and the working parameters of the sparse array radar to construct the angle observation matrix includes:
利用角度测量字典矩阵与稀布阵雷达的布阵位置、子阵个数构建对应的角度观测矩阵,角度观测矩阵表示为:Using the angle measurement dictionary matrix, the array position and the number of sub-arrays of the sparse array radar to construct the corresponding angle observation matrix, the angle observation matrix is expressed as:
A=[Φ1 Φ2 … ΦP]T (3)A=[Φ 1 Φ 2 ... Φ P ] T (3)
式(3)中,[·]T表示矩阵的转置,Φp表示角度测量字典矩阵Φ的第p个行向量,p=1,…,P,P<N。In formula (3), [·] T represents the transposition of the matrix, Φ p represents the pth row vector of the angle measurement dictionary matrix Φ, p=1,...,P, P<N.
优选地,所述的基于所述目标的角度观测向量和所述角度观测矩阵,利用优化算法对不同目标的相位延迟向量进行重构成像,实现所述目标回波信号的多目标超分辨与角度测量,包括:Preferably, based on the angle observation vector of the target and the angle observation matrix, an optimization algorithm is used to reconstruct and image the phase delay vectors of different targets, so as to realize the multi-target super-resolution and angle Measurements, including:
利用稀布阵雷达的目标角度观测向量和角度观测矩阵,采用最优化算法对不同目标的相位延迟向量在观测矩阵基向量下的投影系数进行重构求解;Using the target angle observation vector and angle observation matrix of the sparse array radar, the optimization algorithm is used to reconstruct the projection coefficient of the phase delay vector of different targets under the base vector of the observation matrix;
以最小l1范数为准则,根据稀布阵雷达的角度观测矩阵A和角度观测向量S(k),利用最优化算法求解下述最优化模型:Taking the minimum l 1 norm as the criterion, according to the angle observation matrix A and angle observation vector S(k) of the sparse array radar, the optimization algorithm is used to solve the following optimization model:
式(4)中,为空间各角度单元的目标后向散射系数估计值,argmin为取函数最小值,根据求得的最优解/>将第k个距离单元处的各目标进行分辨,得到各目标的角度。In formula (4), is the estimated value of the backscattering coefficient of the target in each angular unit of the space, argmin is the minimum value of the function, according to the obtained optimal solution /> Each target at the kth distance unit is distinguished to obtain the angle of each target.
由上述本发明的实施例提供的技术方案可以看出,本发明实施例提供了一种基于稀布阵列雷达的密集多目标分辨和测角方法,能够从雷达单脉冲回波数据中对波束内距离和角度相近的多个目标实现角度分辨与精确测量。It can be seen from the technical solutions provided by the above-mentioned embodiments of the present invention that the embodiments of the present invention provide a dense multi-target resolution and angle measurement method based on a sparse array radar, which can detect objects within a beam from the radar single-pulse echo data. Multiple targets with similar distances and angles achieve angular resolution and precise measurement.
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。Additional aspects and advantages of the invention will be set forth in part in the description which follows, and will become apparent from the description, or may be learned by practice of the invention.
附图说明Description of drawings
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following will briefly introduce the accompanying drawings that need to be used in the description of the embodiments. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For Those of ordinary skill in the art can also obtain other drawings based on these drawings without making creative efforts.
图1为本发明实施例提供的一种基于稀布阵雷达的多目标超分辨测角实现方法的原理图;Fig. 1 is a schematic diagram of a multi-target super-resolution angle measurement method based on sparse array radar provided by an embodiment of the present invention;
图2为本发明实施例提供的一种稀布阵雷达对空间密集多目标的观测模型示意图;FIG. 2 is a schematic diagram of an observation model of a sparse array radar for space-dense multi-targets provided by an embodiment of the present invention;
图3为本发明实施例提供的一种波束内目标各阵列回波波程延迟示意图;FIG. 3 is a schematic diagram of wave path delays of echoes from each array of targets within a beam provided by an embodiment of the present invention;
图4为本发明实施例提供的一种基于稀布阵雷达单脉冲多目标回波数据进行脉冲压缩后的一维距离像示意图;Fig. 4 is a schematic diagram of a one-dimensional range image after pulse compression based on the sparse array radar single-pulse multi-target echo data provided by the embodiment of the present invention;
图5为本发明实施例提供的一种利用稀布阵雷达系统工作参数构建的角度测量字典矩阵数据示意图;Fig. 5 is a schematic diagram of angle measurement dictionary matrix data constructed by utilizing the operating parameters of the sparse array radar system provided by the embodiment of the present invention;
图6为本发明实施例提供的一种利用稀布阵雷达系统工作参数构建的角度观测矩阵数据示意图;Fig. 6 is a schematic diagram of angle observation matrix data constructed by using the working parameters of the sparse array radar system provided by the embodiment of the present invention;
图7为本发明实施例提供的一种检测该算法的单目标稀布阵雷达角度观测数据和角度重构结果图;左图为稀布阵雷达距离单元内单目标的角度观测数据;右图为利用优化算法对单目标进行角度重构成像的结果示意图;Figure 7 is a diagram of the angle observation data and angle reconstruction results of a single-target sparse array radar for detecting the algorithm provided by the embodiment of the present invention; the left figure is the angle observation data of a single target in the range unit of the sparse array radar; the right figure Schematic diagram of the results of angle reconstruction imaging of a single target using the optimization algorithm;
图8为本发明实施例提供的一种检测该算法的距离相同的双目标稀布阵雷达角度观测数据和角度重构结果图;左图为稀布阵雷达距离单元内双目标的角度观测数据;右图为利用优化算法对两目标进行角度重构成像的结果示意图。Figure 8 is a diagram of the angle observation data and angle reconstruction results of a dual-target sparse array radar with the same distance to detect the algorithm provided by the embodiment of the present invention; the left figure is the angle observation data of dual targets in the range unit of the sparse array radar ; The right figure is a schematic diagram of the angle reconstruction imaging results of two targets using the optimization algorithm.
具体实施方式Detailed ways
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。Embodiments of the present invention are described in detail below, examples of which are shown in the drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the figures are exemplary only for explaining the present invention and should not be construed as limiting the present invention.
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。Those skilled in the art will understand that unless otherwise stated, the singular forms "a", "an", "said" and "the" used herein may also include plural forms. It should be further understood that the word "comprising" used in the description of the present invention refers to the presence of said features, integers, steps, operations, elements and/or components, but does not exclude the presence or addition of one or more other features, Integers, steps, operations, elements, components, and/or groups thereof.
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。Those skilled in the art can understand that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It should also be understood that terms such as those defined in commonly used dictionaries should be understood to have a meaning consistent with the meaning in the context of the prior art, and unless defined as herein, are not to be interpreted in an idealized or overly formal sense explain.
为便于对本发明实施例的理解,下面将结合附图以几个具体实施例为例做进一步的解释说明,且各个实施例并不构成对本发明实施例的限定。In order to facilitate the understanding of the embodiments of the present invention, several specific embodiments will be taken as examples for further explanation below in conjunction with the accompanying drawings, and each embodiment does not constitute a limitation to the embodiments of the present invention.
本发明实施例提供的一种基于稀布阵雷达的单脉冲多目标超分辨测角实现方法的实现原理图如图1所示,包括如下的处理步骤:The implementation schematic diagram of a single-pulse multi-target super-resolution angle measurement method based on sparse array radar provided by the embodiment of the present invention is shown in Figure 1, including the following processing steps:
步骤101:对目标回波进行一维距离像恒虚警检测,利用目标所在距离单元的各阵列信号,构建目标的角度观测向量。Step 101: Perform one-dimensional range image constant false alarm detection on the target echo, and construct an angular observation vector of the target by using the array signals of the range unit where the target is located.
一维距离像是空间距雷达不同距离处的目标反射回波强度的距离投影。对目标回波进行一维距离像CFAR检测,可以将空间距离不同的强反射目标分辨出来。同时,利用CFAR检测到的目标所在距离单元的各阵列数据,构建该目标的角度观测向量。The one-dimensional distance is like the distance projection of the target reflection echo intensity at different distances from the radar in space. The one-dimensional range image CFAR detection of the target echo can distinguish the strong reflective targets with different spatial distances. At the same time, the angular observation vector of the target is constructed by using the array data of the range unit where the target is detected by CFAR.
图2为本发明实施例一提供的一种稀布阵雷达对空间密集多目标的观测模型示意图。首先,对单脉冲稀布阵雷达多个通道的目标回波信号进行数字波束合成,提高检测数据的信噪比;对数字波束合成后的回波数据进行匹配滤波,得到目标一维距离像数据;然后,对目标一维距离像数据中的每个距离单元的回波数据进行一维距离CFAR(Constant FalseAlarm Rate,恒虚警率)检测,检测出不同距离的目标,得到各个目标所在距离单元的阵列信号。FIG. 2 is a schematic diagram of an observation model of a sparse array radar for space-dense multi-targets provided by
利用检测到的目标所在距离单元的各阵列信号,构建目标的角度观测向量,目标的角度观测向量可以表示为:Using the detected array signals of the distance unit where the target is located, the angle observation vector of the target is constructed, and the angle observation vector of the target can be expressed as:
S(k)=[s1(k) s2(k) … sP(k)] (5)S(k)=[s 1 (k) s 2 (k) … s P (k)] (5)
式(1)中,sp(k)为一维距离像CFAR检测出的第k个目标在稀布阵雷达第p个通道的距离单元复信号数据,p=1,…,P,P<N,N为阵列雷达全阵通道数。In formula (1), sp (k) is the range unit complex signal data of the kth target detected by the one-dimensional range image CFAR in the pth channel of the sparse array radar, p=1,...,P, P< N, N is the number of channels in the array radar array.
步骤102:利用稀布阵雷达系统的工作参数构建角度测量字典矩阵,根据角度测量字典矩阵与稀布阵雷达的工作参数构建角度观测矩阵。Step 102: Construct an angle measurement dictionary matrix by using the operating parameters of the sparse array radar system, and construct an angle measurement matrix according to the angle measurement dictionary matrix and the operating parameters of the sparse array radar.
根据稀布阵雷达系统对应的全阵雷达阵列排布、工作波长、全阵子阵个数、相邻子阵间距及测角分辨间隔构建基于阵列雷达的角度测量字典矩阵。该角度测量字典矩阵可以表示为一N×M的矩阵:The angle measurement dictionary matrix based on the array radar is constructed according to the arrangement of the full array radar array corresponding to the sparse array radar system, the working wavelength, the number of full array sub-arrays, the distance between adjacent sub-arrays and the angular resolution interval. The angle measurement dictionary matrix can be expressed as a matrix of N×M:
式(2)中,为角度θm的导向矢量,θm=(m-N/2)Δθ,M为空间角度划分个数,Δθ为空间角度划分间隔,d为阵列雷达相邻两子阵间距,λ为雷达工作频率。In formula (2), is the steering vector of angle θ m , θ m = (mN/2)Δθ, M is the number of space angle divisions, Δθ is the space angle division interval, d is the distance between two adjacent sub-arrays of the array radar, and λ is the radar operating frequency.
然后,利用角度测量字典矩阵与稀布阵雷达的布阵位置、子阵个数构建角度观测矩阵,角度观测矩阵可表示为:Then, the angle measurement matrix is constructed by using the angle measurement dictionary matrix, the array position and the number of sub-arrays of the sparse array radar, and the angle observation matrix can be expressed as:
A=[Φ1 Φ2 … ΦP]T (7)A=[Φ 1 Φ 2 ... Φ P ] T (7)
式(3)中,[·]T表示矩阵的转置,Φp表示角度测量字典矩阵Φ的第p个行向量,p=1,…,P,P<N。In formula (3), [·] T represents the transposition of the matrix, Φ p represents the pth row vector of the angle measurement dictionary matrix Φ, p=1,...,P, P<N.
步骤103:基于目标角度观测向量和角度观测矩阵,利用优化算法对目标角度进行重构成像,实现多目标的超分辨与角度测量。Step 103: Based on the target angle observation vector and the angle observation matrix, use an optimization algorithm to reconstruct the image of the target angle, and realize multi-target super-resolution and angle measurement.
图3为本发明实施例一提供的一种波束内目标各阵列回波波程延迟示意图。当波束内多个目标距离相同,空间角度相近时,多个目标的雷达一维距离像在距离单元上出现重叠,此时,一维距离像无法分辨距离相同的多个目标。同时,由于多个目标的角度间隔小于雷达波束宽度,依靠雷达波束宽度也无法分辨多个目标。在此情况下,由于不同角度的目标回波到雷达各阵列的相位延迟不同,因此,不同目标的回波在各阵列中的相位差在系统测量字典矩阵基向量下的投影也不同。基于上述性质,利用稀布阵雷达的目标角度观测向量和角度观测矩阵,采用最优化算法对不同目标的相位延迟向量在观测矩阵基向量下的投影系数进行重构求解,就可将距离相同、角度相近的多个目标进行超分辨,并可实现各目标角度的精确测量。FIG. 3 is a schematic diagram of wave path delays of echoes from each array of targets within a beam provided by
具体的,以最小l1范数为准则,根据稀布阵雷达的角度观测矩阵A和角度观测向量S(k),利用最优化算法求解下述最优化模型:Specifically, with the minimum l 1 norm as the criterion, according to the angle observation matrix A and the angle observation vector S(k) of the sparse array radar, the optimization algorithm is used to solve the following optimization model:
式(4)中,为空间各角度单元的目标后向散射系数估计值,arg min为取函数最小值,根据求得的最优解/>就可以将第k个距离单元处的各目标进行分辨,并可得到各目标的精确角度。In formula (4), is the estimated value of the backscatter coefficient of the target in each angular unit of the space, arg min is the minimum value of the function, according to the obtained optimal solution /> Then the targets at the kth distance unit can be distinguished, and the precise angle of each target can be obtained.
下面结合具体实施例对本发明的技术方案作进一步详细说明。The technical solutions of the present invention will be further described in detail below in conjunction with specific embodiments.
实施例1Example 1
图4为检验本发明实施例的算法的稀布阵雷达一维距离像数据示意图。其中,稀布阵雷达工作频率为9.2GHz,发射信号带宽3MHz,波束宽度为0.45°,波束内3个目标距离雷达分别为:[10、10、15]km,目标方位角度分别为:[1.0、1.3、0.9]°。由图中可以看出,第1,2个目标由于距离相同,一维距离像无法分辨,而第3个目标从一维距离像中可以与前两个目标进行分辨。Fig. 4 is a schematic diagram of one-dimensional range image data of sparse array radar for testing the algorithm of the embodiment of the present invention. Among them, the operating frequency of the sparse array radar is 9.2GHz, the transmission signal bandwidth is 3MHz, and the beam width is 0.45°. The distances of the three targets in the beam to the radar are: [10, 10, 15]km, and the azimuth angles of the targets are: [1.0 , 1.3, 0.9]°. It can be seen from the figure that the 1st and 2nd targets are indistinguishable from the one-dimensional range image due to the same distance, while the third target can be distinguished from the first two targets from the one-dimensional range image.
图5为对利用全阵列雷达系统参数构建的角度测量字典矩阵的数据示意图,其中全阵雷达子阵个数为128,字典矩阵角度范围为[-5,5]°,角度划分间隔为0.1°。Figure 5 is a data schematic diagram of the angle measurement dictionary matrix constructed by using the parameters of the full-array radar system, wherein the number of sub-arrays of the full-array radar is 128, the angle range of the dictionary matrix is [-5,5]°, and the angle division interval is 0.1° .
图6为利用稀布阵雷达系统阵列排布构建的角度观测矩阵,其中,稀布阵雷达子阵个数为16个,为全阵雷达子阵数的1/8。Figure 6 shows the angular observation matrix constructed by using the array arrangement of the sparse array radar system. The number of sparse array radar sub-arrays is 16, which is 1/8 of the number of full-array radar sub-arrays.
图7为利用第3个目标所在距离单元的稀布阵雷达角度测量数据和重构求得的目标角度图像示意图。其中,左图为全阵雷达和稀布阵雷达对该目标的角度观测向量数据图,从图7中可以看出稀布阵雷达为全阵雷达观测数据量的1/8。右图为重构求得的该目标角度图像数据,从图7中可以看出,由于该距离单元处只有一个目标,角度重构图像中只有一个峰值,峰值位置为该目标的角度精确值。Fig. 7 is a schematic diagram of the target angle image obtained by using the sparse array radar angle measurement data and reconstruction of the distance unit where the third target is located. Among them, the left figure is the angular observation vector data diagram of the target by the full-array radar and the sparse-array radar. From Figure 7, it can be seen that the sparse-array radar is 1/8 of the observation data of the full-array radar. The figure on the right shows the angle image data of the target obtained by reconstruction. It can be seen from Figure 7 that since there is only one target at this distance unit, there is only one peak in the angle reconstruction image, and the peak position is the precise angle value of the target.
图8为利用第1、2个目标所在距离单元的稀布阵雷达角度测量数据和重构求得的目标角度图像示意图。其中,左图为全阵雷达和稀布阵雷达该距离单元的角度观测向量数据图,右图为重构求得的该距离单元处各目标的角度数据图。从图8中可以看出,由于该距离单元包含第1、2共两个目标,重构角度图像中出现了两个峰值,可将两个目标进行清晰分辨,且各峰值位置为两目标对应的角度准确值。稀布阵雷达仅利用全阵雷达1/8的数据量,实现了距离相同、角度小于波束宽度的两个目标的角度超分辨和精确测量。Fig. 8 is a schematic diagram of the target angle image obtained by using the sparse array radar angle measurement data and reconstruction of the first and second target distance units. Among them, the left figure is the angle observation vector data diagram of the range unit of the full-array radar and the sparse array radar, and the right figure is the angle data diagram of each target at the range unit obtained by reconstruction. It can be seen from Figure 8 that since the distance unit contains two
综上所述,本发明实施例提供了一种基于稀布阵列雷达的密集多目标分辨和测角方法,能够从雷达单脉冲回波数据中对波束内距离和角度相近的多个目标进行角度分辨与精确测量。In summary, the embodiment of the present invention provides a dense multi-target resolution and angle measurement method based on sparse array radar, which can measure the angles of multiple targets with similar distances and angles in the beam from the radar single pulse echo data. Resolution and precise measurement.
通过本发明实施例的方法,不仅可极大降低阵列雷达系统数据量,还可对波束内多个距离相同角度相近的多个目标进行超分辨和精确角度测量。Through the method of the embodiment of the present invention, not only the amount of data in the array radar system can be greatly reduced, but also super-resolution and precise angle measurement can be performed on multiple targets within the beam at the same distance and close to the angle.
本领域普通技术人员可以理解:附图只是一个实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。Those skilled in the art can understand that the accompanying drawing is only a schematic diagram of an embodiment, and the modules or processes in the accompanying drawing are not necessarily necessary for implementing the present invention.
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。It can be known from the above description of the implementation manners that those skilled in the art can clearly understand that the present invention can be implemented by means of software plus a necessary general hardware platform. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art can be embodied in the form of software products, and the computer software products can be stored in storage media, such as ROM/RAM, disk , CD, etc., including several instructions to make a computer device (which may be a personal computer, server, or network device, etc.) execute the methods described in various embodiments or some parts of the embodiments of the present invention.
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。Each embodiment in this specification is described in a progressive manner, the same and similar parts of each embodiment can be referred to each other, and each embodiment focuses on the differences from other embodiments.
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。The above is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any person skilled in the art within the technical scope disclosed in the present invention can easily think of changes or Replacement should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be determined by the protection scope of the claims.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110038771.XA CN112882007B (en) | 2021-01-12 | 2021-01-12 | Realization method of single-pulse multi-target super-resolution angle measurement based on sparse array radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110038771.XA CN112882007B (en) | 2021-01-12 | 2021-01-12 | Realization method of single-pulse multi-target super-resolution angle measurement based on sparse array radar |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112882007A CN112882007A (en) | 2021-06-01 |
CN112882007B true CN112882007B (en) | 2023-06-23 |
Family
ID=76044727
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110038771.XA Active CN112882007B (en) | 2021-01-12 | 2021-01-12 | Realization method of single-pulse multi-target super-resolution angle measurement based on sparse array radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112882007B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114779164B (en) * | 2022-06-17 | 2022-09-16 | 中国电子科技集团公司第五十四研究所 | Wide-area random sparse array multi-target positioning method based on energy separation |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010286403A (en) * | 2009-06-12 | 2010-12-24 | Mitsubishi Electric Corp | Angle measuring device, monopulse angle measuring device, monopulse radar, multi-static radar |
CN107991659A (en) * | 2017-11-17 | 2018-05-04 | 西安电子科技大学 | The low elevation angle target of metre wave radar based on dictionary learning surveys high method |
CN109061622A (en) * | 2018-06-22 | 2018-12-21 | 电子科技大学 | A kind of concealed target multipath detection method based on millimetre-wave radar |
CN109212526A (en) * | 2018-10-17 | 2019-01-15 | 哈尔滨工业大学 | Distributive array target angle measurement method for high-frequency ground wave radar |
CN110488255A (en) * | 2019-09-09 | 2019-11-22 | 上海无线电设备研究所 | A kind of phased-array radar pulse high-resolution angle measuring system and method |
CN111175745A (en) * | 2020-01-19 | 2020-05-19 | 北京航空航天大学 | A method of three-dimensional imaging of moving target based on state-space balance method |
CN111175753A (en) * | 2020-02-24 | 2020-05-19 | 四川九洲电器集团有限责任公司 | Vehicle-mounted anti-collision radar wave-splitting target detection method |
-
2021
- 2021-01-12 CN CN202110038771.XA patent/CN112882007B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010286403A (en) * | 2009-06-12 | 2010-12-24 | Mitsubishi Electric Corp | Angle measuring device, monopulse angle measuring device, monopulse radar, multi-static radar |
CN107991659A (en) * | 2017-11-17 | 2018-05-04 | 西安电子科技大学 | The low elevation angle target of metre wave radar based on dictionary learning surveys high method |
CN109061622A (en) * | 2018-06-22 | 2018-12-21 | 电子科技大学 | A kind of concealed target multipath detection method based on millimetre-wave radar |
CN109212526A (en) * | 2018-10-17 | 2019-01-15 | 哈尔滨工业大学 | Distributive array target angle measurement method for high-frequency ground wave radar |
CN110488255A (en) * | 2019-09-09 | 2019-11-22 | 上海无线电设备研究所 | A kind of phased-array radar pulse high-resolution angle measuring system and method |
CN111175745A (en) * | 2020-01-19 | 2020-05-19 | 北京航空航天大学 | A method of three-dimensional imaging of moving target based on state-space balance method |
CN111175753A (en) * | 2020-02-24 | 2020-05-19 | 四川九洲电器集团有限责任公司 | Vehicle-mounted anti-collision radar wave-splitting target detection method |
Non-Patent Citations (1)
Title |
---|
Subarray-Based Frequency Diverse Array for Target Range-Angle Localization With Monopulse Processing;Shengbin Luo Wang et al.;《IEEE SENSORS JOURNAL》;20180715;第18卷(第14期);第5937-5947页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112882007A (en) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8193967B2 (en) | Method and system for forming very low noise imagery using pixel classification | |
CN109581352B (en) | Super-resolution angle measurement system based on millimeter wave radar | |
CN110780298B (en) | Multi-base ISAR fusion imaging method based on variational Bayes learning | |
US11860270B2 (en) | Sparse linear array approach in automotive radars using matrix completion | |
CN110806577B (en) | Focusing imaging method and device of synthetic aperture radar, equipment and storage medium | |
CN110879391B (en) | Radar image dataset creation method based on electromagnetic simulation and missile echo simulation | |
CN116500626B (en) | Multi-receiving array element data receiving-transmitting combination conversion method | |
CN107278271A (en) | Clutter recognition in ultrasonic image-forming system | |
CA2787425A1 (en) | Super-resolution imaging radar | |
CN108845318B (en) | Satellite-borne high-resolution wide-range imaging method based on Relax algorithm | |
CN108020834B (en) | Moving target detection method, device and electronic equipment based on improved EDPCA | |
CN113391286A (en) | Virtual aperture MIMO radar target detection method based on two-dimensional block sparse recovery | |
CN118011392A (en) | Point cloud high-resolution imaging method based on 4D millimeter wave radar | |
CN112882007B (en) | Realization method of single-pulse multi-target super-resolution angle measurement based on sparse array radar | |
CN113534077B (en) | Radar radiation source power inversion method and device and electronic equipment | |
CN112882008B (en) | Beam internal dense multi-target super-resolution implementation method based on digital array radar | |
Feng et al. | Constained adaptive monopulse algorithm based on sub-array | |
CN109061645A (en) | The virtual related sparse imaging method and system of synthetic aperture radiometer | |
CN115932852A (en) | SAR-GMTI distance ambiguity suppression method based on non-uniform sub-pulse coding | |
CN110045374B (en) | A multi-angle ultra-long spotlight SAR super-resolution imaging algorithm based on Chirp Scaling | |
Jouadé et al. | High resolution radar focusing using spectral estimation methods in wide-band and near-field configurations: Application to millimeter-wave near-range imaging | |
CN113419240A (en) | Moving target detection method based on dual-channel SAR, dual-channel SAR and storage medium | |
Yu et al. | CSAR imaging with data extrapolation and approximate GLRT techniques | |
Brandfass et al. | Real-time video SAR on curved trajectories making use of graphic core processing | |
Hosseiny et al. | Evaluation of a Signal Processing Algorithm in a Ground-Based SAR System in Simulated Environment |
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 |