Simultaneous measurement method of wideband signal multi-domain parameters based on space-time joint processing

Abstract

The application discloses a kind of wideband signal multi-domain parameter simultaneous measurement method based on space-time joint processing, belong to electronic reconnaissance field, first to radar signal AD sampling, after polyphase filtering, filter data are obtained;Filter data are respectively carried out time domain and space domain FFT transformation, and the spectrum sequence of space domain focusing is obtained;Direction correction is carried out to the problem of the spectrum sequence of space domain focusing exists direction deviation, and space-time two-dimensional spectrum matrix is obtained;Space-time two-dimensional spectrum matrix is projected in space domain, and through threshold setting and peak value screening, the angle of arrival measurement value is obtained;Center frequency and bandwidth measurement value are obtained by threshold design and peak value screening using angle of arrival measurement value to extract specific angle of arrival sequence.

Description

TECHNICAL FIELD

[0001] The present application belongs to the field of electronic reconnaissance, and particularly relates to a wideband signal multi-domain parameter simultaneous measurement method based on space-time joint processing. BACKGROUND

[0002] In order to realize effective jamming in a complex combat environment, a wideband receiver of a ground electromagnetic defense system needs to have a rapid sensing capability for target parameter characteristics, and can dynamically and accurately measure multi-domain parameter values of a target. Therefore, the wideband receiver of the ground electromagnetic defense system needs to have a strong signal analysis capability.

[0003] In a conventional parameter characteristic sensing method, parameters of each domain are usually completed independently, and multi-dimensional power focusing of a target is not realized, which seriously affects the characteristic sensing capability of a weak wideband signal. Gong Bin proposed a wideband array signal detection algorithm based on power focusing in the paper of Wideband Array Signal Detection Algorithm Based on Power Focusing. The method uses phase change information of a uniform linear array received signal, focuses power of each array element received signal to a signal arrival angle through an improved double-FFT algorithm, constructs a constant false alarm detection threshold, and realizes effective detection of a weak target. The method only realizes a detection process of a wideband array signal, and does not consider parameter measurement of the signal, which is insufficient to provide dynamic and reliable information support for emitter identification and positioning tracking. Therefore, it is important to study a space-time two-dimensional spectrum graph after power focusing, select a suitable image processing method, realize a multi-domain parameter simultaneous measurement process of a wideband signal based on power focusing, and improve the environmental sensing capability and countermeasure capability of a ground electromagnetic defense system. SUMMARY

[0004] The present application proposes a wideband signal multi-domain parameter simultaneous measurement method based on space-time joint processing, realizes space domain power focusing of signals with different incident directions, improves an environmental signal-to-noise ratio of parameter measurement, and improves the environmental sensing capability of a defense system. The method realizes space domain power focusing of a full reconnaissance frequency band and a full reconnaissance angle, successfully constructs a space-time two-dimensional spectrum amplitude graph of a combat scene, uses linear integral operation of different mapping angles for straight line detection, and realizes simultaneous and accurate measurement of multi-class targets and multi-domain parameters.

[0005] The technical solution of the present application is as follows: a wideband signal multi-domain parameter simultaneous measurement method based on space-time joint processing, comprising the following steps:

[0006] Step 1: after AD sampling and multi-phase filtering of a radar signal, filtered data are obtained, and the step 2 is entered.

[0007] Step 2: time domain FFT calculation is performed on the filtered data, a wideband decomposition spectrum matrix is obtained, and the step 3 is entered.

[0008] Step 3, spatial domain FFT calculation is performed on each row of the broadband decomposed frequency spectrum matrix to obtain a spatially focused frequency spectrum sequence, and step 4 is entered.

[0009] Step 4, the position of each element of the spatially focused frequency spectrum sequence is shifted and the element is intercepted by using the periodicity of the discrete Fourier transform in the frequency domain to obtain a space-time two-dimensional frequency spectrum matrix, and step 5 is entered.

[0010] Step 5, summing each column element of the space-time two-dimensional frequency spectrum matrix to obtain an angle of arrival projection sequence, setting an angle of arrival detection threshold, finding a peak value element greater than the detection threshold, and calculating an angle of arrival measurement value sequence according to the position corresponding to the peak value element of the angle of arrival projection sequence, and step 6 is entered.

[0011] Step 6, according to the angle of arrival peak position sequence, combining the elements of the corresponding column of the space-time two-dimensional frequency spectrum matrix in step 4, calculating a center frequency measurement value sequence and a bandwidth measurement value sequence.

[0012] Compared with the prior art, the present application has the following advantages:

[0013] 1) The present application realizes multi-dimensional power focusing, constructs a space-time two-dimensional frequency spectrum graph, and improves weak target sensing performance. 2) It realizes straight line detection of any mapping angle of the frequency spectrum graph, realizes multi-dimensional analysis ability of the signal; 3) It realizes simultaneous measurement of multi-domain parameters of multiple broadband signals, and improves the cognitive ability in complex combat environment. BRIEF DESCRIPTION OF DRAWINGS

[0014] Figure 1 The flow chart of the method for simultaneously measuring multi-domain parameters of a broadband signal based on space-time joint processing.

[0015] Figure 2 It is a broadband decomposition schematic diagram.

[0016] Figure 3 It is a space-time two-dimensional frequency spectrum matrix schematic diagram.

[0017] Figure 4 It is an angle of arrival projection sequence schematic diagram. DETAILED DESCRIPTION

[0018] The technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0019] In addition, the technical solutions among various embodiments of the present application can be combined with each other, but the combination of the technical solutions should be considered not to exist and not within the protection scope required by the present application on the basis that the combination of the technical solutions can be realized by the ordinary skilled in the art, and the combination of the technical solutions cannot be realized when the combination of the technical solutions appears contradictory.

[0020] The technical difficulties and points of the present application will be further introduced in the specific embodiments combined with the design examples.

[0021] Combined Figure 1 A wideband signal multi-domain parameter simultaneous measurement method based on space-time joint processing, the steps are as follows:

[0022] Step 1, after AD sampling and polyphase filtering of the radar signal, filtered data is obtained. Specifically, the following steps are included:

[0023] Step 11, AD sampling is performed on the radar signal received by each sub-array element in the antenna array to obtain a sampling signal sequence x k (n), wherein the sub-array element serial number k=1, 2, …, K, K represents the number of sub-array elements, the input time variable before filtering n=1, 2, …, N, N represents the sampling length. The polyphase filtering is performed on the sampling signal sequence x k (n) to obtain sub-array element filtered data y k (nM) as follows:

[0024]

[0025] Wherein, the decimation rate M=32, then the output sampling frequency after filtering is 500 / 32=15.625MHz. The channel number Q=32, ω q (q=0, 1, …, Q-1) is the center frequency of each channel, h(n) is a low-pass FIR filter, and j represents an imaginary part.

[0026] Step 12, the K sub-array element filtered data is combined into a matrix to obtain filtered data D=[y k (nM) N×K .

[0027] Step 2, time domain FFT calculation is performed on the filtered data to obtain a wideband decomposition frequency spectrum matrix (as shown in Figure 2 ). Specifically, the following steps are included:

[0028] Step 21, independent wideband decomposition is performed on the sub-array element filtered data y k (nM), that is, FFT operation with a length of N is performed on each column in the filtered data D to obtain a wideband decomposition frequency spectrum matrix XF, and the specific form is as follows:

[0029]

[0030] wherein XF(i, k) represents the spectral element of the i-th frequency bin and the k-th element in XF, the columns of the wideband decomposed spectral matrix XF represent the spectral sequence of the radar signal received by each element, i (i = 1, 2, …, I) represents the frequency bin number, I is the total number of frequency bins, and satisfies I = N, if the frequency range of the filtered data D is [F min ,F max ], then the width of each frequency bin is ΔF = (F max -F min ) / N.

[0031] Step 22, if the received radar signal is at an angle of γ, then the phase difference caused by the time delay between the adjacent column data of the i-th frequency bin in the wideband decomposed spectral matrix XF can be represented as

[0032]

[0033] wherein λ i represents the wavelength value corresponding to the i-th frequency bin, and d represents the distance between adjacent elements.

[0034] Step 23, according to the relationship between the spectral sequences of the radar signals between adjacent elements, the wideband decomposed spectral matrix XF can be represented as

[0035]

[0036] In the wideband decomposed spectral matrix of the above formula, the column data reflects the spectral characteristics of the received signal sequence of the corresponding sub-array element in the corresponding frequency bin, wherein X(i) represents the spectral amplitude of the signal in the i-th frequency bin.

[0037] Step 3, performing spatial FFT calculation on each row of the wideband decomposed spectral matrix to obtain the spatially focused spectral sequence. This step realizes spatial power focusing of signals with different incident directions, improves the signal-to-noise ratio of the parameter measurement environment, and improves the environmental perception ability of the defense system. Specifically, the following steps are included:

[0038] Step 31, taking the frequency bin corresponding to the wavelength λ ref = 2d as the reference frequency bin, i.e. i = i ref , and setting the FFT operation length of the reference frequency bin i ref L ref , then the FFT operation length L i of the i-th frequency bin can be represented as

[0039]

[0040] wherein round(*) represents the integer closest to the number *, and satisfies λ​i / 2d≥1, so the FFT operation length satisfies L i ≥L ref .

[0041] Step 32, spatial domain FFT calculation is performed on each row of the broadband decomposed spectrum matrix XF, and the FFT operation length is L i , to obtain the i-th frequency band spatially focused spectrum sequence BX i , and the sequence length is L i , and the specific form is as follows:

[0042] BX i =BX i (1), BX i (2), …, BX i (L i -1), BX i (L i ).

[0043] Step 4, by using the periodicity of the discrete Fourier transform in the frequency domain, the position of each element of the spatially focused spectrum sequence is shifted and the element is cut off, to obtain a space-time two-dimensional spectrum matrix (as shown in Figure 3 ). This process realizes spatial power focusing of the full detection frequency band and full detection angle, and successfully constructs a space-time two-dimensional spectrum amplitude graph of the combat scene. Specifically, the following steps are included:

[0044] Step 41, initialize the frequency band sequence number i = 1, and go to step 42.

[0045] Step 42, shift the last L i / 2 elements of the i-th frequency band spatially focused spectrum sequence BX ref to the front of BX i , to obtain a shifted spectrum sequence:

[0046] CX i =BX i (L i -L ref / 2), …, BX i (L i ), BX i (1), …, BX i (L i -L ref / 2-1), and go to step 43.

[0047] Step 43, cut off the first L ref elements of the shifted spectrum sequence to form a cut-off spectrum sequence YX i , let i = i + 1, and go to step 44.

[0048] Step 44: Determine if i = N + 1 is true. If not, return to step 42. If true, then obtain N truncated spectrum sequences YX. i YX i The space-time two-dimensional spectrum matrix is ​​obtained by combining the frequency band numbers.

[0049] Step 5: Summate the elements of each column of the spacetime two-dimensional spectrum matrix to obtain the angle-of-arrival projection sequence (e.g., ...). Figure 4 As shown in the diagram, an angle of arrival detection threshold is set, and peak elements exceeding the threshold are identified. Based on the positions corresponding to the peak elements in the angle of arrival projection sequence, the measured angle of arrival sequence is calculated. This step utilizes linear integration operations at different mapping angles for line detection, enabling simultaneous and accurate measurement of spatial parameters for multiple target types. Specifically, the steps include:

[0050] Step 51: Sum the elements of each column of the spacetime two-dimensional spectrum matrix YX to obtain the angle-of-arrival projection sequence Y, in the following form:

[0051]

[0052] Among them, Y l This represents the l-th element of the angle-of-arrival projection sequence Y, with column indices l = 1, 2, ..., L of the two-dimensional spectrum matrix. ref L ref The FFT operation length of the reference frequency band is given by N, the sampling length is given by i, the frequency band number is given by YX(i,l) and the element in the i-th row and l-th column of the space-time two-dimensional spectrum matrix YX is given by YX(i,l).

[0053] Step 52: Set the angle of arrival detection threshold Th based on the angle of arrival projection sequence Y. DOA .

[0054] Step 52_1: Sort the angle of arrival projection sequence Y from smallest to largest to obtain the angle of arrival projection ascending sequence.

[0055] Step 52_2: Based on the ascending sequence of the angle of arrival projection Y′, calculate the median Y of the angle of arrival projection sequence Y. mid for And the median Y mid As the angle of arrival detection threshold, i.e., Th DOA =Y mid .

[0056] Step 53: Compare each element value of the angle of arrival projection sequence Y with the angle of arrival detection threshold Th. DOA Compare and label the results, i.e., when Y l Th greater than the angle of arrival detection threshold DOALab_DOA l is 1, otherwise is 0, and the expression is as follows:

[0057]

[0058] Step 54, find all peak elements in the DOA projection sequence which are greater than the DOA detection threshold using the DOA label sequence Lab_DOA, and confirm the corresponding positions of the peak elements to form a DOA peak position sequence MaxInd_DOA.

[0059] Step 54_1, extract all element numbers corresponding to the elements with value 1 in the DOA label sequence Lab_DOA and form a DOA index sequence Ind_DOA, and the length of the DOA index sequence is L DOA is less than L ref , and go to Step 54_2.

[0060] Step 54_2, perform a difference operation on the DOA index sequence to obtain a DOA index difference sequence DiffInd_DOA, and the specific form is as follows:

[0061]

[0062] wherein Ind_DOA(l DOA ) represents the l DOA th element of Ind_DOA, 1 DOA < L DOA . The length of the DOA index difference sequence is represents the element number in the DOA index difference sequence DiffInd_DOA, It can be known that the DOA index difference sequence satisfies DiffInd_DOA≥1, and go to Step 54_3.

[0063] Step 54_3, initialize the element number of the DOA index difference sequence initialize the DOA peak position sequence MaxInd_DOA=0, and the element number of MaxInd_DOA is initialize the DOA temporary sequence Temp_DOA=0, and the element number of Temp_DOA is go to Step 54_4.

[0064] Step 54_4, if let the l th element value of the DOA temporary sequence Temp_DOA be go to Step 54_5.

[0065] Step 54_5, if Repeat step 54_4; otherwise, let the first position of the arrival angle temporary sequence be... element values Calculate the maximum value of the temporary sequence Temp_DOA (Angle of Arrival), and derive the position corresponding to the peak element, in the following form:

[0066]

[0067] Let the peak position sequence number of the arrival angle Element number of the corner index difference sequence Proceed to step 54_6.

[0068] Step 54_6, if and Then let Let the peak position sequence number of the arrival angle Element number of the corner index difference sequence Proceed to step 54_7.

[0069] Step 54_7, if and Let the element indices of the angle-of-arrival index difference sequence be... Proceed to step 54_8.

[0070] Step 54_8, Judgment With L DOA If they are not equal, return to step 54_4. The resulting peak arrival position sequence is MaxInd_DOA, meaning the sequence length of MaxInd_DOA satisfies...

[0071] Step 55: Based on the peak position sequence of angle of arrival MaxInd_DOA, calculate the sequence of measured angle of arrival DOA = arcsin(MaxInd_DOA), |DOA|≤90°.

[0072] Step 6: Based on the peak position sequence of the angle of arrival, extract the elements of the corresponding columns of the spatiotemporal two-dimensional spectrum matrix to calculate the center frequency measurement sequence and bandwidth measurement sequence. This step utilizes linear integration operations at different mapping angles for line detection, enabling simultaneous and accurate measurement of frequency domain parameters for multiple target types. Specifically, it includes the following steps:

[0073] Step 61: Initialize the element index of the arrival angle peak position sequence MaxInd_DOA as follows: Let the center frequency position sequence Ind_f c sequence number i f =1, the sequence number i of the bandwidth position sequence Ind_BW BW= 1.

[0074] Step 62, extract the column elements of the space-time two-dimensional spectrum matrix YX, compose the specific angle of arrival spectrum sequence YX DOA = YX DOA (i) N , where YX DOA (i) represents the specific angle of arrival spectrum sequence value of the i-th frequency band, and the sequence length of YX DOA is N.

[0075] Step 63, set a frequency domain detection threshold, find the peak value elements in the specific angle of arrival spectrum sequence YX DOA greater than the detection threshold, derive the frequency domain peak position sequence MaxInd_Fre and the frequency domain peak sequence Max_Fre.

[0076] Step 63_1, calculate the minimum value of the specific angle of arrival spectrum sequence YX DOA Set the threshold factor as w, and w>1, calculate the frequency domain detection threshold

[0077] Step 63_2, compare each element value of the specific angle of arrival spectrum sequence YX DOA with the frequency domain detection threshold Th Fre , and mark the result, when the i-th frequency band YX DOA of the specific angle of arrival spectrum sequence YX DOA (i) is greater than the angle of arrival detection threshold Th Fre , set the corresponding element value Lab_Fre i of the frequency domain label sequence Lab_Fre to 1, otherwise to 0, as follows:

[0078]

[0079] Step 63_3, find all peak value elements in the specific angle of arrival spectrum sequence YX DOA greater than the frequency domain detection threshold using the frequency domain label sequence Lab_Fre, and confirm the corresponding positions of the peak value elements, compose the frequency domain peak position sequence MaxInd_Fre, and the corresponding frequency domain peak sequence is Max_Fre=YX DOA (MaxInd_Fre), the sequence length of the frequency domain peak sequence is

[0080] Step 64, according to the frequency domain peak position sequence MaxInd_Fre and the frequency domain peak sequence Max_Fre, calculate the center frequency position sequence Ind_f c and the bandwidth position sequence Ind_BW. ​​

[0081] Step 64_1: Initialize the sequence number of the frequency domain peak sequence. Proceed to step 64_2.

[0082] Step 64_2: Set the bandwidth threshold Th BW This is at the 3dB position of the frequency domain peak sequence Max_Fre, i.e. Proceed to step 64_3.

[0083] Step 64_3: Calculate the spectral magnitude sequence for a specific angle of arrival as |YX DOA | Searching for | YX DOA |the The element whose left side is closest to Th BW (i BW The element index of ) is And the rightmost one is closest to Th BW (i BW The element index of ) is Proceed to step 64_4.

[0084] Step 64_4: Calculate the bandwidth location sequence Let i BW =i BW +1, proceed to step 64_5.

[0085] Step 64_5: Calculate the bandwidth location sequence Let i f =i f +1, proceed to step 64_6.

[0086] Step 64_6, if and Then proceed to step 65, if and make Execute step 62, if and make Perform step 64_2.

[0087] Step 65: Based on the center frequency position sequence Ind_f c The center frequency measurement sequence and bandwidth measurement sequence are calculated using the bandwidth location sequence Ind_BW.

[0088] Step 65_1: Calculate the center frequency measurement sequence f using the bandwidth ΔF. c The specific form is as follows:

[0089] f c =F min +ΔF·Ind_f

[0090] Step 65_2, using the frequency band width ΔF, a sequence of bandwidth measurements BW is calculated, in the form

[0091] BW = ΔF - Ind_BW.

Claims

1. A method for simultaneous multi-domain parameter measurement of wideband signals based on space-time joint processing, characterized in that, The method comprises the following steps: Step 1, after AD sampling and polyphase filtering of the radar signal, filtered data are obtained, and step 2 is entered; Step 2, time domain FFT calculation is performed on the filtered data to obtain a wideband decomposition spectrum matrix, and step 3 is entered; Step 3, space domain FFT calculation is performed on each row of the wideband decomposition spectrum matrix to obtain a space domain focused spectrum sequence, and step 4 is entered; Step 4, by using the periodicity in the frequency domain of the discrete Fourier transform, position shifting and element intercepting are performed on each row element of the space domain focused spectrum sequence to obtain a space-time two-dimensional spectrum matrix, and step 5 is entered; Step 5, summation is performed on each column element of the space-time two-dimensional spectrum matrix to obtain an angle of arrival projection sequence, a detection threshold of the angle of arrival is set, peak value elements greater than the detection threshold are found, and an angle of arrival measurement value sequence is calculated according to the positions corresponding to the peak value elements of the angle of arrival projection sequence, and step 6 is entered; Step 6, according to the angle of arrival peak value position sequence, the elements of the corresponding columns of the space-time two-dimensional spectrum matrix in step 4 are combined to calculate a center frequency measurement value sequence and a bandwidth measurement value sequence.

2. The method according to claim 1, wherein, In step 5, summation is performed on each column element of the space-time two-dimensional spectrum matrix to obtain an angle of arrival projection sequence, a detection threshold of the angle of arrival is set, peak value elements greater than the detection threshold are found, and an angle of arrival measurement value sequence is calculated according to the positions corresponding to the peak value elements of the angle of arrival projection sequence, and the specific process is as follows: Step 51, summation is performed on each column element of the space-time two-dimensional spectrum matrix YX to obtain an angle of arrival projection sequence Y, and the specific form is as follows: wherein Y l represents the lth element of the angle of arrival projection sequence Y, the column number of the two-dimensional frequency spectrum matrix l = 1, 2, …, L ref , L ref is the FFT operation length of the reference frequency band, N represents the sampling length, i represents the frequency band number, and YX(i, l) represents the lth column element of the ith frequency band of the space-time two-dimensional frequency spectrum matrix YX. Step 52, set the angle of arrival detection threshold Th according to the angle of arrival projection sequence Y DOA ; Step 53: Compare each element value of the angle of arrival projection sequence Y with the angle of arrival detection threshold Th. DOA Compare and label the results, i.e., when Y l Th greater than the angle of arrival detection threshold DOA At that time, let the element value corresponding to the reached corner label sequence Lab_DOA be Lab_DOA. l The expression is as follows: 1 for positive and 0 for negative. Step 54, by using an angle of arrival label sequence Lab_DOA, all peak value elements greater than the detection threshold of the angle of arrival in the angle of arrival projection sequence are found, and the positions corresponding to the peak value elements are confirmed to form an angle of arrival peak value position sequence MaxInd_DOA; Step 55, according to the angle of arrival peak value position sequence MaxInd_DOA, an angle of arrival measurement value sequence DOA = arcsin(MaxInd_DOA), |DOA| ≤ 90° is calculated.

3. The method of claim 2, wherein the method is characterized by, In step 52, according to the angle of arrival projection sequence Y, an angle of arrival detection threshold Th is set DOA and specifically includes the following steps: Step 52_1, sort the angle of arrival projection sequence Y in ascending order to obtain an ascending angle of arrival projection sequence Step 52_2, find the median Y of the AOA projected sequence Y according to the ascending sequence Y' mid for Y mid Lref / 2 and take the median Y mid as the AOA detection threshold, i.e. Th DOA = Y mid .​ 4. The method according to claim 3, wherein, In step 54, by using an angle of arrival label sequence Lab_DOA, all peak value elements greater than the detection threshold of the angle of arrival in the angle of arrival projection sequence are found, and the positions corresponding to the peak value elements are confirmed to form an angle of arrival peak value position sequence MaxInd_DOA, and the specific process includes the following steps: Step 54_1, extract the element serial number corresponding to the element with value 1 in the angle of arrival label sequence Lab_DOA and form an angle of arrival index sequence Ind_DOA, the length of the angle of arrival index sequence is L DOA less than L ref , turn to step 54_2; Step 54_2, difference operation is performed on the angle of arrival index sequence to obtain an angle of arrival index difference sequence DiffInd_DOA, and the specific form is as follows: wherein Ind_DOA(l DOA ) denotes the l DOA th element of Ind_DOA, 1 < l DOA < L DOA ; the length of the angle of arrival index difference sequence is denotes the element number in the angle of arrival index difference sequence DiffInd_DOA, It can be known that the angle of arrival index difference sequence satisfies DiffInd_DOA≥1, and turn to step 54_3; Step 54_3, initialize element index of the difference sequence of the angle of arrival Initialize the sequence of the peak position of the angle of arrival MaxInd_DOA = 0, and the element index of MaxInd_DOA is Initialize the temporary sequence of the angle of arrival Temp_DOA = 0, and the element index of Temp_DOA is Go to Step 54_4; Step 54_4, if Let the value of the i-th element of the temporary sequence Temp_DOA be Temp_DOA(i) Go to Step 54_5; Step 54_5, if Repeat step 54_4; otherwise, let the first position of the arrival angle temporary sequence be... element values Calculate the maximum value of the temporary sequence Temp_DOA (Angle of Arrival), and derive the position corresponding to the peak element, in the following form: Element number of the sequence of angle of arrival peak position Element number of the sequence of angle of arrival index difference Go to step 54_6; Step 54_6, if and then let let the sequence of DOA peak position sequence index element index of the sequence of DOA index difference go to Step 54_7; Step 54_7, if and Let the element number of the difference sequence of the angle of arrival index Go to Step 54_8; Step 54_8, judging with L DOA not equal, return to step 54_4, if the sequence of peak position of the angle of arrival MaxInd_DOA is obtained, that is, the sequence length of MaxInd_DOA satisfies 5. The method of claim 4, wherein the method is characterized by, In step 6, according to the angle of arrival peak value position sequence, the elements of the corresponding columns of the space-time two-dimensional spectrum matrix in step 4 are combined to calculate a center frequency measurement value sequence and a bandwidth measurement value sequence, and the specific process includes the following steps: Step 61: Initialize the element index of the arrival angle peak position sequence MaxInd_DOA as follows: Let the center frequency position sequence Ind_f c sequence number i f =1, the sequence number i of the bandwidth position sequence Ind_BW BW =1; Step 62, extract the i-th column element of the space-time two-dimensional spectrum matrix YX, compose a specific angle of arrival spectrum sequence YX DOA DOA (i) N where YX DOA (i) denotes the specific angle of arrival spectrum sequence value of the i-th frequency band, and the sequence length of YX DOA is N;​​ Step 63: Set the frequency domain detection threshold and search for a specific angle-of-arrival spectral sequence YX. DOA For peak elements larger than the detection threshold, derive the frequency domain peak position sequence MaxInd_Fre and the frequency domain peak sequence Max_Fre; Step 64, calculating the center frequency position sequence Ind_f and the bandwidth position sequence Ind_BW according to the frequency domain peak value position sequence MaxInd_Fre and the frequency domain peak value sequence Max_Fre. c and the bandwidth position sequence Ind_BW; Step 65, compute the center frequency measurement sequence and the bandwidth measurement sequence from the center frequency position sequence Ind_f and the bandwidth position sequence Ind_BW. c and the bandwidth position sequence Ind_BW.

6. The method according to claim 5, wherein, In step 63, a frequency domain detection threshold is set to find peak elements in the specific angle spectrum sequence YX DOA greater than the detection threshold, and derive a frequency domain peak position sequence MaxInd_Fre and a frequency domain peak sequence Max_Fre, as follows: Step 63_1, compute the minimum of the specific angle of arrival spectrum sequence YX DOA Set a threshold factor w, and w > 1, compute the frequency domain detection threshold ​ Step 63_2, compare each element value of the specific angle of arrival spectrum sequence YX DOA with the frequency domain detection threshold Th Fre , and mark the result, when the i-th frequency band YX DOA of the specific angle of arrival spectrum sequence YX DOA (i) is greater than the angle of arrival detection threshold Th Fre , the corresponding element value Lab_Fre i of the frequency domain label sequence Lab_Fre is 1, otherwise 0, as follows: Step 63_3, finding the specific angle of arrival spectrum sequence YX with the frequency domain label sequence Lab_Fre DOA All peak elements greater than the frequency domain detection threshold in the middle are confirmed, and the corresponding positions of the peak elements are confirmed to form the frequency domain peak position sequence MaxInd_Fre, and the corresponding frequency domain peak sequence is Max_Fre=YX DOA (MaxInd_Fre), the sequence length of the frequency domain peak sequence is 7. The method of claim 5, wherein the method is characterized by, In step 64, the center frequency position sequence Ind_f and the bandwidth position sequence Ind_BW are calculated according to the frequency domain peak position sequence MaxInd_Fre and the frequency domain peak sequence Max_Fre, and specifically include the following steps: c and the bandwidth position sequence Ind_BW, and specifically include the following steps: Step 64_1, initialize sequence number of sequence of frequency domain peak values Go to Step 64_2; Step 64_2, setting a bandwidth threshold Th BW For the 3dB position of the frequency domain peak sequence Max_Fre, i.e. Go to Step 64_3; Step 64_3: Calculate the spectral magnitude sequence for a specific angle of arrival as |YX DOA | Searching for | YX DOA |the The element whose left side is closest to Th BW (i BW The element index of ) is And the rightmost one is closest to Th BW (i BW The element index of ) is Proceed to step 64_4; Step 64_4, compute the bandwidth position sequence Let i BW = i BW + 1, go to Step 64_5; Step 64_5, compute the bandwidth position sequence Let i f = i f + 1, go to Step 64_6; Step 64_6, if and Step 65 is performed, if and Let Step 62 is performed, if and Let Step 64_2 is performed.

Images