CN106680837A - Interference suppression algorithm for satellite navigation - Google Patents

Interference suppression algorithm for satellite navigation Download PDF

Info

Publication number
CN106680837A
CN106680837A CN201611149124.1A CN201611149124A CN106680837A CN 106680837 A CN106680837 A CN 106680837A CN 201611149124 A CN201611149124 A CN 201611149124A CN 106680837 A CN106680837 A CN 106680837A
Authority
CN
China
Prior art keywords
weight vector
signal
vector
covariance matrix
interference suppression
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
Application number
CN201611149124.1A
Other languages
Chinese (zh)
Other versions
CN106680837B (en
Inventor
王晓宇
谢斌斌
金燕
张骅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC 20 Research Institute
Original Assignee
CETC 20 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CETC 20 Research Institute filed Critical CETC 20 Research Institute
Priority to CN201611149124.1A priority Critical patent/CN106680837B/en
Publication of CN106680837A publication Critical patent/CN106680837A/en
Application granted granted Critical
Publication of CN106680837B publication Critical patent/CN106680837B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Abstract

The invention provides an interference suppression algorithm for satellite navigation. The interference suppression algorithm is characterized in that the estimation value of a received-data covariance matrix is calculated, eigen-decomposition operation is performed on the estimation value of the received-data covariance matrix to obtain eigenvalues and corresponding eigenvectors, reciprocal operation is performed on the minimal eigenvalue, the obtained reciprocal is multiplied with the corresponding eigenvector to obtain an unnormalized weight vector, normalization process is performed, and the weight vector is used to perform weighted stacking on a received medium-frequency digital complex signal to obtain an output signal after interference suppression. By the interference suppression algorithm, strong interference signals and weak interference signals can be effectively inhibited.

Description

A kind of satellite navigation interference suppression algorithm
Technical field
It is a kind of pressing type satellite navigation interference signal to be realized to suppress the invention belongs to the anti-interference field of satellite navigation Algorithm.
Background technology
In the signal processing of anti-interference antenna of satellite navigation system, interference suppression algorithm is at whole digital signal The part of core in reason.Interference suppression algorithm can be divided into adaptive nulling class algorithm and optimum digital multiple beam class algorithm two is big Class.
The Typical Representative of adaptive nulling class algorithm is power inversion (Power Inversion, PI) algorithm, and the algorithm is not The prior informations such as satellite position information, null depth is needed to automatically adjust with the intensity of jamming power, for high reject signal has Preferable inhibition, thus be widely used in the engineer applied in the anti-interference field of satellite navigation.But, it lacks Point is poor for the inhibition of weak jamming signal so that it is in weak jamming environment or while there is strong jamming and weak jamming Application in environment is restricted.
Optimum digital multiple beam algorithm adopts digital multiplex technology, while multiple optimum numerals are formed in visual-field space connecing Receive wave beam.Typical optimum digital beam is the undistorted response response of minimum variance (Minimum Variance Distorted Response, MVDR) wave beam.The main lobe of each optimum digital beam points to an aeronautical satellite, while adaptively in interference Direction forms null.Optimum digital multiple beam algorithm can be with the noise of lift system output signal while to AF panel Than.But the algorithm needs the prior informations such as satellite position, array attitude to aid in, algorithm structure is complicated, amount of calculation is larger, it is difficult to It is used widely in engineer applied.
The content of the invention
In order to overcome the deficiencies in the prior art, the present invention provides a kind of satellite navigation interference suppression algorithm, using normalization The noise subspace characteristic vector of receiving data covariance matrix afterwards is used as in the docking collection of letters number of adaptive nulling weight vector Interference components are suppressed, and the weak powered interferer component and high power interference that can be docked simultaneously in the collection of letters number is carried out effectively Ground suppresses, it is easy to Project Realization.
The technical solution adopted for the present invention to solve the technical problems is comprised the following steps：
The first step, calculates receiving data covariance matrix RxEstimated valueWherein, x (n) is by M The intermediate frequency digital complex signal vector that the array antenna of individual array element is arrived in n receptions, dimension is M × 1, M >=2；Dimension Size is M × M；Subscript H is conjugate transpose operator；N is to calculate the signal that the estimated value of receiving data covariance matrix needs Fast umber of beats, N >=4M；
Second step, the estimated value to receiving data covariance matrix carries out feature decomposition computing using Jacobi algorithms, obtains Obtain eigenvalue λmWith corresponding characteristic vector qm, m=1,2 ..., M；
3rd step：By minimal eigenvalue λminSeek derivative action and corresponding characteristic vector is multiplied and obtains non-normalizing The weight vector of change Dimension size be M × 1；
4th step, to weight vectorIt is normalized, obtains normalized weight vectorWherein,ForFirst element of weight vector；
5th step, using weight vector woptIntermediate frequency digital complex signal x (n) to receiving in the first step is weighted folded Plus process, the output signal after AF panel
The invention has the beneficial effects as follows：Weight vector w obtained using the present inventionoptSave corresponding to interference characteristic vector Power component so that noise signal component is only included in the signal power exported after weighting, can simultaneously to high reject signal and weak Interference signal is effectively suppressed.
Description of the drawings
Fig. 1 is the algorithm flow chart of the present invention.
Specific embodiment
With reference to the accompanying drawings and examples the present invention is further described, and the present invention includes but are not limited to following enforcements Example.
Inventive algorithm realizes that step is as follows：
The first step：Calculate receiving data covariance matrix RxEstimated valueCalculation expression is as follows：
Wherein, x (n) be from the individual array elements of M (M >=2) array antenna n receptions to intermediate frequency digital complex signal to Amount, dimension is M × 1；Dimension size be M × M；Subscript H is conjugate transpose operator；N is calculating receiving data covariance The fast umber of beats of signal that the estimated value of matrix needs, need to meet N >=4M.
Second step：Estimated value to receiving data covariance matrix carries out feature decomposition computing, and feature decomposition is adopted The realization of Jacobi algorithms, obtains eigenvalue λm(m=1,2 ..., M) and corresponding characteristic vector qm(m=1,2 ..., M).
3rd step：By minimal eigenvalue λminSeek derivative action and corresponding characteristic vector is multiplied and obtains non-normalizing The weight vector of change
Dimension size be M × 1.
4th step：To the not normalized weight vector obtained in the 3rd stepIt is normalized, obtains normalized Weight vector wopt, calculating process is shown below：
Wherein,ForFirst element of weight vector.
5th step：Using calculated weight vector w in the 4th stepoptTo the intermediate frequency digital plural number received in the first step Signal x (n) is weighted overlap-add procedure, and calculating process is shown below：
Wherein, subscript H is conjugate transpose operator, and y is the output signal after AF panel.
For Power-inversion algorithm, its weight vector is calculated as follows：
Wherein, wPIFor the calculated weight vector of Power-inversion algorithm, dimension size is M × 1；To receive data association side Difference Matrix Estimation valueInverse matrix, dimension size be M × M；a0=[1,0 ..., 0]TFor constrained vector, dimension size be M × 1；K is interference source number, and K ＜ M；δ2For system noise power；wjamAnd wnoiseIt is divided into interference signal and noise signal two parts The corresponding weight vector of difference, dimension size is M × 1.
When there is high reject signal in the external world (power of interference signal be much larger than noise power), interference signal to Eigenvalue will be much larger than the corresponding eigenvalue of noise signal, i.e.,
From formula (6), with the increase of jamming power, the power component in best initial weights corresponding to interference characteristic vector will It is less and less.Therefore, power inversion method can be upwardly formed deeper null in strong jamming side, and capacity of resisting disturbance is stronger.Conversely, If during the power of interference signal less (being slightly larger than noise signal), corresponding to the power component of interference characteristic vector in best initial weights (in formula (5) second equal sign right side Section 1) will be slightly less than weight vector corresponding to noise characteristic component (second in formula (5) Equal sign right side Section 2), the null of formation shoals, and capacity of resisting disturbance is poor.And adopt weight vector w for obtaining of the inventionoptSave Corresponding to the power component of interference characteristic vector so that only include noise signal component in the signal power exported after weighting, can be with High reject signal and weak jamming signal are effectively suppressed simultaneously.
The present invention is applied to GPS, BDS and GLONASS satellite navigation system AF panel is processed.With 4 unit Bs D2-B3 frequency Point anti-interference antenna resists the specific embodiment that the present invention is illustrated as a example by two broadband interferences.
Step 1：Calculate the estimated value of 4 element array receiving data covariance matrixesThe fast umber of beats N=256 of data sampling, meter ObtainIt is Hermit matrixes that dimension is 4 × 4, calculating process is as follows：
Step 2：Feature decomposition computing is carried out to the estimated value for receiving data covariance matrix using Jacobi algorithms, is obtained Eigenvalue λm(m=1,2 ..., 4) with corresponding characteristic vector qm(m=1,2 ..., 4).
Step 3：By the calculated minimal eigenvalue λ of step 2minAsk derivative action and corresponding characteristic vector Multiplication obtains not normalized weight vectorCalculating process is as follows：
Step 4：To the not normalized weight vector obtained in step 3It is normalized, obtains normalized power Vectorial wopt, shown in calculating process such as formula (3).
Step 5：Using calculated weight vector w in step 4optTo the intermediate frequency digital plural number letter received in the first step Number x (n) is weighted overlap-add procedure and exports, shown in calculating process such as formula (4).
So far, the AF panel for completing to receive 4 unit B D2-B3 frequencies signal is processed.

Claims (1)

1. a kind of satellite navigation interference suppression algorithm, it is characterised in that comprise the steps：
The first step, calculates receiving data covariance matrix RxEstimated valueWherein, x (n) is by M battle array The intermediate frequency digital complex signal vector that the array antenna of unit is arrived in n receptions, dimension is M × 1, M >=2；Dimension size For M × M；Subscript H is conjugate transpose operator；N is to calculate the signal snap that the estimated value of receiving data covariance matrix needs Number, N >=4M；
Second step, the estimated value to receiving data covariance matrix carries out feature decomposition computing using Jacobi algorithms, obtains special Value indicative λmWith corresponding characteristic vector qm, m=1,2 ..., M；
3rd step：By minimal eigenvalue λminSeek derivative action and corresponding characteristic vector is multiplied and obtains not normalized Weight vector Dimension size be M × 1；
4th step, to weight vectorIt is normalized, obtains normalized weight vectorWherein, ForFirst element of weight vector；
5th step, using weight vector woptIntermediate frequency digital complex signal x (n) to receiving in the first step is weighted at superposition Reason, the output signal after AF panel
CN201611149124.1A 2016-12-14 2016-12-14 A kind of satellite navigation interference suppression algorithm Active CN106680837B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611149124.1A CN106680837B (en) 2016-12-14 2016-12-14 A kind of satellite navigation interference suppression algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611149124.1A CN106680837B (en) 2016-12-14 2016-12-14 A kind of satellite navigation interference suppression algorithm

Publications (2)

Publication Number Publication Date
CN106680837A true CN106680837A (en) 2017-05-17
CN106680837B CN106680837B (en) 2019-04-19

Family

ID=58868312

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611149124.1A Active CN106680837B (en) 2016-12-14 2016-12-14 A kind of satellite navigation interference suppression algorithm

Country Status (1)

Country Link
CN (1) CN106680837B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108089162A (en) * 2017-12-29 2018-05-29 中国电子科技集团公司第二十研究所 A kind of detection of pulse interference signal and suppressing method
CN108241160A (en) * 2017-12-29 2018-07-03 中国电子科技集团公司第二十研究所 A kind of satellite navigation interference source space spectral peak method of estimation

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103630910A (en) * 2013-12-13 2014-03-12 武汉大学 Anti-interference method of GNSS (global navigation satellite system) receiver equipment
CN104536018A (en) * 2015-01-06 2015-04-22 中国人民解放军国防科学技术大学 GNSS multi-satellite unified capture method using array antenna anti-interference technology
CN105137454A (en) * 2015-07-22 2015-12-09 北京航空航天大学 Anti-interference algorithm FPGA realization method based on covariance matrix characteristic decomposition and realization device thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103630910A (en) * 2013-12-13 2014-03-12 武汉大学 Anti-interference method of GNSS (global navigation satellite system) receiver equipment
CN104536018A (en) * 2015-01-06 2015-04-22 中国人民解放军国防科学技术大学 GNSS multi-satellite unified capture method using array antenna anti-interference technology
CN105137454A (en) * 2015-07-22 2015-12-09 北京航空航天大学 Anti-interference algorithm FPGA realization method based on covariance matrix characteristic decomposition and realization device thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
江城 等: ""一种改进的功率倒置卫星导航抗干扰方法"", 《现代导航》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108089162A (en) * 2017-12-29 2018-05-29 中国电子科技集团公司第二十研究所 A kind of detection of pulse interference signal and suppressing method
CN108241160A (en) * 2017-12-29 2018-07-03 中国电子科技集团公司第二十研究所 A kind of satellite navigation interference source space spectral peak method of estimation

Also Published As

Publication number Publication date
CN106680837B (en) 2019-04-19

Similar Documents

Publication Publication Date Title
CN102830387B (en) Data preprocessing based covariance matrix orthogonalization wave-beam forming method
CA2791593C (en) Method and system for locating interference by frequency sub-band
CN109254261A (en) Coherent signal null based on uniform circular array EPUMA deepens method
CN103837861B (en) The Subarray linear restriction Adaptive beamformer method of feature based subspace
CN105629266B (en) Formula is cheated in satellite navigation and pressing type disturbs the joint suppressing method of blind adaptive
CN105302936A (en) Self-adaptive beam-forming method based on related calculation and clutter covariance matrix reconstruction
De Lamare Adaptive reduced-rank LCMV beamforming algorithms based on joint iterative optimisation of filters
CN106353738B (en) A kind of robust adaptive beamforming method under new DOA mismatch condition
Vesa Direction of arrival estimation using music and root-music algorithm
CN105204008A (en) Adaptive antenna wave beam forming nulling widening method based on covariance matrix extension
CN105022268A (en) Linear constraint virtual antenna beam forming method
Li et al. High performance robust adaptive beamforming in the presence of array imperfections
CN106680837A (en) Interference suppression algorithm for satellite navigation
CN105699988B (en) For the denoising MVDR cheating interference suppressing methods of navigation neceiver
CN113504549B (en) Navigation space-time anti-interference method based on generalized sidelobe canceller
CN110727915A (en) Robust self-adaptive beam forming method based on data correlation constraint
CN108776347A (en) A kind of dual-polarized antenna array broadens the high-dynamic GNSS disturbance restraining method of technology based on null
US8144057B1 (en) Methods and apparatus for adaptively determining angles of arrival of signals
Qian et al. Robust beamforming based on steering vector and covariance matrix estimation
Sarkar et al. Performance analysis of uniform concentric circular antenna array beamformer using different doa estimation technique
CN104459627A (en) Reduced rank beam forming method based on united alternative optimization
CN108241160A (en) A kind of satellite navigation interference source space spectral peak method of estimation
CN106452548A (en) Adaptive robust beamforming method
CN108833038B (en) Signal power estimation method based on oblique projection operator
Cai et al. Low-complexity reduced-dimension space–time adaptive processing for navigation receivers

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