阅读说明：本技术 基于和差天线体制雷达空时自适应处理的测角方法 (Angle measuring method based on sum and difference antenna system radar space-time adaptive processing ) 是由 许京伟 牛亚丽 阚庆云 李军 朱圣棋 王建新 常文豪 于 2021-09-24 设计创作，主要内容包括：本发明公开一种基于和差天线体制雷达空时自适应处理的测角方法。其方案是：获得雷达的回波数据,估计其杂波的协方差矩阵；计算主瓣指向的空时导向矢量；根据杂波协方差矩阵和空时导向矢量求出权值；将权值作用于回波数据进行杂波抑制,得到杂波抑制后和、方位差、俯仰差三个通道的输出；根据杂波抑制后和通道的输出,得到目标距离门、多普勒通道以及幅值；根据目标距离门和多普勒通道,得到方位差通道、俯仰差通道在目标处的幅值；计算方位差、俯仰差通道分别与和通道的比值,计算方位维与俯仰维鉴角曲线的斜率,并根据斜率和比值计算目标的方位角和俯仰角。本发明空间自由度少,实时性高,测角能力强,可用于杂波背景下雷达对运动目标的估计。(The invention discloses an angle measurement method based on sum and difference antenna system radar space-time adaptive processing. The scheme is as follows: acquiring echo data of a radar, and estimating a covariance matrix of clutter of the radar; calculating a space-time guide vector pointed by a main lobe; solving a weight according to the clutter covariance matrix and the space-time guide vector; the weight is acted on echo data to carry out clutter suppression, and the outputs of three channels of sum, azimuth difference and pitching difference after clutter suppression are obtained; obtaining a target range gate, a Doppler channel and an amplitude according to the clutter suppressed sum channel output; obtaining the amplitudes of the azimuth difference channel and the pitch difference channel at the target according to the target distance gate and the Doppler channel; and calculating the ratio of the azimuth difference channel and the pitch difference channel to the sum channel, calculating the slope of the azimuth dimension and pitch dimension angle identification curve, and calculating the azimuth angle and the pitch angle of the target according to the slope and the ratio. The method has the advantages of less spatial freedom, high real-time performance and strong angle measurement capability, and can be used for estimating a moving target by a radar under a clutter background.)

1. A method for angle measurement based on sum and difference antenna system radar space-time adaptive processing is characterized by comprising the following steps:

(1) obtaining echo data NxK and L, wherein N represents the number of sum-difference channels of an airborne/missile-borne radar receiving antenna, K represents the number of coherent pulses of each pulse repetition period, and L represents the number of range gate samples of sampled data;

(2) clutter suppression is carried out on echo data:

(2a) estimating a sampling covariance matrix of clutter

(2b) According to the characteristic that the response of the difference beam antenna directional diagram in the main lobe direction is zero and the difference branch components in the airspace guiding vector are all zero under the ideal condition, the sum-difference antenna system and the channel airspace guiding vector a are set_sΣAzimuth difference channel airspace guide vectorPitching difference channel airspace guide vectorRespectively as follows:

a_sΣ＝[1,0,0]^T，

the upper mark T represents transposition operation, and the output of the usable sum and difference beam antenna directional diagram in the main lobe direction under the phased array system is used as the value of a space domain steering vector;

(2c) according to the result of (2b), calculating a sum channel space-time two-dimensional guide vector a_ΣSpace-time two-dimensional guide vector of azimuth difference channelPitching difference channel space-time two-dimensional guide vector

Wherein the content of the first and second substances,denotes the kronecker operation, a_tA time domain steering vector representing the main lobe direction of the radar receiving antenna;

(2d) calculating and channel-optimal adaptive weight vector w according to the results of (2a) and (2c)_ΣAzimuth difference channel optimal adaptive weight vectorOptimal adaptive weight vector of pitch difference channel

Wherein mu is a scalar constant, and the superscript-1 represents the inversion operation;

(2e) obtaining the output after clutter suppression according to the result of (2 d):

y_Σ＝[y_1Σ,y_2Σ,…,y_lΣ，…,y_LΣ]，l＝1,2,…,L

wherein, y_ΣRepresenting the sum channel output after clutter suppression, y_lΣRepresenting the output of the ith range gate clutter suppressed sum channel in the sampled echo data block,the output of the azimuth difference channel is represented,representing the output of the l-th range gate clutter suppressed back head channel in the sampled echo data block,the output of the pitch difference channel is represented,representing the output of the pitch difference channel after the first range gate clutter suppression in the sampled echo data block;

(3) according to the output after clutter suppression, obtaining a target azimuth angle and a pitch angle:

(3a) output y from the sum channel after clutter suppression_ΣObtaining the range gate l of the target_tDoppler channel m_tAnd a widthValue y_SΣ；

(3b) Distance door l according to target_tAnd Doppler channel m_tOutput from azimuth difference channelObtaining the amplitude y of the target in the azimuth difference channel_θOutput from the differential pitch channelObtaining the amplitude of the target in the pitch difference channel

(3c) According to the amplitude y of the target sum channel_ΣThe amplitude y of the target in the azimuth difference channel_θAnd the amplitude of the target in the pitch difference channelCalculating the ratio q of the channel to the sum channel of the target prescription potential difference_θAnd the ratio of the pitch difference channel to the sum channel

(3d) Calculating the slope k of the azimuth angle-measuring curve according to the parameters of the antenna under the system of sum and difference antennas_θAnd the slope of the pitch angle curve

(3e) Ratio q of channel to sum channel according to azimuth difference_θAnd the slope k of the azimuthal mapping curve_θCalculating the azimuth angle theta of the target according to the ratio of the pitch difference channel to the sum channelAnd the slope of the pitch angle curveCalculating the pitch angle of the targetAnd finishing the estimation of the target azimuth angle and the target pitch angle.

2. The method of claim 1 wherein the clutter covariance matrix is estimated in (2a)The following maximum likelihood formula is adopted:

wherein x is_lAnd representing the echo vector corresponding to the ith range gate in the sampled echo data block, and the superscript H represents the conjugate transpose operation.

3. The method of claim 1, wherein (2c) is a of the time domain steering vector pointed by the main lobe of the radar receiving antenna_tThe expression is as follows:

wherein j represents an imaginary unit symbol, f_dIndicating the doppler frequency and the superscript T indicating the transpose operation.

4. The method of claim 1, wherein the output y of the ith range gate clutter suppressed sum channel in the echo data block is sampled in (2e)_lΣSampling the output of the first range gate clutter suppression rear potential difference channel in the echo data blockSampling output of pitch difference channel after first range gate clutter suppression in echo data blockThe expression is as follows:

wherein x is_lAnd representing the echo vector corresponding to the ith range gate in the sampled echo data block, and the superscript H represents the conjugate transpose operation.

5. The method according to claim 1, wherein the range gate/, of the target obtained in (3a)_tDoppler channel m_tAnd an amplitude y_ΣBy the output y of the sum channel after clutter suppression_lΣAnd performing simulation, and obtaining coordinates displayed at the target in the simulation diagram.

6. The method of claim 1, wherein the target obtained in (3b) is the amplitude of the azimuth difference channelThe output matrix of the azimuth difference channel is searched by MATLAB softwareAt point (l)_t,m_t) A corresponding value is obtained.

7. The method of claim 1, wherein the magnitude of the target in the pitch difference channel obtained in (3b)The output matrix of the pitch difference channel is searched by MATLAB softwareAt point (l)_t,m_t) A corresponding value is obtained.

8. A method according to claim 1, wherein the ratio q of the channel to the sum channel of the target prescribed bit difference in (3c) is calculated_θRatio of pitch difference channel to sum channel at sum targetThe formula is as follows:

wherein, y_ΣThe representation and the amplitude of the channel at the target,representing the magnitude of the azimuth difference channel at the target, representingThe pitch difference channel is the magnitude at the target.

9. The method of claim 1, wherein the computing party in (3d)Slope k of the bit dimension angle curve_θAnd the slope of the pitch angle curveThe formula is as follows:

wherein j represents an imaginary unit symbol, L_cDenotes the length of the antenna, L_rIndicating the height, theta, of the antenna₀Indicating the azimuth angle at which the main lobe of the radar receiving antenna is directed,the pitch angle of the main lobe of the radar receiving antenna is shown, and lambda represents the wavelength.

10. The method of claim 1, wherein the azimuthal angle θ and the elevation angle are calculated in (3e)The formula is as follows:

wherein, theta₀Azimuth angle, q, representing the main lobe orientation of the radar receiving antenna_θRepresenting the ratio of the potential difference channel to the sum channel, k, for the target_θSlope representing azimuth angle curveThe ratio of the total weight of the particles,represents the pitch angle of the main lobe of the radar receiving antenna,representing the ratio of the pitch difference channel to the sum channel at the target,representing the slope of the azimuthal dimension curve.

Technical Field

The invention belongs to the technical field of radar signal processing, and particularly relates to an angle measurement method based on sum and difference antenna system radar space-time adaptive processing, which can be used for parameter estimation of a radar on a moving target under a clutter background.

Background

Airborne/missile-borne radar is one of the most important sensors in modern battlefields, and because it is mounted on an airplane flying at high altitude, its advantages are incomparable with ground-based radar, however, it also faces some problems. Because the radar is in the downward-looking work, the clutter distribution range is wide, the intensity is large, and meanwhile, because of the movement of the carrier, the clutter spectrum is greatly expanded, so that targets are often submerged in the clutter, and the parameter estimation capability of the radar is seriously influenced. Therefore, how to solve the problem of angle estimation in the clutter background is also a problem to be considered.

The patent CN 109599674A of Beijing remote sensing equipment research provides a 'phased array antenna stable angle tracking method based on decoupling'. The method comprises the following steps: decoupling interruption periodically generates projectile body attitude data at a first moment to obtain a conversion matrix from an inertia rectangular coordinate system to a projectile body rectangular coordinate system; converting the inertial coordinate into the coordinate of the rectangular coordinate of the antenna array surface; calculating a beam pointing azimuth angle and a pitch angle according to the coordinates of the antenna array surface, and controlling the pointing of the antenna beam; compensating the single pulse angle measurement error generated periodically by the angle measurement interruption to the direction cosine coordinate of the antenna array surface; the compensated cosine coordinates of the direction of the antenna array surface finish coordinate inverse transformation to obtain inertia coordinates; decoupling and interrupting to generate projectile body attitude data at the moment n >1 to obtain a conversion matrix from an inertia rectangular coordinate to a projectile body rectangular coordinate and a derotation torque matrix; and calculating the beam pointing azimuth angle and the pitch angle, and controlling the antenna beam pointing direction. However, the method has the disadvantages that clutter is not suppressed, so that the target is submerged in the clutter, and the angle measurement capability is seriously affected.

The university of electronic science and technology of west ampere in its CN 106443663 a patent proposes "a method for phased array radar dimension reduction four-channel and difference beam angle measurement". Determining a phased array radar, and performing subarray division on N array elements included in the phased array radar to obtain M subarrays and the number of the array elements included in each subarray; calculating subarray level and weighting after phased array radar optimization, subarray level pitching difference weighting after phased array radar optimization, subarray level azimuth difference weighting after phased array radar optimization, subarray level double difference weighting after phased array radar optimization, subarray element level weighting after phased array radar optimization and an optimal dimensionality reduction matrix of the phased array radar; respectively calculating a directional function of a phased array radar and a wave beam, a directional function of a phased array radar azimuth difference wave beam, a directional function of a phased array radar pitch difference wave beam and a directional function of a phased array radar double difference wave beam; and respectively calculating the final elevation direction angle estimation of the phased array radar target and the final azimuth direction angle estimation of the phased array radar target. Although the method can better estimate the angle of the target, the method is only suitable for the phased array, so that the angle measurement cannot be carried out on the radar with the sum-difference antenna system.

Disclosure of Invention

The invention aims to provide an angle measurement method based on sum and difference antenna system radar space-time adaptive processing aiming at the defects of the prior art so as to improve the angle measurement capability of the sum and difference antenna system radar and the phased array system radar under the clutter background.

The technical scheme for realizing the purpose of the invention is as follows: clutter suppression is carried out through sum-difference space-time two-dimensional self-adaptive processing, then output of a target in a sum channel and an output of a difference channel are obtained, and finally a target angle is obtained by using a sum-difference angle measurement method, wherein the method comprises the following implementation steps:

(1) obtaining echo data NxK and L, wherein N represents the number of sum-difference channels of an airborne/missile-borne radar receiving antenna, K represents the number of coherent pulses of each pulse repetition period, and L represents the number of range gate samples of sampled data;

(2) clutter suppression is carried out on echo data:

(2a) estimating a sampling covariance matrix of clutter

(2b) According to the characteristic that the response of the difference beam antenna directional diagram in the main lobe direction is zero and the difference branch components in the airspace guiding vector are all zero under the ideal condition, the sum-difference antenna system and the channel airspace guiding vector a are set_sΣAzimuth difference channel airspace guide vectorPitching difference channel airspace guide vectorRespectively as follows:

a_sΣ＝[1,0,0]^T，

the upper mark T represents transposition operation, and the output of the usable sum and difference beam antenna directional diagram in the main lobe direction under the phased array system is used as the value of a space domain steering vector;

(2c) according to the result of (2b), calculating a sum channel space-time two-dimensional guide vector a_ΣSpace-time two-dimensional guide vector of azimuth difference channelPitching difference channel space-time two-dimensional guide vector

Wherein the content of the first and second substances,denotes the kronecker operation, a_tA time domain steering vector representing the main lobe direction of the radar receiving antenna;

(2d) calculating and channel-optimal adaptive weight vector w according to the results of (2a) and (2c)_ΣAzimuth difference channel optimal adaptive weight vectorOptimal adaptive weight vector of pitch difference channel

Wherein mu is a scalar constant, and the superscript-1 represents the inversion operation;

(2e) obtaining the output after clutter suppression according to the result of (2 d):

y_Σ＝[y_1Σ,y_2Σ,…,y_lΣ，…,y_LΣ]，l＝1,2,…,L

wherein, y_ΣRepresenting the sum channel output after clutter suppression, y_lΣRepresenting the output of the ith range gate clutter suppressed sum channel in the sampled echo data block,the output of the azimuth difference channel is represented,output of azimuth difference channel after representing first range gate clutter suppression in sampled echo data blockAnd then the mixture is discharged out of the furnace,the output of the pitch difference channel is represented,representing the output of the pitch difference channel after the first range gate clutter suppression in the sampled echo data block;

(3) according to the output result after clutter suppression, obtaining a target azimuth angle and a pitch angle:

(3a) output y from the sum channel after clutter suppression_ΣObtaining the range gate l of the target_tDoppler channel m_tAnd an amplitude y_SΣ；

(3b) Distance door l according to target_tAnd Doppler channel m_tOutput from azimuth difference channelObtaining the amplitude y of the target in the azimuth difference channel_θOutput from the differential pitch channelObtaining the amplitude of the target in the pitch difference channel

(3c) According to the amplitude y of the target sum channel_ΣThe amplitude y of the target in the azimuth difference channel_θAnd the amplitude of the target in the pitch difference channelCalculating the ratio q of the channel to the sum channel of the target prescription potential difference_θAnd the ratio of the pitch difference channel to the sum channel

(3d) Calculating the slope k of the azimuth angle-measuring curve according to the parameters of the antenna under the system of sum and difference antennas_θAnd pitchSlope of angular curve of vehement

(3e) Ratio q of channel to sum channel according to azimuth difference_θAnd the slope k of the azimuthal mapping curve_θCalculating the azimuth angle theta of the target according to the ratio of the elevation channel to the sum channelAnd the slope of the pitch angle curveCalculating the pitch angle of the targetAnd finishing the estimation of the target azimuth angle and the target pitch angle.

Compared with the prior art, the invention has the following advantages:

first, the invention improves the signal-to-noise ratio of the sum and difference channel output and improves the angle measurement capability of the radar under the clutter background because the sum channel and the difference channel of the sum and difference antenna system radar are subjected to self-adaptive clutter suppression.

Secondly, the echo data of the radar of the sum and difference antenna system and the data obtained by the radar of the phased array system after digital sum and difference beam forming have the same form, so the invention is also suitable for the radar of the phased array system.

Thirdly, the invention adopts the sum-difference system antenna, so the space freedom degree is less and the real-time performance is high.

Drawings

FIG. 1 is a flow chart of an implementation of the present invention;

FIG. 2 is a simulation of the output of the present invention after clutter suppression;

FIG. 3 is a graph of the angle profile of the present invention;

FIG. 4 is a plot of mean square error of the present invention;

fig. 5 is a graph of the average output signal to noise ratio of the present invention.

Detailed Description

Embodiments and effects of the present invention will be described in further detail below with reference to the accompanying drawings.

Referring to fig. 1, the implementation steps for this example are as follows:

step 1, echo data are obtained.

And performing down-conversion and matched filtering operation on echoes reflected by the target to obtain echo data NxK, L, wherein N represents the sum-difference channel number of an airborne/missile-borne radar receiving antenna, K represents the coherent pulse number of each pulse repetition period, and L represents the range gate sample number of the sampled data.

And 2, performing clutter suppression on the echo data.

(2.1) estimating clutter covariance matrix using maximum likelihood formula using echo vectors in echo data blocks

Wherein x is_lRepresenting an echo vector corresponding to the ith range gate in the sampled echo data block, and H represents conjugate transpose operation;

(2.2) setting a space domain guide vector:

the airspace guide vector is composed of a sum beam antenna directional diagram, an azimuth difference beam antenna directional diagram and a pitching difference beam antenna directional diagram, according to the characteristic that the response of the normalized sum beam antenna directional diagram in the main lobe direction is 1, the response of the normalized difference beam antenna directional diagram in the main lobe direction is 0 under an ideal condition, and the difference branch components in the airspace guide vector are 0, the sum channel airspace guide vector a under the sum and difference antenna system is set_sΣAzimuth difference channel airspace guide vectorPitching difference channel airspace guide vector

a_sΣ＝[1,0,0]^T，

The superscript T represents transposition operation, and the echo data form obtained by the phased array system radar after digital beam forming is the same as that of the sum and difference antenna system radar, so the value of the airspace steering vector can be obtained by the output of a sum beam antenna directional diagram, a azimuth difference beam antenna directional diagram and a pitch difference beam antenna directional diagram at the main lobe direction;

(2.3) calculating a space-time two-dimensional steering vector:

spatial domain steering vector a according to sum channel_sΣAzimuth difference channel airspace guide vectorPitching difference channel airspace guide vectorComputing sum channel space-time two-dimensional guide vector a_ΣSpace-time two-dimensional guide vector of azimuth difference channelPitching difference channel space-time two-dimensional guide vector

Wherein the content of the first and second substances,denotes the kronecker operation, a_tA time domain steering vector representing the main lobe direction of the radar receiving antenna; a is_tThe specific expression of (a) is as follows:

wherein j represents an imaginary unit symbol, f_dDenotes the doppler frequency, superscript T denotes the transpose operation;

(2.4) calculating an adaptive weight vector:

from clutter covariance matrixAnd channel space-time two-dimensional guide vector a_ΣSpace-time two-dimensional guide vector of azimuth difference channelPitching difference channel space-time two-dimensional guide vectorComputing and channel-optimized adaptive weight vector w_ΣAzimuth difference channel optimal adaptive weight vectorOptimal adaptive weight vector of pitch difference channel

Wherein mu is a scalar constant, and the superscript-1 represents the inversion operation;

(2.5) performing clutter suppression on the echo data to obtain output after clutter suppression;

optimal adaptive weight vector w according to sum channel_ΣAzimuth difference channel optimal adaptive weight vectorOptimal adaptive weight vector of pitch difference channelPerforming clutter suppression on the echo data to obtain output y of the sum channel after clutter suppression_ΣOutput of azimuth difference channelOutput of pitch difference channel

