CN109188347B - Signal polarization mode three-classification method based on MUSIC spectrum - Google Patents
Signal polarization mode three-classification method based on MUSIC spectrum Download PDFInfo
- Publication number
- CN109188347B CN109188347B CN201811183028.8A CN201811183028A CN109188347B CN 109188347 B CN109188347 B CN 109188347B CN 201811183028 A CN201811183028 A CN 201811183028A CN 109188347 B CN109188347 B CN 109188347B
- Authority
- CN
- China
- Prior art keywords
- polarization
- music spectrum
- music
- signal
- calculating
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
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)
- Optical Recording Or Reproduction (AREA)
Abstract
The invention discloses a signal polarization mode three-classification method based on a MUSIC spectrum, which specifically comprises the following steps: calculating the MUSIC spectrum of the received signal of the polarization sensitive array; extracting features based on MUSIC spectrum; and designing a polarization mode three classifier. The invention can realize high-precision angle estimation and polarization mode three-classification of the polarization sensitive array receiving signals on the premise of not increasing the testing workload.
Description
Technical Field
The invention relates to a polarization mode three-classification method for polarization sensitive array received signals, in particular to a polarization mode three-classification method for polarization sensitive array received signals based on MUSIC spectrums.
Background
Recently, the polarization sensitive array is limited by factors such as wide attention of researchers and tested workload, and generally only performs a plane scan test that the polarization sensitive array receives two linear polarization signals of horizontal polarization and vertical polarization in a test stage with a cover, stores amplitude and phase information received by each array element, and then synthesizes horizontal and vertical polarization guide vectors of the polarization sensitive array by using the information.
However, the polarization sensitive array is usually sensitive to the polarization mode of the received signal, i.e. the polarization sensitive array only carries the guide vectors of the horizontal and vertical polarization modes, and for the received signal with the polarization mode matched with the guide vector, the angle estimation with higher precision can be respectively realized through the spatial spectrum estimation based on the MUSIC algorithm; for the received signals with the polarization mode not matched with the guide vector, the angle estimation with higher precision is difficult to realize through the spatial spectrum estimation based on the MUSIC algorithm.
Therefore, on the premise of not increasing the test workload, if the closed loop of the polarization sensitive array angle estimation system is to be completed, namely for the received signals of different polarization modes, the angle estimation with higher precision can be realized through the spatial spectrum estimation (based on the MUSIC algorithm), a polarization mode three-classification method of the polarization sensitive array received signals is designed, and the polarization mode three-classification of the polarization sensitive array received signals is realized, namely, the polarization mode class 1 (assuming H-type) γ,1 ): linear polarization angle γ =0 ° "," polarization mode class 2 (assuming H) γ,2 ): linear polarization angle γ ≠ 0 ° and γ ≠ 90 ° "and" polarization system class 3 "(assuming that H is γ,3 ): the linear polarization angle γ =90 ° ", is a problem to be solved.
Disclosure of Invention
In view of the above technical problems, an object of the present invention is to provide a polarization-sensitive array received signal polarization mode three-classification method based on MUSIC spectra, so as to implement polarization mode three-classification of polarization-sensitive array received signals.
A signal polarization mode three-classification method based on MUSIC spectrum specifically comprises the following steps:
s1, calculating an MUSIC spectrum of a polarization sensitive array receiving signal;
s2, extracting characteristics based on the MUSIC spectrum;
and S3, designing a polarization mode three classifier.
Further, the step S1 specifically includes:
according to the amplitude of each channel obtained by testing the polarization sensitive array under the conditions of horizontal polarization and vertical polarization respectively:
and phase difference:
wherein N is E Is the number of array elements (number of channels) of the polarization sensitive array, theta andrespectively outputting an incident azimuth angle and a pitch angle of a test result;andrespectively indicate an incident direction ofThe horizontally polarized reception amplitude and the vertically polarized reception amplitude of the time array element i,andrespectively indicate an incident direction ofThe horizontal polarization phase difference and the vertical polarization phase difference of the time array element i and the array element 1 (reference array element).
Respectively calculating horizontal polarization guide vectors a of the polarization sensitive array h And a perpendicular polarization steering vector a v Namely:
and
calculating a signal S with a real linear polarization angle gamma received by the polarization sensitive array γ Respectively based on horizontal polarization steering vectorAnd a perpendicular polarization steering vectorMUSIC spectrum P of γh And P γv The calculation formula is as follows:
Further, the step S2 specifically includes the steps of:
s201, according to a signal S with a real linear polarization angle gamma received by the polarization sensitive array γ MUSI of (A)C spectrum P γh And P γv Extracting statistical characteristics based on the MUSIC spectrum;
s202, respectively calculating the ratio of the MUSIC spectrum peak value to the MUSIC spectrum mean value, respectively calculating the ratio of the MUSIC spectrum standard deviation to the MUSIC spectrum mean value, and respectively calculating the ratio of the MUSIC spectrum peak value to the MUSIC spectrum base;
s203, calculating a spectrum P based on MUSIC γh Characteristic X of γh (S γ ) And based on MUSIC spectrum P γv Characteristic X of γv (S γ ) Obtaining a characteristic X (S) based on the MUSIC spectrum γ )。
Further, the step S201 specifically includes:
respectively calculating MUSIC spectral peak value x γh,1 =max(P γh ) And x γv,1 =max(P γv )。
Further, the step S201 specifically includes:
respectively calculating the average value x of the MUSIC spectrums γh,2 =μ(P γh ) And x γv,2 =μ(P γv )。
Further, the step S201 specifically includes:
respectively calculating the standard deviation x of the MUSIC spectrum γh,3 =σ(P γh ) And x γv,3 =σ(P γv )。
Further, the step S201 specifically includes:
calculating the MUSIC spectrum bases respectively:
and
and
wherein theta is H Andupper bound, theta, of search for azimuth and pitch angles, respectively L Andlower bound, θ, of the search for azimuth and pitch angles, respectively S Andthe search step lengths of the azimuth angle and the pitch angle are respectively.
Further, the step S202 specifically includes:
further, the step S203 specifically includes:
X γh (S γ )={x γh,i :i=1,...,7};
X γv (S γ )={x γv,i :i=1,...,7};
X(S γ )=X γh (S γ )∪X γv (S γ )。
further, the step S3 specifically includes:
the polarization mode three classifiers are H γ (S γ )=C[X(S γ )]。
The invention can realize high-precision angle estimation and polarization mode three-classification of the polarization sensitive array receiving signals on the premise of not increasing the test workload.
Drawings
FIG. 1 is a schematic flow chart of a polarization-sensitive array received signal polarization-based three-classification method according to the present invention.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In this embodiment, taking a polarization sensitive array in a certain structural state as an example, as shown in fig. 1, a polarization sensitive array received signal polarization mode three-classification method based on MUSIC spectrum includes the following specific steps:
in the first step, the MUSIC spectrum of the polarization sensitive array receiving signal is calculated.
According to the amplitude of each channel obtained by testing the polarization sensitive array under the conditions of horizontal polarization and vertical polarization respectively:
and phase difference:
wherein N is E Is the number of array elements (number of channels) of the polarization sensitive array, theta andan incident azimuth angle and a pitch angle for outputting a test result are respectively;andrespectively indicate an incident direction ofThe horizontally polarized reception amplitude and the vertically polarized reception amplitude of the time array element i,andrespectively indicate an incident direction ofThe horizontal polarization phase difference and the vertical polarization phase difference of the time array element i and the array element 1 (reference array element).
Respectively calculating horizontal polarization guide vectors a of the polarization sensitive array h And a perpendicular polarization steering vector a v Namely:
and
calculating a signal S with a real linear polarization angle gamma received by the polarization sensitive array γ Respectively based on horizontally polarized steering vectorsAnd a perpendicular polarization steering vectorMUSIC spectrum P of γh And P γv The calculation formula is as follows:
And secondly, extracting features based on the MUSIC spectrum.
According to a signal S with a real linear polarization angle gamma received by a polarization sensitive array γ MUSIC spectrum P of γh And P γv Extracting statistical characteristics based on the MUSIC spectrum,
respectively calculating MUSIC spectrum peak values x γh,1 =max(P γh ) And x γv,1 =max(P γv ) Where the operator max (·) represents the calculated maximum value.
Respectively calculating the average value x of the MUSIC spectrums γh,2 =μ(P γh ) And x γv,2 =μ(P γv ) Where the operator μ (·) represents the calculated average.
Respectively calculating MUSIC spectrum standardDifference x γh,3 =σ(P γh ) And x γv,3 =σ(P γv ) Where the operator σ (·) represents the calculated standard deviation.
Calculating the MUSIC spectrum bases respectively:
and
and
wherein theta is H Andupper bound, theta, of search for azimuth and pitch angles, respectively L Andlower bound, theta, of search for azimuth and pitch angles, respectively S Andthe search step lengths of the azimuth angle and the pitch angle are respectively.
Respectively calculating the ratio of the peak value of the MUSIC spectrum to the mean value of the MUSIC spectrum according to the statistical characteristics based on the MUSIC spectrumAndrespectively calculating the ratio of the standard deviation of the MUSIC spectrum and the mean value of the MUSIC spectrumAndand respectively calculating the ratio of the MUSIC spectrum peak value to the MUSIC spectrum baseAnd
the calculation is based on the MUSIC spectrum P γh Characteristic X of γh (S γ ) And based on MUSIC spectrum P γv Characteristic X of γh (S γ ) Obtaining a characteristic X (S) based on the MUSIC spectrum γ ),
X γh (S γ )={x γh,i :i=1,...,7};
X γv (S γ )={x γv,i :i=1,...,7};
X(S γ )=X γh (S γ )∪X γv (S γ )。
And thirdly, designing a polarization mode three classifier. According to the feature X (S) based on the MUSIC spectrum γ )=X γh (S γ )∪X γv (S γ ) Design polarization mode three classifiers H γ (S γ )=C[X(S γ )]Three categories of polarization modes of the polarization sensitive array receiving signals, namely 'polarization mode class 1' (supposing H) γ,1 ): linear polarization angle γ =0 ° "," polarization mode class 2 (assuming H) γ,2 ): linear polarization angle γ ≠ 0 ° and γ ≠ 90 ° "and" polarization mode class 3 (assume H) γ,3 ): linear polarization angle γ =90 ° ".
It should be understood that the above embodiments are only examples for clarity of description, and are not limiting. Other variations and modifications will be apparent to persons skilled in the art upon reference to the above description. And are neither required nor exhaustive of all embodiments. And obvious variations or modifications therefrom are intended to be within the scope of the invention.
Claims (9)
1. A signal polarization mode three-classification method based on MUSIC spectrum is characterized by comprising the following steps:
s1, calculating an MUSIC spectrum of a polarization sensitive array receiving signal;
s2, extracting characteristics based on the MUSIC spectrum;
s3, designing a polarization mode three classifier;
the step S1 specifically includes:
according to the amplitude of each channel obtained by testing the polarization sensitive array under the conditions of horizontal polarization and vertical polarization respectively:
and phase difference:
wherein N is E Is the array element number of the polarization sensitive array, theta andan incident azimuth angle and a pitch angle for outputting a test result are respectively;andrespectively indicate an incident direction ofThe horizontally polarized reception amplitude and the vertically polarized reception amplitude of the time array element i,andrespectively indicate an incident direction ofThe horizontal polarization phase difference and the vertical polarization phase difference of the time array element i and the array element 1;
respectively calculating horizontal polarization guide vectors a of the polarization sensitive array h And a perpendicular polarization steering vector a v Namely:
and
calculating a signal S with a real linear polarization angle gamma received by the polarization sensitive array γ Respectively based on horizontally polarized steering vectorsAnd a perpendicular polarization steering vectorMUSIC spectrum P of γh And P γv The calculation formula is as follows:
2. The signal polarization mode three-classification method according to claim 1, wherein the step S2 specifically comprises the steps of:
s201, receiving a signal S with a true linear polarization angle gamma according to a polarization sensitive array γ MUSIC spectrum P of γh And P γv Extracting statistical characteristics based on the MUSIC spectrum;
s202, respectively calculating the ratio of the MUSIC spectrum peak value to the MUSIC spectrum mean value, respectively calculating the ratio of the MUSIC spectrum standard deviation to the MUSIC spectrum mean value, and respectively calculating the ratio of the MUSIC spectrum peak value to the MUSIC spectrum base according to the statistical characteristics based on the MUSIC spectrum;
s203, calculating the spectrum P based on the MUSIC γh Characteristic X of γh (S γ ) And based on MUSIC spectrum P γv Characteristic X of γv (S γ ) Obtaining a characteristic X (S) based on the MUSIC spectrum γ )。
3. The method for three classification of signal polarization modes according to claim 2, wherein the step S201 specifically comprises:
respectively calculating MUSIC spectral peak value x γh,1 =max(P γh ) And x γv,1 =max(P γv )。
4. The method for three classification of signal polarization modes according to claim 3, wherein the step S201 specifically comprises:
respectively calculating the average value x of the MUSIC spectrums γh,2 =μ(P γh ) And x γv,2 =μ(P γv )。
5. The method for three classification of signal polarization modes according to claim 4, wherein the step S201 specifically comprises:
respectively calculating the standard deviation x of the MUSIC spectrum γh,3 =σ(P γh ) And x γv,3 =σ(P γv )。
6. The method for three classification of signal polarization modes according to claim 5, wherein the step S201 specifically comprises:
calculating the MUSIC spectrum bases respectively:
and
and
7. The method for three classification of signal polarization modes according to claim 6, wherein the step S202 specifically includes:
8. the method for three classification of signal polarization modes according to claim 7, wherein the step S203 specifically comprises:
X γh (S γ )={x γh,i :i=1,...,7};
X γv (S γ )={x γv,i :i=1,...,7};
X(S γ )=X γh (S γ )∪X γv (S γ )。
9. the method for three classification of signal polarization modes according to claim 8, wherein the step S3 specifically comprises:
the polarization mode three classifiers are H γ (S γ )=C[X(S γ )]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811183028.8A CN109188347B (en) | 2018-10-11 | 2018-10-11 | Signal polarization mode three-classification method based on MUSIC spectrum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811183028.8A CN109188347B (en) | 2018-10-11 | 2018-10-11 | Signal polarization mode three-classification method based on MUSIC spectrum |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109188347A CN109188347A (en) | 2019-01-11 |
CN109188347B true CN109188347B (en) | 2022-11-22 |
Family
ID=64947665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811183028.8A Active CN109188347B (en) | 2018-10-11 | 2018-10-11 | Signal polarization mode three-classification method based on MUSIC spectrum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109188347B (en) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010121975A (en) * | 2008-11-17 | 2010-06-03 | Advanced Telecommunication Research Institute International | Sound-source localizing device |
CN102122352B (en) * | 2011-03-01 | 2012-10-24 | 西安电子科技大学 | Characteristic value distribution statistical property-based polarized SAR image classification method |
CN103941222B (en) * | 2014-03-12 | 2016-02-03 | 西安电子科技大学 | Based on the method for parameter estimation of Rayleigh entropy dimensionality reduction MUSIC algorithm |
CN103941221B (en) * | 2014-03-12 | 2016-06-08 | 西安电子科技大学 | Space stretching Electromagnetic Vector Sensor Array method for parameter estimation |
CN104766093B (en) * | 2015-04-01 | 2018-02-16 | 中国科学院上海微系统与信息技术研究所 | A kind of acoustic target sorting technique based on microphone array |
CN104898085B (en) * | 2015-05-14 | 2017-11-17 | 电子科技大学 | A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation |
CN108459706A (en) * | 2018-01-24 | 2018-08-28 | 重庆邮电大学 | Wi-Fi gesture identification methods based on relative movement orbit tracking |
-
2018
- 2018-10-11 CN CN201811183028.8A patent/CN109188347B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN109188347A (en) | 2019-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dai et al. | A sparse representation method for DOA estimation with unknown mutual coupling | |
Ye et al. | On the resiliency of MUSIC direction finding against antenna sensor coupling | |
EP2628316B1 (en) | Apparatus and method for deriving a directional information and computer program product | |
Zheng et al. | Classification and localization of mixed near-field and far-field sources using mixed-order statistics | |
CN103954950A (en) | Direction-of-arrival estimation method based on sample covariance matrix sparsity | |
CN106646350B (en) | A kind of modification method when each channel amplitude gain of single vector hydrophone is inconsistent | |
Zhang et al. | Compensation for the mutual coupling effect in the ESPRIT direction finding algorithm by using a more effective method | |
CN111308424A (en) | Transformer substation equipment audible sound source positioning method based on summation and MUSIC combined algorithm | |
CN113253194A (en) | Broadband arrival angle and polarization combined measurement method based on sparse representation | |
CN110045322A (en) | A kind of shortwave automatic direction finding method based on high-resolution direction finding sonagram intelligent recognition | |
Liu et al. | Two dimensional direction of arrival estimation for co-prime L-shaped array using sparse reconstruction | |
CN110018438A (en) | A kind of Wave arrival direction estimating method and device | |
CN109188347B (en) | Signal polarization mode three-classification method based on MUSIC spectrum | |
CN108226855A (en) | The not rounded joint parameter estimation method in far and near field in the case of mutual coupling | |
Yang et al. | A convex geometry-based blind source separation method for separating nonnegative sources | |
Ramamonjy et al. | Source localization and identification with a compact array of digital mems microphones | |
Qi et al. | DOA estimation and self-calibration algorithm for multiple subarrays in the presence of mutual coupling | |
CN108594165A (en) | A kind of narrow band signal Wave arrival direction estimating method based on expectation-maximization algorithm | |
US20240012083A1 (en) | Method and apparatus for measuring directions of arrival of multiple sound sources | |
Lekhovytskiy et al. | Rapidly convergent “superresolving” direction finders of noise radiation sources in adaptive arrays | |
Mahata | A subspace algorithm for wideband source localization without narrowband filtering | |
Cai | A normalized spatial spectrum for DOA estimation with uniform linear arrays in the presence of unknown mutual coupling | |
Zhang et al. | Low-complexity gridless 2D harmonic retrieval via decoupled-ANM covariance reconstruction | |
Meng et al. | A time picking method based on spectral multimanifold clustering in microseismic data | |
Wisudawan et al. | Two dimensional angle of arrival estimation using minimum sparse ruler based rectangular array of antennas |
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 |