CN108445442A - A kind of near-field signals source localization method based on the singular value decomposition blocked - Google Patents

A kind of near-field signals source localization method based on the singular value decomposition blocked Download PDF

Info

Publication number
CN108445442A
CN108445442A CN201810168021.2A CN201810168021A CN108445442A CN 108445442 A CN108445442 A CN 108445442A CN 201810168021 A CN201810168021 A CN 201810168021A CN 108445442 A CN108445442 A CN 108445442A
Authority
CN
China
Prior art keywords
value
angle
field signals
signal
indicate
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.)
Granted
Application number
CN201810168021.2A
Other languages
Chinese (zh)
Other versions
CN108445442B (en
Inventor
辛景民
刘文怡
左炜亮
郑南宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201810168021.2A priority Critical patent/CN108445442B/en
Publication of CN108445442A publication Critical patent/CN108445442A/en
Application granted granted Critical
Publication of CN108445442B publication Critical patent/CN108445442B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computing Systems (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention discloses a kind of near-field signals source localization method based on the singular value decomposition blocked, and includes the following steps:1. the electrical angle of estimation of near field signalInitial value, and calculate initial covariance matrix R;2. calculating the oblique projection operator of electrical angle initial value;3. obtaining updated covariance matrix according to oblique projection operator4. according to electrical angle initial value and updated covariance matrix, electric angle angle value is updated;5. if update electric angle angle value and the difference of electrical angle initial value are less than given threshold, the positioning that final electric angle angle value completes near-field signals source is obtained, otherwise using updated electric angle angle value as initial value, step 2 to 5 is repeated and is iterated.The floating error caused by the characteristic value by noise subspace is eliminated, precision of prediction is improved.Effectively eliminate influencing each other between multiple near-field signals, by iteration for several times can within the scope of certain allowable error approaching to reality value.

Description

A kind of near-field signals source localization method based on the singular value decomposition blocked
Technical field
The invention belongs to array signal processing technologies, are related to a kind of localization method of near field multisignal source, specially A kind of near-field signals source localization method based on the singular value decomposition blocked.
Background technology
Signal source positioning is the root problem of field of signal processing, and it is in microphone array, radar, sonar, wireless The fields such as communication, seismology and robot are all widely used.Therefore, educational circles proposes a large amount of algorithm and solves far field scene Under narrow band signal arrival bearing's estimation problem.Wherein, array received to wave be considered as plane wave.However, near field Under scene, the wave for being incident on array has spherical corrugated, this also cause signal that array received arrives should by direction of arrival and Two parameters of distance describe simultaneously.
The algorithm of most of near-field signals source positioning is all based on does approximation with Taylor's second order expension to spherical corrugated.Recently, Also it has been proposed that some are based on high-order statistic, the algorithm of Cyclostationarity and maximal possibility estimation.Compared to above-mentioned These three algorithms, path following algorithm, the algorithm of Hu Tianyu high-precision signals subspace, the algorithm based on tensor and based on rotation The algorithm computation burden for turning the Signal parameter estimation of invariance is smaller, without calculating high-order statistic or multi-dimensional search.
By using the special property of the second-order statistics for the signal observed from symmetrically and evenly linear array, there is Person proposes a kind of algorithm of the linear prediction of weighting.However, in order to take the covariance matrix of signal into account, it is necessary to about incidence The priori of the design feature of signal covariance matrix, and this is only just effective when number of snapshots are larger.So this method The poor performance when number of snapshots are smaller because the error of covariance estimation cause the Mutual coupling near field signal-to-noise ratio it is low and All occurs saturated phenomenon in the case of height.Here " it is saturated " evaluated error for referring to direction of arrival in the case of high s/n ratio still Keep higher value.
Invention content
For problems of the prior art, the present invention provides a kind of near-field signals based on the singular value decomposition blocked Source localization method is utilized the anti-diagonal element of covariance matrix this design feature, passes through the iteration based on oblique projection operator Algorithm, saturated phenomenon is obtained for effective solution in the case of lower signal-to-noise ratio and smaller number of snapshots.
The present invention is to be achieved through the following technical solutions:
A kind of near-field signals source localization method based on the singular value decomposition blocked, includes the following steps:
Step 1, the electrical angle of estimation of near field signalWithInitial value, and calculate initial covariance square Battle array R;
Step 2, the oblique projection operator of electrical angle initial value is calculated;
Step 3, updated covariance matrix is obtained according to oblique projection operator
Step 4, according to electrical angle initial value and updated covariance matrix, electric angle angle value is updated;
Step 5, if update electric angle angle value and the difference of electrical angle initial value are less than given threshold, final electric angle is obtained Angle value completes the positioning in near-field signals source, otherwise using updated electric angle angle value as initial value, repeats step 2 to 5 and changes Generation.
Preferably, which is characterized in that in the step 1, estimated using the singular value decomposition and linear prediction method blocked The initial value of the electrical angle of near-field signals.
Further, which is characterized in that utilize the singular value decomposition and linear prediction method estimation of near field signal blocked The initial value of electrical angle, is as follows:
It is poised for battle and is listed in the near-field signals data that moment n is received It is handled as follows, wherein 2M+1 is array element element number of array, ()TIndicate transposition;
Step 1.1, initial covariance matrix R is calculated;
X (n) is expressed as vector matrix form:
Wherein, K indicates the total number of near-field signals, θkAnd rkIndicate k-th of incoming signal incident angle and it is incident away from From a (θk,rk) it is direction vector, sk(n) sampled value of k-th of the signal of expression in moment n;A is direction matrix, is expressed asS (n) indicates sampled value of the near-field signals vector in moment n of incidence, It is expressed asω (n) indicates additive noise vector in the sampled value of moment n, table It is shown as
R=E { x (n) x (n)H};
Wherein, expectation, () are asked in E { } expressionsHIndicate Ai Er meter Te transposition;
Step 1.2, linear prediction model is built according to the back-diagonal element of R;
P-th of element of the b articles back-diagonal vector of R can be expressed as,
Wherein,Indicate the power of k-th of near-field signals, σ2Indicate the power of noise signal, electrical angle ψkAnd φkRespectively It is defined asD is array element spacing, and λ is near-field signals wavelength, B=0, ± 1 ..., ± 2M, p=-M+b-,…M-b+,δ () indicates Crow Interior gram of function;
And then linear prediction model is obtained,
Wherein,The M-b of the b articles back-diagonal of representing matrix R+- l-1 elements, ()TIndicate transposition,It is the coefficient vector of linear prediction model, Q representative model exponent numbers;
Step 1.3, the coefficient vector of linear prediction model is estimated using the singular value decomposition blocked;
First, it is R to the Matrix-Vector form of linear prediction model(b)a(b)=g(b)Singular value decomposition is carried out to obtain,
Wherein, ()HAi Er meter Te transposition are represented,
Then, coefficient vector is calculated by following formula,
Wherein, the preceding K characteristic value of representation signal subspace is larger, represents the rear q-1-K characteristic value of noise subspace It is smaller, i.e.,
Step 1.4, according to coefficient vector a(b), structure linear prediction multinomial D(b)(z);
Wherein, z=2j (ψk-bφk),
Step 1.5, according to D(b)(z) justify K nearest zero in z-plane parasang to estimate to obtain the electricity of near-field signals Angle ψkAnd φkInitial value.
Preferably, which is characterized in that in the step 2, K oblique projection operator of electrical angle initial value is as follows,
Wherein, Indicate k-th of direction vector,Expression remove direction to AmountDirection matrix afterwards.
Preferably, which is characterized in that in step 3, using oblique projection operator Separation of Coherent Signals and incoherent signal, obtain Updated covariance matrix
Wherein, I2M+1Indicate the unit matrix of 2M+1 dimensions,Indicate the estimated value of initial covariance matrix R,Indicate kth The signal power of a near-field signals, It is RwEstimated value, Indicate the surplus of non-signal part in initial covariance matrix.
Further, which is characterized in that the estimated value of initial covariance matrix RIndicate as follows,
Wherein, { }*Conjugate operation is sought in expression,Indicate the signal power of k-th of near-field signals,Indicate k-th of side To vector,Direction vector is removed in expressionDirection matrix afterwards,N is indicated Total hits,Signal s is removed in representativek(n) near-field signals vector,Representative and AkCorresponding signal covariance square Battle array.
Preferably, which is characterized in that step 5 is as follows,
Judge whether update electric angle angle value and the difference of electrical angle initial value is less than given threshold ε according to the following formula,
Iteration is terminated if setting up, exports final electric angle angle value
Otherwise, incremented circulation index index (i=i+1) continues iteration.
Preferably, according to mapping relations following between actual parameter and electrical angle, final electrical angle is done and is mapped Direction of arrival to multisignal source and incident distance, complete the positioning in near-field signals source;
Wherein, d is array element spacing, and λ is the wavelength in incident near-field signals source.
Compared with prior art, the present invention has technique effect beneficial below:
In the localization method of near-field signals source of the present invention, the singular value decomposition blocked only utilizes and near-field signals son sky Between relevant characteristic value calculate linear predictor coefficient, eliminate the floating error caused by the characteristic value by noise subspace, carry High precision of prediction.The introducing of oblique projection operator, on the one hand so that the facing arrays aperture of each signal increases, improving can profit Information content;On the other hand, due to only considering a signal every time, the mutual shadow between multiple near-field signals is effectively eliminated Ring, by iteration for several times can within the scope of certain allowable error approaching to reality value.So as in the lower feelings of number of snapshots Also superior performance can be kept under condition, therefore robustness is stronger, successfully solve the " full of parameter Estimation in the case of high s/n ratio With " phenomenon.
Description of the drawings
Fig. 1 is to indicate corrected symmetrical uniform linear array according to the ... of the embodiment of the present invention.
Fig. 2 a are the variations for indicating near-field signals Mutual coupling performance with signal-to-noise ratio.
Fig. 2 b are the variations for indicating near-field signals Mutual coupling performance with number of snapshots.
Specific implementation mode
With reference to specific embodiment, the present invention is described in further detail, it is described be explanation of the invention and It is not to limit.
Data model of the present invention and problem are described as follows.
As shown in Figure 1, it is contemplated that K near field narrowband incoherent signal { sk(n) } one is incident on by 2M+1 in moment n On the uniform linear array that the array element that a spacing is d forms, and assume that the array of these array elements composition is calibrated completely.If The center array element for determining array is phase reference point, then the near-field signals x with noise that m-th of array element receivesmIt (n) can be with table It is shown as
Wherein, m=-M ..., M, sk(n) indicate k-th of signal in the sampled value of moment n, ωm(n) be additive noise when Carve the sampled value of n, τmkIt is k-th of signal phase delay caused by time delay between reference array element and m-th of array elementWherein θkAnd rkIndicate k-th of signal direction of arrival angle and Distance, λ are the wavelength of incident near-field signals.When k-th of signal is in Fresnel region, time delay item τmkIt can be by second order Taylor expansion obtains good approximation, i.e.,
Wherein,Since they are actual parameter It is mapped to simpler form, so being called electrical angle.Ignored exponent number is represented to be equal to or more thanItem.So, the near-field signals vector received It can be expressed as again with the form of vector sum matrix in the sampled value of moment n
Wherein, ()TRepresent transposition, s (n) and ω (n) be byWithGiven incoming signal vector sum additive noise vector respectively.A is By the direction matrix of the uniform linear array by calibration, it is expressed as It is the direction vector of array, can be expressed as
In the present invention, it is assumed that the response matrix A full ranks (for example, the order of A is K) of array, incoming signal { sk(n) } it is The stationary stochastic process of the zero-mean of broad sense, additive noise { ωm(n) } be space-time complex field white Gaussian random process, mean value It is zero, variance σ2, and with signal { sk(n) } uncorrelated.The number K of incoming signal is it is known that and meet K≤M.
Near-field signals source of the present invention localization method is described below, is included the following steps:
1) electrical angle of estimation of near field signalWithInitial value, and calculate initial covariance matrix R;
2) the oblique projection operator of electrical angle initial value is calculated;
3) updated covariance matrix is obtained according to oblique projection operator
4) according to electrical angle initial value and updated covariance matrix, electric angle angle value is updated;
5) if updated electric angle angle value and the difference of electrical angle initial value are less than given threshold, final electric angle is obtained Angle value completes the positioning in near-field signals source, otherwise using updated electric angle angle value as initial value, repeats step 2 to 5 and changes Generation.
In the step 1), the electrical angle of the singular value decomposition and linear prediction method estimation of near field signal blocked is utilized Initial value is as follows:The near-field signals data that 2M+1 is clapped It is handled as follows,
A, initial covariance matrix R is calculated
R=E { x (n) x (n)H} (4)
Wherein, expectation, () are asked in E { } expressionsHIndicate Ai Er meter Te transposition.
B, linear prediction model is built according to the back-diagonal element of R
Wherein,Indicate p-th of element of the vector being made of the b row back-diagonal elements of covariance matrix R, K tables Show near-field signals number, b=-M+b-..., M-b+, p=-M+b_... M-b+, δ () indicate Kronecker function.
All elementsConstitute vector r(b).Vector r(b)Before being divided into that L is overlapped and being made of q element To vector, then entryIt can predict to obtain from the linear combination by other vector elements
Wherein,
(6) the Matrix-Vector form of formula is R(b)a(b)=g(b) (7)
Wherein,It is the coefficient vector of linear prediction model,
q Representative model exponent number.
C, the coefficient vector of linear prediction model described above is estimated using the singular value decomposition blocked.Formula (7) is done very Different value decomposes to obtain
Wherein, ()HAi Er meter Te transposition are represented,
Further, coefficient vector can be calculated by following formula
Wherein, the preceding K characteristic value of representation signal subspace is larger, represents the rear q-1-K characteristic value of noise subspace It is smaller, i.e.,
D, according to coefficient vector a(b), structure linear prediction multinomial D(b)(z)。
Wherein, z=2j (ψk-bφk),
The electrical angle ψ of near-field signalskAnd φkφkInitial value can be from D(b)(z) justify nearest K in z-plane parasang A zero is estimated to obtain.
2) K oblique projection operator of electrical angle initial value is calculated;
The estimated value of initial covariance matrix RIt is expressed as following form
Wherein, { }*Conjugate operation is sought in expression,Indicate the signal power of k-th of near-field signals,Indicate k-th of side To vector,Direction vector is removed in expressionDirection matrix afterwards,N is indicated Total hits,Signal s is removed in representativek(n) near-field signals vector,Representative and AkCorresponding signal covariance square Battle array.
In order to eliminate influencing each other between signal and noise, i.e. Section 2, Section 3 and Section 4 in formula (11), count Calculate K new oblique projection operators
Wherein,
Wherein, I2M+1Indicate the unit matrix of 2M+1 dimensions.
3) oblique projection operator Separation of Coherent Signals and incoherent signal, the covariance matrix obtained after update are utilized
To formula (11), premultiplication and the right side multiply projection operator respectively, and wherein Section 2, Section 3 and Section 4 are due to projection operator Characteristic be eliminated, i.e.,
Wherein, I2M+1Indicate the unit matrix of 2M+1 dimensions,Indicate the estimated value of initial covariance matrix R, It is RwEstimated value, indicate non-near in initial covariance matrix The surplus of field signal and noise signal portions.When signal-to-noise ratio is sufficiently large,It is sufficiently small, then It is obvious that at this timeOnly include the ψ of k-th of near-field signalsk-bφkInformation.
4) according to electrical angle initial value and updated covariance matrix, electric angle angle value is updated;
Updated covariance matrix substitution step 1) can be calculated to the estimated value of this iteration.
5) if update electric angle angle value and the difference of electrical angle initial value are less than given threshold, final electric angle angle value is obtained The positioning for completing near-field signals source repeats step 2 to 5 and is iterated otherwise using updated electric angle angle value as initial value.
Specifically, being decided whether to carry out next step iteration according to the difference of electrical angle between iteration twice, if iteration is poor Value is less than a certain small threshold, such as:
Wherein, ε is arbitrary a smaller normal number, terminates iteration, exports final electric angle angle valueOtherwise, incremented circulation index index (i=i+1) continues iteration.Because of actual parameter and electricity Meet certain mapping relations between angle, i.e., So right Electrical angle, which does corresponding mapping, can solve the direction of arrival of multisignal source and incident distance.
Numerical simulation demonstrates the validity of the method for the present invention,
The symmetrical uniform linear array that method proposed in this paper is made of by one array element that 9 spacing are d Carry out simulating, verifying.It is assumed that two incident near-field signals power are equal, wavelength λ, incident angle and distance is respectively (6 °, 2.5 λ), (20 °, 3.5 λ).The threshold value of iterative algorithm is set as 10-6.Fixed signal-to-noise ratio and number of snapshots are respectively set to 10dB and 200.
Fig. 2 a and Fig. 2 b indicate the estimation performance of angle and distance in the positioning of near-field signals source with signal-to-noise ratio and snap respectively Several change curves.As can be seen that in signal-to-noise ratio and larger number of snapshots, method (solid line in figure) proposed in this paper and performance are estimated The theory lower bound (dotted line in figure) of meter closely, illustrates that iterative algorithm effectively eliminates saturated phenomenon, greatly improves Parameter Estimation Precision.

Claims (8)

1. a kind of near-field signals source localization method based on the singular value decomposition blocked, which is characterized in that include the following steps:
Step 1, the electrical angle of estimation of near field signalWithInitial value, and calculate initial covariance matrix R;
Step 2, the oblique projection operator of electrical angle initial value is calculated;
Step 3, updated covariance matrix is obtained according to oblique projection operator
Step 4, according to electrical angle initial value and updated covariance matrix, electric angle angle value is updated;
Step 5, if update electric angle angle value and the difference of electrical angle initial value are less than given threshold, final electric angle angle value is obtained The positioning for completing near-field signals source repeats step 2 to 5 and is iterated otherwise using updated electric angle angle value as initial value.
2. a kind of near-field signals source localization method based on the singular value decomposition blocked according to claim 1, feature Be, in the step 1, using the singular value decomposition and linear prediction method estimation of near field signal blocked electrical angle it is initial Value.
3. a kind of near-field signals source localization method based on the singular value decomposition blocked according to claim 2, feature It is, utilizes the initial value of the electrical angle of the singular value decomposition and linear prediction method estimation of near field signal blocked, specific steps It is as follows:
It is poised for battle and is listed in the near-field signals data that moment n is received It is handled as follows, wherein 2M+1 is array element element number of array, ()TIndicate transposition;
Step 1.1, initial covariance matrix R is calculated;
X (n) is expressed as vector matrix form:
Wherein, K indicates the total number of near-field signals, θkAnd rkIndicate the incident angle of k-th of incoming signal and incident distance, a (θk,rk) it is direction vector, sk(n) sampled value of k-th of the signal of expression in moment n;A is direction matrix, is expressed asS (n) indicates sampled value of the near-field signals vector in moment n of incidence, It is expressed asω (n) indicates additive noise vector in the sampled value of moment n, table It is shown as
R=E { x (n) x (n)H};
Wherein, expectation, () are asked in E { } expressionsHIndicate Ai Er meter Te transposition;
Step 1.2, linear prediction model is built according to the back-diagonal element of R;
P-th of element of the b articles back-diagonal vector of R can be expressed as,
Wherein,Indicate the power of k-th of near-field signals, σ2Indicate the power of noise signal, electrical angle ψkWithkIt is respectively defined asD be array element spacing, λ be near-field signals wavelength, b=0, ± 1 ..., ± 2M, p=-M+b-,…M-b+,δ () indicates Kronecker letter Number;
And then linear prediction model is obtained,
Wherein,The M-b of the b articles back-diagonal of representing matrix R+- l-1 elements, ()TIndicate transposition,It is the coefficient vector of linear prediction model,q Representative model exponent number;
Step 1.3, the coefficient vector of linear prediction model is estimated using the singular value decomposition blocked;
First, it is R to the Matrix-Vector form of linear prediction model(b)a(b)=g(b)Singular value decomposition is carried out to obtain,
R(b)Hg(b)=U(b)Λ(b)V(b)H,
Wherein, ()HAi Er meter Te transposition are represented,
Then, coefficient vector is calculated by following formula,
Wherein, the preceding K characteristic value of representation signal subspace is larger, and the rear q-1-K characteristic value for representing noise subspace is smaller, I.e.
Step 1.4, according to coefficient vector a(b), structure linear prediction multinomial D(b)(z);
Wherein, z=2j (ψk-bφk),
Step 1.5, according to D(b)(z) justify K nearest zero in z-plane parasang to estimate to obtain the electrical angle ψ of near-field signalsk And φkInitial value.
4. a kind of near-field signals source localization method based on the singular value decomposition blocked according to claim 1, feature It is, in the step 2, K oblique projection operator of electrical angle initial value is as follows,
Wherein, Indicate k-th of direction vector,Direction vector is removed in expression Direction matrix afterwards.
5. a kind of near-field signals source localization method based on the singular value decomposition blocked according to claim 1, feature It is, in step 3, using oblique projection operator Separation of Coherent Signals and incoherent signal, obtains updated covariance matrix
Wherein, I2M+1Indicate the unit matrix of 2M+1 dimensions,Indicate the estimated value of initial covariance matrix R,Indicate k-th closely The signal power of field signal, It is RwEstimated value, indicate The surplus of non-signal part in initial covariance matrix.
6. a kind of near-field signals source localization method based on the singular value decomposition blocked according to claim 5, feature It is, the estimated value of initial covariance matrix RIndicate as follows,
Wherein, { }*Conjugate operation is sought in expression,Indicate the signal power of k-th of near-field signals,Indicate k-th of direction to Amount,Direction vector is removed in expressionDirection matrix afterwards,N expressions are always adopted Sample number,Signal s is removed in representativek(n) near-field signals vector,Representative and AkCorresponding signal covariance matrix.
7. a kind of near-field signals source localization method based on the singular value decomposition blocked according to claim 1, feature It is, step 5 is as follows,
Judge whether update electric angle angle value and the difference of electrical angle initial value is less than given threshold ε according to the following formula,
Iteration is terminated if setting up, exports final electric angle angle value
Otherwise, incremented circulation index index (i=i+1) continues iteration.
8. a kind of near-field signals source localization method based on the singular value decomposition blocked according to claim 1, feature It is, according to mapping relations following between actual parameter and electrical angle, mapping is done to final electrical angle and obtains multisignal source Direction of arrival and incident distance, complete the positioning in near-field signals source;
Wherein, d is array element spacing, and λ is the wavelength in incident near-field signals source.
CN201810168021.2A 2018-02-28 2018-02-28 Near-field signal source positioning method based on truncated singular value decomposition Active CN108445442B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810168021.2A CN108445442B (en) 2018-02-28 2018-02-28 Near-field signal source positioning method based on truncated singular value decomposition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810168021.2A CN108445442B (en) 2018-02-28 2018-02-28 Near-field signal source positioning method based on truncated singular value decomposition

