CN107045131A - Tensor resolution satellite navigation anti-interference method - Google Patents
Tensor resolution satellite navigation anti-interference method Download PDFInfo
- Publication number
- CN107045131A CN107045131A CN201710371234.0A CN201710371234A CN107045131A CN 107045131 A CN107045131 A CN 107045131A CN 201710371234 A CN201710371234 A CN 201710371234A CN 107045131 A CN107045131 A CN 107045131A
- Authority
- CN
- China
- Prior art keywords
- signal
- array
- tensor
- represent
- represented
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention provides a kind of tensor resolution satellite navigation anti-interference method, realize that space is structured the formation using vector sensor, Wave beam forming is carried out using tensor resolution technology, the higher-dimension characteristic that make use of tensor has in itself is decomposed to the polarization sensitive array for meeting linearshiftinvariance, some low-dimensional submatrixs are obtained, Wave beam forming is then carried out respectively to each submatrix.The present invention effectively reduces system operations amount, lifting system real-time is conducive to Project Realization relative to global beam formation technology by reducing finding the inverse matrix dimension.
Description
Technical field
The present invention relates to a kind of satellite navigation anti-interference method, belong to field of satellite navigation.
Background technology
In recent years, GNSS is widely used in military and civilian field and plays military work important all the more
With, but due to terminal receiving device apart from satellite farther out, satellite-signal reach array received face when it is extremely faint, easily by
The interference of the unexpected electromagnetic signal in space, while the AF panel nargin of navigation neceiver itself is again very small, therefore, based on each
The self-adaptive restraining technology for planting standard interference obtains extensive research with applying.
The array type of satellite navigation jamproof system can be divided into scalar array and vector array in terms of antenna performance, from
It can be divided into periodic arrangement array and random arrangement array in terms of antenna element mode.The cycle being used current Anti-Jamming Technique more
Property arrangement scalar array, and array element spacing can not be more than satellite signal wave lengths half, if spacing be more than half-wavelength, wave beam side
Graing lobe will be produced to figure, the constraint limits arrangement of the scalar array element in space.
In addition, Wave beam forming is carried out more than current Anti-Jamming Technique by the way of based on Matrix Solving, with array number
The increase of amount, the dimension of matrix is consequently increased, for correspondence large scale array, currently conventional some FPGA or dsp chip
It will be difficult to bear its matrix operation amount.Therefore engineer applied of the large scale array in satellite navigation jamproof system is limited.
The content of the invention
In order to overcome the deficiencies in the prior art, the present invention provides a kind of anti-interference side of satellite navigation based on tensor resolution
Method, realizes that space is structured the formation using vector sensor, breaks through the limitation of scalar array elements spacing.
The technical solution adopted for the present invention to solve the technical problems comprises the following steps：
A) the N roads radio frequency analog signal arrived to antenna array receiver carries out low noise amplification, then downconverted to obtain N roads
Analog intermediate frequency signal；
B) A/D samplings are carried out to N roads analog intermediate frequency signal, obtains N roads digital intermediate frequency signal；
C) using initial aerial array as referential array, its element number of array NsNot less than 1, G translation is carried out to it will
Obtain final spatial domain array；OrderThe position vector of n-th of vector sensor in referential array is represented, is usedThe spatial direction of desired signal is represented, θ and φ are represented respectively
Desired signal come to azimuth and the angle of pitch, the wavelength of desired signal is represented with λ, then spatial domain of the referential array to desired signal
Steering vector
D) represent that the polarization auxiliary angle and polarization phases of desired signal are poor respectively with γ and η, then array is to desired signal
Polarizing field steering vector
In above formulaPolarization selection matrix is represented,Three dipole antenna arrays are represented,Represent
Full electromagnetic vector sensor；
E) useRepresent the g times translation in l-th of spatial translation vector, wherein g=1,2 ..., G, l=1,2 ...,
Ng, then g-th of translation invariant steering vector
N in formulagRepresent after being translated through the g times, be translated the increased multiple of array；
F) Wave beam forming weights are initialized, w is made1=ap, w2=as, wg+2=ag；
G) the signal tensor of the continuous K time domain discrete sampled point of array receivedWherein NpRepresent
Polarization sensitive array element is subjected to polarization components number, then signal tensor1 ..., n ' -1, n '+1 ..., G+2 } mould expansionN '=1,2 ..., G+2 in formula, ×qRepresent the q moulds of tensor
Product；
H) to Xn′Covariance matrix solution is carried out,
I) the covariance matrix R of step h) generations is utilizedn′And step c), d), e) generation steering vector as、apAnd ag
Wave beam forming weights are updated,
J) signal tensorThe output signal obtained after weightedThen signal output work(
Rate p=| | y | |2；
K) calculating is the output signal power before right value update, and calculates the difference of performance number obtained by the performance number and step j)
Value, if the absolute value of the difference is less than setting thresholding ε, continues executing with next step, otherwise returns to step g)；
L) automatic growth control is carried out to output signal y, controls the number of significant digit of output data, filter out unnecessary symbol
Position；
M) navigation information resolving is carried out to the signal that step l) is exported, obtains current location information.
The beneficial effects of the invention are as follows：Wave beam forming is carried out using tensor resolution technology, make use of what tensor had in itself
Higher-dimension characteristic is decomposed to the polarization sensitive array for meeting linearshiftinvariance, obtains some low-dimensional submatrixs, then right
Each submatrix carries out Wave beam forming respectively, and this method has relative to global beam formation technology by reducing finding the inverse matrix dimension
Effect reduces system operations amount, and lifting system real-time is conducive to Project Realization.
Brief description of the drawings
Fig. 1 is the constant polarization sensitive array schematic diagram of linear translation；
Fig. 2 is the satellite navigation anti-interference method flow chart based on tensor resolution；
Fig. 3 is the constant polarization sensitive array schematic diagram of 8 array element linear translations.
Embodiment
The present invention is further described with reference to the accompanying drawings and examples, and the present invention includes but are not limited to following implementations
Example.
The present invention proposes a kind of satellite navigation anti-interference method based on tensor resolution, is realized using vector sensor empty
Between structure the formation, the limitation of scalar array elements spacing is broken through, in addition, in order to reduce system operations amount, the present invention is to high-dimensional big rule
Mode array carries out tensor resolution, improves the real-time of system, to make vector array meet tensor separability, this programme uses battle array
Column space layout meets linear translation invariant feature, as shown in Figure 1.
The present invention comprises the following steps (corresponding flow is as shown in Figure 2)：
A) analog down：The N roads radio frequency analog signal arrived to antenna array receiver carries out low noise amplification, then passes through down
N roads analog intermediate frequency signal is obtained after frequency conversion.
B) A/D samples：A/D samplings are carried out to N roads analog intermediate frequency signal, N roads digital intermediate frequency signal is obtained.
C) referential array spatial domain steering vector is generated：Referential array refers to the original array without linear translation, its array element
Number is not less than 1, and repeatedly translation is carried out to it to obtain final spatial domain array, is made firstRepresent n-th of arrow in referential array
The position vector of quantity sensor, n=1,2 ..., Ns, useThe expression phase
Hope the spatial direction of signal, θ and φ represent respectively desired signal to azimuth and the angle of pitch, represented to expect letter with λ
Number wavelength.Then referential array is represented by the spatial domain steering vector of desired signal
D) array polarizing field steering vector is generated：Represent the polarization auxiliary angle and polarization phase of desired signal respectively with γ and η
Potential difference, then array the polarizing field steering vector of desired signal is represented by
In above formulaPolarization selection matrix is represented,Three dipole antenna arrays are represented,Represent
Full electromagnetic vector sensor.
E) translation invariant steering vector is generated：WithRepresent l-th of the spatial translation vector in the g times translation, wherein g=
1,2 ..., G, l=1,2 ..., Ng, G representation spaces translation number of times, then g-th of translation invariant steering vector be expressed as
N in formulagRepresent after being translated through the g times, be translated the increased multiple of array, such as N3=5, then it represents that array is through the 3rd
After secondary translation, array elements number will be changed into the 3rd time translation before 5 times.
F) Wave beam forming weights are initialized：Make w1=ap, w2=as, wg+2=ag, g=1,2 ..., G.
G) docking collection of letters tensor carries out expansion operation：The signal tensor of the continuous K time domain discrete sampled point of array receivedWherein NpRepresent that polarization sensitive array element is subjected to polarization components number,Represent implication and step
Suddenly e) describe identical, then signal tensor1 ..., n ' -1, n '+1 ..., G+2 } mould expansion be represented by
N '=1,2 ..., G+2 in formula, ×qRepresent the q moulds product of tensor.
H) covariance matrix is generated：To the signal tensor of g) middle generation1 ..., n ' -1, n '+1 ..., G+2 } mould exhibition
Open Xn′Covariance matrix solution is carried out, i.e.,
N '=1,2 in formula ..., G+2.
I) right value update：Utilize the covariance matrix R of h) middle generationn′And steering vector a c), d), e) generateds、apAnd ag
Wave beam forming weights are updated,
J) output signal power is counted：Signal tensorThe output signal type obtained after weighted is represented by
Then output power signal is represented by
P=| | y | |2 (8)
K) difference judges：Output signal power before being performed according to formula (7), (8) calculation formula (6), i.e. right value update
Preceding output signal power, and the performance number and j) difference of middle gained performance number are calculated, used if the absolute value of the difference is less than
Family defines thresholding ε, then continues executing with next step, if more than or equal to ε, returning to the g) step, ε values are general between 1 to 10.
L) automatic growth control：The output signal y obtained to performing formula (7) carries out automatic growth control, control output
The number of significant digit of data, filters out unnecessary sign bit, reaches the effect of limitation output amplitude.
M) navigation information is resolved：Navigation information resolving is carried out to the signal of the l) step output, current location information is obtained.
Embodiments of the invention propose a kind of satellite navigation anti-interference method based on tensor resolution, its anti-interference method
If flow chart as shown in Fig. 2 referential array is a full electromagnetic vector sensor at rectangular coordinate system in space origin, is passed through
Cross the translation of 3 sublinears and obtain complete space array, andSpace array model with gps satellite as shown in figure 3, believed
Exemplified by number, its center frequency point be 1.57542GHz, it is assumed that it is circularly polarised wave, and its space come to for (θ, φ)=(120 °,
90 °), polarization parameter (γ, η)=(30 °, 0 °), tensor resolution satellite navigation anti-interference method is comprised the following steps that：
Step one：Aerial array receives radio frequency (RF) gps satellite signal and interference signal, Mei Gefang in 8, space orientation
Exportable 6 polarization components of full electromagnetic vector sensor of position, then obtain 48 Lu Zhizhengs after LNA and analog down
Frequency is 46.42MHz analog intermediate frequency signal.
Step 2：The 48 road analog if signals received are sampled by AD sampling modules, obtained in 48 railway digitals
Frequency signal.
Step 3：Spatial domain steering vector is generated according to formula (1)
as=[1-1 1-1 0.67-0.75i-0.67+0.75i 0.67-0.75i-0.67+0.75i]T。
Step 4：Polarizing field steering vector a is generated according to formula (2)p=[- 0.75-0.43-0.50 0.43 0.25
-0.87]T。
Step 5：Translation invariant steering vector a is generated according to formula (3)1=[0-1]T, a2=[0 1]T, a3=[0
0.6661 -0.7458i]T。
Step 6：Initialize Wave beam forming weights：Make w1=ap, w2=as, wg+2=ag, g=1,2,3.
Step 7：The data of continuous 1024 time domain discrete sampled points are counted, obtain receiving the tensor model of signalX is calculated according to formula (4)n′, n '=1,2 ..., 5.
Step 8：Solved according to formula (5) and correspond to Xn′Covariance matrix Rn′。
Step 9：Weights are updated according to formula (6).
Step 10：Signal tensor is weighted according to the weights that formula (7) is obtained with step 9, output signal is obtained,
And solve output signal power.
Step 11：ε values are 6, compare the mould of signal power difference before and after right value update, are returned if its value is more than 6
Step 7, the mould of signal power difference is less than 6 before and after right value update.
Step 12：Automatic growth control is carried out to the output signal for meeting step 11, only retains output data 3
Sign bit, data effective width is set to 8.
Step 13：The signal exported by navigation module to step 12 carries out navigation information resolving, obtains current position
Confidence ceases.
Claims (1)
1. a kind of tensor resolution satellite navigation anti-interference method, it is characterised in that comprise the steps：
A) the N roads radio frequency analog signal arrived to antenna array receiver carries out low noise amplification, then downconverted to obtain N roads intermediate frequency
Analog signal；
B) A/D samplings are carried out to N roads analog intermediate frequency signal, obtains N roads digital intermediate frequency signal；
C) using initial aerial array as referential array, its element number of array NsNot less than 1, G translation is carried out to it will obtain most
Whole spatial domain array；OrderThe position vector of n-th of vector sensor in referential array is represented, is usedThe spatial direction of desired signal is represented, θ and φ are represented respectively
Desired signal come to azimuth and the angle of pitch, the wavelength of desired signal is represented with λ, then spatial domain of the referential array to desired signal
Steering vector
D) represent that the polarization auxiliary angle and polarization phases of desired signal are poor respectively with γ and η, then polarization of the array to desired signal
Domain steering vector
In above formulaPolarization selection matrix is represented,Three dipole antenna arrays are represented,Represent full electricity
Magnetic vector sensor；
E) useRepresent l-th of the spatial translation vector in the g times translation, wherein g=1,2 ..., G, l=1,2 ..., Ng, then
G-th of translation invariant steering vector
N in formulagRepresent after being translated through the g times, be translated the increased multiple of array；
F) Wave beam forming weights are initialized, w is made1=ap, w2=as, wg+2=ag；
G) the signal tensor of the continuous K time domain discrete sampled point of array receivedWherein NpRepresent polarization
Sensitive array element is subjected to polarization components number, then signal tensor1 ..., n ' -1, n '+1 ..., G+2 } mould expansionN '=1,2 ..., G+2 in formula, ×qRepresent the q moulds of tensor
Product；
H) to Xn′Covariance matrix solution is carried out,
I) the covariance matrix R of step h) generations is utilizedn′And step c), d), e) generation steering vector as、apAnd agUpdate
Wave beam forming weights,
J) signal tensorThe output signal obtained after weightedThen output power signal p
=| | y | |2；
K) calculating is the output signal power before right value update, and calculates the difference of performance number obtained by the performance number and step j),
If the absolute value of the difference is less than setting thresholding ε, next step is continued executing with, step g) is otherwise returned to；
L) automatic growth control is carried out to output signal y, controls the number of significant digit of output data, filter out unnecessary sign bit；
M) navigation information resolving is carried out to the signal that step l) is exported, obtains current location information.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710371234.0A CN107045131B (en) | 2017-05-24 | 2017-05-24 | Tensor resolution satellite navigation anti-interference method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710371234.0A CN107045131B (en) | 2017-05-24 | 2017-05-24 | Tensor resolution satellite navigation anti-interference method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107045131A true CN107045131A (en) | 2017-08-15 |
| CN107045131B CN107045131B (en) | 2019-04-19 |
Family
ID=59547457
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710371234.0A Active CN107045131B (en) | 2017-05-24 | 2017-05-24 | Tensor resolution satellite navigation anti-interference method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107045131B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108107451A (en) * | 2017-11-30 | 2018-06-01 | 西北工业大学 | The vector array interference positioning method decomposed based on tensor CP |
| CN109802717A (en) * | 2019-03-12 | 2019-05-24 | 西北工业大学 | A kind of polarization based on tensor-airspace Wave beam forming anti-interference method |
| CN110266363A (en) * | 2019-06-26 | 2019-09-20 | 电子科技大学 | A kind of distributed diffusion self-adapting anti-jamming method based on tensor |
| US20230048116A1 (en) * | 2020-12-16 | 2023-02-16 | Zhejiang University | Composite tensor beamforming method for electromagnetic vector coprime planar array |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6166690A (en) * | 1999-07-02 | 2000-12-26 | Sensor Systems, Inc. | Adaptive nulling methods for GPS reception in multiple-interference environments |
| CN101476893A (en) * | 2009-01-20 | 2009-07-08 | 郑州紫辰科技开发有限公司 | Non-blind zone GPS receiving device and receiving method |
| CN103235327A (en) * | 2013-04-07 | 2013-08-07 | 清华大学 | GNSS/MINS (global navigation satellite system/micro-electro-mechanical systems inertial navigation system) super-deep combination navigation method, system and device |
| CN104049262A (en) * | 2014-06-26 | 2014-09-17 | 哈尔滨工程大学 | Beam forming anti-interference method based on vector tracking loop |
-
2017
- 2017-05-24 CN CN201710371234.0A patent/CN107045131B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6166690A (en) * | 1999-07-02 | 2000-12-26 | Sensor Systems, Inc. | Adaptive nulling methods for GPS reception in multiple-interference environments |
| CN101476893A (en) * | 2009-01-20 | 2009-07-08 | 郑州紫辰科技开发有限公司 | Non-blind zone GPS receiving device and receiving method |
| CN103235327A (en) * | 2013-04-07 | 2013-08-07 | 清华大学 | GNSS/MINS (global navigation satellite system/micro-electro-mechanical systems inertial navigation system) super-deep combination navigation method, system and device |
| CN104049262A (en) * | 2014-06-26 | 2014-09-17 | 哈尔滨工程大学 | Beam forming anti-interference method based on vector tracking loop |
Non-Patent Citations (1)
| Title |
|---|
| 吕潇君等: "幅相误差对数字多波束导航抗干扰影响分析", 《导航定位学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108107451A (en) * | 2017-11-30 | 2018-06-01 | 西北工业大学 | The vector array interference positioning method decomposed based on tensor CP |
| CN108107451B (en) * | 2017-11-30 | 2019-05-10 | 西北工业大学 | The vector array interference positioning method decomposed based on tensor CP |
| CN109802717A (en) * | 2019-03-12 | 2019-05-24 | 西北工业大学 | A kind of polarization based on tensor-airspace Wave beam forming anti-interference method |
| CN110266363A (en) * | 2019-06-26 | 2019-09-20 | 电子科技大学 | A kind of distributed diffusion self-adapting anti-jamming method based on tensor |
| CN110266363B (en) * | 2019-06-26 | 2021-09-24 | 电子科技大学 | Tensor-based distributed diffusion self-adaptive anti-interference method |
| US20230048116A1 (en) * | 2020-12-16 | 2023-02-16 | Zhejiang University | Composite tensor beamforming method for electromagnetic vector coprime planar array |
| US11841448B2 (en) * | 2020-12-16 | 2023-12-12 | Zhejiang University | Composite tensor beamforming method for electromagnetic vector coprime planar array |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107045131B (en) | 2019-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107045131B (en) | Tensor resolution satellite navigation anti-interference method | |
| CN100392426C (en) | Single channel phase control array receiving signal reconstruction and space signal treatment method | |
| Sun et al. | Partial spectral search‐based DOA estimation method for co‐prime linear arrays | |
| CN108896954B (en) | Estimation method of angle of arrival based on joint real-value subspace in co-prime matrix | |
| CN102830386B (en) | Estimation method of arbitrary array weak signal source angle under strong interference | |
| CN109061685B (en) | Satellite navigation receiver anti-interference method and system based on switching antenna array | |
| CN101644760B (en) | Rapid and robust method for detecting information source number suitable for high-resolution array | |
| CN109061554A (en) | A kind of target angle of arrival estimation method updated based on spatial spreading grid dynamic | |
| Forster et al. | Threshold performance analysis of maximum likelihood DOA estimation | |
| CN109143275A (en) | A kind of miniaturization array antenna anti-interference realization method based on population | |
| CN112130111A (en) | Single-snapshot two-dimensional DOA estimation method for large-scale uniform cross array | |
| CN104931923A (en) | Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA) | |
| Lu et al. | Optimal reference element for interference suppression in GNSS antenna arrays under channel mismatch | |
| CN101483280B (en) | Weight solving method for stable wave beam synthesizer | |
| Yang et al. | High-order cumulants based sparse array design via fractal geometries—Part I: Structures and DOFs | |
| CN108107451B (en) | The vector array interference positioning method decomposed based on tensor CP | |
| Xu et al. | A two-step beam-forming method based on carrier phases for GNSS adaptive array anti-jamming | |
| Tian et al. | 2D-DOA estimation in arc-array with a DNN based covariance matrix completion strategy | |
| CN113466782A (en) | Deep Learning (DL) -based cross-coupling correction D O A estimation method | |
| CN109116377A (en) | A kind of satellite navigation anti-interference method and device calculated based on time domain submatrix | |
| Chen et al. | Super-resolution DOA estimation using dynamic metasurface antenna | |
| Yu et al. | A novel 2D off-grid DOA estimation method based on compressive sensing and least square optimization | |
| CN114706034A (en) | Target direction finding method for unmanned aerial vehicle carrying intelligent reflecting surface | |
| CN100544233C (en) | Space radiation source DOA estimation method based on nonopiate decomposition | |
| CN104898094A (en) | Space-frequency two-dimensional beam forming method of compressed sampling array |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |