CN103546191A - Method for detecting main synchronization sequence interference cancellation of full-domain coverage multi-beam S-LTE - Google Patents

Abstract

The invention discloses a method for detecting main synchronization sequence interference cancellation of full-domain coverage multi-beam S-LTE. The method comprises the steps of conducting low pass filtering and down sampling on a time domain receiving signal, respectively conducting moving correlation on a sum sequence obtained by adding main synchronization sequences belonging to the same satellite, and a time domain receiving sequence, conducting peak detection to obtain the initial position of the time domain receiving main synchronization sequence and type information of the satellite which a serving cell belongs to, respectively computing the correlation energy ratio of the time domain receiving main synchronization sequence to a sum sequence of the satellite which the serving cell belongs to and the correlation energy ratio of the time domain receiving main synchronization sequence to a sum sequence of an interference satellite, utilizing the sum sequences for conducting frequency offset estimation and compensation, conducting threshold judgment on the correlation energy ratio, if the correlation ratio is larger than the threshold value, conducting the correlation operation respectively on a local time domain main synchronization sequence and the time domain receiving main synchronization sequence which is processed through the frequency offset compensation, if the correlation value is smaller than the threshold value, subtracting an interference correlation value from the correlation value, and conducting the peak detection to obtain the main synchronization sequence. According to the method for detecting the main synchronization sequence interference cancellation of the full-domain coverage multi-beam S-LTE, the larger frequency offset is resisted, computation complexity is low, and the 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 to realize one of necessary means communicating in any place.In recent years, the 4th 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 Generation Partnership Project Long Term Evolution of key technology, 3GPP-LTE) put it into commercial operation successively.The land LTE(Terrestrial LTE of the features such as two-forty, large capacity, spectral efficient, high power efficiency will be there is, 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 is focus and the difficult point of current satellite mobile communication area research.

The multi-beam S-LTE mobile communication system of identical networking adopts large-scale antenna array to produce a plurality of wave beams on satellite, and arrival forms a plurality of communities behind ground.The same with T-LTE, neighbor cell configures different main synchronizing sequences.Yet the multi-beam S-LTE of identical networking exists the inter-beam interference even more serious compared with T-LTE (Inter-Beam Interference, IBI), the overlapping scope of neighbor cell is larger.Multi-beam S-LTE adopts multi-satellite (may be dissimilar satellite) common networking, thereby realize universe, covers.In order to reduce the interference adopting between identical main synchronizing sequence allocating cell, and make the terminal can identification satellite, thereby realize better Cell searching, the main synchronizing sequence that different satellite configuration is different, the signal of same satellite arrives same ground receiving terminal through identical fading channel, and the signal of different satellites arrives same ground receiving terminal through different fading channels.

Traditional main synchronizing sequence detects the cross-correlation test algorithm adopting based on the main synchronizing sequence of local time domain, this algorithm has the following disadvantages: (1) this algorithm need to slide relevant to time domain receiving sequence respectively by all main synchronizing sequences of local time domain, 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 be resisted larger frequency deviation, although can improve by the method for piecemeal cross-correlation anti-deviation capability, need to increase computation complexity as cost.

Summary of the invention

Goal of the invention: cover the configuration feature of the main synchronizing sequence of multi-beam S-LTE system and the deficiency of the traditional cross-correlation test algorithm based on the main synchronizing sequence of local time domain for universe, the invention provides and a kind ofly can resist the main synchronizing sequence Interference Cancellation detection method that larger frequency deviation, 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) respectively the main synchronizing sequence of same satellite is added to obtain relevant to the slip of time domain receiving sequence with sequence, judgement correlation peak location, obtain the information that time domain receives the affiliated satellite of original position and Serving cell of main synchronizing sequence, further, calculate time domain receive main synchronizing sequence respectively with Serving cell under satellite and ratio sequence and jammer satellite and correlation energy value sequence;

(3) utilize time domain to receive carrying out frequency deviation estimation and compensate with sequence of the affiliated satellite of main synchronizing sequence and Serving cell;

(4) ratio according to correlation energy value carries out threshold judgement, the detection method of the main synchronizing sequence that selection Serving cell configures: if be greater than threshold value, respectively the time domain of the main synchronizing sequence of local time domain of satellite under Serving cell and process compensate of frequency deviation received to main synchronizing sequence and carry out related operation; If be less than threshold value, the correlation that time domain of the main synchronizing sequence of local time domain of satellite under Serving cell and process compensate of frequency deviation is received to main synchronizing sequence deducts interference correlation, and judgement correlation peak, obtains the main synchronizing sequence that 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)}(n) be that the root sequence number of p satellite is the main synchronizing sequence of local time domain of u (i), i=0,1 ..., N _u-1, N _ufor root sequence number sum, sequence is total, and n is main synchronizing sequence sampled point sequence number, n=0, and 1 ..., N-1.

Described step (2) the 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)$

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

Under the original position of the main synchronizing sequence of described step (2) 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 the affiliated satellite of original position and Serving cell of main synchronizing sequence; | x| ²represent to ask the energy of x;

expression is got and made x is peaked p, θ value.

Time domain receives main synchronizing sequence and satellite with the computing formula of ratio sequence and jammer satellite and correlation energy value sequence is respectively with under Serving cell in described step (2):

$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 _rratio for correlation energy value;

for time domain receives correlation energy value satellite and sequence under main synchronizing sequence and Serving cell;

for jammer satellite;

for time domain receives correlation energy value main synchronizing sequence and jammer satellite and sequence.

Described step (3) utilizes time domain to receive carrying out frequency deviation estimation formulas with sequence and being of satellite under 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; ∠ { } represents to ask phase angle;

for time domain receives main synchronizing sequence; be the main synchronizing sequence of local time domain of satellite is added the sequence obtaining.

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

$\tilde{r}(n+\widehat{\mathrm{\&theta;}})=r(n+\widehat{\mathrm{\&theta;}})\exp (-\mathrm{<figure-callout\; id="2"\; label="j"\; filenames="HDA0000407053050000021.png"\; state="\{\{state\}\}">j</figure-callout>}2\mathrm{\&pi;n}\widehat{\mathrm{\&epsiv;}}/N)$

Wherein,

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

Under Serving cell, the main synchronizing sequence of local time domain of satellite and time domain through compensate of frequency deviation receive the computing formula that main synchronizing sequence carries out related operation and are in described step (4):

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

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

be

the root sequence number of satellite is the main synchronizing sequence of local time domain of u (i).

In described step (4), disturb the computing formula of correlation to be:

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

Wherein, P _ifor disturbing 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 the estimated value of the main synchronizing sequence root sequence number configuring for Serving cell;

expression is got and made x is peaked u (i) value.

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: the needed N of cross-correlation test algorithm by tradition based on the main synchronizing sequence of local time domain _uthe huge amount of calculation that the main synchronizing sequence of the local time domain of * P bar is relevant with the slip of time domain receiving sequence is reduced to only needs P bar sequence and the time domain receiving sequence relevant amount of calculation of sliding; The second, obtain the affiliated satellite information (multi-satellite common networking) in Serving cell; The 3rd, can resist larger carrier wave frequency deviation; The 4th, improve main synchronizing sequence and detected correct probability; The 5th, existing under strong inter-satellite signal interference, by disturbing correlation 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 while being positioned at center of housing estate, the performance simulation comparison diagram of main synchronizing sequence Interference Cancellation detection method of the present invention and the cross-correlation test method of tradition based on the main synchronizing sequence of local time domain;

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

Embodiment

Below in conjunction with specific embodiment, further illustrate the present invention, should understand these embodiment is only not used in and limits the scope of the invention for the present invention is described, after having read the present invention, those skilled in the art all fall within the application's claims limited range to the modification of the various equivalent form of values of the present invention.

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

With P=2 satellite and every satellite configuration N _uarticle=3, the different main synchronizing sequence of local time domain: 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) length N=64 are example, by reference to the accompanying drawings the specific embodiment of the present invention are described in further detail, and Fig. 2 is described routine system model figure.

(1) the main synchronizing sequence of the local time domain of same satellite is added, every satellite all obtains one and sequence:

${\tilde{s}}_1\left(n\right)=\underset{i=0}{\overset{N_u-1}{\mathrm{\&Sigma;}}}{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{\&Sigma;}}}{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="HDA0000407053050000021.png"\; state="\{\{state\}\}">s</figure-callout>}}_{2,u\left(i\right)}\left(n\right)={\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="HDA0000407053050000021.png"\; state="\{\{state\}\}">s</figure-callout>}}_{2,u=33}\left(n\right)+{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="HDA0000407053050000021.png"\; state="\{\{state\}\}">s</figure-callout>}}_{2,u=34}\left(n\right)+{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="HDA0000407053050000021.png"\; state="\{\{state\}\}">s</figure-callout>}}_{2,u=38}\left(n\right)$ (formula 2)

Fig. 3 has provided 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}\&NotEqual;0}{1-N\&le;d\&le;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{\&Sigma;}}}s\left(n\right)s^{*}\left((n+\tilde{n})\mathrm{mod}N\right)\exp \left(\frac{\mathrm{<figure-callout\; id="2"\; label="j"\; filenames="HDA0000407053050000021.png"\; state="\{\{state\}\}">j</figure-callout>}2\mathrm{\&pi;n\&epsiv;}}{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) respectively every satellite slided relevant with sequence and time domain receiving sequence, judgement peak, obtain the information of satellite under the original position of main synchronizing sequence and Serving cell and calculate time domain receive main synchronizing sequence respectively with Serving cell under satellite and ratio sequence and jammer satellite and correlation energy value sequence:

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

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

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

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

$\widehat{\mathrm{\&epsiv;}}=\frac{1}{\text{\&pi;}}\&angle;\left\{{\left[\underset{n=0}{\overset{\frac{N}{2}-1}{\mathrm{\&Sigma;}}}r(n+\widehat{\mathrm{\&theta;}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\right]}^{*}\right[\underset{n=\frac{N}{2}}{\overset{N-1}{\mathrm{\&Sigma;}}}r(n+\widehat{\mathrm{\&theta;}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\left]\right\}$ (formula 7)

$\tilde{r}(n+\widehat{\mathrm{\&theta;}})=r(n+\widehat{\mathrm{\&theta;}})\exp (-\mathrm{<figure-callout\; id="2"\; label="j"\; filenames="HDA0000407053050000021.png"\; state="\{\{state\}\}">j</figure-callout>}2\mathrm{\&pi;n}\widehat{\mathrm{\&epsiv;}}/N)$ (formula 8)

(4) ratio according to correlation energy value carries out threshold judgement, the detection method of the main synchronizing sequence that selection Serving cell configures: if be greater than threshold value, respectively the time domain of the main synchronizing sequence of local time domain of satellite under Serving cell and process compensate of frequency deviation is received to main synchronizing sequence and carry out related operation, if be less than threshold value, the correlation that time domain of the main synchronizing sequence of local time domain of satellite under Serving cell and process compensate of frequency deviation is received to main synchronizing sequence deducts interference correlation, judgement correlation peak, obtains the main synchronizing sequence that Serving cell configures:

$X\left\{u\left(i\right)\right\}=\underset{n=0}{\overset{N-1}{\mathrm{\&Sigma;}}}\tilde{r}(n+\widehat{\mathrm{\&theta;}})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 are respectively terminal and are positioned at center of housing estate and cell edge place, the performance simulation comparison diagram of main synchronizing sequence Interference Cancellation detection method of the present invention and the cross-correlation test method of tradition based on the main synchronizing sequence of local time domain.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 the main synchronizing sequence of local time domain, under large frequency deviation, the performance gain of main synchronizing sequence Interference Cancellation detection method of the present invention is larger especially.

Table 1 simulation parameter

Table 2 community relative power arranges

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 (10)

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) respectively the main synchronizing sequence that belongs to same satellite is added to obtain relevant to the slip of time domain receiving sequence with sequence, the peak of judgement correlated results, obtain the information that time domain receives satellite under the original position of main synchronizing sequence and Serving cell, calculate time domain receive main synchronizing sequence respectively with Serving cell under satellite and ratio sequence and jammer satellite and correlation energy value sequence;

(3) utilize time domain to receive carrying out frequency deviation estimation and compensate with sequence of the affiliated satellite of main synchronizing sequence and Serving cell;

(4) ratio according to correlation energy value carries out threshold judgement, the detection method of the main synchronizing sequence that selection Serving cell configures: if be greater than threshold value, respectively the time domain of the main synchronizing sequence of local time domain of satellite under Serving cell and process compensate of frequency deviation is received to main synchronizing sequence and carry out related operation, if be less than threshold value, the correlation that time domain of the main synchronizing sequence of local time domain of satellite under Serving cell and process compensate of frequency deviation is received to main synchronizing sequence deducts interference correlation, judgement correlation peak, obtains 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{\&Sigma;}}}{s}_{p,u\left(i\right)}\left(n\right)$

Wherein,

be p satellite and sequence, p=0,1 ..., P-1; s _{p, u (i)}(n) be that the root sequence number of p satellite is the main synchronizing sequence of local time domain of u (i), i=0,1 ..., N _u-1; N is main synchronizing sequence sampled point sequence number, n=0, and 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: the relevant operational formula of sliding in described step (2) is:

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

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

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: under the original position of the main synchronizing sequence of described step (2) and Serving cell, the judgment formula of the information of satellite is:

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

Wherein,

with

be respectively the estimated value of the affiliated satellite of original position and Serving cell of main synchronizing sequence; | x| ²represent to ask the energy of x;

expression is got and made x is peaked p, θ value.

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: in described step (2), time domain receives main synchronizing sequence and satellite with the computing formula of ratio sequence and jammer satellite and correlation energy value sequence is respectively with under Serving cell:

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

Wherein, P _rratio for correlation energy value; for time domain receives correlation energy value satellite and sequence under main synchronizing sequence and Serving cell; for co-channel interference satellite;

for time domain receives correlation energy value main synchronizing sequence and jammer satellite and sequence.

6. universe according to claim 5 covers the main synchronizing sequence Interference Cancellation of multi-beam S-LTE detection method, it is characterized in that: in described step (3), utilize time domain to receive carrying out frequency deviation estimation formulas with sequence and being of satellite under main synchronizing sequence and Serving cell:

$\widehat{\mathrm{\&epsiv;}}=\frac{1}{\text{\&pi;}}\&angle;\left\{{\left[\underset{n=0}{\overset{\frac{N}{2}-1}{\mathrm{\&Sigma;}}}r(n+\widehat{\mathrm{\&theta;}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\right]}^{*}\right[\underset{n=\frac{N}{2}}{\overset{N-1}{\mathrm{\&Sigma;}}}r(n+\widehat{\mathrm{\&theta;}}){\tilde{s}}_{\hat{p}}^{*}\left(n\right)\left]\right\}$

Wherein,

for frequency deviation estimated value; ∠ { } represents to ask phase angle;

for time domain receives main synchronizing sequence;

be

satellite and sequence.

7. universe according to claim 6 covers the main synchronizing sequence Interference Cancellation of multi-beam S-LTE detection method, it is characterized in that: described step (3) frequency deviation compensation formula is:

$\tilde{r}(n+\widehat{\mathrm{\&theta;}})=r(n+\widehat{\mathrm{\&theta;}})\exp (-j2\mathrm{\&pi;n}\widehat{\mathrm{\&epsiv;}}/N)$

Wherein,

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

8. universe according to claim 7 covers the main synchronizing sequence Interference Cancellation of multi-beam S-LTE detection method, it is characterized in that: in described step (4), under Serving cell, the main synchronizing sequence of local time domain of satellite and time domain through compensate of frequency deviation receive 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{\&Sigma;}}}\tilde{r}(n+\widehat{\mathrm{\&theta;}})s_{\hat{p},u\left(i\right)}^{*}\left(n\right)$

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

be

the root sequence number of satellite is the main synchronizing sequence of local time domain of u (i).

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

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

Wherein, P _ifor disturbing correlation;

for jammer satellite and sequence.

10. universe according to claim 9 covers the main synchronizing sequence Interference Cancellation of multi-beam S-LTE detection method, 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

the estimated value of the main synchronizing sequence root sequence number configuring for Serving cell;

expression is got and made x is peaked u (i) value.

Images