CN104898085B - A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation - Google Patents

A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation Download PDF

Info

Publication number
CN104898085B
CN104898085B CN201510245865.9A CN201510245865A CN104898085B CN 104898085 B CN104898085 B CN 104898085B CN 201510245865 A CN201510245865 A CN 201510245865A CN 104898085 B CN104898085 B CN 104898085B
Authority
CN
China
Prior art keywords
mrow
array
polarization
signal
span
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
Application number
CN201510245865.9A
Other languages
Chinese (zh)
Other versions
CN104898085A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201510245865.9A priority Critical patent/CN104898085B/en
Publication of CN104898085A publication Critical patent/CN104898085A/en
Application granted granted Critical
Publication of CN104898085B publication Critical patent/CN104898085B/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
    • G01S3/00Direction-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/02Direction-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/14Systems for determining direction or deviation from predetermined direction
    • G01S3/146Systems for determining direction or deviation from predetermined direction by comparing linear polarisation components

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation, belong to array signal process technique field.The signal correlation matrix that array received data obtain is carried out Eigenvalues Decomposition and obtains noise subspace by the present invention;Direction of arrival angle and polarization information in the steering vector of signal are separated, obtain intermediate variable;By the search function of intermediate variable construction direction of arrival angle, spectrum peak search obtains the estimation of direction of arrival angle;The estimate of direction of arrival angle is substituted into dimensionality reduction MUSIC Power estimations, coordinate obtains polarization estimation as corresponding to spectral peak.Four-dimensional spectrum peak search is reduced to two-dimentional spectrum peak search by the algorithm of the present invention so that the calculating magnitude of algorithm is from o { n4It is reduced to o { n2, the computational complexity of algorithm is greatly reduced, improves the operation efficiency of system.

Description

A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation
Technical field
The invention belongs to array signal process technique field, more particularly to polarization sensitive array Signal parameter estimation algorithm Improve, be concretely to reduce algorithm operation quantity, improve the dimensionality reduction MUSIC algorithms of operation efficiency.
Background technology
The space that polarization sensitive array can not only obtain signal by the difference of the locus of array element reaches angle information, And can enough obtains signal possessed polarization information in itself by the difference of array element quadrature phase, and the output of its array element is containing pole Change the Vector Message of information rather than common scalar information, this characteristic has established the property of further raising polarization sensitive array The basis of energy, while also cause us to be possible to further improve the disposal ability of array signal.Compared to common scalar array, Performance advantage existing for polarization sensitive array, it can be showed in terms of following four:(1) stronger antijamming capability;(2) Sane detectability;(3) higher resolution capability;(4) ability of polarization multiple access.
Signal parameter estimation problem is a basic problem of array signal processing, and since the 1960s, DOA estimates Meter technology obtains development at full speed, has emerged in large numbers substantial amounts of excellent algorithm.It is that Burg estimated in the maximum entropy spectrum proposed in 1967 first Meter method and SchmidtRO indicate that array antenna exists in the minimum variance Power estimation method proposed in 1969, both algorithms Improvement is obtained in terms of resolution capability.Multiple signal classification (Multiple Signal Classification, MUSIC) algorithm It is proposed so that the development of high-resolution direction finding technology has to be broken through tremendously.MUSIC algorithms be substantially using noise subspace and Obtained from the orthogonality principle of signal subspace.MUSIC algorithms have just obtained extensive research and application since being suggested, It is a kind of high-resolution features Structure Method for producing progressive unbiased esti-mator.
Multiple signal classification algorithm is used for the signal source angle of arrival of polarization sensitive array and the thought of polarization estimation Mainly by carrying out the mathematic decomposition such as feature decomposition or singular value decomposition to receiving data, data will be received and be divided into two phases Mutually orthogonal subspace, i.e. signal subspace and noise subspace, then using the relation of array manifold vector subspace come Tectonic syntaxis is composed, so as to obtain signal source angle of arrival and polarization estimation according to spectral peak.Traditional MUSIC combined spectrums method of estimation It is as follows:
1) array output covariance matrix is first obtained;
2) Eigenvalues Decomposition is carried out to covariance matrix, obtains noise subspace span { Un};
3) by noise subspace and array manifold matrixThe subspace opened and noise subspace span {UnOrthogonal principle, construction MUSIC joint Power estimations:
To estimate that the angle of arrival of signal source and polarization parameter need to carry out four-dimensional spectrum peak search it can be seen from formula (1), need Huge operand and storage capacity are paid, is not easy to realize.
The content of the invention
The invention provides a kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation, estimate in no reduction algorithm While counting error performance, four-dimensional spectrum peak search is reduced to two-dimentional spectrum peak search so that the calculating magnitude of algorithm is from o { n4Drop For o { n2(n be hunting zone in points), the computational complexity of algorithm is greatly reduced, improves the operation efficiency of system.
The present invention solve the above problems the technical scheme that uses for:
The dimensionality reduction MUSIC algorithms of a kind of polarization sensitive array parameter Estimation, it is characterised in that comprise the following steps:
Step 1. establishes the mathematical modeling of array received signal
Space has K incoherent signals to incide in array, according to the angle of arrival of incoming signalAnd polarization ginseng Number (γ, η) and noise obtain the expression formula of array received signal model
Wherein, θ ∈ [- pi/2, pi/2) be the incoming signal angle of pitch,For incoming signal azimuth, γ ∈ [0, Pi/2) for polarization explement, η ∈ [- π, π) for polarization phases it is poor.Manifold matrix A=[a of whole array1,a2,...,aK], s (n) =[s1(n),s2(n),...,sK(n)]TFor signal phasor, v (n) is noise.Array scanning vectorIt is that polarization is led To vectorWith steric direction vectorKronecker product, i.e.,
Its steering vector that polarizes is expressed as:
Wherein,Represent the transformation relation matrix between spherical coordinate system and rectangular coordinate system unit vector, Ek(γ, η) represent the Jones vectors of perfact polarization ripple.
It is noted that the difference in step 1 according to demandSix components can not have to all make With.
The n times snapshot data for the array received signal that step 2. is obtained using step 1, obtain estimating signal correlation matrix
Step 3. obtains noise subspace span { U by carrying out Eigenvalues Decomposition to correlation matrix Rn}
In formula, P is to takeComponent number, M is element number of array.Λ=diag { λ12,...,λPM, and And PM-K minimum characteristic value is equal, is equal to array received noise intensity.Characteristic vector structure corresponding to K big characteristic values Into signal subspace<S>=span { Us, characteristic vector corresponding to PM-K small characteristic values forms noise subspace<N>=span {Un}.Array manifold matrix is identical with signal subspace into subspace and orthogonal with noise subspace, is expressed as:
Span { A }=span { Us} (7)
span{A}⊥span{Un} (8)
Step 4. calculates the angle of pitch and azimuth process of signal source
Dimensionality reduction MUSIC algorithms are peeled off direction of arrival angle and polarization parameter in MUSIC Power estimations, according to traditional MUSIC Power estimationsDefined function
Order
By subspace principal, the subspace of array manifold matrix is orthogonal with noise subspace, i.e.,:
span{A}⊥span{Un} (8)
Formula (8) is substituted into formula (9) to obtain
When γ ∈ (0, pi/2), E is understood by formula (4)H(γ, η) is sequency spectrum, whenFor non-full rank, i.e.,Formula (12) is set up, then the angle of pitch and azimuthal estimation
Step 5. calculates the polarization parameter process of signal source
By the required angle of pitch and azimuthIntermediate variable is substituted into, is obtained
Then polarization estimation:
Signal source angle of arrival and polarization parameter are separated search by the present invention it can be seen from formula (13) and formula (15), One four-dimensional spectrum peak search is changed into two two-dimentional spectrum peak searches, by the method for the present invention, estimates incoming signal angle and pole The amount of calculation for changing parameter substantially reduces, and has good correctness and accuracy.
Brief description of the drawings
Fig. 1 is the even linear array to forming by biorthogonal electric dipole.
Embodiment
Even linear array is formed to being arranged along the y-axis by M biorthogonal electric dipole as shown in Figure 1, biorthogonal electric dipole pair Placed respectively along x-axis and y-axis, array element spacing is d, and dimensionality reduction MUSIC algorithms is implemented as follows in this example:
Step 1. establishes array received signal mathematical modeling:
Known spatial has K signal to incide in polarization sensitive array, and incident narrow band signal is orthogonal, θ ∈ [- pi/2, Pi/2) it is the incoming signal angle of pitch,For incoming signal azimuth, γ ∈ [0, pi/2) for polarization explement, η ∈ [- π, π) poor for polarization phases, the receipt signal model x (n) of array can be expressed as:
A=[a in formula1,a2,...aK] be whole array manifold matrix, wherein ak(k=1,2...K) it is polarizations Domain-spatial domain joint steering vector;Signal phasor s (n)=[s1(n),s2(n),...,sK(n)]T;V (n) is that average is zero, variance For σ2White complex gaussian noise;akIt is polarization steering vector (akp) and steric direction vector (aks) Kronecker product, i.e.,:
Wherein,
In formula,λ represents wavelength.
For the electric dipole pair along x-axis and y-axis placement, its steering vector a that polarizeskpTwo components only are used, are expressed as:
Formula (19) is substituted into formula (17) to obtain
Step 2. seeks the autocorrelation matrix of signal
The autocorrelation matrix of array received signal is:
Rx=E [x (n) xH(n)]=ARsAH2I (21)
Wherein, A be whole array manifold matrix, RsFor the auto-correlation function of incoming signal, I is unit battle array.
In Practical Project, the N number of observation sample value x (0) of generally use, x (1) ..., x (N-1) obtain estimating for correlation matrix MeterWithInstead of RxCalculated, wherein:
Step 3. solves noise subspace:
The estimation for the correlation matrix that step 3 is tried to achieveEigenvalues Decomposition is carried out, tries to achieve noise subspace:
Wherein ΛsAnd ΛnValue matrix is characterized, corresponds respectively to K big characteristic value and 2M-K small characteristic values.UsFor K Characteristic vector corresponding to big characteristic value, its into space be signal subspace, UnFor corresponding to 2M-K small characteristic values Characteristic vector, its into space be noise subspace.
Step 4. calculates the angle of pitch and the azimuth of signal source
Defining detection limit first is:
Formula in step 1 (20) is substituted into formula (24) to obtain
Due toSet up, then formula (25) can be expressed as:
By formula (26) it can be seen that the direction of arrival angle information of incoming signal and polarization parameter information are separated and come, now Define intermediate variable:
By subspace principal, the subspace of array manifold matrix is orthogonal with the characteristic vector of noise subspace
span{A}⊥span{Un} (28)
Formula (28) is substituted into formula (26) to obtain
The E when γ ∈ (0, pi/2)H(γ, η) is sequency spectrum, to be causedThenNon- full rank, I.e.K minimum of the spectrum of estimation zero, the position of these minimums are obtained according to optimized search methodThe as estimation of signal source angle of arrival.For convenience's sake, generally its K maximum is searched for the inverse of zero spectrum.It is fixed Adopted dimensionality reduction MUSIC spectrums:
Wherein coordinate corresponding to K spectral peak is the angle-of- arrival estimation of signal source
Step 5. calculates the polarization parameter of signal source:
The signal source angle of arrival that step 4 is obtained is substituted into formula (27), obtains the estimation of intermediate variable
Estimation for polarization parameter, it can search for obtain by following Power estimation:
For the signal source angle of arrival and polarization estimation of polarization sensitive array, traditional joint MUSIC algorithms need Spectrum peak search is carried out simultaneously to four parameters, its amount of calculation is o (n4), wherein n is to be counted in hunting zone, and the calculation of the present invention The angle of arrival of signal source and polarization parameter are carried out separation search by method, and its amount of calculation is o (2n2), it can be seen that greatly reduce Operand, improve the operation efficiency of system.

Claims (1)

1. the dimensionality reduction MUSIC algorithms of a kind of polarization sensitive array parameter Estimation, it is characterised in that comprise the following steps:
Step 1. establishes the mathematical modeling of array received signal
Space has K incoherent signals to incide in array, according to the angle of arrival of incoming signalAnd polarization parameter (γ, η) and noise obtain the expression formula of array received signal model:
Wherein, θ ∈ [- pi/2, pi/2) be the incoming signal angle of pitch,For incoming signal azimuth, γ ∈ [0, pi/2) be Polarize explement, and η ∈ [- π, π) it is that polarization phases are poor;Manifold matrix A=[a of whole array1,a2,...,aK], s (n)=[s1 (n),s2(n),...,sK(n)]TFor signal phasor, v (n) is noise;Array scanning vectorIt is that polarization is oriented to arrow AmountWith steric direction vectorKronecker product, i.e.,:
Its steering vector that polarizesIt is expressed as:
Wherein,Represent the transformation relation matrix between spherical coordinate system and rectangular coordinate system unit vector, Ek(γ, η) table Show the Jones vectors of perfact polarization ripple;
The n times snapshot data for the array received signal that step 2. is obtained using step 1, obtain estimating signal correlation matrix:
<mrow> <mi>R</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <msup> <mi>x</mi> <mi>H</mi> </msup> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
Step 3. obtains noise subspace span { U by carrying out Eigenvalues Decomposition to correlation matrix Rn}
<mrow> <mi>R</mi> <mo>=</mo> <msup> <mi>U&amp;Lambda;U</mi> <mi>H</mi> </msup> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>P</mi> <mi>M</mi> </mrow> </munderover> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <msub> <mi>u</mi> <mi>i</mi> </msub> <msup> <msub> <mi>u</mi> <mi>i</mi> </msub> <mi>H</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
In formula, P is to takeComponent number, M is element number of array, Λ=diag { λ12,...,λPM, and most PM-K small characteristic value is equal, is equal to array received noise intensity;Characteristic vector corresponding to K big characteristic values forms letter Work song space<S>=span { Us, characteristic vector corresponding to PM-K small characteristic values forms noise subspace<N>=span { Un}; Array manifold matrix is identical with signal subspace into subspace and orthogonal with noise subspace, is expressed as:
Span { A }=span { Us} (7)
span{A}⊥span{Un} (8)
Step 4. calculates the angle of pitch and azimuth process of signal source
Direction of arrival angle and polarization parameter in MUSIC Power estimations are peeled off, according to traditional MUSIC Power estimationsDefined function:
Order
By subspace principal, the subspace of array manifold matrix is orthogonal with noise subspace, i.e.,:
span{A}⊥span{Un} (8)
Formula (8) is substituted into formula (9) to obtain
When γ ∈ (0, pi/2), E is known by formula (4)H(γ, η) is sequency spectrum;WhenFor non-full rank, i.e.,Formula (12) is set up;
The then angle of pitch and azimuthal estimation:
Step 5. calculates the polarization parameter process of signal source
By the required angle of pitch and azimuthIntermediate variable is substituted into, is obtained
Then polarization estimation:
CN201510245865.9A 2015-05-14 2015-05-14 A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation Active CN104898085B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510245865.9A CN104898085B (en) 2015-05-14 2015-05-14 A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510245865.9A CN104898085B (en) 2015-05-14 2015-05-14 A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation

Publications (2)

Publication Number Publication Date
CN104898085A CN104898085A (en) 2015-09-09
CN104898085B true CN104898085B (en) 2017-11-17

Family

ID=54030850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510245865.9A Active CN104898085B (en) 2015-05-14 2015-05-14 A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation

Country Status (1)

Country Link
CN (1) CN104898085B (en)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105372635A (en) * 2015-11-19 2016-03-02 哈尔滨工业大学 Improved dimension-reduction space-time adaptive processing-based ship-borne high-frequency ground wave radar sea clutter suppression method
CN105445760A (en) * 2015-11-20 2016-03-30 北京自动化控制设备研究所 Satellite navigation single-antenna broadband interference prevention method
CN106169941A (en) * 2016-08-31 2016-11-30 电子科技大学 A kind of DOA method of estimation of round battle array based on noise power
JP6874127B2 (en) * 2016-09-22 2021-05-19 フラウンホッファー−ゲゼルシャフト ツァ フェルダールング デァ アンゲヴァンテン フォアシュンク エー.ファオ Arrival direction estimation device and arrival direction estimation method
CN106483494B (en) * 2016-10-13 2021-11-02 哈尔滨工业大学(威海) Full polarization interferometer based on space sampling antenna array and parameter estimation method thereof
CN107015213B (en) * 2017-03-31 2019-12-03 长江大学 Bistatic MIMO radar angle evaluation method based on MUSIC algorithm
CN106990386A (en) * 2017-04-13 2017-07-28 哈尔滨工程大学 A kind of single dipole polarization sensitive rotation array DOA and polarization parameter combined estimation method
CN107015191B (en) * 2017-05-18 2019-09-27 哈尔滨工程大学 One kind single dipole polarization sensitization array dimensionality reduction DOA estimation method under multi-path jamming environment
CN107656239B (en) * 2017-08-22 2020-12-04 哈尔滨工程大学 Coherent information source direction finding method based on polarization sensitive array
CN107870314B (en) * 2017-10-31 2021-06-29 西安电子科技大学 Complete electromagnetic component weighting fusion direction-finding optimization method based on polarization sensitive array
CN108414993B (en) * 2018-02-13 2022-03-04 西安电子科技大学 Resampling-based COLD array direction of arrival and polarization parameter joint estimation method
CN109188347B (en) * 2018-10-11 2022-11-22 北京遥感设备研究所 Signal polarization mode three-classification method based on MUSIC spectrum
CN111220953A (en) * 2018-11-26 2020-06-02 北京华航无线电测量研究所 Anti-interference method based on polarization sensitive array
CN109959891B (en) * 2019-04-11 2023-05-23 南京航空航天大学 Dimension-reducing spectral peak searching method for space angle and polarization parameters of electromagnetic vector L array
CN111585631A (en) * 2020-04-09 2020-08-25 哈尔滨工程大学 Beam forming method for improving transmission gain of distributed polarization sensitive array
CN112748407A (en) * 2020-12-15 2021-05-04 杭州电子科技大学 Airspace-polarization domain combined spectrum estimation method based on polarization sensitive area array
CN113253194B (en) * 2021-04-21 2022-07-08 中国电子科技集团公司第二十九研究所 Broadband arrival angle and polarization combined measurement method based on sparse representation

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102721943A (en) * 2012-06-11 2012-10-10 西安交通大学 L type sensor array-based incoherent signal two-dimensional direction-of-arrival tracking method and device
CN103792509A (en) * 2014-02-24 2014-05-14 西安电子科技大学 Two dimensional direction of arrival angle estimation method of electromagnetic signal
CN103926573A (en) * 2014-04-17 2014-07-16 哈尔滨工程大学 Mono-static MIMO radar distribution type target angle estimation method based on fourth-order cumulant
CN103941221A (en) * 2014-03-12 2014-07-23 西安电子科技大学 Method for estimating parameters of space stretching electromagnetic vector sensor array
CN103941222A (en) * 2014-03-12 2014-07-23 西安电子科技大学 Parameter estimation method based on Rayleigh-Ritz entropy dimensionality reduction MUSIC algorithm
CN104360305A (en) * 2014-10-24 2015-02-18 中国航天科技集团公司第五研究院第五一三研究所 Radiation source direction finding positioning method of uniting compressed sensing and signal cycle stationary characteristics
CN104375115A (en) * 2014-11-04 2015-02-25 中国航天科工集团第三研究院第八三五七研究所 Polarization sensitive array based non-circular signal DOA and polarization parameter joint estimation method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001305202A (en) * 2000-04-24 2001-10-31 Toyota Central Res & Dev Lab Inc Music spectrum computation method, and its device and medium
JP4294634B2 (en) * 2005-10-31 2009-07-15 日本電業工作株式会社 Calibration table creation method and arrival direction estimation method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102721943A (en) * 2012-06-11 2012-10-10 西安交通大学 L type sensor array-based incoherent signal two-dimensional direction-of-arrival tracking method and device
CN103792509A (en) * 2014-02-24 2014-05-14 西安电子科技大学 Two dimensional direction of arrival angle estimation method of electromagnetic signal
CN103941221A (en) * 2014-03-12 2014-07-23 西安电子科技大学 Method for estimating parameters of space stretching electromagnetic vector sensor array
CN103941222A (en) * 2014-03-12 2014-07-23 西安电子科技大学 Parameter estimation method based on Rayleigh-Ritz entropy dimensionality reduction MUSIC algorithm
CN103926573A (en) * 2014-04-17 2014-07-16 哈尔滨工程大学 Mono-static MIMO radar distribution type target angle estimation method based on fourth-order cumulant
CN104360305A (en) * 2014-10-24 2015-02-18 中国航天科技集团公司第五研究院第五一三研究所 Radiation source direction finding positioning method of uniting compressed sensing and signal cycle stationary characteristics
CN104375115A (en) * 2014-11-04 2015-02-25 中国航天科工集团第三研究院第八三五七研究所 Polarization sensitive array based non-circular signal DOA and polarization parameter joint estimation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
平面极化天线阵列的DOA及极化参数降维估计方法;司伟建等;《通信学报》;20141231;第35卷(第12期);28-35 *

Also Published As

Publication number Publication date
CN104898085A (en) 2015-09-09

Similar Documents

Publication Publication Date Title
CN104898085B (en) A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation
CN105510874B (en) A kind of dimensionality reduction MUSIC methods of Near-field sources multi-parameter Combined estimator
CN103941221B (en) Space stretching Electromagnetic Vector Sensor Array method for parameter estimation
CN103091671B (en) Bidimensional wave arrival direction estimation method based on non-concentric electromagnetic vector array radar
CN103941222B (en) Based on the method for parameter estimation of Rayleigh entropy dimensionality reduction MUSIC algorithm
CN102540138B (en) Multi-base-line phase searching type two-dimensional spatial spectrum direction-measuring method
CN104515969B (en) Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method
CN108120967B (en) Plane array DOA estimation method and equipment
CN108896954A (en) A kind of direction of arrival estimation method based on joint real value subspace in relatively prime battle array
CN110197209A (en) A kind of Emitter Recognition based on multi-feature fusion
CN105335615B (en) A kind of two dimension angular and polarization parameter combined estimation method of low complex degree
CN104375115A (en) Polarization sensitive array based non-circular signal DOA and polarization parameter joint estimation method
CN106019214A (en) DOA estimation method for broadband coherent signal source
CN104933290B (en) Multi-parameter Combined estimator quaternary counting method of the double L-shaped orthogonal galvanic couple of stretching to array
CN102662158B (en) Quick processing method for sensor antenna array received signals
CN111610486A (en) High-resolution accurate two-dimensional direction of arrival estimation method based on planar co-prime array virtual domain tensor space spectrum search
WO2021068495A1 (en) Degree of freedom-enhanced spatial spectrum estimation method based on block sampled tensor signal construction using planar co-prime array
Zheng et al. Two-dimensional DOA estimation for coprime planar array: A coarray tensor-based solution
CN101252382B (en) Wide frequency range signal polarizing and DOA estimating method and apparatus
CN107493106A (en) A kind of method of frequency and angle Combined estimator based on compressed sensing
Yan et al. Computationally efficient direction finding using polynomial rooting with reduced-order and real-valued computations
CN109581277B (en) A kind of four-dimensional antenna array DOA estimation method based on compressive sensing theory
CN106872936A (en) Near field sources L-type acoustic vector-sensor array row ambiguity solution Multiple Parameter Estimation Methods
CN106507952B (en) A kind of quick spatial spectrum computational methods based on circle battle array
CN109507634B (en) Blind far-field signal direction-of-arrival estimation method based on propagation operator under any sensor array

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant