An ultra-wideband directional antenna uniform circular array deblurring method, system and medium

Abstract

An ultra-wideband directional antenna uniform circular array deblurring method, system and medium are disclosed.The application discloses an ultra-wideband directional antenna uniform circular array deblurring method, and the application is improved on the basis of the stereoscopic baseline method in view of the problems that the stereoscopic baseline method lacks flexibility in baseline pair selection, and when the method is applied to the directional antenna uniform circular array, part of the antenna receives weak signals, which leads to difficulty in deblurring and a large search range of candidate angles.The application utilizes the directivity of the antenna to participate in the deblurring process, so that the positive and negative of the slopes of the two baselines need not be opposite in the selection of the baseline pair, and only the antenna with the largest amplitude in the antenna array and three or more antennas near the antenna are needed to form a baseline pair, and then deblurring can be performed.The search range of the direction of the incoming wave is limited by the antenna with the largest amplitude, the number of candidate angle values is reduced, the calculation amount and error of the subsequent comparison process are reduced, and the application can be widely applied to the technical field of signal processing.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of signal processing, and in particular to an ultra-wideband directional antenna uniform circular array deblurring method, system and medium. BACKGROUND

[0002] A phase interferometer is a device that uses the phase difference measurement results between different array elements, in combination with the spatial position relationship of the array elements, to obtain the incident direction of a signal through a phase comparison method. In order to obtain high direction-finding accuracy, the spacing between the two array elements on the baseline for obtaining phase information needs to be large enough. However, when the spacing between the array elements is greater than half a wavelength, a phase ambiguity problem occurs, which results in a contradiction between the direction-finding accuracy and the phase ambiguity.

[0003] In order to solve the phase ambiguity problem, some deblurring methods have been proposed, such as the long-short baseline method, the remainder theorem method, the virtual baseline method, etc. These methods have problems such as single array form, physical size limitation, etc., and cannot be used on a uniform circular array. Although the stereoscopic baseline method (arbitrary baseline method) does not have requirements for the array form, it still has the following problems. First, it needs to additionally solve a mirror image ambiguity called azimuth angle, and in order to solve the mirror image ambiguity, the positive and negative of the slopes of the two baselines of the selected baseline pair must be opposite, which makes the baseline pair selection lack flexibility, and some antennas receive weak signals, making it difficult to deblur. Therefore, the stereoscopic baseline method cannot be directly used on a directional antenna array. Second, since the stereoscopic baseline method uses omnidirectional antennas, it needs to search for the angle in all directions during the deblurring stage, which makes the number of candidate angle values large when the direction of arrival is finally obtained through the comparison method, resulting in a large amount of calculation. In addition, since the candidate angle values are in all directions, the probability of incorrect deblurring is large at high frequencies. SUMMARY

[0004] Therefore, the embodiments of the present application provide an ultra-wideband directional antenna uniform circular array deblurring method, which can efficiently deblur and effectively reduce the amount of calculation and errors.

[0005] In one aspect, the embodiments of the present application provide an ultra-wideband directional antenna uniform circular array deblurring method, comprising:

[0006] obtaining amplitude coarse guidance information, determining a target antenna and a range of the direction of arrival of the incoming wave according to the amplitude coarse guidance information;

[0007] determining a deblurring baseline group according to the target antenna and the adjacent antennas of the target antenna;

[0008] determining the range of the phase difference ambiguity value of each baseline pair of the deblurring baseline group according to the range of the direction of arrival of the incoming wave and in combination with the redundant angle value;

[0009] According to the phase difference ambiguity number value range, a fuzzy multi-value solution set is obtained; wherein, the fuzzy multi-value solution set includes fuzzy multi-value solutions of each baseline pair in the ambiguity resolution baseline group;

[0010] According to the fuzzy multi-value solution set, an ambiguity resolution result is obtained through clustering processing.

[0011] Optionally, the amplitude coarse guidance information is acquired, and the target antenna and the wave direction angle range are determined according to the amplitude coarse guidance information, including:

[0012] The amplitude coarse guidance information is acquired, and the target antenna is determined according to the amplitude coarse guidance information;

[0013] The wave direction angle range is determined according to the beam intersection point of the target antenna and adjacent antennas of the target antenna.

[0014] Optionally, the ambiguity resolution baseline group is determined according to the target antenna and adjacent antennas of the target antenna, including:

[0015] The ambiguity resolution baseline group of four baseline pairs is formed according to the target antenna and four adjacent antennas of the target antenna.

[0016] Optionally, the phase difference ambiguity number value range of each baseline pair in the ambiguity resolution baseline group is determined according to the wave direction angle range and in combination with the redundant angle value, including:

[0017] The beam range of each baseline pair in the ambiguity resolution baseline group is determined according to the wave direction angle range and in combination with the redundant angle value;

[0018] The phase difference ambiguity number value range of each baseline pair in the ambiguity resolution baseline group is determined according to the beam range.

[0019] Optionally, the fuzzy multi-value solution set is obtained according to the phase difference ambiguity number value range, including:

[0020] The incident signal azimuth tangent value ambiguity solution and the incident signal elevation cosine value ambiguity solution are determined according to the phase difference ambiguity number value range;

[0021] The incident signal azimuth sine value ambiguity solution and the incident signal azimuth cosine value ambiguity solution are determined according to the incident signal azimuth tangent value and in combination with the quadrant range to which the beam range belongs;

[0022] The fuzzy multi-value solution includes the incident signal azimuth sine value ambiguity solution, the incident signal azimuth cosine value ambiguity solution and the incident signal elevation cosine value ambiguity solution, and the beam range is determined through the wave direction angle range and in combination with the redundant angle value.

[0023] Optionally, the method further comprises:

[0024] According to the fuzzy multi-value solution set, a group of fuzzy multi-value solutions is selected as a reference group;

[0025] According to the fuzzy multi-value solutions other than the reference group in the fuzzy multi-value solution set, a target estimated value of the fuzzy multi-value solution set is determined by calculating the Euclidean distance with the reference group;

[0026] A mean operation is performed on the target estimated value to obtain a deblurring result.

[0027] Optionally, the method further comprises:

[0028] According to a preset condition, a simulation direction finding process is performed to obtain a deblurring simulation result;

[0029] The deblurring simulation result includes a successful deblurring probability and a direction finding error.

[0030] In another aspect, an embodiment of the present application provides a deblurring system for a uniform circular array of an ultra-wideband directional antenna, comprising:

[0031] A first module is configured to obtain amplitude coarse guidance information, and determine a target antenna and a wave direction range according to the amplitude coarse guidance information;

[0032] A second module is configured to determine a deblurring baseline group according to the target antenna and a neighboring antenna of the target antenna;

[0033] A third module is configured to determine a phase difference ambiguity number value range of each baseline pair of the deblurring baseline group according to the wave direction range and a redundant angle value;

[0034] A fourth module is configured to obtain a fuzzy multi-value solution set according to the phase difference ambiguity number value range; the fuzzy multi-value solution set includes fuzzy multi-value solutions of each baseline pair of the deblurring baseline group;

[0035] A fifth module is configured to obtain a deblurring result by clustering processing according to the fuzzy multi-value solution set.

[0036] In another aspect, an embodiment of the present application provides an electronic device, comprising a processor and a memory;

[0037] The memory is configured to store a program;

[0038] The processor is configured to execute the program to implement the method according to the first aspect of the present application.

[0039] In another aspect, an embodiment of the present application provides a computer readable storage medium storing a program, the program being executed by a processor to implement the method as described above.

[0040] The embodiment of the present application also discloses a computer program product or computer program, which comprises computer instructions stored in a computer readable storage medium. A processor of a computer device can read the computer instructions from the computer readable storage medium, and the processor executes the computer instructions to enable the computer device to perform the method as described above.

[0041] The embodiment of the present application first acquires amplitude coarse guidance information, determines a target antenna and a range of azimuth angles of a wave based on the amplitude coarse guidance information, determines an ambiguity resolving baseline group based on the target antenna and adjacent antennas of the target antenna, determines a range of values of phase difference ambiguities of each baseline pair of the ambiguity resolving baseline group based on the range of azimuth angles of the wave and in combination with redundant angle values, obtains an ambiguity multi-value solution set based on the range of values of the phase difference ambiguities, wherein the ambiguity multi-value solution set comprises ambiguity multi-value solutions of each baseline pair of the ambiguity resolving baseline group, and obtains an ambiguity resolving result through clustering processing based on the ambiguity multi-value solution set. The embodiment of the present application determines a target antenna based on coarse guidance information, and forms baseline pairs based on the target antenna and adjacent antennas of the target antenna, so that ambiguity resolving can be performed. Moreover, the embodiment of the present application limits the search range of the direction of a wave based on the target antenna, reduces the number of candidate angle values, and thus reduces the calculation amount and error of a subsequent comparison process. BRIEF DESCRIPTION OF DRAWINGS

[0042] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings needed in the embodiment description. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort based on these drawings.

[0043] Figure 1 A flowchart of the super-wideband directional antenna uniform circular array ambiguity resolving method provided by the embodiment of the present application is shown in the figure.

[0044] Figure 2 A two-dimensional phase interferometer direction finding working model is shown in the figure.

[0045] Figure 3 A multi-value ambiguity resolving flowchart provided by the embodiment of the present application is shown in the figure.

[0046] Figure 4 An antenna geometric model provided by the embodiment of the present application is shown in the figure. DETAILED DESCRIPTION

[0047] In order to make the objects, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and not intended to limit the present application.

[0048] In one aspect, with reference to Figure 1 Embodiments of the present application provide a method for resolving ambiguity of a uniform circular array of an ultra-wideband directional antenna, comprising:

[0049] obtaining amplitude coarse guidance information, determining a target antenna and a range of azimuth angles of incoming waves according to the amplitude coarse guidance information;

[0050] determining a baseline group for resolving ambiguity according to the target antenna and adjacent antennas of the target antenna;

[0051] determining a range of values of phase difference ambiguity numbers of each baseline pair of the baseline group for resolving ambiguity according to the range of azimuth angles of incoming waves and in combination with redundant angle values;

[0052] obtaining a multi-value ambiguity solution set according to the range of values of phase difference ambiguity numbers; wherein the multi-value ambiguity solution set comprises multi-value ambiguity solutions of each baseline pair of the baseline group for resolving ambiguity;

[0053] obtaining a result of resolving ambiguity through clustering processing according to the multi-value ambiguity solution set.

[0054] Optionally, obtaining amplitude coarse guidance information, determining a target antenna and a range of azimuth angles of incoming waves according to the amplitude coarse guidance information, comprises:

[0055] obtaining amplitude coarse guidance information, determining a target antenna according to the amplitude coarse guidance information;

[0056] determining the range of azimuth angles of incoming waves according to beam intersection points of the target antenna and adjacent antennas of the target antenna.

[0057] Optionally, determining a baseline group for resolving ambiguity according to the target antenna and adjacent antennas of the target antenna, comprises:

[0058] determining a baseline group for resolving ambiguity of four baseline pairs according to the target antenna and four adjacent antennas of the target antenna.

[0059] Optionally, determining a range of values of phase difference ambiguity numbers of each baseline pair of the baseline group for resolving ambiguity according to the range of azimuth angles of incoming waves and in combination with redundant angle values, comprises:

[0060] determining a beam range of each baseline pair of the baseline group for resolving ambiguity according to the range of azimuth angles of incoming waves and in combination with redundant angle values;

[0061] determining a range of values of phase difference ambiguity numbers of each baseline pair of the baseline group for resolving ambiguity according to the beam range.

[0062] Optionally, according to the phase difference ambiguity number value range, a fuzzy multi-value solution set is obtained, including:

[0063] According to the phase difference ambiguity number value range, a tangent value ambiguity solution of the incident signal azimuth angle and a cosine value ambiguity solution of the incident signal elevation angle are determined.

[0064] According to the tangent value of the incident signal azimuth angle, in combination with the quadrant range to which the beam range belongs, a sine value ambiguity solution of the incident signal azimuth angle and a cosine value ambiguity solution of the incident signal azimuth angle are determined.

[0065] The fuzzy multi-value solution includes the sine value ambiguity solution of the incident signal azimuth angle, the cosine value ambiguity solution of the incident signal azimuth angle, and the cosine value ambiguity solution of the incident signal elevation angle, and the beam range is determined by the wave direction angle range in combination with the redundancy angle value.

[0066] Optionally, according to the fuzzy multi-value solution set, a solution ambiguity result is obtained through clustering processing, including:

[0067] According to the fuzzy multi-value solution set, a group of fuzzy multi-value solutions is selected as a reference group.

[0068] According to the fuzzy multi-value solutions of the fuzzy multi-value solution set except the reference group, a target estimation value of the fuzzy multi-value solution set is determined by calculating the Euclidean distance with the reference group.

[0069] The target estimation value is subjected to mean value operation to obtain the solution ambiguity result.

[0070] Optionally, the method further includes:

[0071] According to a preset condition, a simulation direction finding processing is performed to obtain a solution ambiguity simulation result.

[0072] The solution ambiguity simulation result includes a successful solution ambiguity probability and a direction finding error.

[0073] On the other hand, an embodiment of the present application provides a uniform circular array solution ambiguity system of an ultra-wideband directional antenna, including:

[0074] A first module is configured to acquire amplitude coarse guidance information, and determine a target antenna and a wave direction angle range according to the amplitude coarse guidance information.

[0075] A second module is configured to determine a solution ambiguity baseline group according to the target antenna and a neighboring antenna of the target antenna.

[0076] A third module is configured to determine a phase difference ambiguity number value range of each baseline pair of the solution ambiguity baseline group according to the wave direction angle range in combination with a redundancy angle value.

[0077] A fourth module is configured to obtain a fuzzy multi-value solution set according to the phase difference ambiguity number value range; wherein the fuzzy multi-value solution set includes fuzzy multi-value solutions of each baseline pair of the solution ambiguity baseline group.

[0078] The fifth module is used for obtaining a solution ambiguity result through a clustering process according to the ambiguity multi-value solution set.

[0079] The contents of the method embodiments of the present application are applicable to the system embodiments of the present application, the system embodiments of the present application specifically implement the same functions as the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0080] Another aspect of the embodiment of the present application further provides an electronic device, comprising a processor and a memory;

[0081] The memory is used for storing a program.

[0082] The processor executes the program to realize the method as described above.

[0083] The contents of the method embodiments of the present application are applicable to the electronic device embodiments of the present application, the electronic device embodiments of the present application specifically implement the same functions as the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0084] Another aspect of the embodiment of the present application further provides a computer readable storage medium, the storage medium stores a program, and the program is executed by a processor to realize the method as described above.

[0085] The contents of the method embodiments of the present application are applicable to the computer readable storage medium embodiments of the present application, the computer readable storage medium embodiments of the present application specifically implement the same functions as the above method embodiments, and achieve the same beneficial effects as the above method embodiments.

[0086] The embodiment of the present application further discloses a computer program product or a computer program, the computer program product or the computer program comprises computer instructions stored in a computer readable storage medium. A processor of a computer device can read the computer instructions from the computer readable storage medium, and the processor executes the computer instructions, so that the computer device executes the method as described above.

[0087] The present application will be further described in detail below in combination with some specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application, and are not used to limit the present application.

[0088] First of all, it needs to be pointed out that the basic principle of direction finding of the two-dimensional phase interferometer method and the multi-value ambiguity problem refer to Figure 2 :

[0089] Let be the signal incident direction, and α and β are the azimuth angle and the elevation angle respectively, the origin O is the reference antenna, and point B i (xi , y i , z i ) is the spatial position of the i-th antenna. The straight line B i passes through point B i A is perpendicular to and intersects with at point A(x A , y A , z A ), OA is O, B i The wave path difference of the two antennas, the phase difference between the two antennas is The relationship between the phase difference and the wave path difference OA is:

[0090]

[0091] where λ is the wavelength of the incident signal. Thus we can get:

[0092]

[0093] And the other antenna can be solved together with α and β.

[0094] In the actual direction finding process, the phase difference range of the two antennas is given as (-π, π), and when the distance between the antennas is greater than the wavelength of the signal, the actual phase difference will differ from the given phase difference by 2π times k (k = 0, ±1, ±2), so there is a multi-value ambiguity problem.

[0095] Consider the case where a certain fixed antenna is not used as a reference antenna, place A, B, C three antennas in XOY plane, then the coordinates of the three antennas are A(x A , y A , 0), B(x B , y B , 0) and C(x C , y C , 0). Considering the phase difference multi-value ambiguity, the phase difference relationship between antennas is written as:

[0096]

[0097]

[0098]

[0099] where α and β are the azimuth and elevation angles of the incident signal respectively, k1 is the ambiguity number of the phase difference φ AB , k2 is the ambiguity number of the phase difference φ AC , k3 is the ambiguity number of the phase difference φ BC , all of which are integers.

[0100] Solving equations (1) and (2) together gives:

[0101]

[0102] Thus, the azimuth angle is obtained:

[0103]

[0104] Or,

[0105]

[0106] Thus, the elevation angle is obtained:

[0107]

[0108] In a wideband direction finding system, the wavelength of the high frequency signal is generally short, and in the actual system, the antenna spacing cannot be less than the signal wavelength due to the physical size of the antenna itself, and a smaller antenna spacing is also difficult to meet the requirement of direction finding accuracy.When the antenna spacing is greater than the signal half wavelength, k1, k2, and k3 are unknown numbers, and the shorter the incident signal wavelength, the more possible values of k1, k2, and k3, so in a wideband direction finding system, the value range of k1, k2, and k3 corresponding to different frequencies is different.

[0109] In view of the multi-value ambiguity problem, and in view of the lack of flexibility in the selection of the baseline pair in the stereo baseline method, and in view of the problem that when the directional antenna uniform circular array is used, part of the antennas receive weak signals, which leads to difficulty in ambiguity resolution, and the search range of the candidate angle is large, an ultra-wideband directional antenna uniform circular array ambiguity resolution method is provided in the embodiment of the application, and the process of resolving multi-value ambiguity of the embodiment of the application is as shown in the figure. Figure 3 According to the amplitude of the directional antenna, the range of the incoming wave direction and the quadrant to which the incoming wave direction belongs are preliminarily determined, so as to limit the range of k1, k2, and k3. All possible values of k1, k2, and k3 are reserved under the corresponding frequency, and any two of (1) to (3) are solved to obtain a group of solutions sinα, cosα, and cosβ containing ambiguous multi-values. Then, the method of simultaneously solving multiple groups of antennas and then comparing is used to resolve ambiguity. In the solutions of the multiple groups of antennas, the true value should be common to all groups of solutions and have a small difference, so the group of solutions that is common to all groups of solutions and has a small difference in each group of solutions is the true value.

[0110] After the range of the direction of arrival and the quadrant to which the direction of arrival belongs are preliminarily determined, since it is determined which quadrant the azimuth angle α belongs to, only the antenna with the largest amplitude and its adjacent three or more antennas are selected to form a baseline pair. Since the range of the direction of arrival is limited, the range of k1, k2 and k3 is limited, and the problem of azimuth angle mirror ambiguity is avoided, and only the problem of multi-value ambiguity needs to be solved. Moreover, the obtained angle candidate values only appear in the limited range of the direction of arrival, and the number of angle candidate values is reduced. Meanwhile, even if the error of the solution ambiguity occurs, the result is in the limited range of the direction of arrival.

[0111] The improvement process of the application includes five steps of limiting the range of the direction of arrival of the azimuth angle according to the amplitude coarse guidance information, determining the antenna group used for solution ambiguity, determining the maximum value and the minimum value of the ambiguity number, solving the azimuth angle and the elevation angle and clustering the solution to obtain the final result. In the step of solving the azimuth angle and the elevation angle, the method of the application does not need to solve the problem of azimuth angle mirror ambiguity.

[0112] The detailed steps are as follows:

[0113] 1. Limiting the range of the direction of arrival of the azimuth angle according to the amplitude coarse guidance information

[0114] First, according to the amplitude coarse guidance information, the range of the direction of arrival of the azimuth angle of the incident signal is determined as [α min ,α max ]. The antenna with the largest amplitude of the received signal is determined, and the intersection of the antenna with the largest amplitude of the received signal and its adjacent antennas is used to determine the range of the direction of arrival of the azimuth angle as [α min ,α max ]. Through the range, the range of the quadrant to which the direction of arrival of the azimuth angle belongs can be determined.

[0115] 2. Determining the baseline group used for solution ambiguity

[0116] Four groups of baseline pairs are used for solution ambiguity, which are [n-1, n, n+1], [n-2, n, n+1], [n-1, n, n+2] and [n-2, n, n+2], wherein n is the number of the antenna with the largest amplitude of the received signal. For the uniform circular array of directional antennas, when n=1, n-1 is recorded as N, and n-2 is recorded as N-1. When n=2, n-2 is recorded as N. When n=N-1, n+2 is recorded as 1. When n=N, n+1 is recorded as 1, and n+2 is recorded as 2, wherein the antennas of the uniform circular array of directional antennas are numbered from 1 to N in turn, and N represents the total number of antennas.

[0117] 3. Determining the maximum value and the minimum value of the phase difference ambiguity number

[0118] In the directional antenna array, since the gain of the antenna is close at the beam intersection point, and the influence of the noise is added, the maximum beam may be misjudged, so the determined beam range is added and subtracted αo improve reliability, α o is a redundant angle value, mainly to reduce the impact of the maximum amplitude beam caused by the noise when the incoming wave direction is near the intersection of two beams. It can be considered that the maximum amplitude beam error only exists near the intersection of two beams, and the setting of the redundant angle can still include the correct incoming wave direction angle in the search range when the maximum beam is wrong. Specifically, α o needs to be determined according to the angle range where the antenna gain at the intersection of two beams is similar. Therefore, the determined beam range is [α min -α o , α max +α o ].

[0119] Taking formula (1) as an example, the process of finding the maximum and minimum values of the phase difference ambiguity number is explained. The determined beam range [α min -α o , α max +α o ] is substituted into the formula. The maximum value of which can be calculated as (φ AB +2k1π) max , and the minimum value is (φ AB +2k1π) m,in , from which the maximum value of the ambiguity number k 1max =mod((φ AB +2k1π) max -φ AB ,2π) and the minimum value of the ambiguity number k 1min =mod((φ AB +2k1π) mmin -φ AB ,2π) can be calculated.

[0120] 4. Solve the azimuth angle and the elevation angle

[0121] Taking the simultaneous equations (1) and (2) as an example, the process of solving the candidate values is explained.

[0122] Firstly, through step 3, the maximum value k 1max and the minimum value k 1min of the ambiguity number in equation (1), and the maximum value k 2max and the minimum value k 2min of the ambiguity number in equation (2) can be calculated. By combining equations (1) and (2), we get:

[0123]

[0124] where k1∈[k 1min , k1min ],k2∈[k 2min k 2min After calculating tanα, combine it with the beam range [α] determined in step 3. min -α o α max +α o The quadrant range to which the azimuth angle α belongs is determined by the current fuzzy numbers k1 and k2. Since the range of azimuth angle α is (0, 2π), the tangent value of the angle in the first and third quadrants is positive, and the tangent value of the angle in the second and fourth quadrants is negative. Combining the initially determined quadrant range and the sign of tanα, the value of azimuth angle α can be solved, and then sinα and cosα can be calculated.

[0125] Then, through equation (1), we can obtain:

[0126]

[0127] Solve for the cosine value of the pitch angle β under the current ambiguity numbers k1 and k2, then remove the meaningless cosβ value introduced by the multivalues ​​(i.e. remove the value of cosβ>1) and its corresponding sinα and cosα, to obtain a set of solutions containing ambiguity multivalues ​​[sinαcosαcosβ].

[0128] 5. Solve for the final result using a clustering algorithm.

[0129] Select other baseline groups and repeat step 4 to obtain several sets of solutions containing fuzzy multivalued values. Then, obtain the final result through clustering. The clustering process can be described as follows: Select a set of solutions from the baseline group as a reference group, and select the first solution [sinα] from it. 11 cosα 11 cosβ 11 ], and then select the other 3 groups in sequence whose Euclidean distance is the same as the solution. The closest solution is obtained, and the corresponding Euclidean distance is calculated. The sum of the three closest Euclidean distances is taken as [sinα]. 11 cosα 11 cosβ 11 The clustering degree value of the reference group is calculated. The clustering degree values ​​of the remaining solutions in the reference group are calculated using the same method, and the solution with the lowest clustering degree value and its associated solutions in the other three groups are selected to form four approximate estimates. Finally, the arcsine and arccosine of these four estimates are calculated, and the average is taken to obtain the estimated value of the signal arrival angle.

[0130] In some embodiments, the simulation results of the above method are as follows:

[0131] Consider a nine-element directional antenna circular array, refer to Figure 4 The coordinates of the 9 antennas are as follows: (unit: mm).

[0132] Antenna 1 (6,-123.8); Antenna 2 (84.2,-91); Antenna 3 (123,-15.6); Antenna 4 (104.2,67.1); Antenna 5 (36.7,118.4); Antenna 6 (-48,114.3); Antenna 7 (-110.2,56.7); Antenna 8 (-120.9,-27.4); Antenna 9 (-75,98.7). The geometric model is shown in Figure 4

[0133] At frequencies 2GHz, 10GHz, 18GHz, respectively, at elevation angles 10°, 30° and 50°, the azimuth angle is simulated from 0° to 359°. Assuming the sampling frequency is 4 times the signal frequency, the signal is sampled for ten cycles. The phase difference information is obtained by sampling the signal, and 100 experiments are performed. The success probability of demodulation, the average error of the direction finding result and the root mean square error of the direction finding result of the stereoscopic baseline method and the method of the present application are counted.

[0134] Assuming the antenna pattern is a Gaussian pattern, the beam width is 80° at frequency 2GHz, and the beam width is 60° at frequencies 10GHz and 18GHz. The simulation is performed under the condition that the signal-to-noise ratio of the antenna with the maximum received signal amplitude is 10dB. When simulating the method of the present application, the range of the incoming wave azimuth angle is limited by the beam range of the antenna with the maximum received signal amplitude in the antenna array.

[0135] 1. Success probability of demodulation

[0136] As shown in Tables 1, 2 and 3, the success probability of demodulation of the method of the present application and the stereoscopic baseline method at different elevation angles at frequencies 2GHz, 10GHz and 18GHz, and signal-to-noise ratio 10dB, respectively.

[0137] Table 1

[0138]

[0139] Table 2

[0140]

[0141] Table 3

[0142]

[0143] ​From the simulation results of the successful deblurring probability, it can be seen that the deblurring probability decreases with the increase of the frequency, because with the increase of the frequency, the number of the blurring numbers increases, the beam width decreases, the gain of the adjacent antennas in the range of the maximum antenna of the received signal amplitude decreases, the signal-to-noise ratio decreases, and the phase difference measurement accuracy decreases. From the comparison of the successful deblurring probability, it can be seen that under the conditions of the frequency of 2GHz, 10GHz and 18GHz, the elevation angle of 50°, 30° and 10°, the successful deblurring probability of the method is higher than that of the stereoscopic baseline method. Since the baseline selected by the stereoscopic baseline method needs to satisfy the condition that the slope is positive and negative, the antennas used for direction finding are relatively fixed. When the signal arrives from the range outside the half-power points of the beam, the gain is small, the signal received by the antenna is weak, and due to the influence of the receiver noise, the signal-to-noise ratio is very small, which leads to difficulty in deblurring. The method of the application selects the antenna group with larger received signal amplitude to form a baseline group through amplitude coarse guidance, and the successful deblurring probability is larger. Since the deblurring probability of the stereoscopic baseline method is lower, the direction finding error is larger, and only the direction finding error of the method of the application is simulated and analyzed.

[0144] 2. Direction finding error

[0145] As shown in Tables 4, 5 and 6, the azimuth angle and elevation angle direction finding errors of the method of the application at different elevation angles when the signal-to-noise ratio is 10dB at frequencies of 2GHz, 10GHz and 18GHz are shown.

[0146] Table 4

[0147]

[0148] Table 5

[0149]

[0150] Table 6

[0151]

[0152] From the above simulation results, it can be seen that the method of the application can effectively improve the successful deblurring probability and the direction finding accuracy by limiting the range of the incoming wave azimuth angle through the amplitude coarse guidance information. Compared with the fixed baseline group selected by the stereoscopic baseline method, the method of the application selects the antenna group composed of the antenna with the largest received amplitude and its four adjacent antennas, the signal-to-noise ratio of the received signal of the selected antenna is larger, and more accurate phase difference information can be obtained, so the successful deblurring probability and the direction finding accuracy are higher.

[0153] In summary, aiming at the problems of lack of flexibility in selection of baseline pairs in stereo baseline method, difficulty in unblurring due to weak signals received by part of antennas when used on directional antenna uniform circular array, and large search range of candidate angle, the present application improves the stereo baseline method. The directivity of the antenna is used to participate in the unblurring process, so that in the selection of baseline pairs, it is not necessary to satisfy that the slopes of the two baselines are opposite in sign, but only to select the antenna with the largest amplitude and 3 or more antennas near it, and to form baseline pairs two by two to unblur. And by limiting the search range of the direction of arrival through the antenna with the largest amplitude, the number of candidate angle values is reduced, thereby reducing the calculation amount and error of the subsequent comparison process. The beneficial effects of the embodiments of the present application include:

[0154] First, in the selection of baseline pairs, the present patent method does not need to select two baselines with opposite slopes, but only needs to select the antenna with the largest amplitude and several antennas adjacent to it to form a baseline, so the selection of baseline pairs is more flexible. Before solving the angle, according to the amplitude of the directional antenna, the range of the direction of arrival and the quadrant to which it belongs are determined, thereby limiting the range of k1, k2, and k3, and avoiding the problem of azimuth angle mirror unblurring, and only the problem of multi-value unblurring needs to be solved.

[0155] Second, in the search range of the final result, the directivity of the antenna is used to preliminarily determine the range of the azimuth angle of the received signal, thereby inferring the range of the ambiguity number of the measured phase difference, and the search range of the direction of arrival is reduced to (360° / N+2×α o ) / 360°, where N is the number of elements of the uniform circular array antenna, the number of candidate angle values is reduced, and since the range of the direction of arrival has been limited, the calculation result is also within this range, even if the calculation result is wrong, the error is also reduced.

[0156] In some alternative embodiments, the functions / operations mentioned in the block diagram can not occur in the order mentioned in the operation diagram. For example, depending on the functions / operations involved, two blocks shown in succession can actually be performed substantially simultaneously or the blocks can sometimes be performed in reverse order. In addition, the embodiments presented and described in the flowcharts of the present application are provided by way of example, with the purpose of providing a more comprehensive understanding of the technology. The disclosed method is not limited to the operations and logical flows presented herein. Alternative embodiments are contemplated in which the order of various operations is changed and in which sub-operations described as part of larger operations are independently performed.

[0157] Furthermore, although the present application is described in the context of functional modules, it is to be understood that one or more of the described functions and / or features can be integrated in a single physical device and / or software module, or one or more functions and / or features can be implemented in separate physical devices or software modules. It will also be appreciated that detailed discussion of the actual implementation of each module is not necessary to an understanding of the application. Rather, the actual implementation of the modules, in combination with their attributes, functions, and internal relationships, are to be taken as being readily understood by one skilled in the art in view of the property disclosed herein. Accordingly, the present application is not limited to the embodiments illustrated in the figures. Rather, one of ordinary skill in the art will be able to make and use the present application without undue experimentation based on the detailed description of the drawings in combination with the property disclosed herein. It is also to be understood that the specific concepts disclosed are merely illustrative and that the scope of the present application is not limited to the specific concepts presented. The scope of the present application is to be limited only by the claims set forth below together with their equivalents.

[0158] If the functions are implemented in the form of software function units and sold or used as independent products, they can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or the part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0159] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a list of executable instructions for implementing logic functions, and can be specifically embodied in any computer readable medium for use by an instruction execution device, device or equipment (such as a computer-based device, a device including a processor, or other devices that can fetch and execute instructions from an instruction execution device, device or equipment), or in conjunction with these instruction execution devices, devices or equipment. For the purpose of this specification, "computer readable medium" can be any device that can contain, store, communicate, propagate or transport programs for use by an instruction execution device, device or equipment or in conjunction with these instruction execution devices, devices or equipment.

[0160] More specific examples (a non-exhaustive list) of the computer-readable medium include the following: an electrical connection (electronic) having one or more wires, a portable computer diskette (magnetic), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, and a portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium can also be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, for example, via optical scanning of the paper or other medium, then compiled, interpreted, or otherwise processed in a suitable manner, if necessary, and then stored in a computer memory.

[0161] It should be understood that aspects of the application can be implemented in hardware, software, firmware or combinations thereof. In the embodiments described above, various steps or methods can be implemented, for example, by software or firmware stored in memory and executed by a suitable instruction execution device. If implemented in hardware, for example, in another embodiment, any of the following techniques can be used in whole or in part: discrete logic circuitry having logic gates for implementing logic functions upon data signals, application specific integrated circuits having appropriate combinational logic gates, programmable gate arrays (PGA), field programmable gate arrays (FPGA), and so forth, as well as or instead.

[0162] In the description of the specification, the description of the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" and the like means that the specific feature, structure, material or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the application. In the specification, illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any suitable manner in one or more embodiments or examples.

[0163] Although embodiments of the application have been shown and described, it would be recognized by those of ordinary skill in the art that various changes, modifications, substitutions, and alterations can be made hereto without departing from the spirit and scope of the application. The scope of the application is limited only by the claims and the equivalents thereof.

[0164] The above is a specific description of the preferred embodiments of the present application, but the present application is not limited to the described embodiments, and those skilled in the art can make various equivalent modifications or replacements without departing from the spirit of the present application, and these equivalent modifications or replacements are all included in the scope defined by the claims of the present application.

Claims

1. A method for deambiguing a uniform circular array of ultra-wideband directional antennas, characterized in that, include: Obtain coarse amplitude guidance information, and determine the target antenna and the range of the azimuth angle of arrival based on the coarse amplitude guidance information; Determine the unambiguous baseline group based on the target antenna and its adjacent antennas; Based on the range of the azimuth angle of the incoming wave, and in combination with the redundant angle value, the range of the phase difference ambiguity number values ​​for each pair of baselines in the unambiguous baseline group is determined. Based on the range of the phase difference ambiguity number, a fuzzy multivalued solution set is obtained; wherein, the fuzzy multivalued solution set includes the fuzzy multivalued solutions of each baseline pair of the unfuzzy baseline group; The step of obtaining the fuzzy multi-valued solution set based on the range of values ​​of the phase difference fuzzy number includes: Based on the range of the phase difference ambiguity number, determine the ambiguity solution of the incident signal azimuth tangent and the ambiguity solution of the incident signal elevation cosine. Based on the tangent of the incident signal azimuth angle, and combined with the quadrant range to which the beam range belongs, determine the fuzzy solution of the sine value and the fuzzy solution of the cosine value of the incident signal azimuth angle. The fuzzy multivalued solution includes the fuzzy solution of the incident signal azimuth angle sine value, the fuzzy solution of the incident signal azimuth angle cosine value, and the fuzzy solution of the incident signal elevation angle cosine value. The beam range is determined by combining the incoming wave azimuth angle range with the redundant angle value. Based on the fuzzy multi-valued solution set, the defuzzification result is obtained through clustering.

2. The method for deambiguation of a uniform circular array of ultra-wideband directional antennas according to claim 1, characterized in that, The step of acquiring coarse amplitude guidance information and determining the target antenna and the range of the azimuth angle of arrival based on the coarse amplitude guidance information includes: Obtain coarse amplitude guidance information, and determine the target antenna based on the coarse amplitude guidance information; The azimuth range of the incoming wave is determined based on the beam intersection point of the target antenna and the adjacent antennas of the target antenna.

3. The method for deambiguation of a uniform circular array of ultra-wideband directional antennas according to claim 1, characterized in that, The step of determining the unambiguous baseline group based on the target antenna and its adjacent antennas includes: Based on the target antenna and the four antennas adjacent to the target antenna, four sets of de-ambiguous baseline groups are formed.

4. The method for deambiguation of a uniform circular array of ultra-wideband directional antennas according to claim 1, characterized in that, The step of determining the range of phase difference ambiguity numbers for each baseline pair in the unambiguous baseline group based on the range of the incoming wave azimuth angle and the redundant angle value includes: Based on the range of the incoming wave azimuth angle, the beam range of each baseline pair in the deblurred baseline group is determined in combination with the redundant angle value. Based on the beam range, determine the range of phase difference ambiguity numbers for each baseline pair in the deambiguity baseline group.

5. The method for deambiguing a uniform circular array of ultra-wideband directional antennas according to claim 1, characterized in that, The step of obtaining the defuzzification result through clustering based on the fuzzy multi-valued solution set includes: Based on the fuzzy multivalued solution set, select a set of fuzzy multivalued solutions as a reference set; Based on the fuzzy multivalued solutions of the fuzzy multivalued solution set excluding the reference group, the target estimated value of the fuzzy multivalued solution set is determined by calculating the Euclidean distance with the reference group. The defuzzing result is obtained by averaging the estimated target value.

6. A method for deambiguation of a uniform circular array of ultra-wideband directional antennas according to any one of claims 1 to 5, characterized in that, Also includes: Simulation direction finding is performed based on preset conditions to obtain unfuzzy simulation results; The unambiguity simulation results include the probability of successful unambiguity resolution and the direction-finding error.

7. A deambiguity resolution system for a uniform circular array of ultra-wideband directional antennas, characterized in that, include: The first module is used to acquire coarse amplitude guidance information and determine the target antenna and the range of the azimuth angle of arrival based on the coarse amplitude guidance information. The second module is used to determine a deambiguity baseline group based on the target antenna and its adjacent antennas; The third module is used to determine the range of phase difference ambiguity numbers for each pair of baselines in the unambiguous baseline group based on the range of the incoming wave azimuth angle and the redundant angle value. The fourth module is used to obtain a fuzzy multivalued solution set based on the range of the phase difference fuzziness number; wherein, the fuzzy multivalued solution set includes the fuzzy multivalued solutions of each pair of baselines in the unfuzzy baseline group; The step of obtaining the fuzzy multi-valued solution set based on the range of values ​​of the phase difference fuzzy number includes: Based on the range of the phase difference ambiguity number, determine the ambiguity solution of the incident signal azimuth tangent and the ambiguity solution of the incident signal elevation cosine. Based on the tangent of the incident signal azimuth angle, and combined with the quadrant range to which the beam range belongs, determine the fuzzy solution of the sine value and the fuzzy solution of the cosine value of the incident signal azimuth angle. The fuzzy multivalued solution includes the fuzzy solution of the incident signal azimuth angle sine value, the fuzzy solution of the incident signal azimuth angle cosine value, and the fuzzy solution of the incident signal elevation angle cosine value. The beam range is determined by combining the incoming wave azimuth angle range with the redundant angle value. The fifth module is used to obtain the defuzzification result by clustering based on the fuzzy multi-valued solution set.

8. An electronic device, characterized in that, Including the processor and memory; The memory is used to store programs; The processor executes the program to implement the method as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The storage medium stores a program that is executed by a processor to implement the method as described in any one of claims 1 to 6.

Images