CN108848446B - Data processing method and device and electronic equipment - Google Patents
Data processing method and device and electronic equipment Download PDFInfo
- Publication number
- CN108848446B CN108848446B CN201810647566.1A CN201810647566A CN108848446B CN 108848446 B CN108848446 B CN 108848446B CN 201810647566 A CN201810647566 A CN 201810647566A CN 108848446 B CN108848446 B CN 108848446B
- Authority
- CN
- China
- Prior art keywords
- phase
- region
- signal
- signal source
- determining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/023—Services making use of location information using mutual or relative location information between multiple location based services [LBS] targets or of distance thresholds
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
Abstract
The invention belongs to the technical field of low-altitude security and protection, and relates to a data processing method and device and electronic equipment. The method comprises the following steps: preprocessing signals received by each directional antenna to obtain the frequency spectrum amplitude and phase of each signal; determining a phase region of a signal source based on the spectrum amplitude of each path of signal and a joint detection rule; carrying out differential processing on the phases of any two adjacent paths of signals; performing sum operation on the differential phases to obtain a phase difference summation result; extending the phase difference summation result to a range specified by the determined phase region; and obtaining the incoming wave direction of the signal source based on the extended phase and a preset phase address table. The phase extension ambiguity resolution method based on the amplitude information can effectively reduce the implementation cost of the ambiguity resolution algorithm in the traditional direction finding system; meanwhile, the method can smooth the phase cycle jump, so that the method can adapt to a phase-address mapping method.
Description
Technical Field
The invention belongs to the technical field of low-altitude security and particularly relates to a data processing method and device and electronic equipment.
Background
Along with the rapid development of civilian unmanned aerial vehicle industry, consumption level unmanned aerial vehicle function is powerful day by day, and the unmanned aerial vehicle quantity that flows into market constantly increases, and this means unmanned aerial vehicle simultaneously and has brought a large amount of hidden dangers in the aspect of public safety and secret privacy. Therefore, the phase direction finding method with the high-precision characteristic is gradually developed from the military field to the unmanned aerial vehicle security field. Unmanned aerial vehicle's functions such as control, location and picture pass all need realize through radio signal. According to the phase direction finding method in the low-altitude field, the phase difference of corresponding signals is measured by accumulating received radio signals, and the incoming wave direction of the signals is solved by utilizing the phase difference, namely the direction of the unmanned aerial vehicle is remotely controlled or the unmanned aerial vehicle is in the same direction. Therefore, a series of countermeasures such as tracking, identification, interference and the like can be further implemented. For example, in a phase direction finding method represented by an interferometer, as shown in fig. 1, a signal direction Φ is estimated from a phase difference introduced by a path difference d between two antennas. If the signal wavelength λ is the phase difference
Because of the periodicity of the phase, when d > lambda exists, the mapping of the phase difference and the wave path difference is no longer unique, namely, phase ambiguity exists.
Disclosure of Invention
In view of the above, the present invention provides a data processing method, an apparatus and an electronic device to effectively solve the above problem.
The embodiment of the invention is realized by the following steps:
in a first aspect, an embodiment of the present invention provides a data processing method, which is applied to an antenna array including at least three directional antennas, where the method includes: preprocessing signals received by each directional antenna to obtain the frequency spectrum amplitude and phase of each signal; determining a phase region of a signal source based on the spectrum amplitude of each path of signal and a joint detection rule; carrying out differential processing on the phases of any two adjacent paths of signals; performing sum operation on the differential phases to obtain a phase difference summation result; extending the phase difference summation result to a range specified by the determined phase region; and obtaining the incoming wave direction of the signal source based on the extended phase and a preset phase address table.
With reference to the first implementation manner of the first aspect, preprocessing signals received by each directional antenna to obtain a spectral amplitude and a phase of each signal, includes: performing windowing Fourier transform on signals received by each directional antenna, and converting each time domain signal into a frequency domain signal; and carrying out complex finger conversion on each path of frequency domain signal to obtain the frequency spectrum amplitude and the phase of each path of signal.
With reference to the second implementation manner of the first aspect, performing windowed fourier transform on each path of signal, and converting each path of time domain signal into a frequency domain signal includes: carrying out N-point windowing on each path of signal, wherein N is the number of data points processed at one time; and carrying out Fourier transform on each channel of signals after windowing, and converting each channel of time domain signals into frequency domain signals.
With reference to the third implementation manner of the first aspect, the spectral amplitude of each path of signal includes: the frequency spectrum amplitudes of the three signals are respectively M1, M2 and M3; determining the phase region of the signal source based on the spectrum amplitude of each path of signal and a joint detection rule, comprising: determining sizes of the M1, the M2, and the M3 with respect to each other; when the M1 is maximum, determining that the phase area of the signal source is in a first area, wherein the first area is [ -4 pi, -2 pi ]; determining that the phase region of the signal source is in a fourth region when the M3 is maximum, the fourth region being [2 pi, 4 pi ] is maximum at the M2, and the M1 is greater than the M3, the phase region of the signal source is in a second region, the second region being [ -2 pi, 0 ]; when the M2 is maximum and the M1 is less than the M3, determining that the phase region of the signal source is in a third region, wherein the third region is [0, 2 pi ]; when the M2 is the largest and the M1 is equal to the M3, determining that the phase region of the signal source is in a fifth region, the fifth region being [ - π, π ].
With reference to the fourth implementation manner of the first aspect, obtaining the incoming wave direction of the signal source based on the extended phase and the preset phase address table includes: multiplying the extended phase by a frequency point coefficient corresponding to the Fourier transform to obtain a corrected phase; and obtaining the incoming wave direction of the signal source based on the corrected phase and a preset phase address table.
In a second aspect, an embodiment of the present invention further provides a data processing apparatus, including: the preprocessing module is used for preprocessing the signals received by each directional antenna to obtain the frequency spectrum amplitude and the phase of each signal; the phase region determining module is used for determining the phase region of the signal source based on the spectrum amplitude of each path of signal and a joint detection rule; the differential processing module is used for carrying out differential processing on the phases of any two adjacent paths of signals; the difference summing module is used for carrying out sum operation on the phase after difference to obtain a phase difference summing result; the phase extension module is used for extending the phase difference summation result to a range appointed by the determined phase region; and the incoming wave direction obtaining module is used for obtaining the incoming wave direction of the signal source based on the extended phase and the preset phase address table.
In combination with the first embodiment of the second aspect, the preprocessing module includes: the windowing Fourier transform unit is used for carrying out windowing Fourier transform on the signals received by each directional antenna and converting each time domain signal into a frequency domain signal; and the complex finger conversion unit is used for carrying out complex finger conversion on each path of frequency domain signal to obtain the frequency spectrum amplitude and the phase of each path of signal.
With reference to the second implementation manner of the second aspect, the windowed fourier transform unit includes: the N-point windowing subunit is used for carrying out N-point windowing on each path of signal, wherein N is the number of data points processed at one time; and the Fourier transform subunit is used for performing Fourier transform on each channel of signals subjected to windowing and converting each channel of time domain signals into frequency domain signals.
With reference to the third implementation manner of the second aspect, the spectral amplitude of each path of signal includes: the frequency spectrum amplitudes of the three signals are respectively M1, M2 and M3; the phase region determination module includes: a size determination unit for determining sizes of the M1, the M2, and the M3 with respect to each other; a first region determining unit, configured to determine that the phase region of the signal source is in a first region when M1 is maximum, where the first region is [ -4 pi, -2 pi ]; a fourth area determination unit, configured to determine that the phase area of the signal source is in a fourth area when M3 is maximum, where the fourth area is [2 pi, 4 pi ]; a second region determining unit for determining that the phase region of the signal source is in a second region when the M2 is maximum and the M1 is greater than the M3, the second region being [ -2 pi, 0 ]; a third region determining unit, configured to determine that the phase region of the signal source is in a third region when the M2 is maximum and the M1 is smaller than the M3, where the third region is [0, 2 pi ]; a fifth region determining unit for determining that the phase region of the signal source is in a fifth region when the M2 is maximum and the M1 is equal to the M3, the fifth region being [ -pi, pi ].
With reference to the fourth implementation manner of the second aspect, the incoming wave direction obtaining module includes: the phase correction unit is used for multiplying the extended phase by a frequency point coefficient corresponding to the Fourier transform to obtain a corrected phase; and the incoming wave direction obtaining unit is used for obtaining the incoming wave direction of the signal source based on the corrected phase and a preset phase address table.
In a third aspect, an embodiment of the present invention further provides an electronic device, including: a memory and a processor; the memory is used for storing programs; the processor is configured to call a program stored in the memory to perform the method of the first aspect.
The data processing method provided by the embodiment of the invention is applied to an antenna array comprising at least three directional antennas, and the method is used for preprocessing signals received by each directional antenna to obtain the frequency spectrum amplitude and the phase of each path of signal; then, determining a phase region of a signal source based on the spectrum amplitude of each path of signal and a joint detection rule; then, carrying out differential processing on the phases of any two adjacent paths of signals; then, carrying out sum operation on the differential phases to obtain a phase difference summation result; then extending the phase difference summation result to a range specified by the determined phase region; and finally, obtaining the incoming wave direction of the signal source based on the extended phase and a preset phase address table. The phase extension ambiguity resolution method based on the amplitude information can effectively reduce the implementation cost of the ambiguity resolution algorithm in the traditional direction finding system; meanwhile, the method can smooth the phase cycle jump, so that the method can adapt to a phase-address mapping method.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by the practice of the embodiments of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts. The above and other objects, features and advantages of the present invention will become more apparent from the accompanying drawings. Like reference numerals refer to like parts throughout the drawings. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention.
Fig. 1 shows a signal reception diagram of a conventional antenna array.
Fig. 2 shows a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Fig. 3 is a schematic diagram illustrating an antenna array provided by an embodiment of the present invention.
Fig. 4 is a flowchart illustrating a data processing method according to an embodiment of the present invention.
Fig. 5 is a schematic processing flow diagram illustrating a data processing method according to an embodiment of the present invention.
Fig. 6 shows a block diagram of a data processing apparatus according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
In the description of the present invention, it should be noted that the terms "first", "second", "third", and the like are used only for distinguishing the description, and are not intended to indicate or imply relative importance.
The process of the present application is being studied by the researchersThe current three direction-finding methods for incoming wave direction all have certain defects, for example, the first phase direction-finding method represented by interferometer estimates the signal direction phi according to the phase difference introduced by the wave path difference d between two antennas, as shown in fig. 1. If the signal wavelength λ is the phase difference
Because of the periodicity of the phase, when d > lambda exists, the mapping of the phase difference and the wave path difference is no longer unique, namely, phase ambiguity exists. And secondly, the mutual correlation operation is carried out on the obtained phase and the stored phase table to obtain an autocorrelation spectrum, then the peak value of the autocorrelation spectrum is searched to obtain the direction of the signal, and the scheme has higher precision. And thirdly, directly solving according to a signal model by using a direct solving method based on an inverse trigonometric function.
In the first method, since there is a phase ambiguity problem, in order to solve the ambiguity problem caused by the excessively long baseline, a deblurring process is required. The phase ambiguity resolution is accomplished by finding the phase difference of the received data of the antennas on a plurality of baselines of different lengths. The phase difference of multiple base lines needs more antenna array elements and more received data, and therefore the implementation cost and the operation time cost of the scheme are increased. The physical errors of the antenna can be accumulated by calculating the difference of the phase values corresponding to the multilevel base lines, the higher the signal frequency is, the shorter the wavelength is, the larger the influence of the received errors is, the higher the requirement on the antenna arrangement is provided, and the method is also one of important reasons for limiting the wide application of the interferometer algorithm in the civil direction finding system.
In the second method, a large number of multiplication operations are required by using a cross-correlation algorithm, and the search for the peak value of the autocorrelation spectrum can only be performed sequentially, which consumes a large amount of system resources and operation time. If the scan length L is 120, at least 240 multiplications and 120 additions and 120 comparison operations are required for L2. The realization cost and the operation time of this scheme are all longer, are difficult to in low cost, have the high real-time in the system application simultaneously, can not adapt to the application occasion to unmanned aerial vehicle.
In the third method, the direct solving method requires at least 1 division, 2 multiplications and two inverse trigonometric function operations, wherein the operations of the division and the inverse trigonometric function require at least 20 iterations of addition. The realization cost and the operation time of this scheme are all longer, are difficult to in low cost, have the high real-time in the system application simultaneously, can not adapt to the application occasion to unmanned aerial vehicle.
It should be noted that the defects existing in the above solutions are the results obtained after the inventor has practiced and studied carefully, and therefore, the discovery process of the above problems and the solutions proposed by the following embodiments of the present invention to the above problems should be the contribution of the inventor to the present invention in the process of the present invention.
As shown in fig. 2, fig. 2 is a block diagram illustrating a structure of an electronic device 100 according to an embodiment of the present invention. The electronic device 100 includes: data processing device 110, memory 120, memory controller 130, and processor 140.
The memory 120, the memory controller 130, and the processor 140 are electrically connected to each other directly or indirectly to realize data transmission or interaction. For example, the components may be electrically connected to each other via one or more communication buses or signal lines. The data processing device 110 includes at least one software function module which can be stored in the memory 120 in the form of software or firmware (firmware) or solidified in an Operating System (OS) of the electronic device 100. The processor 140 is used to execute executable modules stored in the memory 120, such as software functional modules or computer programs included in the data processing apparatus 110.
The Memory 120 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like. The memory 120 is configured to store a program, and the processor 140 executes the program after receiving an execution instruction, and a method executed by the electronic device 100 defined by a flow disclosed in any embodiment of the invention described later may be applied to the processor 140, or implemented by the processor 140.
The processor 140 may be an integrated circuit chip having signal processing capabilities. The processor may be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like; but may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components. The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
Please refer to fig. 3, which is a schematic diagram of an antenna array applicable to the data processing method applied to the electronic device 100 according to an embodiment of the present invention. The antenna array comprises at least three directional antennas, and assuming that the 3dB beam width of each directional antenna is phi, the angle interval between each antenna needs to satisfy theta < phi/3, and the effective direction-finding range is [ -theta, theta ]. Fig. 4 is a schematic diagram of the antenna array when θ is 30 °.
Referring to fig. 4, a data processing method applied to the electronic device 100 according to an embodiment of the present invention is provided. The steps involved will be described below in conjunction with fig. 4.
Step S101: and preprocessing the signals received by each directional antenna to obtain the frequency spectrum amplitude and phase of each path of signal.
For example, the signal received by the antenna 1 is preprocessed to obtain the spectral amplitude and phase of the path of signal, the signal received by the antenna 2 is preprocessed to obtain the spectral amplitude and phase of the path of signal, and the signal received by the antenna 3 is preprocessed to obtain the spectral amplitude and phase of the path of signal.
As an optional implementation manner, the above process includes performing windowed fourier transform on signals received by each directional antenna, and converting each time domain signal into a frequency domain signal; and carrying out complex finger conversion on each path of frequency domain signal to obtain the frequency spectrum amplitude and the phase of each path of signal. Wherein, carrying out windowing Fourier transform on each path of signals comprises: carrying out N-point windowing on each path of signal, wherein N is the number of data points processed at one time; and carrying out Fourier transform on each channel of signals after windowing.
Step S102: and determining the phase region of the signal source based on the spectrum amplitude of each path of signal and a joint detection rule.
And after the frequency spectrum amplitude of each path of signal is obtained, determining the phase region of the signal source based on the frequency spectrum amplitude of each path of signal and a joint detection rule. Wherein, the spectrum amplitude of each path of signal comprises: the spectral amplitudes of the three signals (signal received by antenna 1, signal received by antenna 2, signal received by antenna 3) are M1, M2, and M3, respectively. Then the process includes determining the size of the M1, the M2, and the M3 relative to each other; when the M1 is maximum, determining that the phase area of the signal source is in a first area, wherein the first area is [ -4 pi, -2 pi ]; determining that the phase region of the signal source is in a fourth region when the M3 is maximum, the fourth region being [2 pi, 4 pi ] is maximum at the M2, and the M1 is greater than the M3, the phase region of the signal source is in a second region, the second region being [ -2 pi, 0 ]; when the M2 is maximum and the M1 is less than the M3, determining that the phase region of the signal source is in a third region, wherein the third region is [0, 2 pi ]; when the M2 is the largest and the M1 is equal to the M3, determining that the phase region of the signal source is in a fifth region, the fifth region being [ - π, π ].
Note that M1 is the spectral amplitude of the signal received by the antenna 1, M2 is the spectral amplitude of the signal received by the antenna 2, and M3 is the spectral amplitude of the signal received by the antenna 3.
Wherein, the receiving gain of the directional antenna is influenced by the direction of the incoming wave. If the direct radiation gain is G0Then the signal with the angle of incidence phi (phi epsilon-theta, theta)]) The gain value is approximately equal to
In fig. 3, the reception gain of the three antennas can be approximated as:
therefore, the incident signal can be pre-estimated according to the amplitudes of the three antennas
And
these 5 regions. The phase values assigned to the 5 regions are extended by 2 pi periods, and are respectively constrained to-4 pi, -2 pi]，[－2π，0]，[0，2π]，[2π，4π]And [ - π, π]In each range, the phase value has uniqueness, and the anti-fuzzy performance is 4 times higher than that of a common algorithm under the same condition.
Step S103: and carrying out differential processing on the phases of any two adjacent paths of signals.
After the spectrum amplitude of each path of signal is obtained, the phase of any two adjacent paths of signals is subjected to differential processing. Wherein, the phase place of each said way signal includes: the phases of the three signals (signal received by antenna 1, signal received by antenna 2, signal received by antenna 3) are P1, P2, and P3, respectively. The above procedure includes differential processing of the phase of the signal received by antenna 1 and the phase of the signal received by antenna 2, P1-P2; and the phases of the signals received by the antenna 2 and the phases of the signals received by the antenna 3 are differentially processed P2-P3.
Note that P1 is the phase of the signal received by the antenna 1, P2 is the phase of the signal received by the antenna 2, and P3 is the phase of the signal received by the antenna 3.
Step S104: and carrying out sum operation on the differential phases to obtain a phase difference summation result.
And carrying out sum operation on the differential phases to obtain a phase difference summation result. For example, P1-P2 + P2-P3.
Step S105: and extending the phase difference summation result to the range specified by the determined phase region.
After the phase difference summation results (P1-P2 + P2-P3) are obtained and the phase region of the signal source is determined, the phase difference summation results are extended to the range specified by the determined phase region. For example, the phase difference summation result is extended to the range specified by the fifth region.
Step S106: and obtaining the incoming wave direction of the signal source based on the extended phase and a preset phase address table.
After the extended phase is obtained, the preset phase address table is searched to see the area of the extended phase in the preset phase address table, and then the incoming wave direction of the signal source can be obtained based on the address (the incoming wave direction of the signal source) corresponding to the area.
As an optional implementation manner, after obtaining the extended phase, multiplying the extended phase by a frequency point coefficient corresponding to the fourier transform to correct the phase, so as to obtain a corrected phase; and obtaining the incoming wave direction of the signal source based on the corrected phase and a preset phase address table.
For ease of understanding, the above process may be briefly described with reference to the process flow diagram shown in fig. 5.
(1) Carrying out N-point windowing on signals received by channels corresponding to the antennas;
(2) performing N-point FFT on each path of data, and converting a time domain signal into a frequency domain;
(3) performing complex finger conversion on all complex frequency spectrums, and performing complex finger conversion on the amplitude M of the frequency spectrums1，2，3And phase P1，2，3Carrying out separation;
(4) and carrying out joint detection on the amplitude:
wherein M is1Maximum is in the first area;
M3maximum is in the second area;
M2maximum, and M1＞M3Then in the second area;
M2maximum, and M1＜M3Then in the third area;
in particular, for complete coverage, M2Maximum and M1＝M3Then, in the fifth area;
(5) carrying out differential processing on adjacent channel phases;
(6) carrying out summation operation on the differential phase;
(7) according to the amplitude prejudging result in the step (4), extending the phase after the summation operation to a specified range;
(8) multiplying the extended phase difference by the corresponding frequency point coefficient, and converting the phase difference to a mapping address;
(9) and obtaining a lookup table value corresponding to the mapping address, namely the signal direction.
The present embodiment also provides a data processing apparatus 110, as shown in fig. 6. The data processing apparatus 110 includes: a preprocessing module 111, a phase region determining module 112, a difference processing module 113, a difference summing module 114, a phase extension module 115, and an incoming wave direction obtaining module 116.
The preprocessing module 111 is configured to preprocess the signals received by the directional antennas to obtain a spectrum amplitude and a phase of each path of signal.
And a phase region determining module 112, configured to determine a phase region of the signal source based on the spectral amplitude of each channel of signal and a joint detection rule.
And the difference processing module 113 is configured to perform difference processing on the phases of any two adjacent signals.
And the difference summing module 114 is configured to sum the differential phases to obtain a phase difference summing result.
And a phase extension module 115, configured to extend the phase difference summation result to a range specified by the determined phase region.
And an incoming wave direction obtaining module 116, configured to obtain an incoming wave direction of the signal source based on the extended phase and a preset phase address table.
It should be noted that, in the present specification, the embodiments are all described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments may be referred to each other.
The data processing apparatus 110 according to the embodiment of the present invention has the same implementation principle and technical effect as the foregoing method embodiments, and for the sake of brief description, reference may be made to the corresponding contents in the foregoing method embodiments for the parts of the apparatus embodiments that are not mentioned.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented in other ways. The apparatus embodiments described above are merely illustrative, and for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
In addition, the functional modules in the embodiments of the present invention may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.
It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A data processing method applied to an antenna array comprising at least three directional antennas, the method comprising:
preprocessing signals received by each directional antenna to obtain the frequency spectrum amplitude and phase of each signal;
determining a phase region of a signal source based on the spectrum amplitude of each path of signal and a joint detection rule;
carrying out differential processing on the phases of any two adjacent paths of signals;
performing sum operation on the differential phases to obtain a phase difference summation result;
extending the phase difference summation result to a range specified by the determined phase region;
and obtaining the incoming wave direction of the signal source based on the extended phase and a preset phase address table.
2. The method of claim 1, wherein preprocessing the signals received by each directional antenna to obtain the spectral amplitude and phase of each signal comprises:
performing windowing Fourier transform on signals received by each directional antenna, and converting each time domain signal into a frequency domain signal;
and carrying out complex finger conversion on each path of frequency domain signal to obtain the frequency spectrum amplitude and the phase of each path of signal.
3. The method of claim 2, wherein performing a windowed fourier transform on each of the plurality of signals to convert each of the plurality of time domain signals into a frequency domain signal comprises:
carrying out N-point windowing on each path of signal, wherein N is the number of data points processed at one time;
and carrying out Fourier transform on each channel of signals after windowing, and converting each channel of time domain signals into frequency domain signals.
4. The method according to any one of claims 1-3, wherein the spectral magnitudes of the signals comprise: the frequency spectrum amplitudes of the three signals are respectively M1, M2 and M3; determining the phase region of the signal source based on the spectrum amplitude of each path of signal and a joint detection rule, comprising:
determining sizes of the M1, the M2, and the M3 with respect to each other;
when the M1 is maximum, determining that the phase area of the signal source is in a first area, wherein the first area is [ -4 pi, -2 pi ];
when the M3 is maximum, determining that the phase region of the signal source is in a fourth region, wherein the fourth region is [2 pi, 4 pi ]
When the M2 is maximum and the M1 is greater than the M3, determining that the phase region of the signal source is in a second region, the second region being [ -2 π, 0 ];
when the M2 is maximum and the M1 is less than the M3, determining that the phase region of the signal source is in a third region, wherein the third region is [0, 2 pi ];
when the M2 is the largest and the M1 is equal to the M3, determining that the phase region of the signal source is in a fifth region, the fifth region being [ - π, π ].
5. The method of claim 3, wherein obtaining the incoming wave direction of the signal source based on the extended phase and the preset phase address table comprises:
multiplying the extended phase by a frequency point coefficient corresponding to the Fourier transform to obtain a corrected phase;
and obtaining the incoming wave direction of the signal source based on the corrected phase and a preset phase address table.
6. A data processing apparatus, comprising:
the preprocessing module is used for preprocessing the signals received by each directional antenna to obtain the frequency spectrum amplitude and the phase of each signal;
the phase region determining module is used for determining the phase region of the signal source based on the spectrum amplitude of each path of signal and a joint detection rule;
the differential processing module is used for carrying out differential processing on the phases of any two adjacent paths of signals;
the difference summing module is used for carrying out sum operation on the phase after difference to obtain a phase difference summing result;
the phase extension module is used for extending the phase difference summation result to a range appointed by the determined phase region;
and the incoming wave direction obtaining module is used for obtaining the incoming wave direction of the signal source based on the extended phase and the preset phase address table.
7. The apparatus of claim 6, wherein the pre-processing module comprises:
the windowing Fourier transform unit is used for carrying out windowing Fourier transform on the signals received by each directional antenna and converting each time domain signal into a frequency domain signal;
and the complex finger conversion unit is used for carrying out complex finger conversion on each path of frequency domain signal to obtain the frequency spectrum amplitude and the phase of each path of signal.
8. The apparatus of claim 7, wherein the windowed Fourier transform unit comprises:
the N-point windowing subunit is used for carrying out N-point windowing on each path of signal, wherein N is the number of data points processed at one time;
and the Fourier transform subunit is used for performing Fourier transform on each channel of signals subjected to windowing and converting each channel of time domain signals into frequency domain signals.
9. The apparatus according to any one of claims 6-8, wherein the spectral amplitude of each signal comprises: the frequency spectrum amplitudes of the three signals are respectively M1, M2 and M3; the phase region determination module includes:
a size determination unit for determining sizes of the M1, the M2, and the M3 with respect to each other;
a first region determining unit, configured to determine that the phase region of the signal source is in a first region when M1 is maximum, where the first region is [ -4 pi, -2 pi ];
a fourth area determination unit, configured to determine that the phase area of the signal source is in a fourth area when M3 is maximum, where the fourth area is [2 pi, 4 pi ]
A second region determining unit for determining that the phase region of the signal source is in a second region when the M2 is maximum and the M1 is greater than the M3, the second region being [ -2 pi, 0 ];
a third region determining unit, configured to determine that the phase region of the signal source is in a third region when the M2 is maximum and the M1 is smaller than the M3, where the third region is [0, 2 pi ];
a fifth region determining unit for determining that the phase region of the signal source is in a fifth region when the M2 is maximum and the M1 is equal to the M3, the fifth region being [ -pi, pi ].
10. An electronic device, comprising: a memory and a processor;
the memory is used for storing programs;
the processor is configured to invoke a program stored in the memory to perform the method of any of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810647566.1A CN108848446B (en) | 2018-06-21 | 2018-06-21 | Data processing method and device and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810647566.1A CN108848446B (en) | 2018-06-21 | 2018-06-21 | Data processing method and device and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108848446A CN108848446A (en) | 2018-11-20 |
| CN108848446B true CN108848446B (en) | 2020-10-30 |
Family
ID=64203476
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810647566.1A Active CN108848446B (en) | 2018-06-21 | 2018-06-21 | Data processing method and device and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108848446B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112649678A (en) * | 2020-12-24 | 2021-04-13 | 广州山锋测控技术有限公司 | Antenna feeder measuring method and device, antenna feeder measuring device and tester |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1697360A (en) * | 2005-06-29 | 2005-11-16 | 北京邮电大学 | Method for testing aerrays system in use for multiple inputs and multiple outputs |
| CN105979485A (en) * | 2016-05-11 | 2016-09-28 | 南京邮电大学 | Personnel detection method in indoor environment based on channel state information (CSI) |
| CN106324559A (en) * | 2016-08-30 | 2017-01-11 | 西安空间无线电技术研究所 | Large-baseline four-element array broadband signal direction finding system and method |
| CN108051772A (en) * | 2017-10-25 | 2018-05-18 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Width be combined measurement incoming wave azimuth information method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8786488B2 (en) * | 2010-06-17 | 2014-07-22 | Amtech Systems, LLC | System and method for microwave ranging to a target in presence of clutter and multi-path effects |
-
2018
- 2018-06-21 CN CN201810647566.1A patent/CN108848446B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1697360A (en) * | 2005-06-29 | 2005-11-16 | 北京邮电大学 | Method for testing aerrays system in use for multiple inputs and multiple outputs |
| CN105979485A (en) * | 2016-05-11 | 2016-09-28 | 南京邮电大学 | Personnel detection method in indoor environment based on channel state information (CSI) |
| CN106324559A (en) * | 2016-08-30 | 2017-01-11 | 西安空间无线电技术研究所 | Large-baseline four-element array broadband signal direction finding system and method |
| CN108051772A (en) * | 2017-10-25 | 2018-05-18 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Width be combined measurement incoming wave azimuth information method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108848446A (en) | 2018-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2017351742B2 (en) | Direction of arrival estimation | |
| EP2453258B1 (en) | Radar device | |
| Zheng et al. | Classification and localization of mixed near-field and far-field sources using mixed-order statistics | |
| CN109959892B (en) | Uniform circular array dual-channel interferometer direction finding method, device and system | |
| CN102435978B (en) | Direction finding device for phase interferometer and phase spectrum interferometer based on multiple baselines | |
| Vasylyshyn | Removing the outliers in root-MUSIC via pseudo-noise resampling and conventional beamformer | |
| CA2266275C (en) | Direction finder and device for processing measurement results for the same | |
| CN110031793B (en) | Interferometer direction finding method, device and system | |
| AU2011373542B2 (en) | Wireless localisation system | |
| CN108848446B (en) | Data processing method and device and electronic equipment | |
| US20110264249A1 (en) | Method of, and apparatus for, planar audio tracking | |
| CN106980105A (en) | Electromagnetic Vector Sensor Array Space Rotating decorrelation LMS direction-finding method | |
| JP3808431B2 (en) | Direction finding device | |
| CA2948539C (en) | Calibration device | |
| Hsieh et al. | Efficient acquisition algorithm for long pseudorandom sequence | |
| US20060119514A1 (en) | Radio signal direction finder | |
| US10416277B2 (en) | Positioning device | |
| KR20190134893A (en) | Method of estimating direction of arrival of radar signal based on antenna array extrapolation and apparatus for the same | |
| Ning et al. | Velocity-independent and low-complexity method for 1D DOA estimation using an arbitrary cross-linear array | |
| Jiang et al. | Direction-of-arrival estimation for uniform circular arrays under small sample size | |
| ZA200603665B (en) | Methods and device for the radio determination of a number of spectrally overlapping radio stations | |
| CN113238184A (en) | Two-dimensional DOA estimation method based on non-circular signals | |
| CN114460531A (en) | Uniform linear array MUSIC spatial spectrum estimation method | |
| Swartling et al. | Calibration errors of uniform linear sensor arrays for DOA estimation: an analysis with SRP-PHAT | |
| CN110728022B (en) | Direction finding method, device and equipment based on uniform five-element circular array |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |