CN104716982A - Robust anti-interference processing method and device of spread frequency system - Google Patents

Abstract

The invention provides a robust anti-interference processing method and device of a spread frequency system. The method includes the steps that N sections of intermediate frequency AD sampling signals of which the length is K are accumulated for M receiving array elements of a space frequency two-dimensional processor respectively, and K-point fast Fourier transformation is performed on the intermediate frequency AD sampling signals to obtain a frequency domain array; 102, the number B of frequency points, corresponding a whole bandwidth, of a satellite signal bandwidth is determined; 103, covariance matrixes Rb of the B frequency points in the signal bandwidth are calculated; 104, characteristic decomposition is performed on the covariance matrix Rb of each frequency point, so that the characteristic decomposition expression of each covariance matrix Rb is obtained; 105, according to the characteristic decomposition expression of each covariance matrix, the weight vector of the corresponding frequency point under the constraint of a diagonal loading matrix is determined; 106, the weight vectors wb of all the frequency points are used for performing spatial filtering processing on the frequency domain array; 107, K-point fast Fourier inversion transformation is performed on frequency domain data after spatial filtering, so that time-domain intermediate frequency data are obtained.

Description

A kind of sane anti-interference processing method of spread spectrum system and device

Technical field

The present invention relates to the communications field, particularly relate to the sane anti-interference processing method of a kind of spread spectrum system and device.

Background technology

The arrowband of high power density and broad-band interference have become destroys the topmost factor of military navigation system, array antenna zeroing technology improves the main method of satellite navigation receiver antijamming capability, it is when suppressing space interference, on the arrival direction of interference, zero point is formed, to eliminate space interference by the renewal of weight vector.But, simple airspace filter is subject to the restriction of array freedom number, can not meet the application under complicated interference and multi-path environment, and the complexity that the Space-time domain Beamforming Method adding the time domain degree of freedom calculates self adaptation power improves, and causes AF panel real-time to be deteriorated.In practical application, due to the Rapid Variable Design of the vibration of antenna receiving platform or motion, interference position and relatively too slow etc. the reason of the renewal speed of adaptive weight, the angle of arrival of interference can produce slow change in time during Weight Training, mismatch phenomenon is there is between these data just causing the data of Weight Training and weights to be applied, interference is probably shifted out zero and is fallen into position thus can not effectively be offseted, conventional method complete failure under serious conditions.

Summary of the invention

The invention provides the sane anti-interference processing method of a kind of spread spectrum system and device, the technical problem that solve is the problem of the AF panel when data mismatch.

For solving the problems of the technologies described above, the invention provides following technical scheme:

The sane anti-interference processing method of a kind of spread spectrum system, comprising:

101, receive array element to M of sky frequency two-dimensional process device, accumulate the intermediate frequency AD sampled signal that N segment length is K respectively, respectively the conversion of K point quick Fourier is carried out to described intermediate frequency AD sampled signal and obtain frequency domain array:

X ₁₁,X ₁₂,…X _1K；X ₂₁,X ₂₂,…X _2K；…X _M1,X _M2,…X _MK；

The sequence of each frequency-region signal to be all length be N in described frequency domain array;

102, the frequency points B of the corresponding whole bandwidth of satellite-signal bandwidth is determined;

103, according to described frequency domain array, the covariance matrix R of B frequency in signal bandwidth is calculated _b, wherein R _b=E [X _bx _b ^h], wherein X _b=[X _1k, X _2b... X _mB] ^t, b=1,2 ..., B;

104, respectively to the described covariance matrix R of each frequency _bcarry out feature decomposition, obtain each covariance matrix R _bfeature decomposition expression formula;

105, respectively according to each covariance matrix R _bfeature decomposition expression formula, determine that diagonal angle loads the weight vector w of this frequency under matrix restraint _b;

106, the described weight vector w of each frequency is utilized _bairspace filter process is carried out to described frequency domain array;

107, the inverse transformation of K point quick Fourier is carried out to the frequency domain data after airspace filter, obtain time domain intermediate frequency data.