Publications (2)

Publication Number Publication Date
CN108445442A true CN108445442A (en) 2018-08-24
CN108445442B CN108445442B (en) 2020-08-18

Family

ID=63193036

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810168021.2A Active CN108445442B (en) 2018-02-28 2018-02-28 Near-field signal source positioning method based on truncated singular value decomposition

Country Status (1)

Country Link
CN (1) CN108445442B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109740117A (en) * 2019-01-31 2019-05-10 广东工业大学 A kind of robust and quick magnetic orientation algorithm

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102736063A (en) * 2012-07-06 2012-10-17 绍兴文理学院 Near-field sound source positioning method
CN103344939A (en) * 2013-06-14 2013-10-09 西安交通大学 Two-dimension direction of arrival estimating method of an incoherent and coherent blended signal
CN103954931A (en) * 2014-04-28 2014-07-30 西安交通大学 Method for locating far field and near field mixed signal sources
CN104598732A (en) * 2015-01-19 2015-05-06 西安交通大学 Semi-parameter number estimation method for coherent and incoherent mixed signals
CN107064901A (en) * 2017-04-27 2017-08-18 哈尔滨工业大学 A kind of method for estimating target azimuth of carrier-borne High frequency ground wave over-the-horizon aadar
CN107064892A (en) * 2016-11-11 2017-08-18 长江大学 MIMO radar angle estimation algorithm based on tensor subspace and invariable rotary

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102736063A (en) * 2012-07-06 2012-10-17 绍兴文理学院 Near-field sound source positioning method
CN103344939A (en) * 2013-06-14 2013-10-09 西安交通大学 Two-dimension direction of arrival estimating method of an incoherent and coherent blended signal
CN103954931A (en) * 2014-04-28 2014-07-30 西安交通大学 Method for locating far field and near field mixed signal sources
CN104598732A (en) * 2015-01-19 2015-05-06 西安交通大学 Semi-parameter number estimation method for coherent and incoherent mixed signals
CN107064892A (en) * 2016-11-11 2017-08-18 长江大学 MIMO radar angle estimation algorithm based on tensor subspace and invariable rotary
CN107064901A (en) * 2017-04-27 2017-08-18 哈尔滨工业大学 A kind of method for estimating target azimuth of carrier-borne High frequency ground wave over-the-horizon aadar

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WEILIANG ZUO等: ""A COMPUTATIONALLY EFFICIENT SOURCE LOCALIZATION METHOD FOR A MIXTURE OF NEAR-FIELD AND FAR-FIELD NARROWBAND SIGNALS"", 《2014 IEEE INTERNATIONAL CONFERENCE ON ACOUSTIC, SPEECH AND SIGNAL PROCESSING (ICASSP)》 *
张莜等: ""基于截断奇异值分解法(TSVD)的在线相干到达方向(DOA)估计算法"", 《科学技术与工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109740117A (en) * 2019-01-31 2019-05-10 广东工业大学 A kind of robust and quick magnetic orientation algorithm
CN109740117B (en) * 2019-01-31 2021-03-23 广东工业大学 Robust and fast magnetic positioning algorithm

Also Published As

Publication number Publication date
CN108445442B (en) 2020-08-18

Similar Documents

Publication Publication Date Title
Wu et al. Direction of arrival estimation for off-grid signals based on sparse Bayesian learning
Zhou et al. Direction-of-arrival estimation for coprime array via virtual array interpolation
CN107817465A (en) DOA estimation method based on non-grid compressed sensing under super-Gaussian noise background
CN110045323B (en) Matrix filling-based co-prime matrix robust adaptive beamforming algorithm
CN112565122B (en) Super-large-scale MIMO channel estimation method based on Newton-orthogonal matching pursuit
CN109375153B (en) Compact multipath signal angle estimation method based on impulse response compressed sensing
CN113438682B (en) SAGE-BEM5G wireless channel parameter extraction method based on beam forming
Yan et al. MUSIC-like direction of arrival estimation based on virtual array transformation
CN114884841A (en) Underdetermined parameter joint estimation method based on high-order statistics and non-uniform array
Boujemâa Extension of COMET algorithm to multiple diffuse source localization in azimuth and elevation
CN108445442A (en) A kind of near-field signals source localization method based on the singular value decomposition blocked
Yang et al. A correlation-aware sparse Bayesian perspective for DOA estimation with off-grid sources
CN113030849B (en) Near field source positioning method based on self-encoder and parallel network
CN114460535A (en) Multi-array non-circular source direct positioning method based on weighted Euler ESPRIT
Zhang et al. A parameter estimation algorithm for multiple frequency‐hopping signals based on sparse Bayesian method
Guo et al. Deep Learning Based DOA Estimation With Trainable-Step-Size LMS Algorithm
Wang et al. Off-grid DOA Estimation for Temporally Correlated Source via Robust Block-SBL in Mutual Coupling
Li et al. 2-D parameter estimation of near field source based on robust propagator method
Tayem et al. Capon root-MUSIC-like direction of arrival estimation based on real data
Tang et al. A Direction-Of-Arrival Estimation Algorithm Based On Spatial Translation for Coherent Signals
Lee Improve method on DOA estimation accuracy of signal subspace sensitivity based on mutual coupling matrix
CN117478186B (en) Beam forming method based on covariance matrix estimation and subspace construction
Qian et al. Direct Position Determination of Non-Gaussian Sources for Sensor Arrays via Improved Rooting Subspace Data Fusion Method
Liu et al. A direction of arrival estimation method based on deep learning
Tian et al. Near-field Source Location Estimation Algorithm Based on Coprime 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