CN104122533A - Joint parameter estimation method based on distributed polarization sensitive array - Google Patents
Joint parameter estimation method based on distributed polarization sensitive array Download PDFInfo
- Publication number
- CN104122533A CN104122533A CN201410366085.5A CN201410366085A CN104122533A CN 104122533 A CN104122533 A CN 104122533A CN 201410366085 A CN201410366085 A CN 201410366085A CN 104122533 A CN104122533 A CN 104122533A
- Authority
- CN
- China
- Prior art keywords
- matrix
- array
- sub
- polarization sensitive
- polarization
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 135
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000011159 matrix material Substances 0.000 claims abstract description 110
- 239000013598 vector Substances 0.000 claims description 27
- 238000003491 array Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 230000005684 electric field Effects 0.000 description 6
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 238000002330 electrospray ionisation mass spectrometry Methods 0.000 description 3
- 230000001808 coupling effect Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention provides a joint parameter estimation method based on a distributed polarization sensitive array. The joint parameter estimation method includes dispersedly setting array element components of a polarization sensitive array consisting of electric dipole pairs in the space so as to form the distributed polarization sensitive array, dividing the distributed polarization sensitive array into a first sub array and a second sub array according to sequence of array elements, acquiring a covariance matrix of incident signals after the distributed polarization sensitive array receives the incident signals, establishing a signal sub space according to the covariance matrix, dividing the signal sub space consisting of the covariance matrix into a first sub matrix and a second sub matrix according to sequence of line numbers of the matrix, and acquiring estimation to angle of arrival and polarization parameter of the incident signals according to rotational invariance of the first sub matrix and the second sub matrix.
Description
Technical Field
The invention relates to the field of array signal processing, in particular to a parameter joint estimation method based on a distributed polarization sensitive array.
Background
A complete electromagnetic vector sensor is composed of 3 electric dipoles and 3 magnetic dipoles which are arranged in space and are mutually orthogonal in the same point in space, so that a polarization sensitive array is formed, and can receive all electric field components and magnetic field components of incident electromagnetic waves, so that the polarization sensitive array can receive more information of incident signals compared with a traditional scalar array. And the polarization sensitive array can sense the polarization information of the incident signal, thereby obtaining the polarization parameter of the incident electromagnetic signal. However, conventional scalar arrays cannot obtain polarization parameters of an incident electromagnetic signal because they cannot sense polarization information of the incident signal. And the polarization sensitive array can simultaneously induce polarization information and airspace information of incident electromagnetic waves. Thus, the polarization sensitive array, whether used for polarization parameter estimation or adaptive beam forming, has superior system performance over conventional scalar arrays.
In the application of the polarization sensitive array, by utilizing the vector relation among an electric field, a magnetic field and a poynting vector, when a single complete electromagnetic vector sensor is arranged in a space, the electromagnetic vector sensor can be utilized to simultaneously obtain the DOA (DOA) and the polarization parameter estimation of at most 5 uncorrelated signals, so that the method has important significance in the occasion of limited space physical aperture.
However, for the signal processing of the polarization sensitive array, it is mostly assumed that each array element is composed of 2 to 6 mutually orthogonal electric dipoles or magnetic dipoles which are co-located, and therefore, the co-located spatial mutual coupling of the dipoles inevitably has a severe mutual coupling effect, which degrades the performance of the antenna system.
The mutual coupling phenomenon among the array elements is inevitable, and in order to effectively reduce the mutual influence of the mutual coupling among the common-point channels of the array elements, the prior art provides a distributed polarization sensitive array, wherein the distributed polarization sensitive array is formed by dispersedly placing the common-point components of the array elements of the polarization sensitive array in space, so that the mutual coupling effect among the array elements is greatly reduced, and meanwhile, the electric field information and the polarization information of incident electromagnetic waves can be induced. Most of the existing parameter estimation methods for distributed polarization sensitive arrays aim at complete electromagnetic vector sensors, namely 3 electric dipoles and 3 magnetic dipoles are placed in a spatially dispersed manner, and then the parameter estimation is completed by using an improved vector cross product method. However, in practice, since the space electric field and the magnetic field are time-varying, the time-varying electric field can generate a magnetic field, and the time-varying magnetic field can generate an electric field, and at the same time, a polarization sensitive array is formed by using electric dipoles and magnetic dipoles, and a certain redundancy relationship exists.
Disclosure of Invention
In order to overcome the defects, the invention provides a parameter joint estimation method based on a distributed polarization sensitive array, which comprises the following steps:
dispersing and placing each array element component in a polarization sensitive array formed by electric dipole pairs in a space to form a distributed polarization sensitive array;
dividing the distributed polarization sensitive array into a first sub-array and a second sub-array according to the front and back sequence of the array elements;
after the distributed polarization sensitive array receives an incident signal, obtaining a covariance matrix of the incident signal;
constructing a signal subspace based on the covariance matrix;
dividing the signal subspace formed by the matrix into a first sub-matrix and a second sub-matrix according to the front and back sequence of the row number in the matrix;
and obtaining the estimation of the arrival angle and the polarization parameter of the incident signal by utilizing the rotation invariance between the first sub-matrix and the second sub-matrix.
Further, the dispersing and placing of each array element component in the polarization sensitive array formed by the electric dipole pair in the space specifically includes:
in the established coordinate system comprising a horizontal axis and a vertical axis, the array elements in the polarization sensitive array are distributed on the vertical axis at the same intervals, and the electric dipoles are alternately arranged in the direction parallel to the horizontal axis and the direction parallel to the vertical axis.
Further, when the distributed polarization sensitive array includes M array elements, the distributed polarization sensitive array is divided into a first sub-array and a second sub-array according to the sequence of the array elements, specifically:
and dividing the front M-2 array elements of the distributed polarization sensitive array into a first sub-array, and dividing the rear M-2 array elements of the distributed polarization sensitive array into a second sub-array.
Further, constructing a signal subspace based on the covariance matrix specifically includes:
and carrying out eigenvalue decomposition on the covariance matrix to construct a signal subspace.
Further, according to the front-back order of the number of rows in the matrix, the signal subspace formed by the matrix is divided into a first sub-matrix and a second sub-matrix, specifically:
the first M-2 rows of the matrix in the partitioned signal subspace constitute a first sub-matrix and the last M-2 rows of the matrix in the partitioned signal subspace constitute a second sub-matrix.
Further, obtaining the estimation of the angle of arrival and the polarization parameter of the incident signal by using the rotation invariance between the first sub-matrix and the second sub-matrix specifically includes:
obtaining rotation invariant characteristic parameters between the first sub-matrix and the second sub-matrix based on a total least square method rotation invariant method;
performing characteristic value decomposition on the rotation invariant characteristic parameters to obtain characteristic values and characteristic vectors corresponding to the characteristic values;
based on the eigenvalues and the eigenvectors, estimates of the angle of arrival and polarization parameters of the incident signal are obtained.
Further, obtaining an estimation of an angle of arrival and a polarization parameter of the incident signal based on the eigenvalue and the eigenvector specifically includes:
obtaining an angle of arrival of the incident signal based on the feature value;
obtaining an estimate of a first director sub-matrix corresponding to the first sub-matrix based on a matrix of eigenvalues comprised of eigenvalues and a matrix of eigenvectors comprised of eigenvectors;
based on the estimate of the first steering vector submatrix, an estimate of a polarization phase angle and an estimate of a polarization phase difference of the incident signal are obtained.
The invention has the beneficial effects that: the distributed polarization sensitive array is formed by dispersedly placing the array element components in the polarization sensitive array formed by the electric dipole pairs in the space, so that the influence of mutual coupling among the array elements can be reduced, the performance of the system is improved, information redundancy caused by forming the polarization sensitive array by simultaneously utilizing the electric dipole and the magnetic dipole is avoided, more incident signal electromagnetic information is obtained, the hardware cost of the system can be effectively reduced, and the complexity of realizing parameter joint estimation is reduced;
by dividing the signal subspace into a first sub-matrix and a second sub-matrix and utilizing the rotation invariance between the first sub-matrix and the second sub-matrix, the estimation of the arrival angle and the polarization parameter of the incident signal is obtained, the arrival angle and the polarization parameter can be automatically paired, an additional parameter pairing process is not needed, the complexity of parameter joint estimation is further reduced, and the calculation amount is reduced.
Drawings
FIG. 1 is a flow chart of a parameter joint estimation method based on a distributed polarization sensitive array according to the present invention;
FIG. 2 is a schematic diagram of the location of a distributed polarization sensitive array of the present invention;
FIG. 3 is a comparison graph of the performance of the present invention and the prior art co-point polarization sensitive array based on the estimation of the angle of arrival with the same number of channels;
FIG. 4 is a graph comparing the performance of the present invention and a prior art scalar array based on polarization phase angle estimation for the same array element number.
Detailed Description
The technical scheme of the invention is described in detail in the following with reference to the accompanying drawings.
The parameter joint estimation method based on the distributed polarization sensitive array comprises the following steps: dispersing and placing each array element component in a polarization sensitive array formed by electric dipole pairs in a space to form a distributed polarization sensitive array; dividing the distributed polarization sensitive array into a first sub-array and a second sub-array according to the front and back sequence of the array elements; after the distributed polarization sensitive array receives an incident signal, obtaining a covariance matrix of the incident signal; constructing a signal subspace based on the covariance matrix; dividing the signal subspace formed by the matrix into a first sub-matrix and a second sub-matrix according to the front and back sequence of the row number in the matrix; and obtaining the estimation of the arrival angle and the polarization parameter of the incident signal by utilizing the rotation invariance between the first sub-matrix and the second sub-matrix.
The distributed polarization sensitive array is formed by dispersedly placing the array element components in the polarization sensitive array formed by the electric dipole pairs in the space, so that the influence of mutual coupling between the array elements can be reduced, the performance of the system is improved, information redundancy caused by forming the polarization sensitive array by simultaneously utilizing the electric dipole and the magnetic dipole is avoided, more incident signal electromagnetic information is obtained, the hardware cost of the system can be effectively reduced, and the complexity of realizing parameter joint estimation is reduced. By dividing the signal subspace into a first sub-matrix and a second sub-matrix and utilizing the rotation invariance between the first sub-matrix and the second sub-matrix, the estimation of the arrival angle and the polarization parameter of the incident signal is obtained, the arrival angle and the polarization parameter can be automatically paired, an additional parameter pairing process is not needed, the complexity of parameter joint estimation is further reduced, and the calculation amount is reduced.
In the present application, a parameter joint estimation method based on a distributed polarization sensitive array, as shown in fig. 1, includes:
step 101: and dispersing and placing each array element component in the polarization sensitive array formed by the electric dipole pairs in the space, thereby forming a distributed polarization sensitive array.
Specifically, in the established coordinate system comprising a horizontal axis x axis and a vertical axis y axis, the array elements in the polarization sensitive array are distributed on the vertical axis at the same intervals, and the electric dipoles are alternately arranged in the direction parallel to the horizontal axis and the direction parallel to the vertical axis.
In the specific implementation process, the polarization sensitive array is a uniform linear array, each array element is distributed and placed on the y axis, the distance between the array elements is d, d is lambda/2, and lambda is the wavelength of an incident signal. After splitting, the electric dipole pairs are placed in a mode of alternating a first direction and a second direction, the first direction is a direction parallel to the x axis, the second direction is a direction parallel to the y axis, and finally, the formed distributed polarization sensitive array is as shown in fig. 2, wherein the polarization sensitive array is formed by the electric dipole pairs, so that the number M of array elements of the distributed polarization sensitive array is an even number.
After step 101 is completed, the present application performs step 102: and dividing the distributed polarization sensitive array into a first sub-array and a second sub-array according to the front-back sequence of the array elements.
Specifically, the first M-2 array elements of the distributed polarization sensitive matrix are divided into a first sub-array, and the last M-2 array elements are divided into a second sub-array.
Further, an incident signal is received by a distributed polarization sensitive array comprising a first sub-array and a second sub-array, and step 103 is performed:
and after the distributed polarization sensitive array receives the incident signal, obtaining a covariance matrix of the incident signal.
In the specific implementation process, N times of snapshot sampling is generally used firstly, so that the distributed polarization sensitive array of the M array elements receives data of N snapshots, then the time average is used for replacing the statistical average, and the covariance matrix of the incident signals is estimatedNamely:
where N is the number of sampled fast beats, H represents the conjugate transpose, x (t) is the model of the received incident signal of the distributed polarization sensitive array, t is time, and the incident signal x (t) can be represented as:
where K is the number of incident signals, akSteering vector, s, for the joint polarization-space domain of the k-th incident signalk(t) is the complex envelope of the incident signal, N (t) is complex white Gaussian noise, A is the steering vector matrix of the polarization domain-space domain combination, and S (t) is the incident signal matrix.
It should be noted that when the electric dipole pairs in the polarization sensitive array formed by the electric dipole pairs are disposed in a concurrent manner, the polarization steering vector of each array element can be expressed as:
where j is the complex imaginary unit, px、pyPolarization response of electric dipoles parallel to the x-axis direction and the y-axis direction, respectively, gamma and eta are the polarization phase angle and the polarization phase difference of incident electromagnetic waves, respectively, and satisfy the conditions that gamma is greater than or equal to 0 and less than or equal to pi/2, and-pi is greater than or equal to pi and less than or equal to pi, theta and etaRespectively, the pitch angle and the azimuth angle of the incident signal. When is fixed withAnd when theta belongs to (-pi/2, pi/2), the polarization guide vector of each array element can be simplified as follows:
in this application, however, the composition ofThe electric dipoles of the same type have the same polarization response to the space electromagnetic wave, so that the M × 1-dimensional polarization steering vector a of the distributed polarization sensitive array shown in FIG. 2 can be obtainedpComprises the following steps:
ap＝[px,py,px,…,py]T
and, when the orthogonal electric dipole pairs are separated and distributed in an alternating manner as shown in fig. 2, each array element component has a constant polarization response to the incident electromagnetic wave. However, due to the fact that the components of each array element have different spatial responses caused by the spatial phase delay of the incident electromagnetic wave reaching each array element, the spatial guide vector a of the distributed polarization sensitive arraysIs represented as follows:
wherein, in the formula,
and, according to the polarization steering vector apAnd a space vector asAnd solving the Hadamard product of the two to obtain a combined polarization domain-space domain steering vector a, namely:
further, for K electromagnetic signals incident in space, the joint polarization domain-space domain steering vector matrix a can be expressed as: [ a ] A1,a2,…,ak]
After step 103 is completed, step 104 is performed: and constructing a signal subspace based on the covariance matrix.
In the specific implementation process, the eigenvalue decomposition is performed on the covariance matrix, and the following results are obtained:
wherein λ isiIs a characteristic value, uiIs a feature vector. Further, the eigenvalues are sorted according to the magnitude of the value, and a subspace formed by eigenvectors corresponding to the first K large eigenvalues is a signal subspace ESThe subspace formed by the eigenvectors corresponding to the remaining M-K eigenvalues is the noise subspace ENSpecifically, the method comprises the following steps:
ES＝[u1,u2,…uK]；
EN＝[uK+1,…uM]。
after completing step 104, step 105 is performed:
and dividing the signal subspace formed by the matrix into a first sub-matrix and a second sub-matrix according to the front-back sequence of the row number in the matrix.
Specifically, matrix ESIs divided into a first sub-matrix ES1Will matrix ESIs divided into a second sub-matrix ES2. Due to ESIs expanded into the same subspace as the column vectors of the steering vector matrix A of the incident signal, i.e. span { E }SThe sub-array may be divided into sub-arraysTo obtain ES2＝ES1ψ, wherein ψ reflects a rotation invariant characteristic between signal subspaces of reception signals of two sub arrays as a rotation invariant characteristic parameter.
After completing step 105, step 106 is performed: and obtaining the estimation of the arrival angle and the polarization parameter of the incident signal by utilizing the rotation invariance between the first sub-matrix and the second sub-matrix.
Specifically, step 106 includes:
obtaining rotation invariant characteristic parameters between the first sub-matrix and the second sub-matrix based on a total least square method rotation invariant method;
performing characteristic value decomposition on the rotation invariant characteristic parameters to obtain characteristic values and characteristic vectors corresponding to the characteristic values;
based on the eigenvalues and the eigenvectors, estimates of the angle of arrival and polarization parameters of the incident signal are obtained.
Further, obtaining estimates of the angle of arrival and polarization parameters of the incident signal based on the eigenvalues and the eigenvectors, comprising:
obtaining an angle of arrival of the incident signal based on the feature value;
obtaining an estimate of a first director sub-matrix corresponding to the first sub-matrix based on a matrix of eigenvalues comprised of eigenvalues and a matrix of eigenvectors comprised of eigenvectors;
based on the estimate of the first steering vector submatrix, an estimate of a polarization phase angle and an estimate of a polarization phase difference of the incident signal are obtained.
In a specific implementation process, a total least square method rotation invariant method (TLS-ESPRIT) method is used to obtain a rotation invariant characteristic parameter psi, and then the psi is subjected to eigenvalue decomposition to obtain K eigenvalues and corresponding eigenvectors, wherein the eigenvalue is specifically the eigenvalueThe matrix of eigenvalues is denoted as phi and the matrix of eigenvectors is denoted as V, since the eigenvalues are related to the angle of arrival of the incident signal, use is made ofAn estimate of the angle of arrival of the incident signal can be obtained:
wherein, angle () is the operation of phase angle, arcsin () is the operation of inverse sine.
Further, V, ESAnd A has the following relationship: eS＝AV-1In the same way, the following can be obtained: eS1＝A1V-1，ES2＝A2V-1Wherein A is1And A2Is a director quantum matrix divided based on the director matrix a in the same division as the sub-arrays. Wherein, in addition, A2＝A1φ
Then, using phi and V obtained by decomposing phi by eigenvalue, a first director quantum matrix A can be obtained1Namely:
specifically, the method comprises the following steps:
then, a guide vector model of the distributed polarization sensitive array is utilized, so thatAll odd rows of (A) form a matrixByAll even rows of (a) form a matrixFurther, willEach element of (1) andis divided by the corresponding elements to obtain a matrixSpecifically, the method comprises the following steps:
according to a matrixAn estimate of the polarization parameters can be obtained:
wherein,for an estimate of the polarization phase angle of the kth incident signal,for an estimate of the polarization phase difference of the kth incoming signal,is a matrixThe absolute value of the k-th column element of (c),is a matrixThe phase angle value of the kth column element of (1). Because the eigenvalues and the corresponding eigenvectors are in one-to-one correspondence in the estimation of the polarization parameters, the arrival angle of the incident signal and the polarization parameters are automatically paired, an additional parameter pairing process is not needed, and the complexity and the calculation amount of parameter joint estimation are reduced.
To improve the accuracy of the polarization parameter estimation, each polarization parameter may be determinedAnd then, obtaining an average value of the polarization parameter estimation values.
In the present application, under the condition that the number of channels is the same and M is 12, comparing the performance of the distributed polarization sensitive array-based parameter joint estimation method of the present invention with the performance of the common-point time polarization sensitive array in the prior art, obtaining a variation curve of the root mean square error with the signal-to-noise ratio about the angle of arrival estimation as shown in fig. 3, and giving CRB as a comparison, it can be seen that the performance of the distributed polarization sensitive array-based parameter joint estimation method of the present invention is significantly higher than that of the common-point time polarization sensitive array-based parameter estimation method in the prior art.
In the present application, when the number of array elements is guaranteed to be the same and M is 12, the performance of the scalar array in the prior art is compared with the performance of the distributed polarization sensitive array based parameter joint estimation method of the present invention, so as to obtain a variation curve graph of the performance of the angle of arrival estimation along with the signal to noise ratio as shown in fig. 4, and CRB is given as a comparison.
Claims (7)
1. A parameter joint estimation method based on a distributed polarization sensitive array is characterized by comprising the following steps:
dispersing and placing each array element component in a polarization sensitive array formed by electric dipole pairs in a space to form a distributed polarization sensitive array;
dividing the distributed polarization sensitive array into a first sub-array and a second sub-array according to the front and back sequence of the array elements;
after the distributed polarization sensitive array receives an incident signal, obtaining a covariance matrix of the incident signal;
constructing a signal subspace based on the covariance matrix;
dividing the signal subspace formed by the matrix into a first sub-matrix and a second sub-matrix according to the front and back sequence of the row number in the matrix;
and obtaining the estimation of the arrival angle and the polarization parameter of the incident signal by utilizing the rotation invariance between the first sub-matrix and the second sub-matrix.
2. The method for jointly estimating parameters based on the distributed polarization sensitive array according to claim 1, wherein the components of each array element in the polarization sensitive array formed by electric dipole pairs are dispersedly placed in space, specifically:
in the established coordinate system comprising a horizontal axis and a vertical axis, the array elements in the polarization sensitive array are distributed on the vertical axis at the same intervals, and the electric dipoles are alternately arranged in the direction parallel to the horizontal axis and the direction parallel to the vertical axis.
3. The parameter joint estimation method based on the distributed polarization sensitive array as claimed in claim 1, wherein when the distributed polarization sensitive array includes M array elements, the distributed polarization sensitive array is divided into a first sub-array and a second sub-array according to the front and back order of the array elements, specifically:
and dividing the front M-2 array elements of the distributed polarization sensitive array into a first sub-array, and dividing the rear M-2 array elements of the distributed polarization sensitive array into a second sub-array.
4. The method for jointly estimating parameters based on a distributed polarization sensitive array according to claim 1, wherein the constructing a signal subspace based on the covariance matrix specifically comprises:
and carrying out eigenvalue decomposition on the covariance matrix to construct a signal subspace.
5. The method for jointly estimating parameters based on distributed polarization-sensitive arrays according to claim 3, wherein the signal subspace formed by the matrices is divided into a first sub-matrix and a second sub-matrix according to the sequence of the number of rows in the matrices, specifically:
the first M-2 rows of the matrix in the partitioned signal subspace constitute a first sub-matrix and the last M-2 rows of the matrix in the partitioned signal subspace constitute a second sub-matrix.
6. The method for jointly estimating parameters based on a distributed polarization sensitive array according to claim 3, wherein the obtaining of the estimation of the angle of arrival and the polarization parameters of the incident signal by using the rotational invariance between the first sub-matrix and the second sub-matrix specifically comprises:
obtaining rotation invariant characteristic parameters between the first sub-matrix and the second sub-matrix based on a total least square method rotation invariant method;
performing characteristic value decomposition on the rotation invariant characteristic parameters to obtain characteristic values and characteristic vectors corresponding to the characteristic values;
based on the eigenvalues and the eigenvectors, estimates of the angle of arrival and polarization parameters of the incident signal are obtained.
7. The method of claim 6, wherein the obtaining of the estimation of the angle of arrival and the polarization parameter of the incident signal based on the eigenvalue and the eigenvector specifically comprises:
obtaining an angle of arrival of the incident signal based on the feature value;
obtaining an estimate of a first director sub-matrix corresponding to the first sub-matrix based on a matrix of eigenvalues comprised of eigenvalues and a matrix of eigenvectors comprised of eigenvectors;
based on the estimate of the first steering vector submatrix, an estimate of a polarization phase angle and an estimate of a polarization phase difference of the incident signal are obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410366085.5A CN104122533A (en) | 2014-07-29 | 2014-07-29 | Joint parameter estimation method based on distributed polarization sensitive array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410366085.5A CN104122533A (en) | 2014-07-29 | 2014-07-29 | Joint parameter estimation method based on distributed polarization sensitive array |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104122533A true CN104122533A (en) | 2014-10-29 |
Family
ID=51768022
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410366085.5A Pending CN104122533A (en) | 2014-07-29 | 2014-07-29 | Joint parameter estimation method based on distributed polarization sensitive array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104122533A (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104375115A (en) * | 2014-11-04 | 2015-02-25 | 中国航天科工集团第三研究院第八三五七研究所 | Polarization sensitive array based non-circular signal DOA and polarization parameter joint estimation method |
| CN104537248A (en) * | 2015-01-04 | 2015-04-22 | 哈尔滨工程大学 | Information source number estimation method for polarization sensitive array |
| CN105334489A (en) * | 2015-12-15 | 2016-02-17 | 吉林大学 | Method for multi-parameter joint estimation of distributed type electromagnetic vector sensor array |
| CN105425205A (en) * | 2015-11-03 | 2016-03-23 | 天津津航计算技术研究所 | High-resolution circular and non-circular signal mixed incidence DOA estimation method |
| CN105572644A (en) * | 2016-01-13 | 2016-05-11 | 电子科技大学 | Polarization sensitive FDA radar, and wave beam forming and apparatus method of the same |
| CN106872933A (en) * | 2017-01-22 | 2017-06-20 | 南京航空航天大学 | A kind of electromagnetic vector sensor array element posture position 6 DOF error self-correcting method |
| CN109983359A (en) * | 2016-09-22 | 2019-07-05 | 弗劳恩霍夫应用研究促进协会 | For estimating the device and correlation method of arrival direction |
| CN110320489A (en) * | 2019-06-26 | 2019-10-11 | 中国电子科技集团公司第三十八研究所 | A kind of distribution isomorphism face battle array combines the balancing method and system for estimating angular accuracy |
| CN111698014A (en) * | 2020-06-11 | 2020-09-22 | 上海磐启微电子有限公司 | Antenna array system |
| CN114113808A (en) * | 2021-11-22 | 2022-03-01 | 杭州电子科技大学 | DOA-polarization information joint estimation method based on incomplete electric vector sensor |
| CN116995451A (en) * | 2023-09-27 | 2023-11-03 | 成都金支点科技有限公司 | Polarization sensitive array arranging method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007225544A (en) * | 2006-02-27 | 2007-09-06 | Advanced Telecommunication Research Institute International | Arrival direction estimating device |
| CN102279387A (en) * | 2011-07-18 | 2011-12-14 | 西安电子科技大学 | Method for estimating target arrival angle of multiple input multiple output (MIMO) radar |
-
2014
- 2014-07-29 CN CN201410366085.5A patent/CN104122533A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007225544A (en) * | 2006-02-27 | 2007-09-06 | Advanced Telecommunication Research Institute International | Arrival direction estimating device |
| CN102279387A (en) * | 2011-07-18 | 2011-12-14 | 西安电子科技大学 | Method for estimating target arrival angle of multiple input multiple output (MIMO) radar |
Non-Patent Citations (2)
| Title |
|---|
| 郑桂妹 等: "三正交分离式极化敏感阵列的波达方向估计", 《电子与信息学报》, vol. 36, no. 5, 31 May 2014 (2014-05-31) * |
| 钟耀霞: "基于极化敏感阵列的空间到达角和极化参数的联合估计", 《中国优秀硕士学位论文全文数据库信息科技辑》, no. 12, 15 December 2011 (2011-12-15) * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104375115A (en) * | 2014-11-04 | 2015-02-25 | 中国航天科工集团第三研究院第八三五七研究所 | Polarization sensitive array based non-circular signal DOA and polarization parameter joint estimation method |
| CN104537248A (en) * | 2015-01-04 | 2015-04-22 | 哈尔滨工程大学 | Information source number estimation method for polarization sensitive array |
| CN105425205B (en) * | 2015-11-03 | 2018-09-28 | 天津津航计算技术研究所 | The circle signal of high resolution mixes incident DOA estimation method with not rounded signal |
| CN105425205A (en) * | 2015-11-03 | 2016-03-23 | 天津津航计算技术研究所 | High-resolution circular and non-circular signal mixed incidence DOA estimation method |
| CN105334489A (en) * | 2015-12-15 | 2016-02-17 | 吉林大学 | Method for multi-parameter joint estimation of distributed type electromagnetic vector sensor array |
| CN105572644A (en) * | 2016-01-13 | 2016-05-11 | 电子科技大学 | Polarization sensitive FDA radar, and wave beam forming and apparatus method of the same |
| CN105572644B (en) * | 2016-01-13 | 2018-02-06 | 电子科技大学 | A kind of polarization sensitive FDA radars and its Beamforming Method and device |
| CN109983359A (en) * | 2016-09-22 | 2019-07-05 | 弗劳恩霍夫应用研究促进协会 | For estimating the device and correlation method of arrival direction |
| CN109983359B (en) * | 2016-09-22 | 2023-08-08 | 弗劳恩霍夫应用研究促进协会 | Apparatus for estimating direction of arrival and corresponding method |
| CN106872933A (en) * | 2017-01-22 | 2017-06-20 | 南京航空航天大学 | A kind of electromagnetic vector sensor array element posture position 6 DOF error self-correcting method |
| CN106872933B (en) * | 2017-01-22 | 2019-08-20 | 南京航空航天大学 | A kind of electromagnetic vector sensor array element posture position 6 DOF error self-correcting method |
| CN110320489A (en) * | 2019-06-26 | 2019-10-11 | 中国电子科技集团公司第三十八研究所 | A kind of distribution isomorphism face battle array combines the balancing method and system for estimating angular accuracy |
| CN111698014A (en) * | 2020-06-11 | 2020-09-22 | 上海磐启微电子有限公司 | Antenna array system |
| CN111698014B (en) * | 2020-06-11 | 2023-09-26 | 上海磐启微电子有限公司 | Antenna array system |
| CN114113808A (en) * | 2021-11-22 | 2022-03-01 | 杭州电子科技大学 | DOA-polarization information joint estimation method based on incomplete electric vector sensor |
| CN114113808B (en) * | 2021-11-22 | 2023-08-15 | 杭州电子科技大学 | DOA-polarization information joint estimation method based on incomplete electric vector sensor |
| CN116995451A (en) * | 2023-09-27 | 2023-11-03 | 成都金支点科技有限公司 | Polarization sensitive array arranging method |
| CN116995451B (en) * | 2023-09-27 | 2023-12-15 | 成都金支点科技有限公司 | Polarization sensitive array arranging method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104122533A (en) | Joint parameter estimation method based on distributed polarization sensitive array | |
| Jinli et al. | Angle estimation using ESPRIT without pairing in MIMO radar | |
| Wen et al. | Angle and mutual coupling estimation in bistatic MIMO radar based on PARAFAC decomposition | |
| CN104698433B (en) | DOA Estimation in Coherent Signal method based on single snapshot data | |
| CN103941221B (en) | Space stretching Electromagnetic Vector Sensor Array method for parameter estimation | |
| CN104849694B (en) | Quaternary number ESPRIT method for parameter estimation of the electromagnetic dipole to array | |
| CN106483493B (en) | A kind of sparse double parallel linear array and estimating two-dimensional direction-of-arrival method | |
| CN104933290B (en) | Multi-parameter Combined estimator quaternary counting method of the double L-shaped orthogonal galvanic couple of stretching to array | |
| CN108120967B (en) | Plane array DOA estimation method and equipment | |
| CN104035069B (en) | Arrowband based on partial correction linear array symmetrically and evenly near-field signals source location method | |
| CN105335615B (en) | A kind of two dimension angular and polarization parameter combined estimation method of low complex degree | |
| CN109254272B (en) | Two-dimensional angle estimation method of concurrent polarization MIMO radar | |
| CN102707258A (en) | Joint estimation method for azimuth angle and elevation angle of signal on basis of L-type sensor array | |
| CN109959891A (en) | The dimensionality reduction spectrum peak search method of Space Angle and polarization parameter in L gusts of electromagnetic vector | |
| CN109738853A (en) | A kind of angle and polarization estimate method of the relatively prime battle array of electromagnetic vector based on rotational invariance | |
| CN104991236A (en) | Monostatic MIMO radar non-circular signal coherent source DOA (Direction Of Arrival) estimation method | |
| CN110412537A (en) | A kind of bistatic MIMO radar angle estimating method | |
| CN112731275B (en) | Zero-change interpolation-based mutual mass array partial polarization signal parameter estimation method | |
| CN108363048B (en) | Block-sparse-based angle estimation method for polarization MIMO radar | |
| CN105203989A (en) | Method for estimating L-shaped array coherent signal direction of arrival | |
| Zhang et al. | DOA estimation in MIMO radar with broken sensors by difference co-array processing | |
| Zheng et al. | Angle and polarization estimation using ESPRIT with polarimetric MIMO radar | |
| CN106877918B (en) | Robust adaptive beam forming method under mutual coupling condition | |
| CN104849693B (en) | Parameter estimation quaternion method for non-concurrent point electromagnetic vector sensor array | |
| Gao et al. | A two-dimensional arrival angles estimation for L-shaped array based on tensor decomposition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20141029 |