y_Σ＝[y_1Σ,y_2Σ,...,y_lΣ，...,y_LΣ]，l＝1,2,…,L

Wherein, y_lΣRepresenting the output of the ith range gate clutter suppressed sum channel in the sampled echo data block,representing the output of the l-th range gate clutter suppressed back head channel in the sampled echo data block,representing the output, y, of the ith range gate clutter suppressed back pitch difference channel in a block of sampled echo data_lΣ，The specific expression of (a) is as follows:

wherein x is_lAnd representing the echo vector corresponding to the ith range gate in the sampled echo data block, and the superscript H represents the conjugate transpose operation.

And 3, acquiring an azimuth angle and a pitch angle according to the output after clutter suppression.

(3.1) obtaining a range gate, a Doppler channel and an amplitude of the target in the sum channel according to the output after clutter suppression:

the results of simulation of the outputs after clutter suppression are shown in fig. 2, where:

FIG. 2(a) is a three-dimensional simulation diagram of the output of the sum channel after clutter suppression, wherein the X coordinate axis is a Doppler channel output by the sum channel, the Y coordinate axis is a range gate output by the sum channel, and the Z coordinate axis is an amplitude output by the sum channel;

FIG. 2(b) is a three-dimensional simulation diagram of the output of the clutter suppression rear potential difference channel, wherein the X coordinate axis is a Doppler channel output by the difference channel, the Y coordinate axis is a range gate output by the difference channel, and the Z coordinate axis is an amplitude output by the difference channel;

FIG. 2(c) is a diagram showing a three-dimensional simulation of the output of the pitch difference channel after clutter suppression, wherein the X-axis is a Doppler channel output by the pitch difference channel, the Y-axis is a range gate output by the pitch difference channel, and the Z-axis is an amplitude output by the pitch difference channel;

obtaining a target distance gate l from an output simulation diagram of the clutter suppressed sum channel_tDoppler channel m_tAnd an amplitude y_SΣThat is, from the coordinate values at the box in FIG. 2(a), the range gate l of the target is obtained_t339, Doppler channel m_t188 amplitude y_SΣIs 0.0005181.

(3.2) obtaining the amplitude y of the target in the azimuth difference channel_θAnd the amplitude of the target in the pitch difference channel

Searching output matrix of clutter suppression rear potential difference channel through MATLAB softwareAt point (l)_t,m_t) Obtaining the amplitude y of the target in the azimuth difference channel according to the corresponding value_θSearching an output matrix of the pitch difference channel after clutter suppression through MATLAB softwareAt point (l)_t,m_t) Corresponding value obtains the amplitude of the target in the pitch difference channel

(3.3) calculating the ratio q of the potential difference channel to the sum channel in the target prescription_θAnd the ratio of the pitch difference channel to the sum channel

According to the amplitude y of the target sum channel_SΣThe amplitude y of the target in the azimuth difference channel_θAnd the amplitude of the target in the pitch difference channelCalculating the ratio q of the channel to the sum channel of the target prescription potential difference_θAnd the ratio of the pitch difference channel to the sum channel

Wherein, y_SΣRepresenting the amplitude, y, of the sum channel at the target_θRepresenting the magnitude of the azimuth difference channel at the target, representingThe magnitude of the pitch difference channel at the target;

(3.4) calculating the slope of the angle curve:

the angle identification curve is shown in fig. 3, wherein 3(a) is an azimuth dimension angle identification curve two-dimensional graph, the X coordinate axis is an azimuth angle, the unit is degree, and the Y coordinate axis is the amplitude of the azimuth dimension angle identification curve; FIG. 3(b) is a two-dimensional view of a pitch angle profile, where the X-axis is the pitch angle in degrees and the Y-axis is the amplitude of the pitch angle profile;

respectively calculating the slope k of the azimuth angle-identifying curve according to the parameters of the antennas under the sum and difference antenna system_θAnd the slope of the pitch angle curve