Said method also has following features: described each covariance matrix R _bfeature decomposition expression formula be:

$R_b=U_b{\mathrm{\λ}}_b{U_b}^H=\underset{i=0}{\overset{L}{\mathrm{\Σ}}}({\mathrm{\λ}}_i+{\mathrm{\σ}}^2)q_i{q}_i^H{+\mathrm{\σ}}^2\underset{i=L+1}{\overset{M}{\mathrm{\Σ}}}{q}_i{q}_i^H;$

Wherein, λ _bcharacteristic value space, U _bit is characteristic value characteristic of correspondence vector space;

Wherein, M is the sum receiving array element, and L is the sum at interference incident direction angle, σ ²noise variance, λ _icharacteristic value, q _iit is characteristic value characteristic of correspondence vector.

Said method also has following features: the expression formula that described diagonal angle loads matrix is:

${\left[T\right]}_{k,l}=\exp \{-{0.5\mathrm{\σ}}_{\max }^2{[(k-l)\mathrm{\π}/180]}^2\};$

Wherein, k represents that diagonal angle loads the line number of matrix, and l represents that diagonal angle loads matrix column number, it is the maximum of the variance at interference incident direction angle.

Said method also has following features: the weight vector w under under described diagonal angle loading matrix restraint _bfor:

Wherein, I _m*Mrepresent the unit matrix of M dimension, s ₀it is column vector ° to be multiplied for matrix corresponding element.

Said method also has following features: the output signal obtained after filtering is:

y _k(n)＝w _1kX _1k(n)+w _2kX _2k(n)+...+w _MkX _Mk(n)；

Wherein, y _k(n) for the output of a kth frequency airspace filter, n be the sequence number in moment.

The sane anti-interference process device of a kind of spread spectrum system, comprising:

Conversion module, for M the reception array element to sky frequency two-dimensional process device, accumulates the intermediate frequency AD sampled signal that N segment length is K respectively, carries out the conversion of K point quick Fourier respectively obtain frequency domain array to described intermediate frequency AD sampled signal:

X ₁₁,X ₁₂,…X _1K；X ₂₁,X ₂₂,…X _2K；…X _M1,X _M2,…X _MK；

The sequence of each frequency-region signal to be all length be N in described frequency domain array;

Determination module, for determining the frequency points B of the corresponding whole bandwidth of satellite-signal bandwidth;

Matrix computations module, for according to described frequency domain array, calculates the covariance matrix R of B frequency in signal bandwidth _b;

Wherein R _b=E [X _bx _b ^h], wherein X _b=[X _1k, X _2b... X _mB] ^t, b=1,2 ..., B;

Matrix decomposition module, for respectively to the described covariance matrix R of each frequency _bcarry out feature decomposition, obtain each covariance matrix R _bfeature decomposition expression formula;

Weight vector determination module, for respectively according to each covariance matrix R _bfeature decomposition expression formula, determine that diagonal angle loads the weight vector w of this frequency under matrix restraint _b;

Airspace filter module, for utilizing the described weight vector w of each frequency _bairspace filter process is carried out to described frequency domain array;

Inverse transform module, for carrying out the inverse transformation of K point quick Fourier to the frequency domain data after airspace filter, obtains time domain intermediate frequency data.

Said apparatus also has following features: described each covariance matrix R _bfeature decomposition expression formula be:

$R_b=U_b{\mathrm{\λ}}_b{U_b}^H=\underset{i=0}{\overset{L}{\mathrm{\Σ}}}({\mathrm{\λ}}_i+{\mathrm{\σ}}^2)q_i{q}_i^H{+\mathrm{\σ}}^2\underset{i=L+1}{\overset{M}{\mathrm{\Σ}}}{q}_i{q}_i^H;$

Wherein, λ _bcharacteristic value space, U _bit is characteristic value characteristic of correspondence vector space;

Wherein, M represents the sum receiving array element, and L represents the sum at interference incident direction angle, σ ²represent noise variance, λ _irepresentation feature value, q _irepresentation feature value characteristic of correspondence vector.

Said apparatus also has following features: the expression formula that described diagonal angle loads matrix is:

${\left[T\right]}_{k,l}=\exp \{-{0.5\mathrm{\σ}}_{\max }^2{[(k-l)\mathrm{\π}/180]}^2\};$

Wherein, k represents that diagonal angle loads the line number of matrix, and l represents that diagonal angle loads matrix column number, represent the maximum of the variance at interference incident direction angle.

Said apparatus also has following features: the weight vector w under under described diagonal angle loading matrix restraint _bfor:

Wherein, I _m*Mrepresent the unit matrix of M dimension, s ₀represent column vector ° representing matrix corresponding element is multiplied.

Said apparatus also has following features: the output signal obtained after filtering is:

y _k(n)＝w _1kX _1k(n)+w _2kX _2k(n)+...+w _MkX _Mk(n)；

Wherein, y _k(n) for the output of a kth frequency airspace filter, n be the sequence number in moment.

The invention provides embodiment, the present invention utilizes empty self-adaptive processing structure frequently, by carrying out feature decomposition to the covariance matrix of each frequency, make filter weights Fast Convergent, simultaneously covariance matrix and diagonal angle load matrix and are weighted, and calculate the weight vector of each frequency, ensure that the width of disturbance null, still effectively can suppress broad-band interference when data mismatch, enhance the robustness of spatial domain filter algorithms.

Accompanying drawing explanation

Fig. 1 is the flow chart of the sane anti-interference processing method of spread spectrum system provided by the invention;

Fig. 2 is the structure chart of the sane anti-interference process device of spread spectrum system provided by the invention.

Embodiment

For making the object, technical solutions and advantages of the present invention clearly, the present invention is described in further detail below in conjunction with the accompanying drawings and the specific embodiments.It should be noted that, when not conflicting, the embodiment in the application and the feature in embodiment can combination in any mutually.

Fig. 1 is the flow chart of the sane anti-interference processing method of spread spectrum system provided by the invention.Method shown in Fig. 1, comprising:

101, receive array element to M of sky frequency two-dimensional process device, accumulate the intermediate frequency AD sampled signal that N segment length is K respectively, respectively the conversion of K point quick Fourier is carried out to described intermediate frequency AD sampled signal and obtain frequency domain array:

X ₁₁,X ₁₂,…X _1K；X ₂₁,X ₂₂,…X _2K；…X _M1,X _M2,…X _MK；

The sequence of each frequency-region signal to be all length be N in described frequency domain array;

Wherein X ₁₁refer to the 1st frequency-region signal receiving the 1st sampled point in array element, X _mKrefer to M the frequency-region signal receiving K sampled point in array element, other by that analogy; The sequence of each frequency-region signal to be all length be N in described frequency domain array;

102, the frequency points B of the corresponding whole bandwidth of satellite-signal bandwidth is determined;

103, according to described frequency domain array, the covariance matrix R of B frequency in signal bandwidth is calculated _b, wherein R _b=E [X _bx _b ^h], wherein X _b=[X _1k, X _2b... X _mB] ^t, b=1,2 ..., B;

104, respectively to the described covariance matrix R of each frequency _bcarry out feature decomposition, obtain each covariance matrix R _bfeature decomposition expression formula;

105, respectively according to each covariance matrix R _bfeature decomposition expression formula, determine that diagonal angle loads the weight vector w of this frequency under matrix restraint _b;

106, the described weight vector w of each frequency is utilized _bairspace filter process is carried out to described frequency domain array;

107, the inverse transformation of K point quick Fourier is carried out to the frequency domain data after airspace filter, obtain time domain intermediate frequency data.

In an embodiment of the present embodiment, described each covariance matrix R _bfeature decomposition expression formula be:

$R_b=U_b{\mathrm{\λ}}_b{U_b}^H=\underset{i=0}{\overset{L}{\mathrm{\Σ}}}({\mathrm{\λ}}_i+{\mathrm{\σ}}^2)q_i{q}_i^H{+\mathrm{\σ}}^2\underset{i=L+1}{\overset{M}{\mathrm{\Σ}}}{q}_i{q}_i^H;$

Wherein, λ _brepresentation feature value space, U _brepresentation feature value characteristic of correspondence vector space;

Wherein, M represents the sum receiving array element, and L represents the sum at interference incident direction angle, σ ²represent noise variance, λ _irepresentation feature value, q _irepresentation feature value characteristic of correspondence vector.

This shows, the present invention utilizes the expression formula after covariance matrix feature decomposition, calculates the weight vector of each frequency, substituted for iterative computation mode of the prior art, makes filter weights Fast Convergent.

In an embodiment of the present embodiment, the expression formula that described diagonal angle loads matrix is:

${\left[T\right]}_{k,l}=\exp \{-{0.5\mathrm{\σ}}_{\max }^2{[(k-l)\mathrm{\π}/180]}^2\};$

Wherein, k represents that diagonal angle loads the line number of matrix, and l represents that diagonal angle loads matrix column number, represent the maximum of the variance at interference incident direction angle.

Below the acquisition of above-mentioned diagonal angle loading matrix is described:

The present invention utilizes the interference noise covariance matrix of the half-wavelength linear array of M array element, obtains the diagonal angle loading matrix that zero of interference falls into width.Be described as follows:

The expression formula of the interference noise covariance matrix of the half-wavelength linear array of M array element is:

$R=\underset{p=1}{\overset{P}{\mathrm{\Σ}}}{r}_{p}a\left({\mathrm{\θ}}_p\right)a{\left({\mathrm{\θ}}_p\right)}^H+{\mathrm{\σ}}^{2}I;$

Wherein, P is the sum of interference signal, r _pbe the frequency of p interference signal, θ _pa (θ _p) representation feature value characteristic of correspondence vector, σ ²represent noise variance, I represents the unit matrix that M ties up.

Under interference incidence angle adaptation conditions, can be expressed as:

wherein △ θ _pobeying average is 0, and variance is σ ²normal distribution, have probability theory knowledge to obtain: ${\widehat{\mathrm{\σ}}}_p={\mathrm{\σ}}_p{\mathrm{\π}}^2\cos \left({\mathrm{\θ}}_p\right)/180,$ Namely $\mathrm{\π}\sin \left({\widehat{\mathrm{\θ}}}_p\right)\∈N(\mathrm{\π}\sin \left({\mathrm{\θ}}_p\right),{\widehat{\mathrm{\σ}}}_p^2),$ ? be designated as v _p, can covariance matrix be constructed as follows:

${\left[\hat{R}\right]}_{k,l}=\underset{p=1}{\overset{P}{\mathrm{\Σ}}}{r}_p\∫f\left(v_p\right)a\left({\mathrm{\θ}}_p\right)a{\left({\mathrm{\θ}}_p\right)}^H{\mathrm{dv}}_p+{\mathrm{\σ}}^{2}I$

(k, l) individual element be

${\left[\hat{R}\right]}_{k,l}=\underset{p=1}{\overset{P}{\mathrm{\Σ}}}{r}_p\underset{-\∞}{\overset{+\∞}{\∫}}\frac{1}{{\widehat{\mathrm{\σ}}}_p\sqrt{2\mathrm{\π}}}\exp [-{(v_p-\mathrm{\π}\sin {\mathrm{\θ}}_p)}^2/2{\widehat{\mathrm{\σ}}}_p^2]\exp [-j(k-l)v_p]{\mathrm{dv}}_p+{\mathrm{\σ}}^2\mathrm{\δ}(k,l)$

Wherein $\mathrm{\δ}(k,l)=\left\{\begin{array}{c}1,k=l\\ 0,l\≠k\end{array}\right.,$ Due to $\underset{-\∞}{\overset{+\∞}{\∫}}\exp \left({-x}^2\right)\mathrm{dx}=\sqrt{\mathrm{\π}},$ Can obtain

${\left[\hat{R}\right]}_{k,l}=\underset{p=1}{\overset{P}{\mathrm{\Σ}}}{r}_p\exp [-j(k-l)\mathrm{\π}\sin \left({\mathrm{\θ}}_p\right)]*\exp [-0.5{\widehat{\mathrm{\σ}}}_p^2{(k-l)}^2]+{\mathrm{\σ}}^2\mathrm{\δ}(k,l),$

? matrix for

${\left[T({\mathrm{\θ}}_P,{\mathrm{\σ}}_p^2)\right]}_{k,l}=\exp \{-{0.5\mathrm{\σ}}_p^2{[(k-l)\mathrm{\π}\cos \left({\mathrm{\θ}}_P\right)/180]}^2\},$

Under the less condition of interference incident direction angle change, θ _p→ 0, cos (θ _p) → 1, finally can obtain

The expression formula that zero of interference falls into the diagonal angle loading matrix of width is as follows:

${\left[T\right]}_{k,l}=\exp \{-{0.5\mathrm{\σ}}_{\max }^2{[(k-l)\mathrm{\π}/180]}^2\}$

In an embodiment of the present embodiment, the weight vector w under under described diagonal angle loading matrix restraint _bfor:

Wherein, I _m*Mrepresent the unit matrix of M dimension, s ₀represent column vector ° to be multiplied for matrix corresponding element.

Below the acquisition of above-mentioned weight vector is described:

After the expression formula of feature decomposition obtaining covariance matrix, calculate the weights solution formula of rectangular projection;

Wherein the expression formula of the weights solution formula of rectangular projection is

Wherein s ₀it is column vector

Recycle the diagonal angle got above and load matrix, weights formula can be obtained as follows:

As seen from the above, when calculating weight vector, loaded the result of matrix weights by calculating covariance matrix and diagonal angle, obtain the weight vector of each frequency, ensure that the width of disturbance null, still effectively can suppress broad-band interference when data mismatch, enhance the robustness of spatial domain filter algorithms.

Wherein, the output signal obtained after filtering is:

y _k(n)＝w _1kX _1k(n)+w _2kX _2k(n)+...+w _MkX _Mk(n)；

Wherein, y _k(n) for the output of a kth frequency airspace filter, n be the sequence number in moment, w _1k, w _2k..., w _mkfor the weight vector w of a kth frequency _kmiddle comprised a M weight coefficient.

Wherein, the concrete steps realizing airspace filter process in step 106 are as follows:

(1) each array element receives AD data after identical cascade trap process, converts through K point FFT.Represent that frequency domain array is with X:

X＝[X ₁₁,X ₁₂,…X _1K；X ₂₁,X ₂₂,…X _2K；…X _M1,X _M2,…X _MK] ()

(2) according to frequency domain array, covariance matrix R is calculated _k=E [X _kx _k ^h] (M × M dimension);

X _k＝[X _1k,X _2k,…X _Mk] ^Tk＝1,2,...,K ()

(3) solve the minimum variance equation of linear restriction, comprising:

1) by linearly constrained minimum variance, this processor can be described as following optimization problem:

$\left\{\begin{array}{cc}\mathrm{Min}& E\left\{{\left|y_k\left(n\right)\right|}^2\right\}={w_k}^H{R}_k{w}_k\\ s.t.& a^H{w}_k=1\end{array}\right.k=\mathrm{1,2},...,K$

Wherein, steric direction vector a is expressed as

The solution utilizing method of Lagrange multipliers can derive multiple constraint minimum variance processor is:

w _k＝R _k ^-1a

(4) best initial weights w is utilized _kcarry out airspace filter process to frequency domain array, expression formula is:

y _k(n)＝w _1kX _1k(n)+w _2kX _2k(n)+...+w _MkX _Mk(n)

Y _kn () is the output of a kth frequency airspace filter.

The invention provides embodiment of the method, the present invention utilizes empty self-adaptive processing structure frequently, by carrying out feature decomposition to the covariance matrix of each frequency, make filter weights Fast Convergent, simultaneously covariance matrix and diagonal angle load matrix and are weighted, and calculate the weight vector of each frequency, ensure that the width of disturbance null, still effectively can suppress broad-band interference when data mismatch, enhance the robustness of spatial domain filter algorithms.

Fig. 2 is the structure chart of the sane anti-interference process device of spread spectrum system provided by the invention.Fig. 2 shown device comprises:

Conversion module, for M the reception array element to sky frequency two-dimensional process device, accumulates the intermediate frequency AD sampled signal that N segment length is K respectively, carries out the conversion of K point quick Fourier respectively obtain frequency domain array to described intermediate frequency AD sampled signal:

X ₁₁,X ₁₂,…X _1K；X ₂₁,X ₂₂,…X _2K；…X _M1,X _M2,…X _MK；

The sequence of each frequency-region signal to be all length be N in described frequency domain array;

Determination module, for determining the frequency points B of the corresponding whole bandwidth of satellite-signal bandwidth;

Matrix computations module, for according to described frequency domain array, calculates the covariance matrix R of B frequency in signal bandwidth _b;

Wherein R _b=E [X _bx _b ^h], wherein X _b=[X _1k, X _2b... X _mB] ^t, b=1,2 ..., B;

Matrix decomposition module, for respectively to the described covariance matrix R of each frequency _bcarry out feature decomposition, obtain each covariance matrix R _bfeature decomposition expression formula;

Weight vector determination module, for respectively according to each covariance matrix R _bfeature decomposition expression formula, determine that diagonal angle loads the weight vector w of this frequency under matrix restraint _b;

Airspace filter module, for utilizing the described weight vector w of each frequency _bairspace filter process is carried out to described frequency domain array;

Inverse transform module, for carrying out the inverse transformation of K point quick Fourier to the frequency domain data after airspace filter, obtains time domain intermediate frequency data.

In an embodiment of the present embodiment, described each covariance matrix R _bfeature decomposition expression formula be:

$R_b=U_b{\mathrm{\λ}}_b{U_b}^H=\underset{i=0}{\overset{L}{\mathrm{\Σ}}}({\mathrm{\λ}}_i+{\mathrm{\σ}}^2)q_i{q}_i^H{+\mathrm{\σ}}^2\underset{i=L+1}{\overset{M}{\mathrm{\Σ}}}{q}_i{q}_i^H;$

Wherein, λ _bcharacteristic value space, U _bit is characteristic value characteristic of correspondence vector space;

Wherein, M represents the sum receiving array element, and L represents the sum at interference incident direction angle, σ ²represent noise variance, λ _irepresentation feature value, q _irepresentation feature value characteristic of correspondence vector.

In an embodiment of the present embodiment, the expression formula that described diagonal angle loads matrix is:

${\left[T\right]}_{k,l}=\exp \{-{0.5\mathrm{\σ}}_{\max }^2{[(k-l)\mathrm{\π}/180]}^2\};$

Wherein, k represents that diagonal angle loads the line number of matrix, and l represents that diagonal angle loads matrix column number, represent the maximum of the variance at interference incident direction angle.

In an embodiment of the present embodiment, the weight vector w under under described diagonal angle loading matrix restraint _bfor:

Wherein, I _m*Mrepresent the unit matrix of M dimension, s ₀represent column vector ° representing matrix corresponding element is multiplied.

In an embodiment of the present embodiment, the output signal obtained after filtering is:

y _k(n)＝w _1kX _1k(n)+w _2kX _2k(n)+...+w _MkX _Mk(n)；

Wherein, y _k(n) for the output of a kth frequency airspace filter, n be the sequence number in moment.

The invention provides device embodiment, the present invention utilizes empty self-adaptive processing structure frequently, by carrying out feature decomposition to the covariance matrix of each frequency, make filter weights Fast Convergent, simultaneously covariance matrix and diagonal angle load matrix and are weighted, and calculate the weight vector of each frequency, ensure that the width of disturbance null, still effectively can suppress broad-band interference when data mismatch, enhance the robustness of spatial domain filter algorithms.

One of ordinary skill in the art will appreciate that all or part of step of above-described embodiment can use computer program flow process to realize, described computer program can be stored in a computer-readable recording medium, described computer program (as system, unit, device etc.) on corresponding hardware platform performs, when performing, step comprising embodiment of the method one or a combination set of.

Alternatively, all or part of step of above-described embodiment also can use integrated circuit to realize, and these steps can be made into integrated circuit modules one by one respectively, or the multiple module in them or step are made into single integrated circuit module to realize.Like this, the present invention is not restricted to any specific hardware and software combination.

Each device/functional module/functional unit in above-described embodiment can adopt general calculation element to realize, and they can concentrate on single calculation element, also can be distributed on network that multiple calculation element forms.

Each device/functional module/functional unit in above-described embodiment using the form of software function module realize and as independently production marketing or use time, can be stored in a computer read/write memory medium.The above-mentioned computer read/write memory medium mentioned can be read-only memory, disk or CD etc.

The above; be only the specific embodiment of the present invention, but protection scope of the present invention is not limited thereto, is anyly familiar with those skilled in the art in the technical scope that the present invention discloses; change can be expected easily or replace, all should be encompassed within protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection range described in claim.

Claims (10)

1. the sane anti-interference processing method of spread spectrum system, is characterized in that, comprising:

101, receive array element to M of sky frequency two-dimensional process device, accumulate the intermediate frequency AD sampled signal that N segment length is K respectively, respectively the conversion of K point quick Fourier is carried out to described intermediate frequency AD sampled signal and obtain frequency domain array:

X ₁₁,X ₁₂,…X _1K；X ₂₁,X ₂₂,…X _2K；…X _M1,X _M2,…X _MK；

The sequence of each frequency-region signal to be all length be N in described frequency domain array;

102, the frequency points B of the corresponding whole bandwidth of satellite-signal bandwidth is determined;

103, according to described frequency domain array, the covariance matrix R of B frequency in signal bandwidth is calculated _b, wherein R _b=E [X _bx _b ^h], wherein X _b=[X _1k, X _2b... X _mB] ^t, b=1,2 ..., B;

104, respectively to the described covariance matrix R of each frequency _bcarry out feature decomposition, obtain each covariance matrix R _bfeature decomposition expression formula;

105, respectively according to each covariance matrix R _bfeature decomposition expression formula, determine that diagonal angle loads the weight vector w of this frequency under matrix restraint _b;

106, the described weight vector w of each frequency is utilized _bairspace filter process is carried out to described frequency domain array;

107, the inverse transformation of K point quick Fourier is carried out to the frequency domain data after airspace filter, obtain time domain intermediate frequency data.

2. method according to claim 1, is characterized in that:

Described each covariance matrix R _bfeature decomposition expression formula be:

$R_b=U_b{\mathrm{\λ}}_n{U_b}^H=\underset{i=0}{\overset{L}{\mathrm{\Σ}}}({\mathrm{\λ}}_i+{\mathrm{\σ}}^2)q_i{q}_i^H+{\mathrm{\σ}}^2\underset{i=L+1}{\overset{M}{\mathrm{\Σ}}}{q}_i{q}_i^H;$

Wherein, λ _bcharacteristic value space, U _bit is characteristic value characteristic of correspondence vector space;

Wherein, M is the sum receiving array element, and L is the sum at interference incident direction angle, σ ²noise variance, λ _icharacteristic value, q _iit is characteristic value characteristic of correspondence vector.

3. method according to claim 1, is characterized in that:

The expression formula that described diagonal angle loads matrix is:

${\left[T\right]}_{k,l}=\exp \{-0.5{\mathrm{\σ}}_{\max }^2{[(k-1)\mathrm{\π}/180]}^2\};$

Wherein, k represents that diagonal angle loads the line number of matrix, and l represents that diagonal angle loads matrix column number, it is the maximum of the variance at interference incident direction angle.

4. method as claimed in claim 3, is characterized in that, the weight vector w under under described diagonal angle loading matrix restraint _bfor:

Wherein, I _m*Mrepresent the unit matrix of M dimension, s ₀it is column vector ° to be multiplied for matrix corresponding element.

5. the method according to any one of Claims 1 to 4, is characterized in that, the output signal obtained after filtering is:

y _k(n)＝w _1kX _1k(n)+w _2kX _2k(n)+...+w _MkX _Mk(n)；

Wherein, y _k(n) for the output of a kth frequency airspace filter, n be the sequence number in moment.

6. the sane anti-interference process device of spread spectrum system, is characterized in that, comprising:

Conversion module, for M the reception array element to sky frequency two-dimensional process device, accumulates the intermediate frequency AD sampled signal that N segment length is K respectively, carries out the conversion of K point quick Fourier respectively obtain frequency domain array to described intermediate frequency AD sampled signal:

X ₁₁,X ₁₂,…X _1K；X ₂₁,X ₂₂,…X _2K；…X _M1,X _M2,…X _MK；

The sequence of each frequency-region signal to be all length be N in described frequency domain array;

Determination module, for determining the frequency points B of the corresponding whole bandwidth of satellite-signal bandwidth;

Matrix computations module, for according to described frequency domain array, calculates the covariance matrix R of B frequency in signal bandwidth _b;

Wherein R _b=E [X _bx _b ^h], wherein X _b=[X _1k, X _2b... X _mB] ^t, b=1,2 ..., B;

Matrix decomposition module, for respectively to the described covariance matrix R of each frequency _bcarry out feature decomposition, obtain each covariance matrix R _bfeature decomposition expression formula;

Weight vector determination module, for respectively according to each covariance matrix R _bfeature decomposition expression formula, determine that diagonal angle loads the weight vector w of this frequency under matrix restraint _b;

Airspace filter module, for utilizing the described weight vector w of each frequency _bairspace filter process is carried out to described frequency domain array;

Inverse transform module, for carrying out the inverse transformation of K point quick Fourier to the frequency domain data after airspace filter, obtains time domain intermediate frequency data.

7. device according to claim 6, is characterized in that:

Described each covariance matrix R _bfeature decomposition expression formula be:

$R_b=U_b{\mathrm{\λ}}_n{U_b}^H=\underset{i=0}{\overset{L}{\mathrm{\Σ}}}({\mathrm{\λ}}_i+{\mathrm{\σ}}^2)q_i{q}_i^H+{\mathrm{\σ}}^2\underset{i=L+1}{\overset{M}{\mathrm{\Σ}}}{q}_i{q}_i^H;$

Wherein, λ _bcharacteristic value space, U _bit is characteristic value characteristic of correspondence vector space;

Wherein, M represents the sum receiving array element, and L represents the sum at interference incident direction angle, σ ²represent noise variance, λ _irepresentation feature value, q _irepresentation feature value characteristic of correspondence vector.

8. device according to claim 6, is characterized in that:

The expression formula that described diagonal angle loads matrix is:

${\left[T\right]}_{k,l}=\exp \{-0.5{\mathrm{\σ}}_{\max }^2{[(k-1)\mathrm{\π}/180]}^2\};$

Wherein, k represents that diagonal angle loads the line number of matrix, and l represents that diagonal angle loads matrix column number, represent the maximum of the variance at interference incident direction angle.

9. device as claimed in claim 8, is characterized in that, the weight vector w under under described diagonal angle loading matrix restraint _bfor:

Wherein, I _m*Mrepresent the unit matrix of M dimension, s ₀represent column vector ° representing matrix corresponding element is multiplied.

10. the device according to any one of claim 6 ~ 9, is characterized in that, the output signal obtained after filtering is:

y _k(n)＝w _1kX _1k(n)+w _2kX _2k(n)+...+w _MkX _Mk(n)；

Wherein, y _k(n) for the output of a kth frequency airspace filter, n be the sequence number in moment.