Broadband radar hypersonic acceleration target rapid detection and parameter estimation method

Abstract

The invention discloses a broadband radar hypersonic acceleration target rapid detection and parameter estimation method, relates to the technical field of radar detection, and aims to solve the problem of low detection accuracy of uniform acceleration target detection in the prior art. According to the method, two matrixes O and E are separated, and then based on the result of the TRT, the speed characteristics of the uniform acceleration target can be accurately identified through the VPRT by using the information of the speed correlation matrix O and the targeted search of the speed parameters, so that the detection accuracy of the uniform acceleration target is enhanced.

Description

Technical Field

[0001] The present invention relates to the technical field of radar detection, and specifically to a method for rapid detection and parameter estimation of hypersonic acceleration targets by a broadband radar. Background Art

[0002] In the early stage, the range resolution of radars was relatively low. Therefore, within the coherent processing interval, it could be assumed that the target was located in a single range resolution cell, and the echoes of its different pulses only had phase differences without envelope changes, that is, there was no so-called "range walk". Based on this assumption, the earliest moving target detection (MTD) method was proposed. However, in practical applications, range walk of the target might also occur.

[0003] To solve the problem of joint processing of multiple pulses in the case of target range walk, the Hough transform (HT) was initially proposed. This is a non-coherent accumulation method that only utilizes the amplitude information of the target echo and does not consider the phase information. Subsequently, the Keystone transform (KT) was introduced as a means of coherent accumulation, which can utilize both amplitude and phase information and effectively improves the detection performance of the radar. However, KT requires complex interpolation operations and faces the problem of Doppler ambiguity. For this reason, the Radon-Fourier transform (RFT) based on the Radon transform and the Fourier transform was proposed. By jointly searching for the range and velocity parameters of the target and performing coherent accumulation along the corresponding trajectory, RFT can effectively handle uniformly moving targets. Nevertheless, when the target acceleration is not zero, the performance of KT and RFT will be affected to a certain extent. Therefore, in the prior art, there is a problem of low detection accuracy for the detection of uniformly accelerating targets. Summary of the Invention

[0004] The object of the present invention is to provide a method for rapid detection and parameter estimation of hypersonic acceleration targets by a broadband radar, aiming at the problem of low detection accuracy in the prior art for the detection of uniformly accelerating targets.

[0005] The technical solution adopted by the present invention to solve the above technical problems is as follows:

[0006] A method for rapid detection and parameter estimation of hypersonic acceleration targets by a broadband radar, comprising the following steps:

[0007] Step 1: Obtain a plurality of radar pulse echo sampling data, and use the radar pulse echo sampling data to construct a time-domain echo matrix Z. The time-domain echo matrix Z is a matrix of MxN, where M represents the number of pulses and N represents the number of sampling points for each pulse;

[0008] Step 2: Obtain the target distance range, speed range, and acceleration range, and equally divide the target distance range, speed range, and acceleration value range. Then sort the distances, speeds, and accelerations obtained from the equal division in ascending order respectively;

[0009] Step 3: Perform a fast Fourier transform on each row of the time-domain echo matrix Z to convert the time-domain echo matrix Z into a frequency-domain echo matrix Z f ;

[0010] Step 4: Perform a time reversal transformation on the frequency-domain echo matrix Z f to obtain a reference matrix

[0011] Step 5: Construct a matrix f and a matrix according to the frequency-domain echo matrix Z and the reference matrix

[0012] where, represents the conjugate of the matrix , ⊙ represents element-wise multiplication of matrices, O represents a matrix that only depends on the speed parameter, and E represents a matrix that only depends on the distance and acceleration parameters;

[0013] Step 6: Perform a fast Fourier transform on each column of the matrix O to obtain the matrix O1;

[0014] Step 7: Calculate the VPRT output for each speed in the sorting according to the matrix O1;

[0015] Step 8: According to the result of Step 7, select the maximum value of the VPRT output, and determine whether this maximum value exceeds the threshold. If it exceeds, there is a target, and use this maximum value as the speed estimation value of the target Continue to Step 9, otherwise there is no target and end;

[0016] Step 9: Use the speed estimation value of the target to construct a speed matching filter matrix

[0017] Step 10: Use the matrix E and the matrix to construct the matrix E1, and the matrix E1 is expressed as

[0018] Step 11: Obtain the DMT output for each distance in the distance sorting obtained in Step 2 according to the matrix E1, and the distance corresponding to the maximum value of the DMT output is the distance estimation value of the target

[0019] Among the acceleration rankings obtained in step 2 according to matrix E1, the DMT output of each acceleration is obtained, and the acceleration corresponding to the maximum DMT output is the estimated value of the target acceleration.

[0020] Furthermore, the time-domain echo matrix Z is expressed as:

[0021] Z = αS θ + N

[0022] where α represents the target echo amplitude, S θ represents the target echo matrix, N represents the received noise matrix, θ represents the motion parameters of the target, θ = [r, v, a], r represents the initial distance of the target, v represents the velocity of the target, and a represents the acceleration of the target.

[0023] Furthermore, the frequency-domain echo matrix Z f is expressed as:

[0024] Z f = FFT r (Z)

[0025] where FFT r (·) represents performing a fast Fourier transform on each row of the matrix.

[0026] Furthermore, the reference matrix is expressed as:

[0027]

[0028] where flip represents flipping the matrix up and down.

[0029] Furthermore, the matrix O1 is expressed as:

[0030] O1 = FFT c (O)

[0031] where FFT c (·) represents performing a fast Fourier transform on each column of the matrix.

[0032] Furthermore, the VPRT output is expressed as:

[0033]

[0034] ρ xΔv = 1 - λ xΔv

[0035] where Δf represents the fast-time frequency interval, Δf = f s / N, f srepresents the sampling frequency, N represents the number of sampling points per pulse, Δv represents the velocity search interval, x represents the x-th velocity value, u xΔv , l xΔv , λ xΔv and ρ xΔv represent intermediate variables, f c represents the radar carrier frequency, B represents the radar bandwidth, c represents the speed of light, [·] represents the rounding operation, |·| represents the modulus of a complex number, o n,x The subscript x of o represents the x-th velocity value, n represents the n-th fast-time frequency, corresponding to the n-th column of matrix O1. Different columns of matrix O1 correspond to different fast-time frequencies. The fast-time frequency corresponding to the n-th column is represented by f n From the 1st column to the Nth column, the corresponding fast-time frequencies are:

[0036]

[0037] Different rows of matrix O1 correspond to different slow-time frequencies. The fast-time frequency corresponding to the m-th row is represented by f m From the 1st row to the Mth row, the corresponding slow-time frequencies are:

[0038]

[0039] where M represents the number of pulses, f r represents the pulse repetition frequency;

[0040] Calculate the following linear equation according to the x-th velocity value:

[0041]

[0042] where, respectively represent the slow-time frequency and the fast-time frequency on the trajectory;

[0043] o n,x The value-taking rule is as follows:

[0044] Let f n be the fast-time frequency corresponding to the n-th column of matrix O1, and calculate If the slow-time frequency closest to is the m-th slow-time frequency f m , then o n,x is equal to the element in the m-th row and n-th column of matrix O1, that is, o n,x = O1(m,n), and O1(m,n) represents the element in the m-th row and n-th column of matrix O1.

[0045] Further, the velocity estimate value of the target is expressed as:

[0046]

[0047] Among them, represents the estimated value of x.

[0048] Furthermore, the velocity matching filter matrix is expressed as:

[0049]

[0050] Among them, represents the element in the m-th row and n-th column of matrix , rect(·) is the rectangular window function, exp(·) is the exponential function, t m represents the slow time corresponding to the m-th pulse, T c represents the coherent processing time of the radar, represents the imaginary unit, π represents the pi, K represents the radar chirp rate, and represents the intermediate variable.

[0051] Furthermore, the DMT output is obtained through the following steps:

[0052] Acceleration DMT output:

[0053] Step 1.1: In the acceleration sorting, select the first acceleration;

[0054] Step 1.2: According to the selected acceleration, calculate the matching acceleration output of each column of matrix E1, and record the result as DMT(y,n), where y represents the y-th acceleration and n represents the n-th column of matrix E1;

