CN103546191B - Universe covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE - Google Patents

Abstract

The invention discloses the main synchronizing sequence Interference Cancellation detection method that a kind of universe covers multi-beam S-LTE, to time-domain received signal low-pass filtering and down-sampling; Respectively by belong to same main synchronizing sequence of satellite be added obtain slide relevant with sequence to time domain receiving sequence, peakvalue's checking, obtain the type information that time domain receives satellite belonging to the original position of main synchronizing sequence and Serving cell, calculate time domain receive main synchronizing sequence respectively with satellite belonging to Serving cell with sequence and jammer satellite with the ratio of the correlation energy value of sequence; Utilization and sequence are carried out frequency deviation and are estimated and compensate; Ratio according to correlation energy value carries out threshold judgement, if be greater than threshold value, then respectively local time-domain primary synchronization sequence is received main synchronizing sequence with the time domain through compensate of frequency deviation and carry out related operation, if be less than threshold value, then above-mentioned correlation is deducted interference correlation, peakvalue's checking, obtains main synchronizing sequence.The present invention resists larger frequency deviation, computation complexity is low and precision is high.

Description

Universe covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE

Technical field

The present invention relates to the main synchronizing sequence Interference Cancellation detection method that a kind of universe covers multi-beam S-LTE, belong to broadband wireless communication technique field.

Background technology

Satellite mobile communication is one of necessary means realizing carrying out in any place communicating.In recent years, forth generation (The Fourth Generation Mobile Communication Systems, 4G) land honeycomb Mobile Communication Development reaches its maturity, with OFDM (Orthogonal Frequency Division Multiplexing, OFDM), multiple-input and multiple-output (Multiple-Input Multiple-Output, MIMO) and identical networking etc. as third generation partner program Long Term Evolution (the The 3rd GenerationPartnership Project Long Term Evolution of key technology, 3GPP-LTE) put it into commercial operation successively.To there is the land LTE(Terrestrial LTE of the features such as two-forty, Large Copacity, spectral efficient, high power efficiency, T-LTE) be applied in satellite mobile communication, set up multi-beam S-LTE(Satellite LTE, the S-LTE of identical networking) mobile communication system be present satellites moving communicating field research focus and difficult point.

The multi-beam S-LTE mobile communication system of identical networking adopts large-scale antenna array to produce multiple wave beam on satellite, forms multiple community after arriving ground.The same with T-LTE, neighbor cell configures different main synchronizing sequences.But the multi-beam S-LTE of identical networking also exists the inter-beam interference (Inter-Beam Interference, IBI) even more serious compared with T-LTE, and the overlapping ranges of neighbor cell is larger.Multi-beam S-LTE adopts multi-satellite (may be dissimilar satellite) common networking, thus realizes universe covering.In order to reduce the interference adopted between identical main synchronizing sequence allocating cell, and make terminal can identification satellite, thus realize Cell searching better, the main synchronizing sequence that different satellite configuration is different, the signal of same satellite arrives same ground receiver end through identical fading channel, and the signal of different satellite then arrives same ground receiver end through different fading channels.

Traditional main synchronizing sequence detects the cross-correlation test algorithm adopted based on local time-domain primary synchronization sequence, this algorithm has the following disadvantages: (1) this algorithm needs all local time-domain primary synchronization sequences to slide relevant to time domain receiving sequence respectively, the main synchronizing sequence quantity that the multi-beam S-LTE that universe covers configures is larger, therefore, computation complexity is high; (2) this algorithm cannot resist larger frequency deviation, although can improve anti-deviation capability by the method for piecemeal cross-correlation, needs to increase computation complexity as cost.

Summary of the invention

Goal of the invention: cover the configuration feature of multi-beam S-LTE system main synchronizing sequence and the deficiency of traditional cross-correlation test algorithm based on local time-domain primary synchronization sequence for universe, the invention provides and a kind ofly can resist larger frequency deviation, main synchronizing sequence Interference Cancellation detection method that low computation complexity and high-precision universe cover multi-beam S-LTE.

Technical scheme: a kind of universe covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, comprises the steps:

(1) to time domain receiving sequence low-pass filtering and down-sampling;

(2) what main for same satellite synchronizing sequence addition obtained respectively slides relevant with sequence to time domain receiving sequence, judge correlation peak location, obtain the information that time domain receives satellite belonging to the original position of main synchronizing sequence and Serving cell, further, calculate time domain receive main synchronizing sequence respectively with satellite belonging to Serving cell with sequence and jammer satellite with the ratio of the correlation energy value of sequence;

(3) time domain is utilized to receive the carrying out frequency deviation estimation with sequence and compensating of satellite belonging to main synchronizing sequence and Serving cell;

(4) ratio according to correlation energy value carries out threshold judgement, select the detection method of main synchronizing sequence that Serving cell configures: if be greater than threshold value, then respectively the local time-domain primary synchronization sequence of satellite belonging to Serving cell is received main synchronizing sequence with the time domain through compensate of frequency deviation and carry out related operation; If be less than threshold value, then the correlation that the local time-domain primary synchronization sequence of satellite belonging to Serving cell and the time domain through compensate of frequency deviation receive main synchronizing sequence deducted interference correlation, judge correlation peak, the main synchronizing sequence that acquisition Serving cell configures.

In described step (2), main synchronizing sequence sum operation formula is:

${\tilde{s}}_p\left(n\right)=\underset{i=0}{\overset{N_u-1}{\mathrm{\Σ}}}{s}_{p,u\left(i\right)}\left(n\right)$

Wherein, be p satellite and sequence, p=0,1 ..., P-1; s _{p, u (i)}the local time-domain primary synchronization sequence of (n) to be the root sequence number of p satellite be u (i), i=0,1 ..., N _u-1, N _ufor root sequence number sum, i.e. sequence sum, n is main synchronizing sequence sampled point sequence number, n=0,1 ..., N-1.

The operational formula that described step (2) slip is relevant is:

$M_p\left(\mathrm{\θ}\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}r(n+\mathrm{\θ}){\tilde{s}}_p^{*}\left(n\right)$

Wherein, M _p(θ) be the result of related operation, θ is the original position of sliding window; R (n) is time domain receiving sequence; () ^*represent conjugation.

Belonging to the original position of described step (2) main synchronizing sequence and Serving cell, the judgment formula of the information of satellite is:

$\{\widehat{\mathrm{\θ}},\hat{p}\}=\arg \underset{p,\mathrm{\θ}}{\max }\left\{{\left|M_p\left(\mathrm{\θ}\right)\right|}^2\right\}$

Wherein, with be respectively the estimated value of satellite belonging to the original position of main synchronizing sequence and Serving cell; | x| ²represent the energy asking x; representing to get makes x be the p of maximum, θ value.

In described step (2) time domain receive main synchronizing sequence respectively with satellite belonging to Serving cell with sequence and jammer satellite with the computing formula of the ratio of the correlation energy value of sequence be:

$P_r=\frac{{\left|M_{\hat{p}}\left(\widehat{\mathrm{\θ}}\right)\right|}^2}{{\left|M_{\tilde{p}}\left(\widehat{\mathrm{\θ}}\right)\right|}^2}$

Wherein, P _rfor the ratio of correlation energy value; that receive satellite belonging to main synchronizing sequence and Serving cell for time domain with correlation energy value that is sequence; for jammer satellite; that receive main synchronizing sequence and jammer satellite for time domain with correlation energy value that is sequence.

Described step (3) utilizes time domain to receive carrying out frequency deviation estimation formulas with sequence and being of satellite belonging to main synchronizing sequence and Serving cell:

$\widehat{\mathrm{\ϵ}}=\frac{1}{\text{\π}}\∠\left\{{\left[\underset{n=0}{\overset{\frac{N}{2}-1}{\mathrm{\Σ}}}r(n+\widehat{\mathrm{\θ}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\right]}^{*}\right[\underset{n=\frac{N}{2}}{\overset{N-1}{\mathrm{\Σ}}}r(n+\widehat{\mathrm{\θ}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\left]\right\}$

Wherein, for frequency deviation estimated value; Phase angle is asked in ∠ { } expression; for time domain receives main synchronizing sequence; be the local time-domain primary synchronization sequence of satellite is added the sequence obtained.

Described step (3) compensate of frequency deviation formula is:

$\tilde{r}(n+\widehat{\mathrm{\θ}})=r(n+\widehat{\mathrm{\θ}})\exp (-j2\mathrm{\πn}\widehat{\mathrm{\ϵ}}/N)$

Wherein, for the time domain after compensate of frequency deviation receives main synchronizing sequence.

Local time-domain primary synchronization sequence and the time domain through compensate of frequency deviation of satellite belonging to Serving cell receives the computing formula that main synchronizing sequence carries out related operation and is in described step (4):

$X\left\{u\left(i\right)\right\}=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}\tilde{r}(n+\widehat{\mathrm{\θ}})s_{\hat{p},u\left(i\right)}^{*}\left(n\right)$

Wherein, X{u (i) } belonging to Serving cell, the root sequence number of satellite is the result that time domain after the local time-domain primary synchronization sequence of u (i) and compensate of frequency deviation receives main synchronizing sequence related operation; be the root sequence number of satellite is the local time-domain primary synchronization sequence of u (i).

In described step (4), the computing formula of interference correlation is:

$P_I=\frac{\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}\tilde{r}(n+\widehat{\mathrm{\θ}})s_{\tilde{p}}^{*}\left(n\right)}{\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{\left|s_{\tilde{p}}\left(n\right)\right|}^2}\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{s}_{\hat{p},u\left(i\right)}^{*}\left(n\right)s_{\tilde{p}}\left(n\right)$

Wherein, P _ifor interference correlation; for jammer satellite and sequence.

When in described step (4), correlation energy value is greater than threshold value and is less than threshold value, the judgment formula of the main synchronizing sequence that Serving cell configures is respectively:

$\hat{u}=\arg \underset{u\left(i\right)}{\max }\left\{{\left|X\right\{u\left(i\right)\left\}\right|}^2\right\}$

$\hat{u}=\arg \underset{u\left(i\right)}{\max }\left\{{\left|X\right\{u\left(i\right)\}-P_I|}^2\right\}$

Wherein for the estimated value of the main synchronizing sequence root sequence number that Serving cell configures; representing to get makes x be u (i) value of maximum.

Beneficial effect: compared with prior art, universe provided by the invention covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, and tool has the following advantages:

The first, significantly reduce computation complexity: by tradition based on local time-domain primary synchronization sequence cross-correlation test algorithm required for N _uthe local time-domain primary synchronization sequence of × P bar and time domain receiving sequence relevant huge amount of calculation of sliding is reduced to and only needs P bar sequence and time domain receiving sequence to slide the amount of calculation of being correlated with; The second, obtain the satellite information (multi-satellite common networking) belonging to Serving cell; 3rd, larger carrier wave frequency deviation can be resisted; 4th, improve main synchronizing sequence and detect correct probability; 5th, exist under strong inter-satellite signal disturbing, by correlation will be disturbed to offset, can effectively detect main synchronizing sequence.

Accompanying drawing explanation

Fig. 1 is method flow diagram of the present invention;

Fig. 2 is the system model figure that the embodiment of the present invention adopts;

Fig. 3 is sequence s in the embodiment of the present invention _{p, u (i)}(n) and with frequency deviation susceptibility figure;

Fig. 4 is ground based terminal when being positioned at center of housing estate, main synchronizing sequence Interference Cancellation detection method of the present invention and the performance simulation comparison diagram of tradition based on the cross-correlation test method of local time-domain primary synchronization sequence;

Fig. 5 is ground based terminal when being positioned at cell edge, main synchronizing sequence Interference Cancellation detection method of the present invention and the performance simulation comparison diagram of tradition based on the cross-correlation test method of local time-domain primary synchronization sequence.

Embodiment

Below in conjunction with specific embodiment, illustrate the present invention further, these embodiments should be understood only be not used in for illustration of the present invention and limit the scope of the invention, after having read the present invention, the amendment of those skilled in the art to the various equivalent form of value of the present invention has all fallen within the application's claims limited range.

As shown in Figure 1, universe covers the flow process of the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, introduces it below by concrete example.

With P=2 satellite and every satellite configuration N _uarticle=3, different local time-domain primary synchronization sequence: be respectively u={25,29,30} and u={33,34,38}, in corresponding cell set, ID is respectively and and s _{p, u (i)}n the length N=64 of () is example, be described in further detail by reference to the accompanying drawings to the specific embodiment of the present invention, and Fig. 2 is described routine system model figure.

(1) the local time-domain primary synchronization sequence of same satellite be added, every satellite all obtains one and sequence:

${\tilde{s}}_1\left(n\right)=\underset{i=0}{\overset{N_u-1}{\mathrm{\Σ}}}{s}_{1,u\left(i\right)}\left(n\right)=s_{1,u=25}\left(n\right)+s_{1,u=29}\left(n\right)+s_{1,u=30}\left(n\right)$ (formula 1)

${\tilde{s}}_2\left(n\right)=\underset{i=0}{\overset{N_u-1}{\mathrm{\Σ}}}{s}_{2,u\left(i\right)}\left(n\right)=s_{2,u=33}\left(n\right)+s_{2,u=34}\left(n\right)+s_{2,u=38}\left(n\right)$ (formula 2)

Fig. 3 gives sequence s _{p, u (i)}(n) and with frequency deviation susceptibility, the computing formula of the frequency deviation susceptibility FoS of sequence s (n) is:

$\mathrm{FoS}\left\{s\right\}=10l{\mathrm{og}}_{10}\left(\frac{\underset{\underset{\tilde{n}\≠0}{1-N\≤d\≤N-1}}{\max }\left\{{\tilde{R}}_u\left(\tilde{n}\right)\right\}}{{\tilde{R}}_u\left(0\right)}\right)\left(\mathrm{dB}\right)$ (formula 3)

Wherein, ${\tilde{R}}_u\left(\tilde{n}\right)=\left|\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}s\left(n\right)s^{*}\left((n+\tilde{n})\mathrm{mod}N\right)\exp \left(\frac{j2\mathrm{\πn\ϵ}}{N}\right)\right|;$ S is the vector form of sequence s (n); for adding the time domain autocorrelation value of sequence after carrier wave frequency deviation; ε is normalization frequency deviation.

(2) slide relevant with sequence to time domain receiving sequence respectively by every satellite, judge peak, obtain satellite belonging to the original position of main synchronizing sequence and Serving cell information and calculate time domain receive main synchronizing sequence respectively with satellite belonging to Serving cell with sequence and jammer satellite with the ratio of the correlation energy value of sequence:

$M_p\left(\mathrm{\θ}\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}r(n+\mathrm{\θ}){\tilde{s}}_p^{*}\left(n\right),p=\mathrm{1,2}$ (formula 4)

$\{\widehat{\mathrm{\θ}},\hat{p}\}=\arg \underset{p\∈\left\{\mathrm{1,2}\right\},\mathrm{\θ}}{\max }\left\{{\left|M_p\left(\mathrm{\θ}\right)\right|}^2\right\}$ (formula 5)

$P_r=\frac{{\left|M_{\hat{p}}\left(\widehat{\mathrm{\θ}}\right)\right|}^2}{{\left|M_{\tilde{p}}\left(\widehat{\mathrm{\θ}}\right)\right|}^2}$ (formula 6)

(3) sequence utilizing time domain to receive satellite addition acquisition belonging to main synchronizing sequence and Serving cell is carried out frequency deviation and is estimated and compensate:

$\widehat{\mathrm{\ϵ}}=\frac{1}{\text{\π}}\∠\left\{{\left[\underset{n=0}{\overset{\frac{N}{2}-1}{\mathrm{\Σ}}}r(n+\widehat{\mathrm{\θ}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\right]}^{*}\right[\underset{n=\frac{N}{2}}{\overset{N-1}{\mathrm{\Σ}}}r(n+\widehat{\mathrm{\θ}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\left]\right\}$ (formula 7)

$\tilde{r}(n+\widehat{\mathrm{\θ}})=r(n+\widehat{\mathrm{\θ}})\exp (-j2\mathrm{\πn}\widehat{\mathrm{\ϵ}}/N)$ (formula 8)

(4) ratio according to correlation energy value carries out threshold judgement, the detection method of the main synchronizing sequence selecting Serving cell to configure: if be greater than threshold value, then respectively the local time-domain primary synchronization sequence of satellite belonging to Serving cell is received main synchronizing sequence with the time domain through compensate of frequency deviation and carry out related operation, if be less than threshold value, then the correlation that the local time-domain primary synchronization sequence of satellite belonging to Serving cell and the time domain through compensate of frequency deviation receive main synchronizing sequence is deducted interference correlation, judge correlation peak, obtain the main synchronizing sequence that Serving cell configures:

$X\left\{u\left(i\right)\right\}=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}\tilde{r}(n+\widehat{\mathrm{\θ}})s_{\hat{p},u\left(i\right)}^{*}\left(n\right)$ (formula 9)

$\hat{u}=\arg \underset{u\left(i\right)}{\max }\left\{{\left|X\right\{u\left(i\right)\left\}\right|}^2\right\}$ (formula 10)

$\hat{u}=\arg \underset{u\left(i\right)}{\max }\left\{{\left|X\right\{u\left(i\right)\}-P_I|}^2\right\}$ (formula 11)

Fig. 4 and Fig. 5 is respectively terminal and is positioned at center of housing estate and cell edge place, main synchronizing sequence Interference Cancellation detection method of the present invention and the performance simulation comparison diagram of tradition based on the cross-correlation test method of local time-domain primary synchronization sequence.Concrete simulation parameter in table 1 to table 3.Simulation result shows, the performance of main synchronizing sequence Interference Cancellation detection method of the present invention is better than the cross-correlation test method of tradition based on local time-domain primary synchronization sequence, special under large frequency deviation, the performance gain of main synchronizing sequence Interference Cancellation detection method of the present invention is larger.

Table 1 simulation parameter

Table 2 community relative power is arranged

Table 3 CIC decimation filter design parameter

Sequence number	Parameter	Value
			1.	Extract factor R	32
2.	Delay parameter D	1
			3.	Filter cascade exponent number N	2

Claims (5)

1. universe covers a main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, it is characterized in that, comprises the steps:

(1) to time domain receiving sequence low-pass filtering and down-sampling;

(2) what the main synchronizing sequence addition belonging to same satellite obtained respectively slides relevant with sequence to time domain receiving sequence, judge the peak of correlated results, obtain the information that time domain receives satellite belonging to the original position of main synchronizing sequence and Serving cell, calculate time domain receive main synchronizing sequence respectively with satellite belonging to Serving cell with sequence and jammer satellite with the ratio of the correlation energy value of sequence;

Relevant operational formula of sliding is $M_p\left(\mathrm{\θ}\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}r(n+\mathrm{\θ}){\tilde{s}}_p^{*}\left(n\right),p=\mathrm{1,2};$

Belonging to the original position of main synchronizing sequence and Serving cell, the judgment formula of the information of satellite is: $\{\widehat{\mathrm{\θ}},\hat{p}\}=\arg \underset{p\∈\left\{\mathrm{1,2}\right\},\mathrm{\θ}}{\max }\left\{{\left|M_p\left(\mathrm{\θ}\right)\right|}^2\right\};$

Time domain receive main synchronizing sequence respectively with satellite belonging to Serving cell with sequence and jammer satellite with the computing formula of the ratio of the correlation energy value of sequence be:

Wherein, M _p(θ) be the result of related operation, θ is the original position of sliding window; R (n+ θ) is for original position is in the time domain receiving sequence of θ; () ^*represent conjugation; with be respectively the estimated value of satellite belonging to the original position of main synchronizing sequence and Serving cell; | x| ²represent the energy asking x; representing to get makes x be the p of maximum, θ value; P _rfor the ratio of correlation energy value; that receive satellite belonging to main synchronizing sequence and Serving cell for time domain with correlation energy value that is sequence; for co-channel interference satellite; that receive main synchronizing sequence and jammer satellite for time domain with correlation energy value that is sequence;

(3) time domain is utilized to receive the carrying out frequency deviation estimation with sequence and compensating of satellite belonging to main synchronizing sequence and Serving cell;

Time domain is utilized to receive carrying out frequency deviation estimation formulas with sequence and being of satellite belonging to main synchronizing sequence and Serving cell:

$\widehat{\mathrm{\ϵ}}=\frac{1}{\mathrm{\π}}\∠\left\{{\left[\underset{n=0}{\overset{\frac{N}{2}-1}{\mathrm{\Σ}}}r(n+\widehat{\mathrm{\θ}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\right]}^{*}\right[\underset{n=\frac{N}{2}}{\overset{N-1}{\mathrm{\Σ}}}r(n+\widehat{\mathrm{\θ}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\left]\right\}$

Wherein, for frequency deviation estimated value; Phase angle is asked in ∠ { } expression; for time domain receives main synchronizing sequence; be satellite and sequence;

Compensate of frequency deviation formula is:

$\tilde{r}(n+\widehat{\mathrm{\θ}})=r(n+\widehat{\mathrm{\θ}})\exp (-j2\mathrm{\πn}\widehat{\mathrm{\ϵ}}/N)$

Wherein, for the time domain after compensate of frequency deviation receives main synchronizing sequence;

(4) ratio according to correlation energy value carries out threshold judgement, the detection method of the main synchronizing sequence selecting Serving cell to configure: if be greater than threshold value, then respectively the local time-domain primary synchronization sequence of satellite belonging to Serving cell is received main synchronizing sequence with the time domain through compensate of frequency deviation and carry out related operation, if be less than threshold value, then the correlation that the local time-domain primary synchronization sequence of satellite belonging to Serving cell and the time domain through compensate of frequency deviation receive main synchronizing sequence is deducted interference correlation, judge correlation peak, obtain the main synchronizing sequence that Serving cell configures.

2. universe according to claim 1 covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, it is characterized in that: in described step (2), main synchronizing sequence sum operation formula is:

${\tilde{s}}_p\left(n\right)=\underset{i=0}{\overset{N_u-1}{\mathrm{\Σ}}}{s}_{p,u\left(i\right)}\left(n\right)$

Wherein, be p satellite and sequence, p=0,1 ..., P-1; s _{p, u (i)}the local time-domain primary synchronization sequence of (n) to be the root sequence number of p satellite be u (i), i=0,1 ..., N _u-1; N is main synchronizing sequence sampled point sequence number, n=0,1 ..., N-1.

3. universe according to claim 2 covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, it is characterized in that: in described step (4), local time-domain primary synchronization sequence and the time domain through compensate of frequency deviation of satellite belonging to Serving cell receives the computing formula that main synchronizing sequence carries out related operation and be:

$X\left\{u\left(i\right)\right\}=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}\tilde{r}(n+\widehat{\mathrm{\θ}})s_{\hat{p},u\left(i\right)}^{*}\left(n\right)$

Wherein, X{u (i) } belonging to Serving cell, the root sequence number of satellite is the result that time domain after the local time-domain primary synchronization sequence of u (i) and compensate of frequency deviation receives main synchronizing sequence related operation; be the root sequence number of satellite is the local time-domain primary synchronization sequence of u (i).

4. universe according to claim 3 covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, it is characterized in that: in described step (4), the computing formula of interference correlation is:

$P_I=\frac{\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}\tilde{r}(n+\widehat{\mathrm{\θ}})s_{\tilde{p}}^{*}\left(n\right)}{\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{\left|s_{\tilde{p}}\left(n\right)\right|}^2}\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{s}_{\hat{p},u\left(i\right)}^{*}\left(n\right)s_{\tilde{p}}\left(n\right)$

Wherein, P _ifor interference correlation; for jammer satellite and sequence.

5. universe according to claim 4 covers the main synchronizing sequence Interference Cancellation detection method of multi-beam S-LTE, it is characterized in that: when in described step (4), correlation energy value is greater than threshold value and is less than threshold value, the judgment formula of the main synchronizing sequence that Serving cell configures is respectively:

$\hat{u}=\arg \underset{u\left(i\right)}{\max }\left\{{\left|X\right\{u\left(i\right)\left\}\right|}^2\right\}$

$\hat{u}=\arg \underset{u\left(i\right)}{\max }\left\{{\left|X\right\{u\left(i\right)\}-P_I|}^2\right\}$

Wherein for the estimated value of the main synchronizing sequence root sequence number that Serving cell configures; representing to get makes x be u (i) value of maximum.