CN107728102A - A kind of Wave arrival direction estimating method of sensor array - Google Patents
A kind of Wave arrival direction estimating method of sensor array Download PDFInfo
- Publication number
- CN107728102A CN107728102A CN201710911728.3A CN201710911728A CN107728102A CN 107728102 A CN107728102 A CN 107728102A CN 201710911728 A CN201710911728 A CN 201710911728A CN 107728102 A CN107728102 A CN 107728102A
- Authority
- CN
- China
- Prior art keywords
- sub
- array element
- matrix
- estimating method
- arrival direction
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
Abstract
The present invention relates to sensor technical field, more particularly to a kind of Wave arrival direction estimating method of sensor array, including：3D uniform planar linearly disposed antenna arrays are divided into several 3D sub-antenna array elements, the antenna opening of each 3D sub-antennas array element is evenly distributed in tri- directions of X, Y, Z；To reception signal and postpones signal progress relevant treatment acquisition 3D sub-antenna array element covariance matrix of each 3D sub-antennas array element in search window；Each 3D sub-antennas array element covariance matrix is subjected to singular value decomposition and obtains each corresponding signal subspace vector matrix of 3D sub-antennas array element covariance matrix the corresponding signal subspace vector matrix progress operation of direct product of 3D sub-antenna array element covariance matrixes on tri- directions of X, Y, Z.It this method solve the problem of sensor used in unmanned vehicle, unmanned plane, radio frequency identification etc. easily disturbs, error is big, resolution is low.
Description
Technical field
The present invention relates to sensor technical field, more particularly to a kind of Wave arrival direction estimating method of sensor array.
Background technology
The unmanned auto-control Intelligent internet of things equipment application radar sound such as traditional unmanned vehicle, unmanned plane, radio frequency identification
Wave sensor, to determine the position with surrounding objects, a kind of method is that the comprehensive realization of radar itself scans, another method
It is to use phased-array radar.Both approaches, it is required for radar that there is the sensor of high resolution bearing estimate.
The most of high-resolution DOA estimation algorithms proposed at present be all premised on array manifold is accurately known,
Therefore its is functional.However, due to being intercoupled between temperature, component aging and array element in the engineer applied of reality
Deng influence, various sensor errors are inevitably present.These errors directly results in actual array manifold and ideal
Array manifold between there is a certain degree of deviation, now, these high resolution remote sensing data algorithm performances can reduce, in addition lose
Effect.Therefore, array error is that array signal process technique moves towards a practical bottleneck, and technical staff has to face
A problem.
In summary, it is necessary to design it is a kind of anti-interference, low-power consumption, applied to unmanned vehicle, unmanned plane, radio frequency identification etc. have
There is the sensor array of high resolution bearing estimate ability.
The content of the invention
The present invention in view of the shortcomings of the prior art, proposes a kind of Wave arrival direction estimating method of sensor array, solved
Sensor used in unmanned vehicle, unmanned plane, radio frequency identification etc. easily disturbs, error is big, differentiates the problem of low.
To achieve the above object, the present invention adopts the following technical scheme that：
The present invention proposes a kind of Wave arrival direction estimating method of sensor array, including：
3D uniform planar linearly disposed antenna arrays are divided into several 3D sub-antenna array elements, the antenna of each 3D sub-antennas array element
Mouth is evenly distributed in tri- directions of X, Y, Z；
To reception signal and postpones signal progress relevant treatment acquisition 3D of each 3D sub-antennas array element in search window
Bay covariance matrix；
Each 3D sub-antennas array element covariance matrix is subjected to singular value decomposition and obtains each 3D sub-antennas array element covariance
The corresponding signal subspace vector matrix of matrix；
The corresponding signal subspace vector matrix of 3D sub-antenna array element covariance matrixes on tri- directions of X, Y, Z is entered
Row operation of direct productWherein, θ is direction of arrival angle.
Preferably, the reception signal in the search window and postpones signal are carried out at dimensionality reduction using dimensionality reduction transition matrix
Reason.
Preferably, the dimensionality reduction transition matrix is unitary transformation matrix, and the reception signal in the search window and delay are believed
Number become real domain matrix.
Preferably, the 3D sub-antennas array element covariance matrix can be analyzed to characteristic signal subspace and characteristic noise is empty
Between.
Preferably, the characteristic signal subspace can then carry out real-time update weighting square by using recurrent least square method
Battle array.
Preferably, the characteristic signal subspace can update tracking by subspace iteration method.
Preferably, the characteristic signal subspace can update tracking by using invariable rotary ESPRIT rules.
Preferably, the 3D sub-antennas array element on tri- directions of X, Y, Z according to its character pair value size reset after
Participate in operation of direct product.
Preferably, if the position relationship between two 3D sub-antennas array elements meets that adjacent 3D sub-antennas array element is closed
System, interpolation calculation, the interpolation calculation are carried out to two maximum angles of Wave beam forming gain corresponding to the 3D sub-antennas array element
As a result it is the direction of arrival of the wave beam.
Preferably, estimated according to the instantaneous of the one-step prediction value of the direction of arrival and the noise subspace when coming
The instantaneous measured value in ripple direction.
Beneficial effects of the present invention：A kind of Wave arrival direction estimating method of sensor array of the present invention, solves nobody
Sensor used in car, unmanned plane, radio frequency identification etc. easily disturbs, error is big, differentiates the problem of low.
Brief description of the drawings
With accompanying drawing, the invention will be further described, but the embodiment in accompanying drawing does not form any limitation of the invention.
Fig. 1 is the embodiment schematic flow sheet of Wave arrival direction estimating method one of sensor array of the present invention.
Embodiment
Technical solution of the present invention is further described with embodiment below in conjunction with the accompanying drawings, this is the preferable reality of the present invention
Apply example.A kind of Wave arrival direction estimating method of sensor array provided in an embodiment of the present invention can apply to sensor technology neck
Each scene in domain, including but not limited to 2G GSM, 3G CDMA, 4G LTE/LTE-A, 5G eMBB mobile communication, collection
Group communication, satellite communication, laser communication, fiber optic communication, DTV, radio frequency identification, power carrier, unmanned vehicle, unmanned plane, thing
The systems such as networking, radar, the embodiment of the present invention are not particularly limited to this.
As shown in figure 1, the present invention considers the linear antenna array situation of array element spacing d=λ/2 (λ is signal wavelength)
Under representative instance, and consider two incoherent narrow band signal far field scenes when centre frequency fc be 2.5G, propose one kind
The Wave arrival direction estimating method of sensor array：
Step S101,3D uniform planar linearly disposed antenna arrays are divided into several 3D sub-antenna array elements, each 3D sub-antennas
The antenna opening of array element is evenly distributed in tri- directions of X, Y, Z；
Step S102, to reception signal and postpones signal progress relevant treatment of each 3D sub-antennas array element in search window
Obtain 3D sub-antenna array element covariance matrixes；
Step S103, each 3D sub-antennas array element covariance matrix is subjected to singular value decomposition and obtains each 3D sub-antennas battle array
First corresponding signal subspace vector matrix of covariance matrix；
Step S104, the corresponding signal subspace of 3D sub-antenna array element covariance matrixes on tri- directions of X, Y, Z
Vector matrix carries out operation of direct productWherein, θ is direction of arrival angle.
Usually, sender uses frame format structural transmission data, because the data volume of transmission is bigger, if profit completely
Handled with all data, generally require to consume substantial amounts of hardware resource, delay is also larger, thus using dimensionality reduction turns
Change matrix and dimension-reduction treatment is carried out to the reception signal in search window and postpones signal.
Specifically, dimensionality reduction transition matrix in the present embodiment is unitary transformation matrix, by the reception signal in search window and
Postpones signal becomes real domain matrix.It should be noted that it can also implement to drop using the equivalent other matrixes of unitary transformation matrix
Dimension conversion, the present invention are not particularly limited to this.
The 3D sub-antenna array element covariance matrixes that data after dimensionality reduction conversion are constructed can be analyzed to characteristic signal subspace
With characteristic noise subspace, and meet characteristic signal subspace and characteristic noise Orthogonal Subspaces.The process of orthogonalization, can
To be completed using Schmidt Schmidt or amendment Schmidt Schmidt rules, other correlation methods are not particularly limited.
Due to the abominable of wireless transmission environments condition, the signal of transmission causes often by various external interference
Characteristic signal subspace can produce corresponding fluctuation, thus need to be tracked renewal to characteristic signal subspace.The present embodiment
In characteristic signal subspace can then carry out real-time update weighting matrix by using recurrent least square method.
Specifically, characteristic signal subspace can update tracking by subspace iteration method.Wherein use invariable rotary
ESPRIT rules track to update, and are a kind of conventional and effective methods.Similarly, phase can also be used for characteristic noise subspace
Method with principle is updated tracking, and other correlation method present invention are not particularly limited.
Cross-polarized arrays with vertical and horizontal polarization tend to lead to the sendaisle of good separation, and this is for multithread
Data send attractive.From this view, the +/- 45 degree of cross-polarized arrays typically used come from vertical polarization because sending
With two not same polarizations mixing in both horizontal polarizations.Thus, the planar array antenna in practical application can be analyzed to X, Y, Z
Multiple 3D sub-antennas array elements on three directions, the 3D sub-antennas array element on tri- directions of X, Y, Z according to its character pair
It is worth after size is reset and participates in operation of direct product.
Further, if the position relationship between two 3D sub-antenna array elements meets adjacent 3D sub-antennas array element relation,
Two maximum angles of Wave beam forming gain corresponding to 3D sub-antenna array elements are carried out with interpolation calculation, the interpolation calculation result is should
The direction of arrival of wave beam.
In order to further provide for estimated accuracy, can according to the one-step prediction value of direction of arrival and noise subspace it is instantaneous come
Estimate the instantaneous measured value of current arrival bearing.
Finally it should be noted that：Various embodiments above is merely illustrative of the technical solution of the present invention, rather than its limitations；To the greatest extent
The present invention is described in detail with reference to foregoing embodiments for pipe, it will be understood by those within the art that：Its according to
The technical scheme described in foregoing embodiments can so be modified, either which part or all technical characteristic are entered
Row equivalent substitution；And these modifications or replacement, the essence of appropriate technical solution is departed from various embodiments of the present invention technology
The scope of scheme.
Claims (10)
- A kind of 1. Wave arrival direction estimating method of sensor array, it is characterised in that including：3D uniform planar linearly disposed antenna arrays are divided into several 3D sub-antenna array elements, the antenna opening of each 3D sub-antennas array element exists X, tri- directions of Y, Z are evenly distributed；To reception signal and postpones signal progress relevant treatment acquisition 3D sub-antenna of each 3D sub-antennas array element in search window Array element covariance matrix；Each 3D sub-antennas array element covariance matrix is subjected to singular value decomposition and obtains each 3D sub-antennas array element covariance matrix Corresponding signal subspace vector matrix；The corresponding signal subspace vector matrix of 3D sub-antenna array element covariance matrixes on tri- directions of X, Y, Z is carried out straight Product computingWherein, θ is direction of arrival angle.
- 2. the Wave arrival direction estimating method of sensor array according to claim 1, it is characterised in that changed using dimensionality reduction Matrix carries out dimension-reduction treatment to the reception signal in the search window and postpones signal.
- 3. the Wave arrival direction estimating method of sensor array according to claim 2, it is characterised in that the dimensionality reduction conversion Matrix is unitary transformation matrix, and the reception signal in the search window and postpones signal are become into real domain matrix.
- 4. according to the Wave arrival direction estimating method of any described sensor arrays of claim 1-3, it is characterised in that the 3D Sub-antenna array element covariance matrix can be analyzed to characteristic signal subspace and characteristic noise subspace.
- 5. the Wave arrival direction estimating method of sensor array according to claim 4, it is characterised in that the characteristic signal Subspace can then carry out real-time update weighting matrix by using recurrent least square method.
- 6. according to the Wave arrival direction estimating method of any described sensor arrays of claim 1-3, it is characterised in that the spy Sign signal subspace can update tracking by subspace iteration method.
- 7. the Wave arrival direction estimating method of sensor array according to claim 6, it is characterised in that the characteristic signal Subspace can update tracking by using invariable rotary ESPRIT rules.
- 8. the Wave arrival direction estimating method of sensor array according to claim 1, it is characterised in that the 3D sub-antennas Array element participates in operation of direct product after being reset according to its character pair value size on tri- directions of X, Y, Z.
- 9. according to the Wave arrival direction estimating method of any described sensor array in claim 1,2,3,5,7 or 8, its feature exists In if the position relationship between two 3D sub-antennas array elements meets adjacent 3D sub-antennas array element relation, to 3D Two maximum angles of Wave beam forming gain corresponding to bay carry out interpolation calculation, and the interpolation calculation result is the wave beam Direction of arrival.
- 10. the Wave arrival direction estimating method of sensor array according to claim 9, it is characterised in that according to the ripple The instantaneous of one-step prediction value and the noise subspace up to direction estimates the instantaneous measured value of current arrival bearing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710911728.3A CN107728102A (en) | 2017-09-29 | 2017-09-29 | A kind of Wave arrival direction estimating method of sensor array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710911728.3A CN107728102A (en) | 2017-09-29 | 2017-09-29 | A kind of Wave arrival direction estimating method of sensor array |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107728102A true CN107728102A (en) | 2018-02-23 |
Family
ID=61209337
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710911728.3A Pending CN107728102A (en) | 2017-09-29 | 2017-09-29 | A kind of Wave arrival direction estimating method of sensor array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107728102A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109507634A (en) * | 2018-11-08 | 2019-03-22 | 中国电子科技集团公司第二十八研究所 | A kind of blind far-field signal Wave arrival direction estimating method based on sensing operator under any sensor array |
| CN110187304A (en) * | 2019-05-21 | 2019-08-30 | 泰凌微电子（上海）有限公司 | A kind of Direction-of-arrival method and device |
| CN112363108A (en) * | 2020-10-28 | 2021-02-12 | 西安电子科技大学 | Signal subspace weighted super-resolution direction-of-arrival detection method and system |
| CN112740058A (en) * | 2019-12-26 | 2021-04-30 | 华为技术有限公司 | Method and device for measuring angle of arrival |
-
2017
- 2017-09-29 CN CN201710911728.3A patent/CN107728102A/en active Pending
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109507634A (en) * | 2018-11-08 | 2019-03-22 | 中国电子科技集团公司第二十八研究所 | A kind of blind far-field signal Wave arrival direction estimating method based on sensing operator under any sensor array |
| CN109507634B (en) * | 2018-11-08 | 2020-08-11 | 中国电子科技集团公司第二十八研究所 | Blind far-field signal direction-of-arrival estimation method based on propagation operator under any sensor array |
| CN110187304A (en) * | 2019-05-21 | 2019-08-30 | 泰凌微电子（上海）有限公司 | A kind of Direction-of-arrival method and device |
| CN110187304B (en) * | 2019-05-21 | 2021-05-04 | 泰凌微电子(上海)股份有限公司 | Signal arrival angle estimation method and device |
| CN112740058A (en) * | 2019-12-26 | 2021-04-30 | 华为技术有限公司 | Method and device for measuring angle of arrival |
| CN112740058B (en) * | 2019-12-26 | 2022-04-22 | 华为技术有限公司 | Method and device for measuring angle of arrival |
| CN112363108A (en) * | 2020-10-28 | 2021-02-12 | 西安电子科技大学 | Signal subspace weighted super-resolution direction-of-arrival detection method and system |
| CN112363108B (en) * | 2020-10-28 | 2022-11-01 | 西安电子科技大学 | Signal subspace weighting super-resolution direction-of-arrival detection method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107728102A (en) | A kind of Wave arrival direction estimating method of sensor array | |
| CN107290709B (en) | The relatively prime array Wave arrival direction estimating method decomposed based on vandermonde | |
| Balamurugan et al. | DOA tracking for seamless connectivity in beamformed IoT-based drones | |
| CN106682293B (en) | Sparse beam pattern comprehensive designing method | |
| CN107991659B (en) | Method for measuring height of low-elevation target of meter wave radar based on dictionary learning | |
| CN110113088A (en) | A kind of divergence type numerical model analysis antenna system direction of arrival intelligence estimation method | |
| Qi et al. | Deep learning-based BackCom multiple beamforming for 6G UAV IoT networks | |
| CN113466848A (en) | Angle flicker noise scene-oriented co-location MIMO radar multi-target tracking resource optimal allocation method | |
| Athley et al. | On radar detection and direction finding using sparse arrays | |
| Yi et al. | Sea clutter suppression method for shipborne high‐frequency surface‐wave radar | |
| CN108614234A (en) | Wave arrival direction estimating method based on more sampling relatively prime array received signal inverse fast Fourier transforms of snap | |
| CN109507634B (en) | Blind far-field signal direction-of-arrival estimation method based on propagation operator under any sensor array | |
| Yang et al. | Coprime L‐shaped array connected by a triangular spatially‐spread electromagnetic‐vector‐sensor for two‐dimensional direction of arrival estimation | |
| CN107402371A (en) | A kind of smart antenna array | |
| Liu et al. | Two-dimensional DOA estimation using a co-prime symmetric cross array | |
| Bingbing et al. | An indoor positioning algorithm and its experiment research based on RFID | |
| Xiong et al. | Non‐stationary clutter suppression method for bistatic airborne radar based on adaptive segmentation and space–time compensation | |
| Cui et al. | Hybrid precoding for millimetre wave MIMO systems based on particle swarm optimisation | |
| CN114554398A (en) | Indoor positioning method, first positioning server and indoor positioning system | |
| Li et al. | Robust DOA Estimation in Satellite Systems in Presence of Coherent Signals Subject to Low SNR | |
| CN110398733B (en) | Radar cluster single snapshot collaborative imaging method based on time synchronization | |
| Sun et al. | Through‐the‐wall radar imaging algorithm for moving target under wall parameter uncertainties | |
| Wang et al. | Low elevation and range joint estimation method with metre wave frequency diverse array‐multiple input multiple output radar based on sparse array | |
| Feng et al. | Constained adaptive monopulse algorithm based on sub-array | |
| Xiao et al. | Lstm based multiple beamforming for 5g haps iot networks |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180223 |