[0055] Step 1.3: Repeat Step 1.1 and Step 1.2 until the DMT(y,n) corresponding to all accelerations in the acceleration sorting is obtained;

[0056] DMT(y,n) is expressed as:

[0057]

[0058] Among them, e n ∈C M×1 represents the n-th column of matrix E1, M represents the number of pulses, Δa represents the acceleration search interval, u y,n represents the intermediate variable;

[0059] Range DMT output:

[0060] Step 2.1: In the range sorting, select the first range;

[0061] Step 2.2: Calculate the matching distance output for each column of matrix E1 according to the selected distance, and denote the result as DMT(y,z), where y represents the y-th distance and z represents the z-th distance parameter value;

[0062] Step 2.3: Repeat Step 2.1 and Step 2.2 until DMT(y,z) corresponding to all distances in the distance sorting is obtained;

[0063] DMT(y,z) is expressed as:

[0064] DMT(y,z) = d y v y,z

[0065] where, d y represents the y-th row of matrix DMT(y,n), and v y,z (n) = exp(j2πnΔf(zΔr)), Δr represents the distance search interval, and v y,z represents an intermediate variable.

[0066] Furthermore, the distance estimation value of the target and the acceleration estimation value of the target are expressed as:

[0067]

[0068]

[0069] The beneficial effects of the present invention are:

[0070] In this application, through the time reversal operation of the TRT on the echo signal, different motion parameters are effectively separated, manifested as separating two matrices O and E. Then, based on the results of the TRT, through VPRT, using the information of the velocity-related matrix O, the targeted search for velocity parameters can accurately identify the velocity characteristics of uniformly accelerating targets, thereby enhancing the accuracy of uniformly accelerating target detection. BRIEF DESCRIPTION OF THE DRAWINGS

[0071] Figure 1 is the overall flowchart of this application;

[0072] Figure 2 is a schematic diagram of matrix O1;

[0073] Figure 3 is a schematic diagram of the VPRT output;

[0074] Figure 4 is a schematic diagram of the DMT output;

[0075] Figure 5 is a schematic diagram of the DMT output along the distance profile;

[0076] Figure 6 It is a schematic diagram of the DMT output along the acceleration profile. Specific implementation manners

[0077] It should be particularly noted that, without conflict, the various implementation manners disclosed in this application can be combined with each other.

[0078] Specific implementation manner one: The method for fast detection and parameter estimation of broadband radar hypersonic acceleration targets described in this implementation manner includes the following steps:

[0079] Step 1: Obtain multiple radar pulse echo sampling data, and construct a time-domain echo matrix Z using the radar pulse echo sampling data. The time-domain echo matrix Z is a matrix of MxN, where M represents the number of pulses and N represents the number of sampling points for each pulse;

[0080] Step 2: Obtain the target distance range, speed range, and acceleration range, equally spaced divide the target distance range, speed range, and acceleration value range, and sort the equally spaced divided distances, speeds, and accelerations from smallest to largest;

[0081] Step 3: Perform a fast Fourier transform on each row of the time-domain echo matrix Z to convert the time-domain echo matrix Z into a frequency-domain echo matrix Z f ;

[0082] Step 4: Perform a time reversal transformation on the frequency-domain echo matrix Z f to obtain a reference matrix

[0083] Step 5: Construct matrix O and matrix E according to the frequency-domain echo matrix Z f and the reference matrix Matrix O is expressed as:

[0084]

[0085] Matrix E is expressed as:

[0086]

[0087] where represents the conjugate of matrix ⊙ represents element-by-element multiplication of matrices, O represents a matrix that only depends on the speed parameter, and E represents a matrix that only depends on the distance and acceleration parameters;

[0088] Step 6: Perform a fast Fourier transform on each column of matrix O to obtain matrix O1; as Figure 2 shown.

[0089] Step 7: Calculate the VPRT output for each speed in the sorting according to matrix O1; as Figure 3 shown.

