A circular array interferometer direction finding method, apparatus, system, and storage medium

Abstract

This invention provides a direction finding method, apparatus, system, and storage medium for a circular array interferometer, belonging to the field of electronic reconnaissance. The method includes: obtaining a radio signal from a directional antenna array element; processing the radio signal to obtain the original phase and original amplitude; correcting the original phase and original amplitude to obtain the corrected phase and corrected amplitude; analyzing the corrected amplitude to obtain coarse direction finding angles; sorting the coarse direction finding angles in descending order to obtain a sequence of coarse direction finding angles, and using the corrected phases of the directional antennas corresponding to the first N coarse direction finding angles in the sequence as the target phase. This invention can achieve instantaneous 360° omnidirectional high-sensitivity and high-precision direction finding, reducing equipment size, power consumption, and cost requirements, and has significant application value, while also improving the accuracy of direction finding results.

Description

Technical Field

[0001] This invention mainly relates to the field of electronic reconnaissance technology, specifically to a method, device, system, and storage medium for direction finding using a circular array interferometer. Background Technology

[0002] Currently, most interferometer direction finding systems employ a digital multi-baseline architecture. Long baselines are used to ensure direction finding accuracy, while multiple baselines are used for de-ambiguity. However, the shortest baseline length is affected by the size of the individual antenna elements, leading to high-end de-ambiguity issues; the maximum baseline length is limited by the platform's installation dimensions, making it difficult to improve direction finding accuracy.

[0003] To adapt to instantaneous 360° omnidirectional direction finding, four-sided interferometer direction finding stitching, circular array correlative interferometer direction finding, or circular array amplitude and phase comparison direction finding are commonly used. A typical single-sided interferometer direction finding requires four antenna elements to form a multi-baseline interferometer, while a four-sided array requires at least 16 antenna elements and corresponding processing channels. Circular array correlative interferometer direction finding requires omnidirectional antenna elements; for high-precision direction finding, at least nine antenna elements and corresponding processing channels are needed. Furthermore, due to the gain limitation of omnidirectional antennas, the system sensitivity is relatively low, failing to meet the high sensitivity requirements. Circular array amplitude and phase comparison direction finding, through amplitude comparison-guided phase comparison deambiguity resolution, meets the high amplitude comparison direction finding accuracy in the Ku band, requiring at least 12 antenna elements and corresponding processing channels. To meet the technical requirements of instantaneous 360° omnidirectional high-sensitivity and high-precision direction finding, the above three direction finding methods, to some extent, increase equipment size, power consumption, and hardware costs, making installation difficult and hindering the development of equipment miniaturization. Summary of the Invention

[0004] The technical problem to be solved by the present invention is to provide a method, device, system and storage medium for direction finding of a circular array interferometer, which addresses the shortcomings of the prior art.

[0005] The technical solution of the present invention to solve the above-mentioned technical problems is as follows: A direction finding method for a circular array interferometer, comprising the following steps:

[0006] The radio signals of each directional antenna are obtained from multiple directional antenna array elements, and the radio signals of each directional antenna are processed to obtain the original phase and original amplitude of each directional antenna.

[0007] The original phase and original amplitude of each of the directional antennas are corrected to obtain the corrected phase and corrected amplitude of each of the directional antennas.

[0008] The amplitude-ratio coarse direction finding angle of each directional antenna is analyzed by the correction amplitude of each directional antenna to obtain the amplitude-ratio coarse direction finding angle of each directional antenna.

[0009] The amplitude coarse direction finding angles of all the directional antennas are sorted in descending order to obtain the amplitude coarse direction finding angle sequence, and the corrected phases of the directional antennas corresponding to the first N amplitude coarse direction finding angles in the amplitude coarse direction finding angle sequence are all taken as the target phases.

[0010] The target signal test azimuth angle of all the target phases is analyzed to obtain the interferometer direction finding azimuth angle, and the interferometer direction finding azimuth angle is used as the interferometer direction finding result.

[0011] Another technical solution of the present invention to solve the above-mentioned technical problems is as follows: A direction-finding device for a circular array interferometer, comprising:

[0012] The signal processing module is used to obtain the radio signals of each directional antenna from multiple directional antenna array elements, and to process the radio signals of each directional antenna to obtain the original phase and original amplitude of each directional antenna.

[0013] The correction module is used to correct the original phase and original amplitude of each of the directional antennas to obtain the corrected phase and corrected amplitude of each of the directional antennas.

[0014] An angle analysis module is used to analyze the amplitude of each directional antenna after correction to obtain the amplitude coarse direction finding angle of each directional antenna.

[0015] The sorting module is used to sort the amplitude coarse direction finding angles of all the directional antennas in descending order to obtain the amplitude coarse direction finding angle sequence, and to take the corrected phase of the directional antenna corresponding to the first N amplitude coarse direction finding angles in the amplitude coarse direction finding angle sequence as the target phase.

[0016] The direction finding result acquisition module is used to analyze the target signal test azimuth angle of all the target phases, obtain the interferometer direction finding azimuth angle, and use the interferometer direction finding azimuth angle as the interferometer direction finding result.

[0017] Based on the above-mentioned circular array interferometer direction finding method, the present invention also provides a circular array interferometer direction finding system.

[0018] Another technical solution of the present invention to solve the above-mentioned technical problems is as follows: a circular array interferometer direction finding system, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the circular array interferometer direction finding method as described above is implemented.

[0019] Based on the above-described direction finding method using a circular array interferometer, this invention also provides a computer-readable storage medium.

[0020] Another technical solution of the present invention to solve the above-mentioned technical problems is as follows: a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the circular array interferometer direction finding method as described above.

[0021] The beneficial effects of this invention are as follows: By obtaining radio signals from directional antenna array elements, processing the radio signals to obtain the original phase and original amplitude, correcting the original phase and original amplitude to obtain the corrected phase and corrected amplitude, analyzing the amplitude coarse direction finding angles of the corrected amplitude to obtain the amplitude coarse direction finding angles, sorting the amplitude coarse direction finding angles in descending order to obtain the amplitude coarse direction finding angle sequence, and taking the corrected phases of the directional antennas corresponding to the first N amplitude coarse direction finding angles in the amplitude coarse direction finding angle sequence as the target phases, and analyzing the target signal test azimuth angles of the target phases to obtain the interferometer direction finding results, it is possible to achieve instantaneous 360° all-round high-sensitivity and high-precision direction finding, which can reduce the requirements of equipment size, power consumption and cost, has good application value, and also improves the accuracy of direction finding results. Attached Figure Description

[0022] Figure 1 This is a flowchart illustrating a direction-finding method using a circular interferometer provided in an embodiment of the present invention.

[0023] Figure 2 A schematic diagram of a direction finding method using a circular interferometer provided in an embodiment of the present invention;

[0024] Figure 3 This is a block diagram of a circular array interferometer direction finding device provided in an embodiment of the present invention. Detailed Implementation

[0025] The principles and features of the present invention are described below with reference to the accompanying drawings. The examples given are only for explaining the present invention and are not intended to limit the scope of the present invention.

[0026] Figure 1 This is a flowchart illustrating a direction-finding method using a circular array interferometer, as provided in an embodiment of the present invention.

[0027] like Figure 1 and 2 As shown, a direction-finding method using a circular interferometer array includes the following steps:

[0028] The radio signals of each directional antenna are obtained from multiple directional antenna array elements, and the radio signals of each directional antenna are processed to obtain the original phase and original amplitude of each directional antenna.

[0029] The original phase and original amplitude of each of the directional antennas are corrected to obtain the corrected phase and corrected amplitude of each of the directional antennas.

[0030] The amplitude-ratio coarse direction finding angle of each directional antenna is analyzed by the correction amplitude of each directional antenna to obtain the amplitude-ratio coarse direction finding angle of each directional antenna.

[0031] The amplitude coarse direction finding angles of all the directional antennas are sorted in descending order to obtain the amplitude coarse direction finding angle sequence, and the corrected phases of the directional antennas corresponding to the first N amplitude coarse direction finding angles in the amplitude coarse direction finding angle sequence are all taken as the target phases.

[0032] The target signal test azimuth angle of all the target phases is analyzed to obtain the interferometer direction finding azimuth angle, and the interferometer direction finding azimuth angle is used as the interferometer direction finding result.

[0033] Preferably, N can be 3.

[0034] It should be understood that the radio signals are received instantaneously in a 360° airspace by means of eight directional antenna elements arranged in a uniform circular array.

[0035] Specifically, eight directional antenna elements are arranged in a uniform circular array, with each pair of antenna elements spaced 45° apart, to receive radio signals in a 360° airspace instantaneously.

[0036] It should be understood that amplitude / phase measurements are performed on the received signal (i.e., the radio signal).

[0037] It should be understood that amplitude / phase consistency correction is performed on the signal amplitude / phase (i.e., the original phase and the original amplitude) to reduce the impact of factors such as internal noise, external environment and temperature changes on the amplitude / phase consistency of the receiving channel of the equipment.

[0038] Specifically, the amplitude and phase shift differences of the signals received by the 8 antenna array elements are calculated, and an amplitude and phase shift difference correction table is established with the signal frequency value as the address. When the equipment is working normally, the table is looked up according to the signal frequency value to realize the amplitude and phase consistency correction of the signal amplitude / phase, and reduce the impact of external environment and temperature changes on the amplitude and phase consistency of the equipment receiving channel.

[0039] It should be understood that, based on the magnitude of the signal amplitude received by the antenna array element (i.e., the amplitude-ratio coarse direction finding angle), the antenna array element with the largest signal amplitude (i.e., the amplitude-ratio coarse direction finding angle) is selected, and the antenna array element baseline combination corresponding to the 22.5° region is the optimal signal baseline combination.

[0040] In the above embodiments, radio signals are obtained from directional antenna array elements. The original phase and original amplitude are obtained through signal processing of the radio signals. The corrected phase and corrected amplitude are obtained through correction of the original phase and original amplitude. The coarse direction finding angle of the amplitude is obtained through coarse direction finding angle analysis of the corrected amplitude. The coarse direction finding angle of the amplitude is obtained by sorting the coarse direction finding angle of the amplitude in descending order. The corrected phase of the directional antenna corresponding to the first N coarse direction finding angles of the amplitude in the coarse direction finding angle sequence is taken as the target phase. The interferometer direction finding result is obtained by analyzing the azimuth angle of the target signal test of the target phase. This can achieve instantaneous 360° all-round high sensitivity and high precision direction finding, which can reduce the requirements of equipment size, power consumption and cost, and has good application value. It also improves the accuracy of direction finding results.

[0041] Optionally, as an embodiment of the present invention, the process of performing signal processing on the radio signals of each of the directional antennas to obtain the original phase and original amplitude of each of the directional antennas includes:

[0042] The radio signals of each of the directional antennas are converted from analog to digital to obtain the digital signals of each of the directional antennas.

[0043] The digital signals of each of the directional antennas are converted to obtain the original phase and original amplitude of each directional antenna.

[0044] It should be understood that M directional antenna elements are uniformly distributed on a circular array of radius R. A Cartesian coordinate system is established, with the extension of the line connecting the center O of the circle to the first antenna element of the circular array as the Y-axis, and θ as the target azimuth angle. The spatial coordinate position of the m-th antenna element is...

[0045] Specifically, the received signal model of array element m can be expressed as:

[0046]

[0047] The phase shift of array element m is:

[0048]

[0049] In the formula, f c τ is the signal frequency. m Let θ be the time delay from the signal to the center of the antenna array element m, and θ be the incident azimuth angle of the target signal.

[0050] The amplitude of array element m is:

[0051]

[0052] In the formula, θ is the incident azimuth angle of the target signal, θr It is the angle between the two antenna array elements.

[0053] In the above embodiments, the radio signals of each directional antenna are processed to obtain the original phase and original amplitude of each directional antenna, laying the foundation for subsequent data processing. This reduces the requirements for equipment size, power consumption and cost, and has great application value.

[0054] Optionally, as an embodiment of the present invention, the plurality of directional antennas are arranged sequentially;

[0055] The process of analyzing the amplitude of each directional antenna after correction to obtain the amplitude-comparison coarse direction-finding angle includes:

[0056] The difference between the corrected amplitude of each directional antenna and the corrected amplitude of the previous directional antenna is calculated to obtain the antenna element amplitude difference of each directional antenna.

[0057] Import the antenna element angles of each of the directional antennas, and calculate the amplitude comparison direction finding azimuth angles of each of the directional antennas based on the amplitude difference of each antenna element and the antenna element angles.

[0058] Import the signal frequencies of each of the directional antennas, and obtain the first and second parameters of each of the directional antennas from a preset coarse direction finding table based on the signal frequencies of each of the directional antennas and the amplitude difference of the antenna array elements.

[0059] Based on the first formula, the amplitude-comparison coarse direction-finding angle of each directional antenna is calculated according to the amplitude-comparison direction-finding azimuth angle, the first parameter, and the second parameter, to obtain the amplitude-comparison coarse direction-finding angle of each directional antenna. The first formula is:

[0060] θ 比幅 =k·Δθ+b,

[0061] Where, θ 比幅 Δθ is the amplitude-comparison coarse direction finding angle, k is the first parameter, b is the second parameter, and Δθ is the amplitude-comparison direction finding azimuth angle.

[0062] It should be understood that amplitude comparison and direction finding are performed on the signal (i.e., the corrected amplitude) to optimize the baseline combination.

[0063] Specifically, based on the signal measurement frequency (i.e., the signal frequency) and amplitude difference (i.e., the antenna element amplitude difference), the k value (i.e., the first parameter) and b value (i.e., the second parameter) in the corresponding coarse direction finding table (i.e., the preset coarse direction finding table) are looked up, and then the formula θ is applied. 比幅 =k·Δθ+b, which gives the coarse direction finding angle of the target signal.

[0064] In the above embodiments, the amplitude-ratio coarse direction finding angle of each directional antenna is analyzed by comparing the corrected amplitude, so as to obtain the amplitude-ratio coarse direction finding angle of each directional antenna. This provides data for the subsequent optimal signal baseline combination, realizes coarse direction finding, and improves the accuracy of the direction finding results.

[0065] Optionally, as an embodiment of the present invention, the process of calculating the amplitude comparison direction-finding azimuth angle of each of the directional antennas based on the amplitude difference of the antenna elements and the included angle of the antenna elements includes:

[0066] Based on the second equation, the amplitude-comparison direction-finding azimuth angle of each directional antenna is calculated according to the amplitude difference of the antenna array elements and the included angle of the antenna array elements. The second equation is:

[0067]

[0068] Where Δθ is the amplitude-comparison azimuth angle, ΔR is the amplitude difference between antenna elements, and θ 3dB θ is the beamwidth of the antenna array element. γ The angle between the antenna array elements.

[0069] It should be understood that the amplitude-based direction finding formula is:

[0070] In the formula, Δθ is the amplitude-comparison azimuth angle, ΔR is the amplitude difference between the two antenna elements (i.e., the amplitude difference between the antenna elements), and θ 3dB θ is the beamwidth of the antenna array element. r It is the angle between the two antenna elements (i.e., the angle between the antenna elements).

[0071] In the above embodiments, the amplitude-comparison azimuth angle of each directional antenna is calculated based on the second formula according to the amplitude difference of each directional antenna element and the included angle of each directional antenna element. This provides data for the subsequent optimal signal baseline combination, realizes coarse direction finding, and improves the accuracy of the direction finding results.

[0072] Optionally, as an embodiment of the present invention, N is 3, and all the target phases include a first target phase, a second target phase, and a third target phase arranged in sequence;

[0073] The process of analyzing the target signal test azimuth angle of all the target phases to obtain the interferometer direction finding azimuth angle includes:

[0074] The difference between the phase of the first target and the phase of the second target is calculated to obtain the short baseline phase difference;

[0075] The difference between the phase of the first target and the phase of the third target is calculated to obtain the long baseline phase difference;

[0076] Based on the third equation, the short baseline ambiguity number is calculated according to the short baseline phase difference, resulting in multiple short baseline ambiguity numbers. The third equation is:

[0077]

[0078] in,

[0079] Where, k 8,1 d represents the number of short baseline ambiguities. 8,1 For short baselines, Δφ 8,1 Let λ be the short baseline phase difference, λ be the signal wavelength, and R be the radius;

[0080] The long baseline ambiguity number is calculated based on the long baseline phase difference, the short baseline phase difference, and each of the short baseline ambiguity numbers, to obtain a long baseline ambiguity number that corresponds one-to-one with each of the short baseline ambiguity numbers.

[0081] An analysis of the estimated values ​​of all the aforementioned long baseline ambiguity numbers yields the estimated values;

[0082] Based on the fourth equation, the interferometer direction-finding azimuth angle is calculated according to the estimated value and the long baseline phase difference, and the interferometer direction-finding azimuth angle is obtained. The fourth equation is:

[0083]

[0084] in,

[0085] Where θ is the interferometer direction finding azimuth angle, λ is the signal wavelength, k0 is the estimated value, and d 8,2 For a long baseline, Δφ 8,2 This represents the long baseline phase difference.

[0086] It should be understood that the optimal signal baseline combination is selected for direction finding calculations using a circular array dual-baseline 3-antenna element interferometer.

[0087] Specifically, taking the construction of a dual-baseline three-antenna nonlinear array antenna array with antenna elements (8, 1, 2) as an example, the ratio of the dual baselines of the circular array satisfies The relationship (where p and q are coprime numbers) allows us to obtain the corresponding short baseline d. 8,1 and long baseline d 8,2 The true phase difference is shown in the following formula:

[0088]

[0089]

[0090] Among them, short baseline Long baseline φ 8,1 For the true phase difference of the short baseline, φ 8,2 For the true phase difference of the long baseline, k 8,1 k is the number of short baseline ambiguities. 8,2 This represents the number of long baseline ambiguities.

[0091] Based on baseline d 8,1 d 8,2 Coprime property, we get:

[0092]

[0093] After sorting, we can obtain:

[0094]

[0095] Because |sinθ|≤1, for the fuzzy number k 8,1 Further constraints on the search range yield:

[0096]

[0097] It should be understood that, according to the formula:

[0098] In the formula, θ is the target signal test azimuth angle (i.e., the interferometer direction finding azimuth angle), λ is the signal wavelength, and k is the wavelength. 8,2,0 That is, the estimated value k0.

[0099] In the above embodiments, the interferometer direction finding azimuth angle is obtained by analyzing the target signal test azimuth angle of all target phases. This can achieve instantaneous 360° all-round high sensitivity and high precision direction finding, and also achieve fine direction finding. It can reduce the requirements of equipment size, power consumption and cost, and has great application value, and also improves the accuracy of direction finding results.

[0100] Optionally, as an embodiment of the present invention, the process of calculating the long baseline ambiguity number based on the long baseline phase difference, the short baseline phase difference, and each of the short baseline ambiguity numbers to obtain the long baseline ambiguity number corresponding one-to-one with each of the short baseline ambiguity numbers includes:

[0101] Based on the fifth equation, the long baseline ambiguity number is calculated according to the long baseline phase difference, the short baseline phase difference, and each of the short baseline ambiguity numbers, to obtain the long baseline ambiguity number corresponding one-to-one with each of the short baseline ambiguity numbers. The fifth equation is:

[0102]

[0103] Where, k 8,2Let d be the number of long baseline ambiguities. 8,2 For a long baseline, Δφ 8,2 For long baseline phase difference, d 8,1 For short baselines, Δφ 8,1 For short baseline phase difference, k 8,1 This represents the number of short baseline ambiguities.

[0104] It should be understood that when k 8,1 When the value changes within the specified search range, the corresponding determined k 8,2 Also with k 8,1 The value changes with the change:

[0105] In the above embodiments, based on the fifth formula, the long baseline ambiguity number is calculated according to the long baseline phase difference, the short baseline phase difference, and each short baseline ambiguity number to obtain the long baseline ambiguity number that corresponds one-to-one with each short baseline ambiguity number. This can reduce the requirements for equipment size, power consumption, and cost, and has great application value, while also improving the accuracy of direction finding results.

[0106] Optionally, as an embodiment of the present invention, the process of analyzing the estimated values ​​of all the long baseline ambiguity numbers to obtain the estimated values ​​includes:

[0107] Based on the sixth equation, the unprocessed ambiguity number is calculated for each of the long baseline ambiguity numbers that corresponds one-to-one with each of the short baseline ambiguity numbers, to obtain the unprocessed ambiguity number for each of the long baseline ambiguity numbers. The sixth equation is:

[0108] Δk i =(k 82i -round(k 82i )) 2 ,

[0109] Where, Δk i Let k be the number of ambiguities to be processed for the i-th long baseline ambiguity number. 82i Let i be the number of long baseline ambiguities;

[0110] The minimum value of the unprocessed fuzzy number among all the long baseline fuzzy numbers is selected, and the minimum unprocessed fuzzy number is obtained after selection. The long baseline fuzzy number corresponding to the minimum unprocessed fuzzy number is then used as an estimated value.

[0111] Specifically, for k 81 The minimum value within the specified search range is taken as the first possible fuzzy number, and the result is taken as k. 8,2 The closest integer k 8,2,0 As k 8,2 The estimated value is as follows, according to the relevant function formula:

[0112] Δk i =(k82i -round(k 82i )) 2

[0113] Compare k 81 Given all fuzzy numbers within the specified search range (i.e., the fuzzy numbers to be processed), find the minimum value of the correlation function Δk (i.e., the minimum fuzzy number to be processed), and the corresponding k. 8,2,0 That is, the long baseline d 8,2 The corresponding true fuzzy number (i.e., the estimated value) and k 8,2,0 The corresponding k 81 That is, the short baseline d 8,1 Real fuzzy number k 8,1,0 Therefore, there is no ambiguity phase difference.

[0114] In the above embodiments, the estimated values ​​are obtained by analyzing the estimated values ​​of all long baseline ambiguity numbers, which improves the accuracy of direction finding results and can reduce the requirements for equipment size, power consumption and cost, and has great application value.

[0115] Optionally, as another embodiment of the present invention, the present invention includes eight directional antenna elements arranged in a uniform circular array, with a 45° interval between each pair of antenna elements, to complete the instantaneous reception of radio signals in a 360° airspace; the direction finding divides the 360° omnidirectional range into 16 22.5° regions, firstly by using amplitude comparison coarse direction finding to select the 22.5° regions, then by using a circular array dual-baseline interferometer for fine direction finding according to the optimal baseline combination, and finally by combining the amplitude comparison direction finding results and the circular array dual-baseline interferometer direction finding results, and judging the reliability of the direction finding results according to certain constraints, thereby achieving high-precision direction finding of the target signal.

[0116] Optionally, as another embodiment of the present invention, the present invention optimizes the combination of eight directional antenna elements placed in a uniform circular array, makes full use of the phase relationship between the elements to construct the phase vector, uses amplitude comparison direction finding to achieve coarse direction finding, and uses a circular array dual-baseline interferometer to achieve fine direction finding. It can achieve instantaneous 360° all-round high-sensitivity and high-precision direction finding. The performance of the equipment is comparable to the direction finding of a four-sided interferometer with 16 antenna elements, and is superior to the amplitude and phase comparison circular array direction finding of 12 antenna elements. It can reduce the requirements of equipment size, power consumption and cost, and has good economic value and application prospects.

[0117] Optionally, as another embodiment of the present invention, the receiving antenna array of the present invention consists of eight directional antenna elements arranged in a uniform circular array, with a 45° interval between each pair of antenna elements, receiving radio signals in a 360° instantaneous airspace. By optimizing the combination of directional antenna elements and making full use of the phase relationship between the elements to construct the phase vector, coarse direction finding is achieved by amplitude comparison direction finding, and fine direction finding is achieved by a circular array dual-baseline interferometer. Instantaneous 360° omnidirectional high-sensitivity and high-precision direction finding can be achieved. The performance of the equipment is comparable to that of a 16-antenna-element four-sided array interferometer direction finding, and superior to that of a 12-antenna-element amplitude and phase comparison circular array direction finding. It can reduce the requirements for equipment size, power consumption and cost, and has great application value.

[0118] Optionally, as another embodiment of the present invention, the measurement results of amplitude difference and phase difference of the present invention always have certain errors. It is necessary to combine the amplitude comparison direction finding results and the direction finding results of the circular array dual baseline interferometer, and judge the reliability of the direction finding results according to certain constraints, and output the direction finding results.

[0119] Figure 3 This is a block diagram of a circular array interferometer direction finding device provided in an embodiment of the present invention.

[0120] Alternatively, as another embodiment of the present invention, such as Figure 3 As shown, a circular array interferometer direction finding device includes:

[0121] The signal processing module is used to obtain the radio signals of each directional antenna from multiple directional antenna array elements, and to process the radio signals of each directional antenna to obtain the original phase and original amplitude of each directional antenna.

[0122] The correction module is used to correct the original phase and original amplitude of each of the directional antennas to obtain the corrected phase and corrected amplitude of each of the directional antennas.

[0123] An angle analysis module is used to analyze the amplitude of each directional antenna after correction to obtain the amplitude coarse direction finding angle of each directional antenna.

[0124] The sorting module is used to sort the amplitude coarse direction finding angles of all the directional antennas in descending order to obtain the amplitude coarse direction finding angle sequence, and to take the corrected phase of the directional antenna corresponding to the first N amplitude coarse direction finding angles in the amplitude coarse direction finding angle sequence as the target phase.

[0125] The direction finding result acquisition module is used to analyze the target signal test azimuth angle of all the target phases, obtain the interferometer direction finding azimuth angle, and use the interferometer direction finding azimuth angle as the interferometer direction finding result.

[0126] Optionally, another embodiment of the present invention provides a circular array interferometer direction finding system, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the circular array interferometer direction finding method as described above. This system can be a computer or similar system.

[0127] Optionally, another embodiment of the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the circular array interferometer direction finding method as described above.

[0128] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0129] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working process of the above-described apparatus and unit can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0130] In the several embodiments provided in this application, it should be understood that the disclosed apparatus and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative. For instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed.

[0131] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of the embodiments of the present invention, depending on actual needs.

[0132] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0133] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0134] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A direction-finding method for a circular interferometer array, characterized in that, Includes the following steps: The radio signals of each directional antenna are obtained from multiple directional antenna array elements, and the radio signals of each directional antenna are processed to obtain the original phase and original amplitude of each directional antenna. The original phase and original amplitude of each of the directional antennas are corrected to obtain the corrected phase and corrected amplitude of each of the directional antennas. The amplitude-ratio coarse direction finding angle of each directional antenna is analyzed by the correction amplitude of each directional antenna to obtain the amplitude-ratio coarse direction finding angle of each directional antenna. The amplitude coarse direction finding angles of all the directional antennas are sorted in descending order to obtain the amplitude coarse direction finding angle sequence, and the corrected phases of the directional antennas corresponding to the first N amplitude coarse direction finding angles in the amplitude coarse direction finding angle sequence are all taken as the target phases. The target signal test azimuth angle of all the target phases is analyzed to obtain the interferometer direction finding azimuth angle, and the interferometer direction finding azimuth angle is used as the interferometer direction finding result; The N is 3, and all the target phases include a first target phase, a second target phase, and a third target phase arranged in sequence; The process of analyzing the target signal test azimuth angle of all the target phases to obtain the interferometer direction finding azimuth angle includes: The difference between the phase of the first target and the phase of the second target is calculated to obtain the short baseline phase difference; The difference between the phase of the first target and the phase of the third target is calculated to obtain the long baseline phase difference; Based on the third equation, the short baseline ambiguity number is calculated according to the short baseline phase difference, resulting in multiple short baseline ambiguity numbers. The third equation is: ， in, , in, For short baseline ambiguity, For short baselines, For short baseline phase difference, For the signal wavelength, The radius is ; The long baseline ambiguity number is calculated based on the long baseline phase difference, the short baseline phase difference, and each of the short baseline ambiguity numbers, to obtain a long baseline ambiguity number that corresponds one-to-one with each of the short baseline ambiguity numbers. An analysis of the estimated values ​​of all the aforementioned long baseline ambiguity numbers yields the estimated values; Based on the fourth equation, the interferometer direction-finding azimuth angle is calculated according to the estimated value and the long baseline phase difference, and the interferometer direction-finding azimuth angle is obtained. The fourth equation is: ， in, , in, For the azimuth angle of the interferometer direction finding, For the signal wavelength, This is an estimated value. For long baselines, This represents the long baseline phase difference.

2. The direction finding method for a circular interferometer according to claim 1, characterized in that, The process of processing the radio signals of each of the directional antennas to obtain the original phase and original amplitude of each directional antenna includes: The radio signals of each of the directional antennas are converted from analog to digital to obtain the digital signals of each of the directional antennas. The digital signals of each of the directional antennas are converted to obtain the original phase and original amplitude of each directional antenna.

3. The direction finding method for a circular interferometer according to claim 1, characterized in that, The directional antennas are arranged in sequence; The process of analyzing the amplitude of each directional antenna after correction to obtain the amplitude-comparison coarse direction-finding angle includes: The difference between the corrected amplitude of each directional antenna and the corrected amplitude of the previous directional antenna is calculated to obtain the antenna element amplitude difference of each directional antenna. Import the antenna element angles of each of the directional antennas, and calculate the amplitude comparison direction finding azimuth angles of each of the directional antennas based on the amplitude difference of each antenna element and the antenna element angles. Import the signal frequencies of each of the directional antennas, and obtain the first and second parameters of each of the directional antennas from a preset coarse direction finding table based on the signal frequencies of each of the directional antennas and the amplitude difference of the antenna array elements. Based on the first formula, the amplitude-comparison coarse direction-finding angle of each directional antenna is calculated according to the amplitude-comparison direction-finding azimuth angle, the first parameter, and the second parameter, to obtain the amplitude-comparison coarse direction-finding angle of each directional antenna. The first formula is: ， in, For the coarse direction finding angle of the amplitude comparison, As the first parameter, For the second parameter, This is the azimuth angle for amplitude-based direction finding.

4. The direction finding method for a circular interferometer according to claim 3, characterized in that, The process of calculating the amplitude-comparison direction-finding azimuth angle of each directional antenna based on the amplitude difference and the included angle of each antenna element includes: Based on the second equation, the amplitude-comparison direction-finding azimuth angle of each directional antenna is calculated according to the amplitude difference of the antenna array elements and the included angle of the antenna array elements. The second equation is: ， in, For amplitude comparison direction finding azimuth, For the amplitude difference of the antenna array elements, The beamwidth of the antenna array element. The angle between the antenna array elements.

5. The direction finding method for a circular interferometer according to claim 1, characterized in that, The process of calculating the long baseline ambiguity number based on the long baseline phase difference, the short baseline phase difference, and each of the short baseline ambiguity numbers, to obtain the long baseline ambiguity number corresponding one-to-one with each of the short baseline ambiguity numbers includes: Based on the fifth equation, the long baseline ambiguity number is calculated according to the long baseline phase difference, the short baseline phase difference, and each of the short baseline ambiguity numbers, to obtain the long baseline ambiguity number corresponding one-to-one with each of the short baseline ambiguity numbers. The fifth equation is: ， in, For long baseline ambiguity, For long baselines, For long baseline phase difference, For short baselines, For short baseline phase difference, This represents the number of short baseline ambiguities.

6. The direction finding method for a circular interferometer array according to claim 1, characterized in that, The process of estimating all the long baseline ambiguity numbers to obtain the estimated values ​​includes: Based on the sixth equation, the unprocessed ambiguity number is calculated for each of the long baseline ambiguity numbers that corresponds one-to-one with each of the short baseline ambiguity numbers, to obtain the unprocessed ambiguity number for each of the long baseline ambiguity numbers. The sixth equation is: ， in, For the first The number of fuzzy numbers to be processed for a long baseline fuzzy number, For the first Number of long baseline ambiguities; The minimum value of the unprocessed fuzzy number among all the long baseline fuzzy numbers is selected, and the minimum unprocessed fuzzy number is obtained after selection. The long baseline fuzzy number corresponding to the minimum unprocessed fuzzy number is then used as an estimated value.

7. A direction-finding device for a circular array interferometer, characterized in that, include: The signal processing module is used to obtain the radio signals of each directional antenna from multiple directional antenna array elements, and to process the radio signals of each directional antenna to obtain the original phase and original amplitude of each directional antenna. The correction module is used to correct the original phase and original amplitude of each of the directional antennas to obtain the corrected phase and corrected amplitude of each of the directional antennas. An angle analysis module is used to analyze the amplitude of each directional antenna after correction to obtain the amplitude coarse direction finding angle of each directional antenna. The sorting module is used to sort the amplitude coarse direction finding angles of all the directional antennas in descending order to obtain the amplitude coarse direction finding angle sequence, and to take the corrected phase of the directional antenna corresponding to the first N amplitude coarse direction finding angles in the amplitude coarse direction finding angle sequence as the target phase. The direction finding result acquisition module is used to analyze the target signal test azimuth angle of all the target phases, obtain the interferometer direction finding azimuth angle, and use the interferometer direction finding azimuth angle as the interferometer direction finding result; The N is 3, and all the target phases include a first target phase, a second target phase, and a third target phase arranged in sequence; The angle analysis module is specifically used for: The difference between the phase of the first target and the phase of the second target is calculated to obtain the short baseline phase difference; The difference between the phase of the first target and the phase of the third target is calculated to obtain the long baseline phase difference; Based on the third equation, the short baseline ambiguity number is calculated according to the short baseline phase difference, resulting in multiple short baseline ambiguity numbers. The third equation is: ， in, , in, For short baseline ambiguity, For short baselines, For short baseline phase difference, For the signal wavelength, The radius is ; The long baseline ambiguity number is calculated based on the long baseline phase difference, the short baseline phase difference, and each of the short baseline ambiguity numbers, to obtain a long baseline ambiguity number that corresponds one-to-one with each of the short baseline ambiguity numbers. An analysis of the estimated values ​​of all the aforementioned long baseline ambiguity numbers yields the estimated values; Based on the fourth equation, the interferometer direction-finding azimuth angle is calculated according to the estimated value and the long baseline phase difference, and the interferometer direction-finding azimuth angle is obtained. The fourth equation is: ， in, , in, For the azimuth angle of the interferometer direction finding, For the signal wavelength, This is an estimated value. For long baselines, This represents the long baseline phase difference.

8. A circular array interferometer direction finding system, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the circular array interferometer direction finding method as described in any one of claims 1 to 6.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the circular array interferometer direction finding method as described in any one of claims 1 to 6.

Images