CN101325807B - Method for estimating signal wave direction - Google Patents
Method for estimating signal wave direction Download PDFInfo
- Publication number
- CN101325807B CN101325807B CN2008100226553A CN200810022655A CN101325807B CN 101325807 B CN101325807 B CN 101325807B CN 2008100226553 A CN2008100226553 A CN 2008100226553A CN 200810022655 A CN200810022655 A CN 200810022655A CN 101325807 B CN101325807 B CN 101325807B
- Authority
- CN
- China
- Prior art keywords
- signal
- matrix
- covariance
- vector
- noise
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Noise Elimination (AREA)
Abstract
The invention provides a signal DOA estimation method. The method supposes that an unknown noise covariance matrix has a symmetrical Toeplitz matrix character; according to the difference conception of the traditional covariance matrix, an imaginary number j is introduced so that the newly constructed covariance matrix changes into a central Hermitian matrix, thereby, the character decomposition can obtain the same DOA as the incident signal numbers. The method greatly reduces the number of the antenna array element required by estimating the signal wave arrival direction, and a wave crest is only formed in the wave arrival direction of the real signal. The method can reduce cost greatly in the practical application.
Description
Technical field
The present invention relates to a kind of network positions technology, it is background with the cognitive radio technology, utilize smart antenna, the direction of arrival of antenna receiving signal is estimated, thereby can provide technical basis for security service, commerce services, network management and the information service relevant with the position.
Background technology
Along with the continuous development of mobile communication technology, location technology obtains general application in the mobile communications network in recent years.Location technology can be applied to security service, commerce services, network management and the information service relevant with the position.For example: positioning service can be used for the personnel rescue under the situations such as earthquake, battlefield, also can be used for the fields such as communication navigation, target following in city.
According to the demand of different location services to positioning accuracy, the embodiment of navigation system is divided into based on terminal, Network Based, and both such as mix at three kinds of different schemes.Navigation system based on portable terminal (as mobile phone) claims the travelling carriage freedom positioning system again, is representative with GPS assist location etc.What the GPS location utilized is satellite-signal, and locating effect is bad under conditions such as indoor, basement, tunnel, subway.And based on network wireless location system utilizes existing mobile communication network, in the signal of network terminal test travelling carriage generation and the various parameters of estimating signal, as direction of arrival parameter, the time of advent and time difference parameter, the target location is estimated according to these parameters.The enforcement of these class methods is expensive less, instant effect, is the main flow of localization method in the mobile communications network.
Wherein the method that positions according to the direction of arrival of signal is utilized the Energy distribution of signal on all directions of space, to the space different come to signal differentiate.The traditional classical method of estimated signal direction of arrival comprises methods such as multiple signal classification and invariable rotary subspace.Signal in the space when arriving antenna with different time delays, produces correlation because mulitpath is propagated between signal, cause the effectively direction of arrival of estimated signal of said method.Traditional space covariance difference method is from direction of arrival that at first can only estimated signal, to the development that can estimate coherent signal, all there are two significant drawbacks, the one, the requirement bay is counted the twice that M must count q greater than incoming signal, be M>2q, the 2nd, all form the spectrum peak at ± θ place, wherein θ is that the true ripple of signal reaches the angle.Because top two shortcomings cause the cost increase in actual applications of traditional method.
Summary of the invention
Technical problem: the purpose of this invention is to provide a kind of method for estimating signal wave direction, this method is the correlation between ring off signal effectively, thereby can accurately estimate to arrive the direction of arrival of aerial signal, and can under limited antenna sensor number, estimate more signal number.
Technical scheme: the present invention is intended to utilize a kind of new space covariance difference method, and the direction of arrival of estimated signal when signal coherence, is removed because factor affecting such as multipath cause the correlation of signal.The present invention utilizes imaginary number j that traditional method performance is improved, here
Advantage of the present invention is mainly reflected in three aspects: the one, and array number M and incoming signal are counted q and are satisfied M>q, and the 2nd, only reach θ place, angle and form the spectrum peak at signal wave, the 3rd, can eliminate noise effect, with the Bo Dajiao of Search Space Smoothing binding energy estimation coherent signal.
The present invention can solve the direction of arrival estimation problem of signal under the Colored Noise,, combines with Search Space Smoothing and can remove this correlation because factor affecting such as multipath produce when being concerned with when signal, and therefore this invention also can be estimated the direction of arrival of coherent signal.
Method for estimating signal wave direction of the present invention utilizes the space covariance difference method, the direction of arrival of estimated signal, when signal coherence, remove because factor affecting such as multipath cause the correlation of signal, be embodied in three aspects: the one, array number M and incoming signal are counted q and are satisfied M>q, the 2nd, only reach θ place, angle and form the spectrum peak at signal wave, the 3rd, can eliminate noise effect, with the Bo Dajiao of Search Space Smoothing binding energy estimation coherent signal;
The concrete steps of this method are:
1.) receive data vector X (t), according to formula according to aerial array
Construction data covariance matrix R
X, wherein N is a data snap length, n is the discrete time point, ()
HExpression conjugate transpose, X (n) are n reception data constantly,
2.) when signal is irrelevant,, get covariance difference matrix Δ R for eliminating the influence of background noise
XFor
When signal coherence, for forward direction space smoothing covariance matrix R
f, construct new covariance difference matrix
Wherein
J is a transposed matrix, ()
*The expression conjugate operation, Δ R
X fBe called R
fThe covariance difference matrix,
3.) to covariance difference matrix Δ R
XCarry out feature decomposition, wherein M characteristic vector u
1, u
2, L, u
MConstitute matrix U=[u
1, u
2, L, u
M], with M eigenvalue
1, λ
2, L, λ
MBe diagonal entry, all the other elements are zero formation diagonal matrix sigma=diag (λ
1, λ
2, L, λ
M), diag () expression is a diagonal line value with the element in the round parentheses, other element is zero diagonal matrix,
4.) with the characteristic value of covariance difference matrix according to descending, q big characteristic value u
1, u
2, L, u
qPairing characteristic vector constitutes the signal vector matrix U of signal subspace
s=[u
1, u
2, L, u
q]; And M-q little characteristic value u
Q+1, u
Q+2, L, u
MThe characteristic of correspondence vector then constitutes noise subspace, i.e. the noise vector matrix U
N=[u
Q+1, u
Q+2, L, u
M],
5.) the noise vector matrix U
NAnd the steering vector a (θ) of signal=[1, e
-j2 π d (sin θ)/λ, L, e
-j2 π d (M-1) (sin θ)/λ]
TSubstitution
Here, M is an array number, λ is a carrier wavelength, d is the array element interval, P (θ) is called power spectrum, when θ from 0 to 180 degree is got different values, obtains different power spectral value, q the maximum that equates with signal number wherein must be arranged, and the pairing abscissa of this q maximum point is exactly that the ripple of all q signal reaches angle.
Beneficial effect: the present invention proposes a kind of covariance difference method, this method can be estimated the direction of arrival of signal under the Colored Noise, suppose that promptly unknown coloured noise covariance matrix has the Toeplitz matrix characteristic of symmetry, according to this characteristic, by constructing the difference of two covariance matrixes, signal wave is reached the angle estimation effect thereby eliminate coloured noise.Introduce imaginary number j and make this method, can significantly reduce the bay number, and only form crest at the true direction of arrival θ of signal place with respect to traditional method.This invention can reduce cost in actual applications greatly.
Description of drawings
Fig. 1 provides even linear array received signal model;
Fig. 2 provides the invention process flow chart.
Embodiment
For unknown Colored Noise, we suppose that noise covariance matrix has the Toeplitz matrix properties of symmetry.Aerial array is a uniform straight line array, and each array element all is omnidirectional antenna, and array number is M, and array element is spaced apart d.Suppose that far field, q arrowband incoherent signal is respectively from direction θ
i, i=1, L, q incides antenna array.Then can be expressed as at moment t array received data vector:
X(t)=AS(t)+n(t)
X (t)=[x wherein
1(t), L, x
M(t)]
TFor receiving data vector, A=[a (θ
1), L, a (θ
q)] be M * q dimension direction matrix, a (θ)=[1, e
-j2 π d (sin θ)/λ, L, e
-j2 π d (M-1) (sin θ)/λ]
TBe the steering vector corresponding to incidence angle θ, λ is a signal wavelength, S (t)=[s
1(t), L, s
q(t)]
TBe q * 1 dimensional signal vector, n (t)=[n
1(t), L n
M(t)]
TBe M * 1 dimension noise vector, ()
TThe computing of expression transposition.Here suppose that signal and array element noise statistics are independent.
Suppose x
i(t) be the zero-mean random process, structure array covariance matrix is:
R
X=E{X(t)X
H(t)}=AR
SA
H+Q
Wherein go up label " H " representing matrix conjugate transpose, R
S=E{S (t) S
H(t) } be q * q dimensional signal covariance matrix, Q is that M * M ties up noise covariance matrix.According to above-mentioned assumed condition as can be known, matrix R
SFull rank, and matrix Q has symmetry Toeplitz character.
For sake of convenience, introduce following three theorems:
Theorem 1: Matrix C is called central Hermitian matrix, if C satisfies following condition:
C=JC
*J, wherein J is a permutation matrix, is defined as:
Theorem 2: if C is symmetry Toeplitz matrix, C so
H, C
TWith JCJ also be the symmetry Toeplitz matrix.
Theorem 3: if n * n rank Matrix C is central Hermitian matrix, then its respective element satisfies:
Owing to the existence of noise, feasible estimation to the signal direction of arrival is had a strong impact in the communication environment.Existing covariance difference method utilizes the structure priori of noise covariance matrix Q, adopts a kind of reasonable manner transform array covariance matrix R
X, in this conversion, noise section is constant, and signal section obtains changing.Poor by the covariance matrix before and after the conversion, remove The noise, thereby obtain new covariance matrix.When tradition covariance difference method is removed noise, but brought two big shortcomings, the one, array number M and number of signals q must satisfy: M>2q, the 2nd, ± the θ place forms the spectrum peak of symmetry.These two shortcomings cause conventional method to run into very big trouble in actual applications, use cost and heighten.We propose new covariance difference method for this reason.
In order to eliminate The noise, the method below traditional covariance calculus of finite differences has been constructed:
ΔR
X=R
X-JR
XJ
The covariance matrix of this method construct is negative antisymmetric matrix, and these matrix characteristics are the always positive and negative paired appearance of characteristic value, and the result causes above-mentioned shortcoming, and in order to solve top shortcoming, we have invented following new method, make new covariance difference matrix be:
Here
The introducing of j makes Δ R
XBecome central Hermitian matrix, make the not positive and negative appearance of classical inverse symmetrical matrix characteristic value for another example of covariance matrix of neotectonics, so this method estimated signals power spectrum has only q spectrum peak, promptly only at real signal incident direction θ
iThe place forms the spectrum peak.Simultaneously, new method requires array number M and number of signals q to satisfy M>q to get final product.Can get after the simplification:
Top method only is fit to signal not height correlation or incoherent situation, when signal was concerned with owing to factors such as multipaths, conventional method lost efficacy, and we utilize the thought of forward direction Search Space Smoothing and this method to combine, produce the direction of arrival estimation problem that is fit under the coherent condition, even also
R wherein
fBe forward direction mean space smoothed covariance matrix.Substitution and reduced equation get:
From top two kinds of new methods as can be seen, noise section is eliminated from the covariance of neotectonics fully through after the difference, thereby can eliminate noise to signal direction of arrival estimation effect.
At last to the structure new covariance matrix Δ R
XOr Δ R
X f, carry out characteristic value decomposition, can get q nonzero eigenvalue and M-q zero eigenvalue, wherein open into noise subspace corresponding to the characteristic vector of M-q zero eigenvalue, be constructed as follows power spectrum spectrum estimation formulas:
U wherein
n=[u
1, L, u
M-q], u
1, L, u
M-qCharacteristic vector for noise subspace.
The present invention is further described below in conjunction with accompanying drawing.
1.) receive data vector X (t) and obtain data covariance matrix R according to aerial array
X
2.), when signal is irrelevant, construct new covariance difference matrix for eliminating The noise
When signal coherence, structure covariance difference matrix
3.) to Δ R
XCarry out feature decomposition,
4.) determine signal subspace U
SWith noise subspace U
N,
5.) according to noise subspace U
N, obtain power spectrum
6.) when θ from 0 to 180 degree is got different values, obtain different power spectral value, the abscissa of q maximum point correspondence is exactly the incident direction of q signal.
Claims (1)
1. method for estimating signal wave direction, it is characterized in that this method utilizes the space covariance difference method, the direction of arrival of estimated signal, when signal coherence, remove because the signal correlation that the multipath factor affecting causes is embodied in three aspects: the one, array number M and incoming signal are counted q and are satisfied M>q, and the 2nd, only reach θ place, angle and form the spectrum peak at signal wave, the 3rd, can eliminate noise effect, with the Bo Dajiao of Search Space Smoothing binding energy estimation coherent signal;
The concrete steps of this method are:
1.) receive data vector X (t), according to formula according to aerial array
Construction data covariance matrix R
X, wherein N is a data snap length, and n is the discrete time point, and H represents conjugate transpose, and X (n) is n reception data constantly,
2.) when signal is irrelevant,, get covariance difference matrix Δ R for eliminating the influence of background noise
XFor
When signal coherence, for forward direction space smoothing covariance matrix R
f, construct new covariance difference matrix
Wherein
J is a permutation matrix, ()
*The expression conjugate operation,
Be called R
fThe covariance difference matrix,
3.) to covariance difference matrix Δ R
XCarry out feature decomposition, wherein M characteristic vector u
1, u
2..., u
MConstitute matrix U=[u
1, u
2..., u
M], with M eigenvalue
1, λ
2..., λ
MBe diagonal entry, all the other elements are zero formation diagonal matrix sigma=diag (λ
1, λ
2..., λ
M), diag () expression is a diagonal line value with the element in the round parentheses, other element is zero diagonal matrix,
4.) with the characteristic value of covariance difference matrix according to descending, q big characteristic value u
1, u
2..., u
qPairing characteristic vector constitutes the signal vector matrix U of signal subspace
S=[u
1, u
2..., u
q]; And M-q little characteristic value u
Q+1, u
Q+2..., u
MThe characteristic of correspondence vector then constitutes noise subspace, i.e. the noise vector matrix U
N=[u
Q+1, u
Q+2..., u
M],
5.) the noise vector matrix U
NAnd the steering vector of signal
Here, M is an array number, and λ is a carrier wavelength, and d is the array element interval, ()
TThe computing of expression transposition, P (θ) is called power spectrum, when θ from 0 to 180 degree is got different values, obtains different power spectral value, q the maximum that equates with signal number wherein must be arranged, and the pairing abscissa of this q maximum point is exactly that the ripple of all q signal reaches angle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100226553A CN101325807B (en) | 2008-07-24 | 2008-07-24 | Method for estimating signal wave direction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100226553A CN101325807B (en) | 2008-07-24 | 2008-07-24 | Method for estimating signal wave direction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101325807A CN101325807A (en) | 2008-12-17 |
CN101325807B true CN101325807B (en) | 2011-01-12 |
Family
ID=40189070
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100226553A Expired - Fee Related CN101325807B (en) | 2008-07-24 | 2008-07-24 | Method for estimating signal wave direction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101325807B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054130A (en) * | 2016-06-06 | 2016-10-26 | 南京工程学院 | Indoor positioning method and device based on MUSIC algorithm |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592721B (en) * | 2009-06-19 | 2011-09-14 | 中国人民解放军理工大学 | Eigenvalue reconstruction based method for estimating angle of arrival of coherent signal |
CN101795150B (en) * | 2010-03-09 | 2013-01-23 | 西安电子科技大学 | Method for estimating direction of arrival and information source number of strong and weak signals |
CN103135091A (en) * | 2011-11-25 | 2013-06-05 | 上海无线电设备研究所 | Adaptive impulsive noise elimination method of DOA (direction of arrival) estimation system |
CN102706385B (en) * | 2012-05-10 | 2014-12-31 | 西安交通大学苏州研究院 | Method for detecting number of signals under condition of mixing of uncorrelated and correlated signals in uniform linear array |
CN103412188B (en) * | 2013-08-27 | 2016-03-23 | 重庆邮电大学 | Based on the SFM signal parameter estimation of Bessel's function and Toeplitz algorithm |
CN103546221B (en) * | 2013-10-25 | 2015-09-16 | 东南大学 | A kind of wide-band coherent signal direction of arrival estimation method |
CN103634903B (en) * | 2013-11-22 | 2016-08-31 | 中国人民解放军理工大学 | Low-power consumption mini-plant anti-noise localization method |
CN104820216B (en) * | 2015-05-08 | 2017-03-08 | 西北工业大学 | Multipath signal direction of arrival estimation method based on array response rotational invariance |
CN105548957B (en) * | 2016-01-18 | 2017-11-17 | 吉林大学 | Multiple target distance field mixing source localization method under a kind of unknown coloured noise |
CN107045117B (en) * | 2016-09-23 | 2019-10-25 | 南京工程学院 | Based on Capon Wave beam forming localization method and device |
CN106569180B (en) * | 2016-11-10 | 2020-01-21 | 中国人民解放军理工大学 | Prony method-based orientation estimation algorithm |
CN107255796B (en) * | 2017-07-25 | 2020-03-13 | 西安交通大学 | Method for positioning narrow-band near-field signal source under non-uniform noise |
CN107515841B (en) * | 2017-08-10 | 2021-02-23 | 中国科学院电子学研究所 | Method for calculating electromagnetic wave arrival angle |
CN109471063B (en) * | 2018-11-06 | 2022-10-28 | 江西师范大学 | Uniform linear array high-resolution direction-of-arrival estimation method based on delayed snapshot |
CN110320491A (en) * | 2019-07-22 | 2019-10-11 | 天津职业技术师范大学(中国职业培训指导教师进修中心) | Calculation method, device and the electronic equipment of direction of arrival |
CN110531311A (en) * | 2019-08-27 | 2019-12-03 | 武汉大学深圳研究院 | A kind of LTE external illuminators-based radar DOA estimation method based on matrix recombination |
CN113126021B (en) * | 2021-04-19 | 2022-03-29 | 电子科技大学 | Single-snapshot two-dimensional DOA estimation method based on three parallel linear arrays |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1746697A (en) * | 2005-10-18 | 2006-03-15 | 电子科技大学 | Multi-signal sorting algorithm with chip realization |
CN1968046A (en) * | 2006-10-19 | 2007-05-23 | 北京交通大学 | Estimating method of reach direction of user signal wave of array antenna MC-CDMA system |
-
2008
- 2008-07-24 CN CN2008100226553A patent/CN101325807B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1746697A (en) * | 2005-10-18 | 2006-03-15 | 电子科技大学 | Multi-signal sorting algorithm with chip realization |
CN1968046A (en) * | 2006-10-19 | 2007-05-23 | 北京交通大学 | Estimating method of reach direction of user signal wave of array antenna MC-CDMA system |
Non-Patent Citations (1)
Title |
---|
赵春晖等.基于无噪协方差矩阵变换的宽带测向算法.《兵工学报》.2008,第29卷(第3期),第291-295页. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054130A (en) * | 2016-06-06 | 2016-10-26 | 南京工程学院 | Indoor positioning method and device based on MUSIC algorithm |
CN106054130B (en) * | 2016-06-06 | 2019-01-22 | 南京工程学院 | A kind of indoor orientation method and device based on MUSIC algorithm |
Also Published As
Publication number | Publication date |
---|---|
CN101325807A (en) | 2008-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101325807B (en) | Method for estimating signal wave direction | |
CN105611627A (en) | Method for estimating AOA of WLAN access point based on double antennas | |
CN101917235B (en) | Angle estimating method of MIMO communication system positioning | |
CN101431354B (en) | Direction of arrival estimation method | |
CN105182293A (en) | Method for estimating DOA and DOD of MIMO radar based on co-prime array | |
CN103546221B (en) | A kind of wide-band coherent signal direction of arrival estimation method | |
CN104515969B (en) | Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method | |
CN102841344A (en) | Method for estimating parameters of near-field broadband signal resources by utilizing less array elements | |
CN103353596A (en) | Wave beam space domain meter wave radar height measurement method based on compressed sensing | |
CN102752849A (en) | Single receiving machine location method based on signal detection probability and wave angle estimation | |
CN102175989A (en) | Method for measuring incoherently distributed signal two-dimensional DOA (direction of arrival) | |
CN102694588A (en) | Arrival direction estimation method based on conjugation expansion | |
CN107589399A (en) | Based on the relatively prime array Wave arrival direction estimating method for sampling virtual signal singular values decomposition more | |
CN107255796A (en) | Arrowband near-field signals source localization method under a kind of non-uniform noise | |
CN103747455A (en) | Channel modeling method and parameter matching method based on non-uniform scatterer distribution | |
CN106896340A (en) | A kind of relatively prime array high accuracy Wave arrival direction estimating method based on compressed sensing | |
CN101908918B (en) | Beam synthesizing method in wireless communication receiver | |
CN107315161A (en) | Not rounded signal wave based on compressed sensing reaches bearing estimate method | |
CN104333407A (en) | Acquiring method and device of channel information | |
CN101232316B (en) | Method for estimating wave arriving based on mobile WIMAX | |
CN106291453A (en) | A kind of two-dimentional ripple based on beam space conversion reaches angle assay method | |
CN106226729A (en) | Relatively prime array direction of arrival angular estimation method based on fourth-order cumulant | |
CN102621539A (en) | Single searching wideband signal direction-of-arrival estimating method | |
CN101098178B (en) | Special channel broad business beam shaping implementation method and system | |
CN103499811B (en) | Antenna number distribution method capable of improving radar target estimation performance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110112 Termination date: 20130724 |