[0090] Step 8: According to the result of Step 7, select the maximum value of the VPRT output, and determine whether the maximum value exceeds the threshold. If it exceeds, there is a target, and take this maximum value as the speed estimate value of the target Continue to Step 9, otherwise there is no target, end;

[0091] Step 9: Use the speed estimate value of the target to construct a speed matching filter matrix

[0092] Step 10: Use matrix E and matrix to construct matrix E1, and matrix E1 is expressed as

[0093] Step 11: Obtain the DMT output for each distance in the distance sorting obtained in Step 2 according to matrix E1, where the maximum value of the DMT output is the distance estimate value of the target as Figure 4 , Figure 5 and Figure 6 shown.

[0094] Step 12: Obtain the DMT output for each acceleration in the sorting according to matrix E1, where the maximum value of the DMT output is the acceleration estimate value of the target

[0095] Step 1: Assume that the target is moving with a constant acceleration, and obtain the time-domain sampling data of multiple radar pulse echoes. Each pulse echo sampling data is arranged as a row vector, and the echo data of different pulses are arranged as an echo matrix Z. The first row of the matrix represents the first pulse echo sampling data, and so on. Assume that the discrete sampling data received by the radar is:

[0096] Z = αS θ + N

[0097] where α represents the target echo amplitude, S θ represents the target echo matrix, N represents the received noise matrix, θ = [r, v, a] represents the motion parameters of the target, r represents the initial distance of the target, v represents the speed of the target, and a represents the acceleration of the target. Make an equally spaced division of the possible distance, speed, and acceleration value parameters of the target to obtain a series of different distance, speed, and acceleration values, and arrange the distance, speed, and acceleration values from small to large. The first distance, speed, and acceleration values correspond to the minimum distance, speed, and acceleration values.

[0098] Step 2: Perform a fast Fourier transform on each row of the echo matrix to convert the time-domain echo matrix into a frequency-domain echo matrix Z f , which is expressed by the formula:

[0099] Z f = FFT r (Z)

[0100] where FFT r (·) represents performing a fast Fourier transform on each row of the matrix.

[0101] Step 3: Perform a time reversal transform (TRT) on matrix Z f , that is, construct a reference matrix according to the frequency-domain echo matrix which is obtained by flipping the frequency-domain echo matrix Z f vertically, that is the first row of corresponds to the last pulse, and so on. It is expressed by the formula:

[0102]

[0103] where flip represents vertically flipping the matrix.

[0104] Step 4: Calculate matrices O and E according to the frequency-domain echo matrix Z f and the reference matrix , which is expressed by the formula

[0105]

[0106] where represents the conjugate of matrix , ⊙ represents element-wise multiplication of matrices, O represents a matrix that only depends on the velocity parameter, and E represents a matrix that only depends on the distance and acceleration parameters. Through this step, the target parameters can be separated, thus realizing the separate estimation of parameters and avoiding coupling.

[0107] Step 5: Perform a fast Fourier transform on each column of matrix O, and the result is denoted as matrix O1. It is expressed by the formula

[0108] O1 = FFT c (O)

[0109] where FFT c (·) represents performing a fast Fourier transform on each column of the matrix.

[0110] Step 6: Calculate the VPRT output of each velocity value according to matrix O1, which is specifically divided into the following sub-steps:

[0111] (1) Select the first velocity parameter value;

[0112] (2) Based on the selected speed parameter, calculate the VPRT output of the speed value;

[0113] (3) Select the next speed parameter value in turn and repeat step (2) to obtain the VPRT output of all parameter values;

[0114] The VPRT output for the xth velocity value is expressed as:

[0115]

[0116] Where Δf = f s / N represents the fast time frequency interval, f s represents the sampling frequency, N represents the number of sampling points for each pulse. Δv represents the speed search interval, and the search speed v s It can be expressed as v s =xΔv. f c represents the radar carrier frequency, ρ xΔv =1-λ xΔv , B represents the radar bandwidth, c represents the speed of light, o n,x Represents the corresponding search speed discrete trajectory The data on m represents the fast time frequency corresponding to the mth row, the symbol [·] represents the rounding operation, and |·| represents the modulus value of the complex number.

[0117] Step 7: Based on the result of step 6, calculate the maximum value of the VPRT output, compare the maximum value with the threshold, if it exceeds the threshold, the corresponding parameter point is judged as the target, and record the parameter value at the same time, and use the parameter value as the target speed estimation value Continue to step 8; otherwise, it is determined that there is no target and the algorithm ends. The estimated result of the target speed is expressed as follows:

[0118]

[0119] in Indicates the estimated value of speed. Whether the target exists can be determined by the following formula:

[0120]

[0121] Where T represents the detection threshold.

[0122] Step 8: Estimated value of speed parameters according to step 7 Constructing the Matrix The formula is as follows

[0123]

[0124] in Represents a matrix The element in the m-th row and n-th column, rect(·) is the rectangular window function, exp(·) is the exponential function, t m Represents the slow time corresponding to the m-th pulse, T c Represents the radar coherent processing time Represents the imaginary unit, π represents the pi, K represents the radar frequency modulation rate

[0125] Step Nine: According to matrix E in Step Four and matrix Calculate matrix E1, which is expressed by the formula:

[0126]

[0127] Step Ten: Calculate the DMT output of each distance and acceleration value according to matrix E1, which is specifically divided into the following sub-steps:

[0128] Acceleration matching step:

[0129] (1-1) Select the first acceleration parameter value;

[0130] (1-2) According to the selected velocity parameter value, calculate the matching acceleration output of different columns of matrix E1, and record the result as DMT(y,n), where y represents the y-th acceleration parameter value, here y = 1, n represents the n-th column of matrix E1;

[0131] (1-3) Select the next acceleration parameter value in turn, and repeat step (1-2) until all acceleration parameter values are calculated.

[0132] Among them, DMT(y,n) is expressed by the formula as follows:

[0133]

[0134] Among them, e n ∈C M×1 Represents the n-th column of matrix E1, M represents the number of pulses, Δa is the acceleration search interval, and the search acceleration can be expressed as a s =yΔa.

[0135] Distance matching step:

[0136] (2-1) Select the first distance parameter value;

[0137] (2-2) Calculate the matching distance outputs of different rows of DMT(y,n) according to the selected distance parameter value, and record the result as DMT(y,z), where y represents the y-th row of DMT(y,n), corresponding to y acceleration parameter values, and z represents the z-th distance parameter value, where z = 1;

[0138] (2-3) Select the next distance parameter value in sequence, and repeat step (2-2) until all distance parameter values are calculated.

[0139] Among them, DMT(y,z) is expressed by the formula as follows:

[0140] DMT(y,z) = d y v y,z

[0141] Among them, d y represents the y-th row of the matrix DMT(y,n), and v y,z (n) = exp(j2πnΔf(zΔr)), Δr represents the distance search interval, and the search distance can be expressed as r s = zΔr.

[0142] Step Eleven: According to the result of Step Ten, calculate the maximum value of the output of DMT(y,z), and estimate the distance and acceleration of the target using the parameter value corresponding to the maximum value. It is expressed by the formula as

[0143]

[0144] Among them represents the acceleration estimation value, represents the distance estimation value.

[0145] Among them, Δf represents the fast-time frequency interval, Δf = f s / N, f s represents the sampling frequency, N represents the number of sampling points per pulse, Δv represents the velocity search interval, x represents the x-th velocity value, u xΔv , l xΔv , λ xΔv and ρ xΔv are intermediate variables, f c represents the radar carrier frequency, B represents the radar bandwidth, c represents the speed of light, [·] represents the rounding operation, and |·| represents the modulus value of a complex number.

[0146] o n,x The subscript x of represents the x-th velocity value, and n represents the n-th fast-time frequency, corresponding to the n-th column of the matrix O1. Different columns of the matrix O1 correspond to different fast-time frequencies, and the fast-time frequency corresponding to the n-th column is represented by f n , and the fast-time frequencies corresponding to the columns from the 1st column to the Nth column are respectively:

[0147]

[0148] Different rows of matrix O1 correspond to different slow-time frequencies, and the fast-time frequency corresponding to the nth column is represented by f m The slow-time frequencies corresponding to the first row to the Mth row are respectively:

[0149]

[0150] where M represents the number of pulses, and f r represents the pulse repetition frequency.

[0151] Calculate the following linear equation according to the xth speed value:

[0152]

[0153] where respectively represent the slow-time frequency and the fast-time frequency on the trajectory. o n,x The value rule of o is as follows: Let f n be the fast-time frequency corresponding to the nth column of matrix O1, calculate Judge the slow-time frequency closest to Suppose it is the mth slow-time frequency f m , then o n,x is equal to the element in the mth row and nth column of matrix O1, that is, o n,x = O1(m, n), and O1(m, n) represents the element in the mth row and nth column of matrix O1.

[0154] Use DMT (Distance and Acceleration Double Matching Transform) to estimate the target distance and acceleration: After completing the speed estimation, this step focuses on the accurate estimation of distance and acceleration. DMT is divided into two sub-steps:

[0155] (a) Match Acceleration (MA): First, perform matching processing on the acceleration of the target to ensure the accuracy of the acceleration parameters.

[0156] (b) Match Distance (MR): Subsequently, perform matching processing on the distance parameters of the target to achieve high-precision distance estimation.

[0157] This step-by-step processing method not only ensures the independence and accuracy of each parameter estimation, effectively solves the problems of high computational complexity and difficulty in real-time processing in traditional methods, but also improves the computational efficiency by accurately processing different types of parameters in stages.

[0158] It should be noted that the specific implementation manners are only explanations and illustrations of the technical solutions of the present invention, and the scope of the claimed protection cannot be limited thereby. Any changes that are merely partial based on the claims and the description of the present invention should still fall within the protection scope of the present invention.

Claims

1. A broadband radar hypersonic acceleration target rapid detection and parameter estimation method, characterized in that The following steps are involved: Step 1: Acquire multiple radar pulse echo sampling data, and use the radar pulse echo sampling data to construct a time domain echo matrix Z, wherein the time domain echo matrix Z is an MxN matrix, where M represents the number of pulses and N represents the number of sampling points for each pulse; Step 2: Obtain the target distance range, speed range, and acceleration range, and divide the target distance range, speed range, and acceleration value range into equal intervals, and sort the distances, speeds, and accelerations obtained by the equal interval divisions from small to large; Step 3: Perform a fast Fourier transform on each row of the time domain echo matrix Z to convert the time domain echo matrix Z into the frequency domain echo matrix Z f ; Step 4: Frequency domain echo matrix Z f Perform a time flip transformation to obtain the reference matrix Step 5: According to the frequency domain echo matrix Z f and the reference matrix Constructing the Matrix and matrix in, Representation Matrix The conjugate of , ⊙ represents the element-by-element multiplication of matrices, O represents the matrix that depends only on the velocity parameter, and E represents the matrix that depends only on the distance and acceleration parameters; Step 6: Perform fast Fourier transform on each column of matrix O to obtain matrix O1; Step 7: Calculate the VPRT output of each speed in the sorting according to the matrix O1; Step 8: According to the result of step 7, select the maximum value of VPRT output and determine whether the maximum value exceeds the threshold. If it exceeds the threshold, there is a target and the maximum value is used as the target speed estimate. Continue to step 9, otherwise there is no target, end; Step 9: Using the target's velocity estimate Construct velocity matched filter matrix Step 10: Using the matrix E and the matrix Construct matrix E1, which is expressed as Step 11: Obtain the DMT output of each distance in the distance sorting obtained in step 2 according to matrix E1, where the distance corresponding to the maximum DMT output is the estimated distance value of the target. According to matrix E1, the DMT output of each acceleration in the acceleration sorting obtained in step 2 is obtained, where the acceleration corresponding to the maximum DMT output is the estimated acceleration value of the target.

2. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 1 is characterized in that The time domain echo matrix Z is expressed as: Z=αS θ +N Among them, α represents the target echo amplitude, S θ represents the target echo matrix, N represents the received noise matrix, θ represents the motion parameters of the target, θ=[r,v,a], r represents the initial distance of the target, v represents the speed of the target, and a represents the acceleration of the target.

3. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 2 is characterized in that The frequency domain echo matrix Z f It is expressed as: Z f =FFT r (Z) Among them, FFT r (·) represents fast Fourier transform of each row of the matrix.

4. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 3 is characterized in that The reference matrix It is expressed as: Among them, flip means flipping the matrix upside down.

5. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 4 is characterized in that The matrix O1 is expressed as: O1=FFT c (O) Among them, FFT c (·) represents fast Fourier transform of each column of the matrix.

6. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 5 is characterized in that The VPRT output is expressed as: r xΔv =1-λ xΔv Where Δf represents the fast time frequency interval, Δf = f s / N,f s represents the sampling frequency, N represents the number of sampling points for each pulse, Δv represents the speed search interval, x represents the xth speed value, u xΔv , l xΔv , xΔv and ρ xΔv represents the intermediate variable, f c represents the radar carrier frequency, B represents the radar bandwidth, c represents the speed of light, [·] represents the rounding operation, |·| represents the modulus value of the complex number, and o n,x The subscript x represents the xth velocity value, n represents the nth fast time frequency, corresponding to the nth column of the matrix O1. Different columns of the matrix O1 correspond to different fast time frequencies. The fast time frequency corresponding to the nth column is represented by f n Indicates that the fast time frequencies corresponding to the 1st column to the Nth column are: Different rows of the matrix O1 correspond to different slow time frequencies, and the fast time frequency corresponding to the mth row is represented by f m Indicates that the slow time frequencies corresponding to the 1st row to the Mth row are: Where M represents the number of pulses, f r represents the pulse repetition frequency; Calculate the following straight line equation based on the x-th speed value: in, represent the slow time frequency and fast time frequency on the trajectory respectively; o n,x The value rules are as follows: make f n is the fast time frequency corresponding to the nth column of matrix O1, and calculates If The closest slow time frequency is the mth slow time frequency f m , then n,x is equal to the element in the mth row and nth column of matrix O1, that is, o n,x =O1(m,n), O1(m,n) represents the element in the mth row and nth column of the matrix O1.

7. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 6 is characterized in that The target's estimated velocity It is expressed as: in, represents the estimated value of x.

8. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 7 is characterized in that The velocity matching filter matrix It is expressed as: in, Representation Matrix The element in the mth row and nth column, rect(·) is the rectangular window function, exp(·) is the exponential function, t m represents the slow time corresponding to the mth pulse, T c represents the radar coherent processing time, represents the imaginary unit, π represents the ratio of circumference to circumference, K represents the radar modulation frequency, and Represents an intermediate variable.

9. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 8 is characterized in that The DMT output is obtained by the following steps: Acceleration DMT output: Step 1.1: In the acceleration sorting, select the first acceleration; Step 1.2: According to the selected acceleration, calculate the matching acceleration output of each column of the matrix E1, and record the result as DMT(y,n), where y represents the yth acceleration and n represents the nth column of the matrix E1; Step 1.3: Repeat steps 1.1 and 1.2 until the DMT (y, n) corresponding to all accelerations in the acceleration sorting is obtained; DMT(y,n) is expressed as: in, represents the nth column of matrix E1, M represents the number of pulses, Δa represents the acceleration search interval, u y,n represents an intermediate variable; Distance DMT output: Step 2.1: In the distance sorting, select the first distance; Step 2.2: According to the selected distance, calculate the matching distance output of each column of the matrix E1, and record the result as DMT(y,z), where y represents the yth distance and z represents the zth distance parameter value; Step 2.3: Repeat steps 2.1 and 2.2 until the DMT (y, z) corresponding to all distances in the distance sorting is obtained; DMT(y,z) is expressed as: DMT(y,z)=d y v y,z Among them, d y represents the y-th row of the matrix DMT(y,n), v y,z (n) = exp(j2πnΔf(zΔr)), Δr represents the distance search interval, v y,z Represents an intermediate variable.

10. The broadband radar hypersonic acceleration target rapid detection and parameter estimation method according to claim 9 is characterized in that The estimated distance to the target and the target's acceleration estimate It is expressed as:

Images