CN105652234B - A kind of cyclic space spectrum direction-finding method - Google Patents
A kind of cyclic space spectrum direction-finding method Download PDFInfo
- Publication number
- CN105652234B CN105652234B CN201610100784.4A CN201610100784A CN105652234B CN 105652234 B CN105652234 B CN 105652234B CN 201610100784 A CN201610100784 A CN 201610100784A CN 105652234 B CN105652234 B CN 105652234B
- Authority
- CN
- China
- Prior art keywords
- signal
- cyclic
- frequency
- cycle
- autocorrelation matrix
- 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
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/74—Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention discloses a kind of cyclic space spectrum direction-finding method, including:A) cycle frequency or Conjugate Cyclic frequency of the signal acquisition signal received according to direction-finder antenna array；B) the cycle autocorrelation matrix or Conjugate Cyclic autocorrelation matrix of cycle frequency either Conjugate Cyclic frequency acquisition signal are utilized；C) characteristic vector of signal space and spatial noise is calculated using the cycle autocorrelation matrix or the Conjugate Cyclic autocorrelation matrix；The cyclic space that signal is calculated using the spatial noise characteristic vector is composed；The peak value searching for carrying out the cyclic space spectrum determines the arrival direction of signal.According to the cycle frequency of the signal received, either Conjugate Cyclic frequency recycles autocorrelation matrix or Conjugate Cyclic autocorrelation matrix to energy of the invention to calculate, to realize that the cyclic space composes direction finding, it solves the problems, such as the accuracy decline of low signal interference ratio signal direction-finding under common-frequency multi-signal, the close signal direction-finding in incoming wave orientation, direction finding can be carried out to the same frequency multiple target more than number of antennas.
Description
Technical field
The present invention relates to a kind of cyclic spaces to compose direction-finding method.
Background technology
With the development of electronic technology and the communication technology, application of the direction-finding station technology in civil and military field is more next
It is more extensive, the DF and location of the radiation sources such as radio station may be implemented, Practical significance is very big.Spatial spectrum direction finding technology is to develop for nearly 30 years
A kind of emerging spatial domain signal processing technology got up, mainly grows up on the basis of Wave beam forming and time domain Power estimation
's.Since array aperture limits the Rayleigh of arrival bearing in its Power estimation of breaching over, there is higher direction finding precision,
Through being widely used in practical direction detecting positioning system.
Although existing space composes direction finding technology also can carry out direction finding to multiple homogenous frequency signals, in actual use, especially
When being that the larger or arrival bearing of multiple homogenous frequency signal amplitudes differences is close, it is limited to the low letter of Traditional Space spectrum Direction Finding Algorithm
It does than direction finding ability and with frequency spatial resolving power, direction finding precision declines, and direction finding result is severely impacted.Therefore, it is necessary to
A kind of new direction-finding method is provided to solve the above problems.
Invention content
The purpose of the present invention is to provide a kind of cycle that can utilize signal cycle frequency measurement direction of arrival of signal is empty
Between compose direction-finding method.
To achieve the goals above, the technical solution adopted in the present invention is as follows：
A kind of cyclic space spectrum direction-finding method, includes the following steps:
A) cycle frequency or Conjugate Cyclic frequency of the signal acquisition signal received according to direction-finder antenna array；
B) utilize cycle frequency either the cycle autocorrelation matrix of Conjugate Cyclic frequency acquisition signal or Conjugate Cyclic from
Correlation matrix；
C) signal space and noise are calculated using the cycle autocorrelation matrix or the Conjugate Cyclic autocorrelation matrix
The characteristic vector in space；The cyclic space that signal is calculated using the spatial noise characteristic vector is composed；Carry out the cyclic space
The peak value searching of spectrum determines the arrival direction of signal.
Preferably, the cycle frequency of signal or the acquisition methods of Conjugate Cyclic frequency include the following steps in step a)：
A1 the signal sequence that direction-finder antenna array receives) is denoted as x (n), the conjugation of the signal sequence is denoted as
B1) according to x (n) andProduct sequence calculate related square of spectrum X (f) and related square of conjugation is composed
C1) to X (f) andCarry out spectral line sharpen pretreatment, then by find spectral peak determine signal cycle frequency or
Person's Conjugate Cyclic frequency.
Preferably, the acquisition methods of the cycle autocorrelation matrix of signal are in step b)：Cycle frequency is substituted into auto-correlation
Matrix obtains cycle autocorrelation matrix.
Preferably, the acquisition methods of the Conjugate Cyclic autocorrelation matrix of signal are in step b)：By Conjugate Cyclic frequency generation
Enter autocorrelation matrix and obtains Conjugate Cyclic autocorrelation matrix.
Compared with prior art, a kind of advantageous effect of cyclic space spectrum direction-finding method of the present invention is：The present invention can
According to the cycle frequency of the signal received, either Conjugate Cyclic frequency recycles autocorrelation matrix or Conjugate Cyclic certainly to calculate
Correlation matrix solves low signal interference ratio signal direction-finding under common-frequency multi-signal, incoming wave orientation to realize that the cyclic space composes direction finding
The accuracy decline problem of close signal direction-finding can carry out direction finding to the same frequency multiple target more than number of antennas.
Description of the drawings
Fig. 1 is the step schematic diagram that the cyclic space of the present invention composes direction-finding method；
Fig. 2 is that the cyclic space of the present invention composes direction-finding method and Traditional Space spectrum compares low signal interference ratio signal direction-finding ability
Figure；
Fig. 3 is the spatial resolving power comparison diagram that the cyclic space of the present invention composes direction-finding method and Traditional Space spectrum.
Specific implementation mode
The present invention is further described below with reference to specific embodiment.
Shown in please referring to Fig.1 to Fig.3, the present invention provides a kind of cyclic space spectrum direction-finding method, includes the following steps:
A) cycle frequency or Conjugate Cyclic frequency of the signal acquisition signal received according to direction-finder antenna array；
B) utilize cycle frequency either the cycle autocorrelation matrix of Conjugate Cyclic frequency acquisition signal or Conjugate Cyclic from
Correlation matrix；
C) signal space and noise are calculated using the cycle autocorrelation matrix or the Conjugate Cyclic autocorrelation matrix
The characteristic vector in space；The cyclic space that signal is calculated using the spatial noise characteristic vector is composed；Carry out the cyclic space
The peak value searching of spectrum determines the arrival direction of signal.
Preferably, the cycle frequency of signal or the acquisition methods of Conjugate Cyclic frequency include the following steps in step a)：
A1 the signal sequence that direction-finder antenna array receives) is denoted as x (n), the conjugation of the signal sequence is denoted as
B1) according to x (n) andProduct sequence calculate related square of spectrum X (f) and related square of conjugation is composed
C1) to X (f) andCarry out spectral line sharpen pretreatment, then by find spectral peak determine signal cycle frequency or
Person's Conjugate Cyclic frequency.
Preferably, the acquisition methods of the cycle autocorrelation matrix of signal are in step b)：Cycle frequency is substituted into auto-correlation
Matrix obtains cycle autocorrelation matrix.
Preferably, the acquisition methods of the Conjugate Cyclic autocorrelation matrix of signal are in step b)：By Conjugate Cyclic frequency generation
Enter autocorrelation matrix and obtains Conjugate Cyclic autocorrelation matrix.
Embodiment one：
The aerial array for receiving signal is the L-type array of 9 whip antennas composition, in the hair outside 500 kilometers of receiving point
Exit point arranges 200 watts of short-wave transmitters as target radio station.The selection of target radio station is in the identical transmitting frequency of high power broad-casting station
Signal transmitting is carried out in rate, the signal amplitude for receiving broadcasting station is higher 10dB or more than target radio station, empty to verify cycle
Between raising of the spectrum direction-finding method to the direction finding ability of low signal interference ratio signal.
Target radio station emits amplitude-modulated signal in high power station identical frequency, to the target radio signals of reception and broadcast
The mixed signal of radio signals carries out cycle frequency according to step a) or Conjugate Cyclic frequency calculates, and obtains target radio signals
Cycle frequency or Conjugate Cyclic frequency；Cycle autocorrelation matrix or Conjugate Cyclic auto-correlation square are calculated further according to step b)
Battle array finally utilizes step c) to complete the direction finding to target radio station.The selection of target radio station is enterprising in different high power station frequencies
Row transmitting, repeats the above steps.Typically, in 13.123MHz frequencies, the conjugation for obtaining target radio station amplitude-modulated signal first is followed
Ring frequency is fcycThen Conjugate Cyclic autocorrelation matrix can be obtained in=- 95.36Hz
(wherein N is the data length that direction finding uses, and Fs is signal sampling rate), it is rightCarry out Subspace Decomposition, cyclic space spectrum calculates
And spectrum peak search, it is 264.4 degree that can obtain echo signal arrival direction, with 265 degree of actual signal direction, 0.6 degree of error.
The cyclic space, which can be obtained, using the above process composes direction finding as a result, can be obtained using Traditional Space spectral method simultaneously
Traditional Space composes direction finding result.The error condition such as Fig. 2 institutes for the direction finding result and the realistic objective bearing of station that two methods obtain
Show.It can be seen that Traditional Space composes direction-finding method in low signal interference ratio angle measurement error in 4~12 degree of ranges, and recycle empty
Between spectrum direction-finding method angle measurement error be less than 1 degree.Compared with Traditional Space spectral method, cyclic space spectrum direction-finding method effectively increases
Direction finding precision under low signal interference ratio.
Embodiment two：
Signal array is received equally using the L-type array of 9 whip antennas composition, apart from about 3 kilometers of position cloth of receiving point
Target radio station similar in two orientation is set, true bearing angle is respectively 117 degree and 120 degree.By in series of identical frequency
Emit while two target radio station of upper progress, raising of the test loop spatial spectrum direction-finding method to same frequency spatial resolving power.
Two target radio station emit amplitude-modulated signal and QPSK signals respectively on the same frequency, mixed to two target radio station of reception
It closes signal and cycle frequency or the calculating of Conjugate Cyclic frequency is carried out according to step a), respectively obtain the cycle of two target radio signals
Frequency or Conjugate Cyclic frequency；Cycle autocorrelation matrix or Conjugate Cyclic autocorrelation matrix are calculated further according to step b), finally
The direction finding to two target radio station is completed using step c).Typically, in 16.123MHz frequencies, target radio station 1 is respectively obtained
Amplitude-modulated signal Conjugate Cyclic frequency fcyc1The QPSK signal cycle frequencies f of=- 11.92Hz and target radio station 2cyc2=600Hz, so
After calculate separately corresponding Conjugate Cyclic autocorrelation matrixWith cycle autocorrelation matrix R2, and pass through Subspace Decomposition, cycle
Spatial spectrum calculates and spectrum peak search, and the arrival direction for obtaining two target radio station respectively is 121.0 degree and 117.8 degree.
The cyclic space, which can be obtained, using the above process composes direction finding as a result, can be obtained using Traditional Space spectral method simultaneously
Traditional Space compose direction finding as a result, the obtained direction finding of two methods the results are shown in Figure 3.It can be seen that Traditional Space spectral method is not
The two close target radio station of this space angle can be differentiated, have to a direction finding as a result, and cyclic space spectrum direction-finding method
The direction-finding bearing of the two signals can be respectively obtained, is had preferably with frequency spatial resolving power.
In conclusion compared with prior art, the present invention has the following advantages：
1) direction finding precision of low Signal to Interference plus Noise Ratio signal under common-frequency multi-signal is improved；
2) spatial resolving power of the close signal in orientation under common-frequency multi-signal is improved；
3) the signal-selectivity direction finding according to cycle frequency is realized；
It 4) can be to the same frequency multiple target direction finding more than number of antennas.
Schematically the present invention and embodiments thereof are described above, description is not limiting, institute in attached drawing
What is shown is also one of embodiments of the present invention, and actual structure is not limited to this.So if the common skill of this field
Art personnel are enlightened by it, without departing from the spirit of the invention, are not inventively designed and the technical solution
Similar frame mode and embodiment, are within the scope of protection of the invention.
Claims (3)
1. a kind of cyclic space composes direction-finding method, which is characterized in that include the following steps：
A) cycle frequency or Conjugate Cyclic frequency of the signal acquisition signal received according to direction-finder antenna array；
B) the cycle autocorrelation matrix or Conjugate Cyclic auto-correlation of cycle frequency either Conjugate Cyclic frequency acquisition signal are utilized
Matrix；
C) signal space and spatial noise are calculated using the cycle autocorrelation matrix or the Conjugate Cyclic autocorrelation matrix
Characteristic vector；The cyclic space that signal is calculated using the spatial noise characteristic vector is composed；Carry out the cyclic space spectrum
Peak value searching determines the arrival direction of signal.
It is characterized in that, the cycle frequency of signal or the acquisition methods of Conjugate Cyclic frequency include the following steps in step a)：
A1 direction-finder antenna array) is denoted as x (n) by the signal sequence received, the conjugation of the signal sequence is denoted as
B1) according to x (n) andProduct sequence calculate related square of spectrum X (f) and related square of conjugation is composed
C1) to X (f) andCarry out spectral line sharpen pretreatment, then by find spectral peak determine signal cycle frequency or altogether
Yoke cycle frequency.
2. a kind of cyclic space as described in claim 1 composes direction-finding method, which is characterized in that in step b) cycle of signal from
The acquisition methods of correlation matrix are：Cycle frequency substitution autocorrelation matrix is obtained into cycle autocorrelation matrix.
3. a kind of cyclic space as described in claim 1 composes direction-finding method, which is characterized in that the conjugation of signal is followed in step b)
The acquisition methods of ring autocorrelation matrix are：Conjugate Cyclic frequency substitution autocorrelation matrix is obtained into Conjugate Cyclic autocorrelation matrix.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610100784.4A CN105652234B (en) | 2016-02-24 | 2016-02-24 | A kind of cyclic space spectrum direction-finding method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610100784.4A CN105652234B (en) | 2016-02-24 | 2016-02-24 | A kind of cyclic space spectrum direction-finding method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105652234A CN105652234A (en) | 2016-06-08 |
| CN105652234B true CN105652234B (en) | 2018-07-20 |
Family
ID=56488674
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610100784.4A Active CN105652234B (en) | 2016-02-24 | 2016-02-24 | A kind of cyclic space spectrum direction-finding method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105652234B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107144811B (en) * | 2017-05-12 | 2019-08-02 | 电子科技大学 | A kind of cyclic subspace direction-finding method of single channel receiving array signal |
| CN111220941A (en) * | 2019-10-25 | 2020-06-02 | 北京交通大学 | Distributed irregular array short wave direction finding method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006030834A1 (en) * | 2004-09-14 | 2006-03-23 | National University Corporation Hokkaido University | Signal arrival direction deducing device, signal arrival direction deducing method, and signal arrival direction deducing program |
| CN103634038A (en) * | 2013-12-09 | 2014-03-12 | 南京邮电大学 | Multi-antenna based DOA (direction of arrival) estimation and beam forming combined multipath signal receiving method |
| CN104360305A (en) * | 2014-10-24 | 2015-02-18 | 中国航天科技集团公司第五研究院第五一三研究所 | Radiation source direction finding positioning method of uniting compressed sensing and signal cycle stationary characteristics |
| CN105301354A (en) * | 2015-07-12 | 2016-02-03 | 九江学院 | Harmonic signal frequency estimation method in multiplicative and additive noise |
-
2016
- 2016-02-24 CN CN201610100784.4A patent/CN105652234B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006030834A1 (en) * | 2004-09-14 | 2006-03-23 | National University Corporation Hokkaido University | Signal arrival direction deducing device, signal arrival direction deducing method, and signal arrival direction deducing program |
| CN103634038A (en) * | 2013-12-09 | 2014-03-12 | 南京邮电大学 | Multi-antenna based DOA (direction of arrival) estimation and beam forming combined multipath signal receiving method |
| CN104360305A (en) * | 2014-10-24 | 2015-02-18 | 中国航天科技集团公司第五研究院第五一三研究所 | Radiation source direction finding positioning method of uniting compressed sensing and signal cycle stationary characteristics |
| CN105301354A (en) * | 2015-07-12 | 2016-02-03 | 九江学院 | Harmonic signal frequency estimation method in multiplicative and additive noise |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105652234A (en) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3254133B1 (en) | Direction finding using signal power | |
| CN106526589B (en) | A kind of radar target two-dimensional imaging method based on vortex electromagnetic wave | |
| CN105717493B (en) | A kind of passive radar radiation source polarization identification and direction-finding system based on synthesis polarization method | |
| CN106054165B (en) | A method of realizing distributed passive radar target detection | |
| CN102055512B (en) | Circulating self-correlation-based signal phase difference estimation device and method for antenna array | |
| CN101969329B (en) | Signal phase difference estimation device and method in antenna aligning array based on circulating cross-correlation | |
| CN105891817B (en) | A kind of method of distributed passive radar target detection under the conditions of no direct wave | |
| Cuntz et al. | Field test: jamming the DLR adaptive antenna receiver | |
| CN103002576A (en) | Antenna array single base station positioning method based on pulse amplitude ratio fingerprints | |
| CN106909779A (en) | MIMO radar Cramér-Rao lower bound computational methods based on distributed treatment | |
| CN110515038A (en) | It is a kind of based on the adaptive passive location device of unmanned plane-array and implementation method | |
| CN105652234B (en) | A kind of cyclic space spectrum direction-finding method | |
| CN108683274A (en) | A kind of parallel wireless power transfer method and device of the multiple target of frequency division multiple access | |
| CN106324602B (en) | A kind of MIMO Sonar system | |
| Khallaayoun et al. | A blind iterative calibration method for high resolution DOA estimation | |
| CN108828505A (en) | Angle-of- arrival estimation algorithm research and application based on machine learning | |
| Chen et al. | Compressed sensing-based DOA and DOD estimation in bistatic co-prime MIMO arrays | |
| Ma et al. | Passive ranging by low-directivity antennas with quality estimate | |
| Lagunas et al. | UWB joint TOA and DOA estimation | |
| RU2309425C2 (en) | Method of forming calibration data for radio direction finder/ range finder (versions) | |
| RU2510708C1 (en) | Radio-frequency radiation source direction-finding method | |
| Yang et al. | Wideband LFM interference suppression of GNSS receiver based on fractional Fourier transform and power inversion | |
| Shi et al. | Cramér-Rao lower bounds for joint target parameter estimation in FM-based distributed passive radar network with antenna arrays | |
| Erić et al. | An approach for determination of antenna positions in distributed antenna system used for UWB indoor self-localization: Experimental results | |
| Depold et al. | A light weight 3x3 switched polarity URA antenna and receive system for direction finding |
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 |