CN104330808A - Multiclass satellite navigation interference suppression method based on solution expanding technology - Google Patents

Multiclass satellite navigation interference suppression method based on solution expanding technology Download PDF

Info

Publication number
CN104330808A
CN104330808A CN201410606254.8A CN201410606254A CN104330808A CN 104330808 A CN104330808 A CN 104330808A CN 201410606254 A CN201410606254 A CN 201410606254A CN 104330808 A CN104330808 A CN 104330808A
Authority
CN
China
Prior art keywords
satellite
heavily
solution
signal
deceiving interference
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
CN201410606254.8A
Other languages
Chinese (zh)
Other versions
CN104330808B (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.)
Civil Aviation University of China
Original Assignee
Civil Aviation University of China
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 Civil Aviation University of China filed Critical Civil Aviation University of China
Priority to CN201410606254.8A priority Critical patent/CN104330808B/en
Publication of CN104330808A publication Critical patent/CN104330808A/en
Application granted granted Critical
Publication of CN104330808B publication Critical patent/CN104330808B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/21Interference 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)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Radio Relay Systems (AREA)

Abstract

The invention discloses a multiclass satellite navigation interference suppression method based on solution expanding technology. The method comprises the steps as follows: projecting the signal received by the array antenna to the interference orthogonal complement space for restraining the blanket jamming, calculating the solution expanding weight vector matrix for the projected data, judging whether the matrix is full rank for identifying the deception jamming. If having deception jamming, calculating the solution expanding weight vector dependency, identifying the number of deception interference source, the deception interference source comprises fake satellite signal and true satellite signal. Using the solution expanding weight vector of the deception jamming for forming the rectangular projection matrix of the deception jamming, having the projection restraining deception jamming on signal, using the weight vector of true satellite for forming multi-wave beam and obtaining the signal processing gain. The method can generate high gain multi-wave beam with suppress interference, transmitted deception jamming and generated type deception jamming resisting capability, the multi-wave beam forming capability is not limited by the array element, and there is no need to know the satellite signal and interference information.

Description

Based on the multiclass satellite navigation disturbance restraining method of despreading and respreading technology
Technical field
The present invention relates to the anti-interference field of satellite navigation, particularly relate to a kind of multiclass satellite navigation disturbance restraining method based on despreading and respreading technology.
Background technology
GPS (Global Position System) (GNSS) utilizes artificial satellite to launch navigation signal, thus the information such as high-precision position and time can be provided at any time and any place for the receiver in coverage, therefore satellite navigation system has critical importance.Because GNSS satellite signal is extremely faint, how to guarantee that the safety and reliability of GNSS application facility is particularly important, anti-interference problem has become the matter of utmost importance that GPS (Global Position System) faces.At present, the research of satellite navigation AF panel aspect mainly concentrates on the suppression of suppression jamming.Wherein power inversion method utilizes suppression jamming power much larger than satellite-signal power and this feature of noise power, thus suppression jamming upwards formed zero fall into, but because the power of Deceiving interference is lower than noise level, therefore power inversion method cannot be used for the suppression of Deceiving interference.
Existing Deceiving interference suppressing method be from interference to, catch link, aspect such as tracking link, positioning result and navigation message etc. suppresses Deceiving interference.S.Daneshmand proposes a kind of Deceiving interference suppression technology based on double antenna, the Received signal strength of two antennas is carried out computing cross-correlation and obtain coming to angle of Deceiving interference, and then Deceiving interference is suppressed, but the required known array stream shape of this method, pair array stream shape error sensitive, be only applicable to the suppression of single relay type Deceiving interference, and require that the power of Deceiving interference is higher.The people such as G.B.Moon propose to utilize track loop to follow the tracks of all possible deception formula satellite-signal and real satellite signal, then distinguish deception satellite-signal and real satellite signal according to tracking parameter.This technology needs more tracking channel, and the complexity of receiver and cost are increased greatly.History Vincent proposes a kind of DOA estimated value of the DR position of receiver and Received signal strength that utilizes and carrys out method to suppress Deceiving interference, and the method needs to calculate receiver location, and carries out demodulation to navigation message, and computation complexity is high.C.E.McDowell is formed multi-antenna beam and the zeroing theory of technology and performance are analyzed, utilize GPS correlator parameter and track loop parameter to form multi-beam identify and suppress Deceiving interference, but this technology does not possess ability real satellite signal being formed to multi-beam gain while suppressing Deceiving interference.It should be noted that these methods equal cisco unity malfunction when suppression jamming and Deceiving interference coexist.In order to ensure the safety applications of GNSS, Deceiving interference can be identified in the urgent need to research, then suppress suppression jamming, relay type Deceiving interference and production Deceiving interference simultaneously, and pair array stream shape error and jamming power insensitive multiclass satellite navigation disturbance restraining method, this not only has important meaning in theoretical research, and also has great importance in the application in the fields such as Military Application to propelling GNSS.
Summary of the invention
In order to solve the problem, the object of the present invention is to provide one not need known array stream shape and receiver feedback information, satellite navigation system suppression jamming can be suppressed, identify with suppression relay type Deceiving interference and production Deceiving interference and the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology of gain can be provided satellite-signal.
In order to achieve the above object, the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology provided by the invention comprises the following step carried out in order:
(1) by array antenna received to radiofrequency signal change into digital medium-frequency signal:
(2) according to the repeatedly fast beat of data that array antenna received arrives, calculate the interference orthogonal complement space, signal suppressed suppression jamming to the projection of the interference orthogonal complement space:
(3) get a road to the data after projection to calculate solution and heavily to extend power vector matrix, and judge that solution is heavily extended power vector matrix whether full rank, to identify whether Deceiving interference exists;
(4) if there is Deceiving interference, the correlativity that solution is heavily extended power between vector is calculated, the pseudo satellite, pseudolite asterisk comprised according to the number in correlated results determination Deceiving interference source, each Deceiving interference source and real satellite asterisk;
(5) utilize the solution of Deceiving interference heavily extend power vector structure Deceiving interference orthogonal intersection cast shadow matrix, project to suppress Deceiving interference to above-mentioned Received signal strength;
(6) utilize the solution of real satellite vector of heavily extending power to carry out Multibeam synthesis, obtain signal transacting gain.
In step (3), described a road is got to the data after above-mentioned projection calculate solution and heavily to extend power vector matrix, and judge that solution is heavily extended power vector matrix whether full rank, to identify that method that whether Deceiving interference exists heavily to be extended power vector to the solution that each satellite captured all calculates its correspondence according to despreading and respreading technology, this solution heavily extend power vector include signal come to information, then moment of a vector rank of matrix of heavily extending power is solved, the vector matrix full rank if solution is heavily extended power, then show that satellite-signal is separate each other, from different directions, Deceiving interference is there is not in Received signal strength, otherwise, the vector matrix not full rank if solution is heavily extended power, then show have from unidirectional satellite-signal in Received signal strength, there is Deceiving interference.
In step (4), the correlativity that described calculating solution is heavily extended power between vector, classify according to heavily the extend power vector that the solution of satellite-signal heavily extended power that compares of normalized correlation coefficient and thresholding between vector of solution according to the method for correlated results determination Deceiving interference source number, pseudo satellite, pseudolite asterisk that each Deceiving interference source comprises and real satellite asterisk, if exist from unidirectional different satellite-signal, the correlativity that solution then corresponding to these satellite-signals is heavily extended power between vector is very high, think Deceiving interference, otherwise, think real satellite signal.
In step (5), the described solution of Deceiving interference that utilizes heavily is extended power the orthogonal intersection cast shadow matrix of vector structure Deceiving interference, project to suppress the method for Deceiving interference to be corresponding with Deceiving interference subspace according to heavily the extend power subspace of vector of the solution of Deceiving interference to above-mentioned Received signal strength, recycling orthogonal complement space computing formula tries to achieve the orthogonal intersection cast shadow matrix of Deceiving interference.
In step (6), the described solution of the real satellite vector of heavily extending power that utilizes carries out Multibeam synthesis, the method obtaining signal transacting gain heavily extends power vectors converge in the array manifold of real satellite signal according to the solution of real satellite, can be formed further array beams point to corresponding satellite come to.
Multiclass satellite navigation disturbance restraining method based on despreading and respreading technology provided by the invention can produce multiple wave beams with anti-suppression jamming, relay type Deceiving interference and production Deceiving interference ability simultaneously of high-gain, Multibeam synthesis ability does not limit by array number, without the need to knowing that satellite-signal and interference come to information, pair array stream shape error and Deceiving interference watt level insensitive, do not need the feedback information of navigation neceiver.
Accompanying drawing explanation
Fig. 1 adopts the solution that obtains of the inventive method heavily to extend power the beam pattern of vector matrix.
Fig. 2 is the total beam pattern adopting the inventive method to obtain all satellite-signals.
Fig. 3 be do not take after carrying out suppression jamming rectangular projection the inventive method suppress Deceiving interference catch result.
Fig. 4 catches result after taking the inventive method to suppress suppression jamming and Deceiving interference.
Fig. 5 does not take the inventive method to suppress the two-dimensional search result of catching PRN1 during Deceiving interference after carrying out suppression jamming rectangular projection.
Fig. 6 takes the inventive method to suppress the two-dimensional search result of catching PRN1 during Deceiving interference after carrying out suppression jamming rectangular projection.
Fig. 7 is the process flow diagram of the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology provided by the invention.
Embodiment
With specific embodiment, the multiclass satellite navigation disturbance restraining method that the present invention is based on despreading and respreading technology is described in detail with reference to the accompanying drawings.For convenience's sake, be illustrated for the array antenna of GPS (Global Positioning System) system and N number of array element below, in fact, the method is applicable to General Cell antenna, and does not need to know array antenna stream shape.
As shown in Figure 7, the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology provided by the invention comprises the following step carried out in order:
(1) by array antenna received to radiofrequency signal change into digital medium-frequency signal:
Deceiving interference is made up of multi-satellite signal (being called Pseudolite signal) and launches in the same direction.Wherein, relay type Deceiving interference is by receiving current satellite-signal, postpones and amplification satellite-signal according to deception object, then forwards and realizes interference; Production Deceiving interference is the information utilizing commercial satellite navigation signal form known, directly generates Pseudolite signal and launches.
When suppression jamming, relay type Deceiving interference and production Deceiving interference incide on array antenna simultaneously, the intermediate-freuqncy signal that GPS receives is:
x ( nT s ) = Σ m = 1 M a m s m a ( nT s ) + Σ p = 1 P b p s p rs ( nT s ) + Σ q = 1 Q b q s q gs ( nT s ) + Σ k = 1 K c k j k ( nT s ) + n ( nT s ) - - - ( 1 )
Wherein, M, P, Q and K represent real satellite signal number, relay type Deceiving interference source number, production Deceiving interference source number and suppression jamming number respectively. represent the navigation signal of m real satellite signal, represent p relay type Deceiving interference, represent q production Deceiving interference, subscript a, rs and gs indicate real satellite signal, relay type Deceiving interference and production Deceiving interference respectively.Wherein,
s m a ( nT s ) = p m a D m a ( nT s - τ m a ) c m a ( nT s - τ m a ) e j φ m a + j 2 π f m ′ a n T s - - - ( 2 )
s p rs ( nT s ) = Σ l = 1 L p p l rs D l rs ( nT s - τ l rs ) c l rs ( nT s - τ l rs ) e j φ l rs + j 2 π f l ′ rs n T s - - - ( 3 )
s q gs ( nT s ) = Σ l = 1 L q p l gs D l gs ( nT s - τ l gs ) c l gs ( nT s - τ l gs ) e j φ l gs + j 2 π f l ′ gs n T s - - - ( 4 )
φ, f ', p and τ represent phase place, Doppler frequency and code time delay respectively.D (nT s) and c (nT s) represent navigation message and C/A code respectively.L in formula (3) prepresent the number of the Pseudolite signal that p relay type Deceiving interference source comprises, these pseudo satellite, pseudolites are usually identical with real satellite asterisk; L in formula (4) qrepresent the number of the Pseudolite signal that q production Deceiving interference source comprises, these pseudo satellite, pseudolites are not necessarily identical with real satellite asterisk.J k(nT s) represent a kth suppression jamming, n (nT s) be additive white Gaussian noise, a m, b p, b q, c kcorresponding array manifold, T sfor sampling interval.
(2) according to the repeatedly fast beat of data that array antenna received arrives, calculate the interference orthogonal complement space, signal suppressed suppression jamming to the projection of the interference orthogonal complement space:
Due to the low 20dB of power ratio noise level of satellite-signal, Deceiving interference power is suitable with satellite-signal, and the power of suppression jamming is usually all very strong, generally will far above noise power, so the covariance matrix of Received signal strength determines primarily of the covariance matrix of suppression jamming and noise.That is:
R xx≈R j+R v(5)
Wherein, R xxrepresent the covariance matrix of Received signal strength, R vrepresent noise covariance matrix, R jrepresent suppression jamming covariance matrix.Now, the covariance matrix of Received signal strength is inverted and to be directly proportional to the orthogonal complement space of suppression jamming subspace, and Received signal strength is suppressed suppression jamming to the projection of the interference orthogonal complement space:
y ( nT s ) = R xx - 1 x ( nT s ) - - - ( 6 )
(3) get a road to the data after projection to calculate solution and heavily to extend power vector matrix W, and judge that solution is heavily extended power vector matrix whether full rank, to identify whether Deceiving interference exists:
W=[w 1w 2… w l] (7)
Wherein, w l=[w l1w lN] trepresent that the solution of l satellite is heavily extended power vector.
Despreading and respreading technology refers to and detects that the n-th number of bits transmitted is according to b by spreading codes information n, and utilize b nagain spread spectrum goes out to transmit at time cycle [(n-1) T b, nT b] in signal model, then the signal after heavily expanding just can upgrade the weight vector of Adaptive beamformer as reference signal.Utilize the feature in gps signal spreading code codeword structure known and navigation data cycle much larger than spreading code cycle, heavily the extend power cost function of vector of the solution solving l satellite is:
w l = min w Σ k = 1 K ‾ | w l H y ( k ) - r l ( k ) | 2 - - - ( 8 )
Wherein, r lk () is the reference signal of l satellite, being the length of data block, is the sample number in the spreading code cycle.In formula (8), the more new formula of weight vector is:
w l ( h + 1 ) = [ Y ( i ) Y H ( i ) ] - 1 Yr l * ( h ) - - - ( 9 )
Y ( i ) = [ y ( 1 + i K ‾ ) , . . . , y ( K ‾ + i K ‾ ) ] T - - - ( 10 )
r l ( h ) = [ c l ( 1 + h K ‾ - τ ^ l ) e j 2 π f ^ l ( 1 + h K ‾ - τ ^ l ) , . . . , c l ( K ‾ + h K ‾ - τ ^ l ) e j 2 π f ^ l ( K ‾ + h K ‾ - τ ^ l ) ] T - - - ( 11 )
Wherein, h represents the iterations that weight vector calculates, with the code delay of l satellite and the valuation of carrier frequency of catching.When Deceiving interference exists, the correlated results of catching of same asterisk is not often a peak value, as shown in Figure 4.Adopt the condition of ratio as prize judgment of peak value and average noise, the whole real satellite and the pseudo satellite, pseudolite that capture same asterisk can be ensured, and obtain the valuation of code delay and carrier frequency, to reconstruct the reference signal of whole real satellite and pseudo satellite, pseudolite, thus obtain corresponding solution and heavily to extend power vector.
Algorithm concrete steps are as follows:
Step 1: pair array exports and carries out acquisition and tracking, obtains the c of the satellite of having caught l(k), etc. parameter, this is the process of despreading;
Step 2: the spread-spectrum signal r generating every effective satellite according to step 1 parameters obtained by formula (11) respectively l(h), this is the process heavily expanded;
Step 3: respectively with r lh (), as with reference to signal, utilizes formula (9) to upgrade the weight vector of l satellite;
Step 4: computing array output vector;
Step 5: repeat step 1-4, until algorithm convergence.
Judge that solution is heavily extended power vector matrix whether full rank:
rank W=m<l (12)
If formula (12) is false, then w is described 1..., w luncorrelated mutually between any two, do not exist from unidirectional signal, namely currently to catch, there is not Deceiving interference in tracking satellite.Otherwise, if formula (12) is set up, then illustrate and currently to catch, there is Deceiving interference in tracking satellite.
(4) if there is Deceiving interference, the correlativity that solution is heavily extended power between vector is calculated, the pseudo satellite, pseudolite asterisk comprised according to the number in correlated results determination Deceiving interference source, each Deceiving interference source and real satellite asterisk;
If satellite-signal S 1and S 2solution vector of heavily extending power be respectively w 1and w 2, ρ represents their normalized correlation coefficient,
&rho; = | ( w 1 | w 1 | ) H w 2 | w 2 | | - - - ( 13 )
When heavily the extend power normalized correlation coefficient ρ of vector of the solution of two signals is greater than decision threshold D time, can think that these two signals are from same direction, namely
&rho; = | ( w 1 | w 1 | ) H w 2 | w 2 | | &GreaterEqual; D &Delta; - - - ( 14 )
Due to real satellite signal from difference come to, and the pseudo satellite, pseudolite that Deceiving interference comprises is launched from same Deceiving interference source all respectively, therefore, can judge that the satellite that each solution that normalized correlation coefficient is greater than decision threshold is heavily extended power corresponding to vector is pseudo satellite, pseudolite, all the other are real satellite signal.
The method of this identification Deceiving interference, does not utilize the power information of Deceiving interference and real satellite signal, insensitive to Deceiving interference power.
(5) utilize the solution of Deceiving interference heavily extend power vector structure Deceiving interference orthogonal intersection cast shadow matrix, project to suppress Deceiving interference to above-mentioned Received signal strength;
According to the judged result of step (4), with the solution of Deceiving interference heavily extend power vector structure Deceiving interference subspace, be designated as U=span{w 1..., w p+Q.Then its orthogonal complement space is:
U =I-U(U HU) -1U H(15)
Obtaining Deceiving interference orthogonal intersection cast shadow matrix is:
P =U (U ) H(16)
The data without Deceiving interference after being projected by array received signal are:
z(nT s)=P y(nT s) (17)
(6) utilize the solution of real satellite vector of heavily extending power to carry out Multibeam synthesis, obtain signal transacting gain:
Utilize the solution of real satellite vector of heavily extending power to carry out Multibeam synthesis to the output vector of formula (17), obtaining final array antenna output signal is:
output l ( nT s ) = w l H z ( nT s ) - - - ( 18 )
Wherein, subscript l represents it is output signal to l real satellite.Due to acquisition procedure ask for solution heavily extend power vector process in complete, and obtained the c of satellite l(k), etc. parameter, the output of formula (18) is utilized directly to follow the tracks of l real satellite.
Fig. 1 adopts the solution that obtains of the inventive method heavily to extend power the beam pattern of vector matrix.In the present embodiment, array antenna is ten yuan of even linear arrays, four gps satellite signal PRN1, and PRN2, PRN3, PRN25 are respectively from-5 °, and 5 ° ,-25 °, 15 ° of directions are incided on array antenna, and signal to noise ratio (S/N ratio) is-20dB.Relay type Deceiving interference source 1 is made up of the Pseudolite signal of PRN1, PRN2, PRN3, PRN25, incides array antenna from 25 °.Production Deceiving interference source 2 is made up of the Pseudolite signal of PRN6, PRN14, PRN22, incides array antenna from-40 °.The signal to noise ratio (S/N ratio) of each Pseudolite signal in Deceiving interference source is-18dB.Suppression jamming incides array antenna from 60 °, dryly makes an uproar than being 40dB.Black dotted lines in Fig. 1 represents real satellite direction of signal, solid black lines represent Pseudolite signal to.As can be seen from Figure 1, solution heavily extend power vector aim at whole satellite-signal come to Deceiving interference come to.
Fig. 2 is the total beam pattern adopting the inventive method to obtain all satellite-signals.In the present embodiment, array antenna is ten yuan of even linear arrays, four gps satellite signal PRN1, and PRN2, PRN3, PRN25 are respectively from-5 °, and 5 ° ,-25 °, 15 ° of directions are incided on array antenna, and signal to noise ratio (S/N ratio) is-20dB.Relay type Deceiving interference source 1 is made up of the Pseudolite signal of PRN1, PRN2, PRN3, PRN25, incides array antenna from 25 °.Production Deceiving interference source 2 is made up of the Pseudolite signal of PRN6, PRN14, PRN22, incides array antenna from-40 °.The signal to noise ratio (S/N ratio) of each Pseudolite signal in Deceiving interference source is-18dB.Suppression jamming incides array antenna from 60 °, dryly makes an uproar than being 40dB.Black dotted lines in Fig. 2 represents real satellite direction of signal, solid black lines represent Pseudolite signal to and suppression jamming come to.As can be seen from Figure 2 the inventive method upwards can form beam gain to every satellite at satellite-signal, and comes to fall into formation zero at suppression jamming and Deceiving interference.
Fig. 3 be do not take after carrying out suppression jamming rectangular projection the inventive method suppress Deceiving interference catch result.In the present embodiment, array antenna is ten yuan of even linear arrays, four gps satellite signal PRN1, and PRN2, PRN3, PRN25 are respectively from-5 °, and 5 ° ,-25 °, 15 ° of directions are incided on array, and signal to noise ratio (S/N ratio) is-20dB.Relay type Deceiving interference source 1 is made up of the Pseudolite signal of PRN1, PRN2, PRN3, PRN25, incides array antenna from 25 °.Production Deceiving interference source 2 is made up of the Pseudolite signal of PRN6, PRN14, PRN22, incides array antenna from-40 °.The signal to noise ratio (S/N ratio) of each Pseudolite signal in Deceiving interference source is-18dB.Suppression jamming incides array antenna from 60 °, dryly makes an uproar than being 40dB.As can be seen from Figure 3 receiver acquisition to satellite-signal in both comprised real satellite signal and also comprised Pseudolite signal.
Fig. 4 be take after carrying out suppression jamming rectangular projection the inventive method suppress Deceiving interference catch result.In the present embodiment, array antenna is ten yuan of even linear arrays, four gps satellite signal PRN1, and PRN2, PRN3, PRN25 are respectively from-5 °, and 5 ° ,-25 °, 15 ° of directions are incided on array antenna, and signal to noise ratio (S/N ratio) is-20dB.Relay type Deceiving interference source 1 is made up of the Pseudolite signal of PRN1, PRN2, PRN3, PRN25, incides array antenna from 25 °.Production Deceiving interference source 2 is made up of the Pseudolite signal of PRN6, PRN14, PRN22, incides array antenna from-40 °.The signal to noise ratio (S/N ratio) of each Pseudolite signal in Deceiving interference source is-18dB.Suppression jamming incides array antenna from 60 °, dryly makes an uproar than being 40dB.The inventive method successfully curbs Deceiving interference as can be seen from Figure 4, and receiver acquisition is to four real satellite signals.
Fig. 5 does not take the inventive method to suppress the two-dimensional search result of catching PRN1 during Deceiving interference after carrying out suppression jamming rectangular projection.The correlation peak of Deceiving interference and satellite-signal obviously can be observed from Fig. 5.
Fig. 6 takes the inventive method to suppress the two-dimensional search result of catching PRN1 during Deceiving interference after carrying out suppression jamming rectangular projection.After adopting the inventive method to suppress Deceiving interference as can be seen from Figure 6, the correlation peak of Deceiving interference is effectively suppressed, only there is the correlation peak of satellite-signal.

Claims (5)

1., based on a multiclass satellite navigation disturbance restraining method for despreading and respreading technology, it is characterized in that, it comprises the following step carried out in order:
(1) by array antenna received to radiofrequency signal change into digital medium-frequency signal:;
(2) according to the repeatedly fast beat of data that array antenna received arrives, calculate the interference orthogonal complement space, signal suppressed suppression jamming to the projection of the interference orthogonal complement space:
(3) get a road to the data after projection to calculate solution and heavily to extend power vector matrix, and judge that solution is heavily extended power vector matrix whether full rank, to identify whether Deceiving interference exists;
(4) if there is Deceiving interference, the correlativity that solution is heavily extended power between vector is calculated, the pseudo satellite, pseudolite asterisk comprised according to the number in correlated results determination Deceiving interference source, each Deceiving interference source and real satellite asterisk;
(5) utilize the solution of Deceiving interference heavily extend power vector structure Deceiving interference orthogonal intersection cast shadow matrix, project to suppress Deceiving interference to above-mentioned Received signal strength;
(6) utilize the solution of real satellite vector of heavily extending power to carry out Multibeam synthesis, obtain signal transacting gain.
2. the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology according to claim 1, it is characterized in that, in step (3), described a road is got to the data after above-mentioned projection calculate solution and heavily to extend power vector matrix, and judge that solution is heavily extended power vector matrix whether full rank, to identify that method that whether Deceiving interference exists heavily to be extended power vector to the solution that each satellite captured all calculates its correspondence according to despreading and respreading technology, this solution heavily extend power vector include signal come to information, then moment of a vector rank of matrix of heavily extending power is solved, the vector matrix full rank if solution is heavily extended power, then show that satellite-signal is separate each other, from different directions, Deceiving interference is there is not in Received signal strength, otherwise, the vector matrix not full rank if solution is heavily extended power, then show have from unidirectional satellite-signal in Received signal strength, there is Deceiving interference.
3. the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology according to claim 1, it is characterized in that, in step (4), the correlativity that described calculating solution is heavily extended power between vector, according to correlated results determination Deceiving interference source number, the method of the pseudo satellite, pseudolite asterisk that each Deceiving interference source comprises and real satellite asterisk is classified according to heavily the extend power vector of the solution of satellite-signal heavily being extended power that compares of normalized correlation coefficient and thresholding between vector of solution, if exist from unidirectional different satellite-signal, the correlativity that solution then corresponding to these satellite-signals is heavily extended power between vector is very high, think Deceiving interference, otherwise, think real satellite signal.
4. the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology according to claim 1, it is characterized in that, in step (5), the described solution of Deceiving interference that utilizes heavily is extended power the orthogonal intersection cast shadow matrix of vector structure Deceiving interference, project to suppress the method for Deceiving interference to be corresponding with Deceiving interference subspace according to heavily the extend power subspace of vector of the solution of Deceiving interference to above-mentioned Received signal strength, recycling orthogonal complement space computing formula tries to achieve the orthogonal intersection cast shadow matrix of Deceiving interference.
5. the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology according to claim 1, it is characterized in that, in step (6), the described solution of the real satellite vector of heavily extending power that utilizes carries out Multibeam synthesis, the method obtaining signal transacting gain heavily extends power vectors converge in the array manifold of real satellite signal according to the solution of real satellite, can be formed further array beams point to corresponding satellite come to.
CN201410606254.8A 2014-11-01 2014-11-01 Multiclass satellite navigation interference suppression method based on solution expanding technology Expired - Fee Related CN104330808B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410606254.8A CN104330808B (en) 2014-11-01 2014-11-01 Multiclass satellite navigation interference suppression method based on solution expanding technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410606254.8A CN104330808B (en) 2014-11-01 2014-11-01 Multiclass satellite navigation interference suppression method based on solution expanding technology

Publications (2)

Publication Number Publication Date
CN104330808A true CN104330808A (en) 2015-02-04
CN104330808B CN104330808B (en) 2017-01-25

Family

ID=52405569

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410606254.8A Expired - Fee Related CN104330808B (en) 2014-11-01 2014-11-01 Multiclass satellite navigation interference suppression method based on solution expanding technology

Country Status (1)

Country Link
CN (1) CN104330808B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629266A (en) * 2015-12-30 2016-06-01 中国电子科技集团公司第五十四研究所 Blind-adaptive combined suppression method for satellite navigation cheating jamming and suppressing jamming
CN107179540A (en) * 2017-05-26 2017-09-19 中国民航大学 Based on the GNSS vector receiver anti-interference methods for solving re-spread algorithm
CN109425875A (en) * 2017-09-05 2019-03-05 清华大学 Satellite-signal separating and treating apparatus and method
CN110244323A (en) * 2019-05-24 2019-09-17 中国科学院光电研究院 Micro-, light-duty unmanned plane GNSS anti-spoofing system and curve detection and air navigation aid
CN112255649A (en) * 2020-10-31 2021-01-22 江苏集萃未来城市应用技术研究所有限公司 Subspace projection-based method for weakening near-far effect of pseudolite
CN112904379A (en) * 2021-01-28 2021-06-04 中国民航大学 GNSS induced deception jamming data generation method based on track
CN113534198A (en) * 2021-06-16 2021-10-22 北京遥感设备研究所 Satellite navigation dynamic anti-interference method and system based on covariance matrix reconstruction
CN113608242A (en) * 2021-06-18 2021-11-05 西安空间无线电技术研究所 Navigation signal security enhancement method based on code period spread spectrum code authentication
CN114609651A (en) * 2022-03-28 2022-06-10 电子科技大学 Space domain anti-interference method of satellite navigation receiver based on small sample data
CN115951374A (en) * 2023-03-13 2023-04-11 钛玛科(北京)工业科技有限公司 Signal interference suppression method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477189A (en) * 2009-01-22 2009-07-08 中国民航大学 Multiple kinds of interference suppression method of universal satellite navigation system
CN101799551A (en) * 2010-02-08 2010-08-11 中国民航大学 Inhibition method of space-time blind self-adaptive GPS (Global Positioning System) interference based on despreading and respreading technology
CN101561490B (en) * 2009-06-01 2012-06-27 中国民航大学 Blind self-adaptive GPS interference suppression method based on codeword structure
US20130002477A1 (en) * 2011-06-29 2013-01-03 Sintayehu Dehnie Methods and systems for detecting GPS spoofing attacks

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477189A (en) * 2009-01-22 2009-07-08 中国民航大学 Multiple kinds of interference suppression method of universal satellite navigation system
CN101561490B (en) * 2009-06-01 2012-06-27 中国民航大学 Blind self-adaptive GPS interference suppression method based on codeword structure
CN101799551A (en) * 2010-02-08 2010-08-11 中国民航大学 Inhibition method of space-time blind self-adaptive GPS (Global Positioning System) interference based on despreading and respreading technology
US20130002477A1 (en) * 2011-06-29 2013-01-03 Sintayehu Dehnie Methods and systems for detecting GPS spoofing attacks

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YUHANG WU 等: "A New GPS Adaptive Anti-jamming Algorithm in Multi-interference", 《2009 WORLD CONGRESS ON COMPUTER SCIENCE AND INFORMATION ENGINEERING》 *
卢丹: ""稳健的全球卫星导航系统抗干扰技术研究"", 《中国博士学位论文全文数据库信息科技辑》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629266A (en) * 2015-12-30 2016-06-01 中国电子科技集团公司第五十四研究所 Blind-adaptive combined suppression method for satellite navigation cheating jamming and suppressing jamming
CN107179540A (en) * 2017-05-26 2017-09-19 中国民航大学 Based on the GNSS vector receiver anti-interference methods for solving re-spread algorithm
CN107179540B (en) * 2017-05-26 2020-02-18 中国民航大学 GNSS vector receiver anti-interference method based on despreading algorithm
CN109425875A (en) * 2017-09-05 2019-03-05 清华大学 Satellite-signal separating and treating apparatus and method
CN109425875B (en) * 2017-09-05 2021-05-04 清华大学 Satellite signal separation processing device and method
CN110244323B (en) * 2019-05-24 2021-04-20 中国科学院光电研究院 GNSS anti-spoofing system of micro and light unmanned aerial vehicle and spoofing signal detection and navigation method
CN110244323A (en) * 2019-05-24 2019-09-17 中国科学院光电研究院 Micro-, light-duty unmanned plane GNSS anti-spoofing system and curve detection and air navigation aid
CN112255649A (en) * 2020-10-31 2021-01-22 江苏集萃未来城市应用技术研究所有限公司 Subspace projection-based method for weakening near-far effect of pseudolite
CN112904379A (en) * 2021-01-28 2021-06-04 中国民航大学 GNSS induced deception jamming data generation method based on track
CN113534198A (en) * 2021-06-16 2021-10-22 北京遥感设备研究所 Satellite navigation dynamic anti-interference method and system based on covariance matrix reconstruction
CN113534198B (en) * 2021-06-16 2023-05-23 北京遥感设备研究所 Satellite navigation dynamic anti-interference method and system based on covariance matrix reconstruction
CN113608242A (en) * 2021-06-18 2021-11-05 西安空间无线电技术研究所 Navigation signal security enhancement method based on code period spread spectrum code authentication
CN113608242B (en) * 2021-06-18 2023-08-11 西安空间无线电技术研究所 Navigation signal security enhancement method based on code period spread spectrum code authentication
CN114609651A (en) * 2022-03-28 2022-06-10 电子科技大学 Space domain anti-interference method of satellite navigation receiver based on small sample data
CN115951374A (en) * 2023-03-13 2023-04-11 钛玛科(北京)工业科技有限公司 Signal interference suppression method

Also Published As

Publication number Publication date
CN104330808B (en) 2017-01-25

Similar Documents

Publication Publication Date Title
CN104330808A (en) Multiclass satellite navigation interference suppression method based on solution expanding technology
CN104330809A (en) Multi-information-source estimation based satellite navigation deception jamming inhibition method
US10739466B2 (en) Mitigation of spoofer satellite signals
CN101477189B (en) Multiple kinds of interference suppression method of universal satellite navigation system
CN102393525B (en) Navigation interference suppression and signal amplification method for subspace projection
CN101561490B (en) Blind self-adaptive GPS interference suppression method based on codeword structure
CN104536018B (en) GNSS multi-satellite unified capture method using array antenna anti-interference technology
CN103941267A (en) Satellite navigation deception interference suppression method combined with denoising and DOA estimation
CN105629266B (en) Formula is cheated in satellite navigation and pressing type disturbs the joint suppressing method of blind adaptive
CN104155662A (en) Self-adaptive mutual interference restraining method based on GNSS (global navigation satellite system) related peak value detector
CN102707270B (en) Automatic estimation method for antenna pattern relative to high-frequency ground wave radar
CN104391305A (en) Satellite navigation deception interference suppression method based on deception interference DOA estimation
CN104656104A (en) Satellite navigation deceptive signal identification method and system based on maximum likelihood estimation
CN109143265B (en) GNSS anti-spoofing method based on spatial correlation identification
CN102353947A (en) Method for estimating target echo signal subspaces of passive radars based on CSA-MWF (correlation subtraction algorithm-multistage wiener filter)
CN104181556A (en) BOC modulating signal capturing method based on overlapped difference circulation coherent integration
Xu et al. DOA classification and CCPM-PC based GNSS spoofing detection technique
Elghamrawy et al. High-resolution spectral estimation for continuous wave jamming mitigation of GNSS signals in autonomous vehicles
He et al. Civilian unmanned aerial vehicle vulnerability to gps spoofing attacks
Ortega et al. Theoretical evaluation of the GNSS synchronization performance degradation under interferences
CN101788675B (en) GPS (Global Positioning System) interference inhibition method based on single-channel multi-delay cross-correlation treatment
CN103885031A (en) Moving satellite interference source positioning method based on searching optimization filtering
CN104101868A (en) Jamming subspace reconstruction-based radar multi-false target jamming suppression method
CN101819277B (en) GPS (Global Position System) interference suppression method based on single channel single delay mutual-correlation processing
CN104297723A (en) Moving satellite interference source positioning method based on searching optimization and nonlinear filtering

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170125

Termination date: 20181101