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

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

Multiclass satellite navigation disturbance restraining method based on despreading and respreading technology

Technical field

The present invention relates to the anti-interference field of satellite navigation, more particularly to a kind of led based on the multiclass satellite of despreading and respreading technology Boat disturbance restraining method.

Background technology

GPS (gnss) utilizes artificial satellite to launch navigation signal, so as to for connecing in coverage Receipts machine provides the information such as high-precision position and time at any time and any place, and therefore satellite navigation system has extremely 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 faces.At present, the grinding of satellite navigation AF panel aspect Study carefully the suppression being concentrated mainly on suppression jamming.Wherein power inversion method is to be much larger than satellite-signal using suppression jamming power Power and this feature of noise power, thus to be upwardly formed null in suppression jamming, yet with the work(of Deceiving interference Rate is less than noise level, and therefore power inversion method is not used to the suppression of Deceiving interference.

Existing Deceiving interference suppressing method is to come to, capture link, follow the tracks of link, positioning result and lead from interference The aspects such as avionics literary composition suppress to Deceiving interference.S.daneshmand proposes a kind of Deceiving interference based on double antenna Suppression technology, the receipt signal of two antennas is carried out computing cross-correlation and obtains the coming to angle of Deceiving interference, and then to deception Formula interference is suppressed, but the required known array manifold of this method, to array manifold error sensitive, it is only applicable to single turn The suppression of hairdo Deceiving interference, and require Deceiving interference power higher.G.b.moon et al. proposes using tracking ring All possible deception formula satellite-signal and real satellite signal are followed the tracks of in road, then distinguish deception satellite letter according to tracking parameter Number and real satellite signal.This technology needs more tracking channels, so that the complexity of receiver and cost is greatly increased. History Vincent proposes the doa estimated value of a kind of DR position using receiver and receipt signal and suppresses Deceiving interference side Method, the method needs receiver location is calculated, and navigation message is demodulated, and computation complexity is high. C.e.mcdowell is analyzed to the theory of multi-antenna beam formation and zeroing technology and performance, using gps correlator ginseng Number and track loop parameter form multi-beam to identify and to suppress Deceiving interference, but this technology does not possess in suppression deception formula While interference, real satellite signal is formed with the ability of multi-beam gain.It should be noted that these methods are done in briquettability Disturb equal cisco unity malfunction when coexisting with Deceiving interference.In order to ensure the safety applications of gnss, can know in the urgent need to research Other Deceiving interference, then suppresses suppression jamming, relay type Deceiving interference and production Deceiving interference simultaneously, and right Array manifold error and the insensitive multiclass satellite navigation disturbance restraining method of jamming power, this not only has in terms of theoretical research Important meaning, and the propulsion application in fields such as Military Application for the gnss is also had great importance.

Content of the invention

In order to solve the above problems, it is an object of the invention to provide one kind does not need known array manifold and receiver anti- Feedforward information, can suppress satellite navigation system suppression jamming, identification and suppression relay type Deceiving interference and production deception Formula is disturbed and satellite-signal can be provided with the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology of gain.

In order to achieve the above object, the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology that the present invention provides The following step including carrying out in order:

(1) by array antenna received to radiofrequency signal change into digital medium-frequency signal:

(2) the multiple snapshot data being arrived according to array antenna received, calculates the interference orthogonal complement space, by signal to interference just Hand over complementary space to project to suppress suppression jamming:

(3) take a road to calculate the data after projection and solve re-spread weight vector matrix, and judge that solving re-spread weight vector matrix is No full rank, to identify that Deceiving interference whether there is；

(4) if there is Deceiving interference, calculate the dependency solving between re-spread weight vector, deception is determined according to correlated results Pseudo satellite, pseudolite asterisk and real satellite asterisk that the number of formula interference source, each Deceiving interference source comprise；

(5) utilize the orthogonal intersection cast shadow matrix of the solution re-spread weight vector construction Deceiving interference of Deceiving interference, connect to above-mentioned The collection of letters number is projected to suppress Deceiving interference；

(6) carry out Multibeam synthesis using the re-spread weight vector of solution of real satellite, obtain signal processing gain.

In step (3), a described road that the data after above-mentioned projection is taken calculates the re-spread weight vector matrix of solution, and judges Solve re-spread weight vector matrix whether full rank, to identify that the method that Deceiving interference whether there is is to each according to despreading and respreading technology The satellite capturing all calculates its re-spread weight vector of corresponding solution, and what the re-spread weight vector of this solution included signal comes to information, Then solving re-spread weight vector rank of matrix, if solving re-spread weight vector matrix full rank, showing satellite-signal phase each other Mutually independent, from different directions, in receipt signal, there is not Deceiving interference, whereas if solving re-spread weight vector matrix not Full rank, then show have from unidirectional satellite-signal in receipt signal, there is Deceiving interference.

In step (4), dependency between the described re-spread weight vector of calculating solution, deception formula is determined according to correlated results The method of interference source number, the pseudo satellite, pseudolite asterisk that each Deceiving interference source comprises and real satellite asterisk is according to the re-spread power of solution Between vector, re-spread for the solution of satellite-signal weight vector is classified by normalizated correlation coefficient and the comparison of thresholding, if exist being derived from same side To different satellite-signals, then corresponding to these satellite-signals solve re-spread weight vector between dependency very high it is believed that being to take advantage of Deceive formula interference, conversely, being considered real satellite signal.

In step (5), the described re-spread weight vector of solution using Deceiving interference constructs just trading of Deceiving interference Shadow matrix, is projected the re-spread power of solution to suppress the method for Deceiving interference to be according to Deceiving interference to above-mentioned receipt signal The subspace of vector is corresponding with Deceiving interference subspace, recycles orthogonal complement space computing formula to try to achieve deception formula and does The orthogonal intersection cast shadow matrix disturbed.

In step (6), the described re-spread weight vector of the solution using real satellite carries out Multibeam synthesis, obtains at signal The method of reason gain is to converge on the array manifold of real satellite signal according to the re-spread weight vector of solution of real satellite, further may be used With formed array beamses point to corresponding satellite come to.

The tool that can produce high-gain based on the multiclass satellite navigation disturbance restraining method of despreading and respreading technology that the present invention provides There are multiple wave beams of anti-suppression jamming, relay type Deceiving interference and production Deceiving interference ability simultaneously, multi-beam shape One-tenth ability is not limited by array number, need not know that satellite-signal and interference, to information, array manifold error and deception formula to be done Disturb watt level insensitive it is not necessary to the feedback information of navigation neceiver.

Brief description

Fig. 1 is the beam pattern of the re-spread weight vector matrix of solution being obtained using the inventive method.

Fig. 2 is the total beam pattern all satellite-signals being obtained using the inventive method.

The inventive method of Fig. 3 not taking after being by suppression jamming rectangular projection suppresses the capture knot of Deceiving interference Really.

Fig. 4 is the capture result after taking the inventive method to suppress suppression jamming and Deceiving interference.

Fig. 5 catches to prn1 when not taking the inventive method suppression Deceiving interference after being by suppression jamming rectangular projection The two-dimensional search result obtaining.

Fig. 6 captures to prn1 when taking the inventive method suppression Deceiving interference after being by suppression jamming rectangular projection Two-dimensional search result.

Fig. 7 is the flow chart of the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology that the present invention provides.

Specific embodiment

With reference to the accompanying drawings with specific embodiment to the present invention multiclass satellite navigation AF panel based on despreading and respreading technology Method gives to describe in detail.For convenience's sake, below with gps (global positioning system) system and n battle array It is illustrated as a example the array antenna of unit, in fact, the method is applied to General Cell antenna, and require no knowledge about array antenna Manifold.

As shown in fig. 7, the present invention provide based on the multiclass satellite navigation disturbance restraining method of despreading and respreading technology include by The following step that order is carried out:

(1) by array antenna received to radiofrequency signal change into digital medium-frequency signal:

Deceiving interference is formed and is launched in the same direction by multi-satellite signal (referred to as Pseudolite signal).Wherein, turn Hairdo Deceiving interference is by receiving current satellite-signal, according to deception purpose, satellite-signal is postponed and amplifies, then Forward and to realize disturbing；Production Deceiving interference is using information known to commercial satellite navigation signal form, directly generates Pseudolite signal is simultaneously launched.

When suppression jamming, relay type Deceiving interference and production Deceiving interference are made simultaneously incident on array antenna When, the intermediate-freuqncy signal that gps receiver receives is:

$x\left({\mathrm{nt}}_s\right)=\underset{m=1}{\overset{m}{\mathrm{\σ}}}{a}_m{s}_m^a\left({\mathrm{nt}}_s\right)+\underset{p=1}{\overset{p}{\mathrm{\σ}}}{b}_p{s}_p^{\mathrm{rs}}\left({\mathrm{nt}}_s\right)+\underset{q=1}{\overset{q}{\mathrm{\σ}}}{b}_q{s}_q^{\mathrm{gs}}\left({\mathrm{nt}}_s\right)+\underset{k=1}{\overset{k}{\mathrm{\σ}}}{c}_k{j}_k\left({\mathrm{nt}}_s\right)+n\left({\mathrm{nt}}_s\right)---\left(1\right)$

Wherein, m, p, q and k represent that real satellite signal number, relay type Deceiving interference source number, production are taken advantage of respectively Deceive formula interference source number and suppression jamming number.Represent the navigation signal of m-th real satellite signal, Represent p-th relay type Deceiving interference,Represent q-th production Deceiving interference, subscript a, rs and gs refer to respectively Show real satellite signal, relay type Deceiving interference and production Deceiving interference.Wherein,

$s_m^a\left({\mathrm{nt}}_s\right)=\sqrt{p_m^a}{d}_m^a({\mathrm{nt}}_s-{\mathrm{\τ}}_m^a)c_m^a({\mathrm{nt}}_s-{\mathrm{\τ}}_m^a)e^{j{\mathrm{\φ}}_m^a+j2\mathrm{\π}{f}_m^{\′a}n{t}_s}---\left(2\right)$

$s_p^{\mathrm{rs}}\left({\mathrm{nt}}_s\right)=\underset{l=1}{\overset{l_p}{\mathrm{\σ}}}\sqrt{p_l^{\mathrm{rs}}}{d}_l^{\mathrm{rs}}({\mathrm{nt}}_s-{\mathrm{\τ}}_l^{\mathrm{rs}})c_l^{\mathrm{rs}}({\mathrm{nt}}_s-{\mathrm{\τ}}_l^{\mathrm{rs}})e^{j{\mathrm{\φ}}_l^{\mathrm{rs}}+j2\mathrm{\π}{f}_l^{\′\mathrm{rs}}n{t}_s}---\left(3\right)$

$s_q^{\mathrm{gs}}\left({\mathrm{nt}}_s\right)=\underset{l=1}{\overset{l_q}{\mathrm{\σ}}}\sqrt{p_l^{\mathrm{gs}}}{d}_l^{\mathrm{gs}}({\mathrm{nt}}_s-{\mathrm{\τ}}_l^{\mathrm{gs}})c_l^{\mathrm{gs}}({\mathrm{nt}}_s-{\mathrm{\τ}}_l^{\mathrm{gs}})e^{j{\mathrm{\φ}}_l^{\mathrm{gs}}+j2\mathrm{\π}{f}_l^{\′\mathrm{gs}}n{t}_s}---\left(4\right)$

φ, f ', p and τ represent phase place, Doppler frequency and code time delay respectively.d(nt_s) and c (nt_s) represent navigation respectively Text and c/a code.L in formula (3)_pRepresent the number of the Pseudolite signal that p-th relay type Deceiving interference source is comprised, this A little pseudo satellite, pseudolites are generally identical with real satellite asterisk；L in formula (4)_qRepresent that q-th production Deceiving interference source is comprised The number of Pseudolite signal, these pseudo satellite, pseudolites are not necessarily identical with real satellite asterisk.j_k(nt_s) represent that k-th briquettability is done Disturb, n (nt_s) it is additive white Gaussian noise, a_m、b_p、b_q、c_kIt is corresponding array manifold, t_sFor the sampling interval.

(2) the multiple snapshot data being arrived according to array antenna received, calculates the interference orthogonal complement space, by signal to interference just Hand over complementary space to project to suppress suppression jamming:

Due to the low 20db of power ratio noise level of satellite-signal, Deceiving interference power is suitable with satellite-signal, and presses The power of property processed interference is usually all very strong, noise power typically to be far above, thus the covariance matrix of receipt signal mainly by The covariance matrix of suppression jamming and noise determines.That is:

r_xx≈r_j+r_v(5)

Wherein, r_xxRepresent the covariance matrix of receipt signal, r_vRepresent noise covariance matrix, r_jRepresent suppression jamming Covariance matrix.Now, the covariance matrix of receipt signal is inverted and is directly proportional to the orthogonal complement space of suppression jamming subspace, Receipt signal is projected to suppress suppression jamming to the interference orthogonal complement space:

$y\left({\mathrm{nt}}_s\right)=r_{\mathrm{xx}}^{-1}x\left({\mathrm{nt}}_s\right)---\left(6\right)$

(3) take a road to calculate the data after projection and solve re-spread weight vector matrix w, and judge that solving re-spread weight vector matrix is No full rank, to identify that Deceiving interference whether there is:

W=[w₁w₂… w_l] (7)

Wherein, w_l=[w_l1… w_ln]^tRepresent the re-spread weight vector of solution of l satellite.

Despreading and respreading technology refers to detect n-th number of bits of transmission signal according to b by spreading codes information_n, and utilize b_n Again spread out transmission signal in time cycle [(n-1) t_b,nt_b] in signal model, then re-spread after signal can serve as Reference signal is updating the weight vector of Adaptive beamformer.Using known to gps signal spread-spectrum code codeword structure and navigation data Cycle is much larger than the feature in spreading code cycle, and the cost function solving the re-spread weight vector of solution of l satellite is:

$w_l=\underset{w}{\min }\underset{k=1}{\overset{\stackrel{\&overbar;}{k}}{\mathrm{\σ}}}{|{w_l}^{h}y\left(k\right)-r_l\left(k\right)|}^2---\left(8\right)$

Wherein, r_lK () is the reference signal of l satellite,It is the length of data block, be in the spreading code cycle Sample number.In formula (8), the more new formula of weight vector is:

$w_l(h+1)={\left[y\left(i\right)y^h\left(i\right)\right]}^{-1}{\mathrm{yr}}_l^{*}\left(h\right)---\left(9\right)$

$y\left(i\right)={[y(1+i\stackrel{\&overbar;}{k}),...,y(\stackrel{\&overbar;}{k}+i\stackrel{\&overbar;}{k})]}^t---\left(10\right)$

$r_l\left(h\right)={[c_l(1+h\stackrel{\&overbar;}{k}-{\widehat{\mathrm{\τ}}}_l)e^{j2\mathrm{\π}{\hat{f}}_l(1+h\stackrel{\&overbar;}{k}-{\widehat{\mathrm{\τ}}}_l)},...,c_l(\stackrel{\&overbar;}{k}+h\stackrel{\&overbar;}{k}-{\widehat{\mathrm{\τ}}}_l)e^{j2\mathrm{\π}{\hat{f}}_l(\stackrel{\&overbar;}{k}+h\stackrel{\&overbar;}{k}-{\widehat{\mathrm{\τ}}}_l)}]}^t---\left(11\right)$

Wherein, h represents the iterationses that weight vector calculates,WithIt is the code delay of l satellite and the carrier wave of capture The valuation of frequency.When Deceiving interference exists, the capture correlated results of same asterisk is frequently not a peak value, as shown in Figure 4. Ratio using peak value and average noise is truly defended it is ensured that capturing the whole of same asterisk as the condition of prize judgment Star and pseudo satellite, pseudolite, and obtain the valuation of code delay and carrier frequency, to reconstruct the reference of whole real satellite and pseudo satellite, pseudolite Signal, thus solved re-spread weight vector accordingly.

Algorithm specifically comprises the following steps that

Step 1: acquisition and tracking is carried out to array output, obtains the c of the satellite of capture_l(k)、Etc. parameter, this is solution The process expanding；

Step 2: generate the spread-spectrum signal r of every effective satellite according to step 1 parameters obtained by formula (11) respectively_l(h), this It is re-spread process；

Step 3: respectively with r_lH (), as reference signal, updates the weight vector of l-th satellite using formula (9)；

Step 4: computing array output vector；

Step 5: repeat step 14, until algorithmic statement.

Judge to solve re-spread weight vector matrix whether full rank:

Rank w=m < l (12)

If formula (12) is false, w is described₁,…,w_lOrthogonal between any two, do not exist from unidirectional letter Number, that is, currently capture, there is not Deceiving interference in tracking satellite.If conversely, formula (12) is set up, illustrate currently to capture, with There is Deceiving interference in track satellite.

(4) if there is Deceiving interference, calculate the dependency solving between re-spread weight vector, deception is determined according to correlated results Pseudo satellite, pseudolite asterisk and real satellite asterisk that the number of formula interference source, each Deceiving interference source comprise；

If satellite-signal s₁And s₂The re-spread weight vector of solution be respectively w₁And w₂, ρ represents their normalizated correlation coefficient,

$\mathrm{\&rho;}=\left|{\left(\frac{w_1}{\left|w_1\right|}\right)}^h\frac{w_2}{\left|w_2\right|}\right|---\left(13\right)$

When the normalizated correlation coefficient ρ of the re-spread weight vector of solution of two signals is more than decision threshold d_△When it is believed that this two Individual signal is from same direction, that is,

$\mathrm{\&rho;}=\left|{\left(\frac{w_1}{\left|w_1\right|}\right)}^h\frac{w_2}{\left|w_2\right|}\right|\&greaterequal;d_{\mathrm{\&delta;}}---\left(14\right)$

Because real satellite signal comes to and the pseudo satellite, pseudolite that Deceiving interference comprises all is taken advantage of from same respectively from difference Deceive the transmitting of formula interference source, therefore, can determine whether that normalizated correlation coefficient is more than corresponding to the re-spread weight vector of each solution of decision threshold Satellite is pseudo satellite, pseudolite, and remaining is real satellite signal.

The method of this identification Deceiving interference, does not utilize the power information of Deceiving interference and real satellite signal, right Deceiving interference power is insensitive.

(5) utilize the orthogonal intersection cast shadow matrix of the solution re-spread weight vector construction Deceiving interference of Deceiving interference, connect to above-mentioned The collection of letters number is projected to suppress Deceiving interference；

According to the judged result of step (4), construct Deceiving interference subspace with the re-spread weight vector of the solution of Deceiving interference, It is designated as u=span { w₁,…,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 of the no Deceiving interference after array received signal is projected is:

z(nt_s)=p_⊥y(nt_s) (17)

(6) Multibeam synthesis are carried out using the re-spread weight vector of solution of real satellite, acquisition signal processing gain:

Using the solution of real satellite re-spread weight vector, Multibeam synthesis are carried out to the output vector of formula (17), obtain final Array antenna output signal is:

${\mathrm{output}}_l\left({\mathrm{nt}}_s\right)=w_l^{h}z\left({\mathrm{nt}}_s\right)---\left(18\right)$

Wherein, subscript l represents is output signal to l real satellite.Because acquisition procedure is asking for solution weight Complete during power expansion vector, and obtained the c of satellite_l(k)、Etc. parameter, the output using formula (18) can be directly right The l real satellite is tracked.

Fig. 1 is the beam pattern of the re-spread weight vector matrix of solution being obtained using the inventive method.In the present embodiment, array sky Line be ten yuan of even linear arrays, four gps satellite-signal prn1, prn2, prn3, prn25 respectively from -5 °, 5 °, -25 °, 15 ° of directions Incide on array antenna, signal to noise ratio is -20db.Relay type Deceiving interference source 1 is by prn1, the puppet of prn2, prn3, prn25 Satellite-signal forms, and incides array antenna from 25 °.Production Deceiving interference source 2 is by prn6, the puppet of prn14, prn22 Satellite-signal forms, and incides array antenna from -40 °.The signal to noise ratio of each Pseudolite signal in Deceiving interference source be- 18db.Suppression jamming incides array antenna from 60 °, dry makes an uproar than for 40db.Black dotted lines in Fig. 1 represent real satellite Direction of signal, solid black lines represent Pseudolite signal to.It will be seen from figure 1 that solution re-spread weight vector be aligned whole satellite letter Number come to Deceiving interference come to.

Fig. 2 is the total beam pattern all satellite-signals being obtained using the inventive method.In the present embodiment, array antenna For ten yuan of even linear arrays, four gps satellite-signal prn1, respectively from -5 °, 5 °, -25 °, 15 ° of directions enter for prn2, prn3, prn25 It is mapped on array antenna, signal to noise ratio is -20db.Relay type Deceiving interference source 1 is defended by prn1, the puppet of prn2, prn3, prn25 Star signal forms, and incides array antenna from 25 °.Production Deceiving interference source 2 is defended by prn6, the puppet of prn14, prn22 Star signal forms, and incides array antenna from -40 °.The signal to noise ratio of each Pseudolite signal in Deceiving interference source be- 18db.Suppression jamming incides array antenna from 60 °, dry makes an uproar than for 40db.Black dotted lines in Fig. 2 represent 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 is permissible In satellite-signal, beam gain is upwardly formed to every satellite, and comes to formation zero in suppression jamming and Deceiving interference Fall into.

The inventive method of Fig. 3 not taking after being by suppression jamming rectangular projection suppresses the capture knot of Deceiving interference Really.In the present embodiment, array antenna is ten yuan of even linear arrays, four gps satellite-signal prn1, and prn2, prn3, prn25 divide Not from -5 °, 5 °, -25 °, 15 ° of directions are incided on array, and signal to noise ratio is -20db.Relay type Deceiving interference source 1 by prn1, The Pseudolite signal composition of prn2, prn3, prn25, incides array antenna from 25 °.Production Deceiving interference source 2 by The Pseudolite signal composition of prn6, prn14, prn22, incides array antenna from -40 °.Each in Deceiving interference source The signal to noise ratio of Pseudolite signal is -18db.Suppression jamming incides array antenna from 60 °, dry makes an uproar than for 40db.From Fig. 3 It can be seen that receiver acquisition to satellite-signal in both comprised real satellite signal and comprised Pseudolite signal.

Fig. 4 takes the inventive method to suppress the capture result of Deceiving interference after being by suppression jamming rectangular projection. 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 °, 5 °, -25 °, 15 ° of directions are incided on array antenna, and signal to noise ratio is -20db.Relay type Deceiving interference source 1 by The Pseudolite signal composition of prn1, prn2, prn3, prn25, incides array antenna from 25 °.Production Deceiving interference source 2, by prn6, the Pseudolite signal composition of prn14, prn22, incide array antenna from -40 °.Every in Deceiving interference source The signal to noise ratio of individual Pseudolite signal is -18db.Suppression jamming incides array antenna from 60 °, dry makes an uproar than for 40db.From figure 4 it can be seen that the inventive method successfully curbs Deceiving interference, and receiver acquisition is to four real satellite signals.

Fig. 5 catches to prn1 when not taking the inventive method suppression Deceiving interference after being by suppression jamming rectangular projection The two-dimensional search result obtaining.The correlation peak of Deceiving interference and satellite-signal can substantially be observed from Fig. 5.

Fig. 6 captures to prn1 when taking the inventive method suppression Deceiving interference after being by suppression jamming rectangular projection Two-dimensional search result.After adopting the inventive method suppression Deceiving interference as can be seen from Figure 6, the relevant peaks of Deceiving interference Value is effectively suppressed, only exists the correlation peak of satellite-signal.

Claims (5)

1. a kind of multiclass satellite navigation disturbance restraining method based on despreading and respreading technology is it is characterised in that it includes entering in order The following step of row:

(1) by array antenna received to radiofrequency signal change into digital medium-frequency signal；

(2) the multiple snapshot data being arrived according to array antenna received, calculates the interference orthogonal complement space, by signal to interference orthocomplement, orthogonal complement Space projection is suppressing suppression jamming；

(3) take a road to calculate the data after projection and solve re-spread weight vector matrix, and judge whether completely to solve re-spread weight vector matrix Order, to identify that Deceiving interference whether there is；

(4) if there is Deceiving interference, calculating the dependency solving between re-spread weight vector, determining that deception formula is done according to correlated results Disturb the number in source, pseudo satellite, pseudolite asterisk that each Deceiving interference source comprises and real satellite asterisk；

(5) utilize the orthogonal intersection cast shadow matrix of the solution re-spread weight vector construction Deceiving interference of Deceiving interference, to above-mentioned through interference The signal eliminating suppression jamming after orthogonal complement space projection is projected to suppress Deceiving interference；

(6) carry out Multibeam synthesis using the re-spread weight vector of solution of real satellite, obtain signal processing gain.

2. the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology according to claim 1 it is characterised in that In step (3), a described road that the data after above-mentioned projection is taken calculates the re-spread weight vector matrix of solution, and judges to solve re-spread power Vector matrix whether full rank, to identify that the method that Deceiving interference whether there is is according to despreading and respreading technology, each to be captured Satellite all calculate its corresponding solution re-spread weight vector, the re-spread weight vector of this solution include signal come to information, then solve Re-spread weight vector rank of matrix, if solving re-spread weight vector matrix full rank, shows that satellite-signal is separate each other, comes From different directions, in receipt signal, there is not Deceiving interference, whereas if solving re-spread weight vector matrix not full rank, then table Have from unidirectional satellite-signal in bright receipt signal, there is Deceiving interference.

3. the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology according to claim 1 it is characterised in that In step (4), dependency between the described re-spread weight vector of calculating solution, Deceiving interference source is determined according to correlated results The pseudo satellite, pseudolite asterisk that several, each Deceiving interference source comprises and the method for real satellite asterisk are to return according between the re-spread weight vector of solution Re-spread for the solution of satellite-signal weight vector is classified by the comparison of one change correlation coefficient and thresholding, if exist being derived from unidirectional difference Satellite-signal, then corresponding to these satellite-signals solve re-spread weight vector between dependency very high it is believed that being Deceiving interference, Conversely, being considered real satellite signal.

4. the multiclass satellite navigation disturbance restraining method based on despreading and respreading technology according to claim 1 it is characterised in that In step (5), the re-spread weight vector of solution of described utilization Deceiving interference constructs the orthogonal intersection cast shadow matrix of Deceiving interference, right The signal eliminating suppression jamming after the above-mentioned orthogonal complement space projection through interference is projected to suppress the side of Deceiving interference Method is to solve that the subspace opened of re-spread weight vector is corresponding with Deceiving interference subspace according to Deceiving interference, and recycling is just Complementary space computing formula is handed over to try 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 characterised in that In step (6), the described re-spread weight vector of the solution using real satellite carries out Multibeam synthesis, obtains signal processing gain Method is to converge on the array manifold of real satellite signal according to the re-spread weight vector of solution of real satellite, can form battle array further The corresponding satellite of row beam position come to.