Wherein j represents an imaginary unit symbol, L_cDenotes the length of the antenna, L_rIndicating the height, theta, of the antenna₀Indicating the azimuth angle at which the main lobe of the radar receiving antenna is directed,the pitch angle of the main lobe direction of the radar receiving antenna is represented, and lambda represents the wavelength;

(3.5) calculating the azimuth angle theta and the pitch angle of the target

Ratio q of channel to sum channel according to azimuth difference_θAnd the slope k of the azimuthal mapping curve_θCalculating the azimuth angle theta of the target:

according to the ratio of the elevation channel to the sum channelAnd the slope of the pitch angle curveCalculating the pitch angle of the target

Wherein, theta₀Azimuth angle, q, representing the main lobe orientation of the radar receiving antenna_θRepresenting the ratio of the potential difference channel to the sum channel, k, for the target_θRepresenting the slope of the azimuthal triangulation angle curve,represents the pitch angle of the main lobe of the radar receiving antenna,representing the ratio of the pitch difference channel to the sum channel at the target,representing the slope of the azimuthal dimension curve.

The effect of the present invention is further explained by combining the simulation experiment as follows:

1. simulation experiment conditions are as follows:

the hardware platform of the simulation experiment of the invention is as follows: the processor is Intel (R) core (TM) i7-10700 CPU, the main frequency is 2.90GHz, and the memory is 16 GB.

The software platform of the simulation experiment of the invention is as follows: the Windows 10 operating system and MATLAB R2020 b.

The parameters of the simulation experiment of the invention are set as follows: adopting 16 × 16 half-wavelength equidistant area array, radar working wavelength λ being 0.24m, array element spacing d being 0.12m, coherent pulse number K being 16 in each pulse repetition period, aircraft flying height H being 8Km, aircraft moving speed V being 120m/s, pulse repetition frequency f being 8 m/s_prf2000Hz, sample rate f_sThe noise to noise ratio CNR is 30dB at 5MHz, and the main beam width is 20 °.

2. Simulation content and result analysis thereof:

simulation 1, selecting multiple input snr points within the range of [ -20dB,10dB ] by using the above simulation conditions, performing 100 Monte-Carlo experiments on each snr point, averaging the corresponding angle error values to obtain the angle mean square error value corresponding to each snr point, and connecting the values to obtain the mean square error curve of the present invention, as shown in fig. 4. In fig. 4, the horizontal axis represents the input signal-to-noise ratio in dB, the vertical axis represents the mean square error in degrees, the curve marked by the diamond symbols represents the change of the mean square error of the pitch angle with the input signal-to-noise ratio, and the curve marked by the circle symbols represents the change of the mean square error of the azimuth angle with the input signal-to-noise ratio.

Simulation 2, selecting a plurality of input signal-to-noise ratio points within the range of [ -20dB,10dB ] by using the simulation conditions, and connecting the ratio of the signal output power corresponding to each signal-to-noise ratio point of the present invention to the clutter noise output power to obtain an output signal-to-noise ratio curve of the present invention, as shown in fig. 5. In fig. 5, the horizontal axis represents the input signal-to-noise ratio in dB, the vertical axis represents the output signal-to-noise ratio in dB, the curve marked with diamond symbols represents the curve of the output signal-to-noise ratio of the sum channel varying with the input signal-to-noise ratio, the curve marked with circle symbols represents the curve of the output signal-to-noise ratio of the azimuth difference channel varying with the input signal-to-noise ratio, and the curve marked with dot symbols represents the curve of the output signal-to-noise ratio of the pitch difference channel varying with the input signal-to-noise ratio.

As can be seen from the simulation result of FIG. 4, the input SNR is-13 dB, the output SNR corresponding to the channel and FIG. 5 is 8dB, and the angle measurement accuracy of the present invention reaches 1/10 beam width; the input signal-to-noise ratio is-8 dB, the output signal-to-noise ratio corresponding to the channel and the diagram of FIG. 5 is 12dB, and the angle measurement precision of the invention reaches 1/20 beam width.

The simulation result shows that the method can solve the problems of detection and angle measurement of the moving target under the clutter background, and verifies the correctness, effectiveness and reliability of the method.