CN103916348A - Calculation methods and systems for phase deviant, timing deviation and frequency deviation - Google Patents

Abstract

The embodiment of the invention relates to the technical field of communication and discloses a calculation method and system for a phase deviant, a calculation method and system for a timing deviant and a calculation method and system for a frequency deviant. The calculation method for the phase deviant includes: receiving time-domain signals r<t> transmitted through a channel; performing noise cancelling on the time-domain signals r<t>; performing rapid fourier transformation on the received time-domain signals r<t> so as to obtain frequency-domain signals r<f>; selecting at least one signal r<M>, which adopts a multiple phase shift keying (MPSK) modulation method, in the frequency-domain signals r<f> and determining an M value of the value of a constant A of the MPSK modulation method adopted by the frequency-domain signals r<M>; and performing M times of self-correlation on the frequency-domain signals r<M> so that the phase deviant can be obtained. Through the calculation methods and systems for the phase deviant, the timing deviation and the frequency deviation, a large quantity of local reference signals are not needed to be saved, calculation can be simplified and the calculation speed can be improved.

Description

The computational methods of phase pushing figure, timing offset, frequency departure and system

Technical field

The present invention relates to communication technical field, be specifically related to a kind of computational methods and system of phase pushing figure, a kind of computational methods of timing offset value and system, a kind of computational methods of exemplary frequency deviation values and system.

Background technology

Long Term Evolution (being called for short LTE) technology has two kinds of frame structure types.(be called for short, in 3GPP TS 36.211 (R81), specified two kinds of frame structure types of LTE system at the 3rd generation evolved general land wireless access of partner program technical specification group wireless access network (E-UTRA) (version 8).Frame structure Class1 is applicable to Frequency Division Duplexing (FDD) (the being called for short FDD) pattern of full and half duplex.A radio frames comprises 20 time slots, and sequence number is O to 19, and a subframe definition is two continuous slots, i.e. subframes _ccomprise time slot _2cwith _2c+1.Frame structure type 2 is applicable to time division duplex (being called for short TDD) pattern.Each wireless frame length is made up of two fields, and each field is made up of five subframes.Each subframe _cby 2 time slots _2cwith _2c+1represent, the signal transmission in a time slot can be described with resource lattice, and its size is downlink bandwidth configuration (abbreviation ) be multiplied by resource block size on frequency domain and (be called for short ) the OFDM symbolic number that comprises in individual subcarrier and a descending time slot (is called for short ) individual OFDM (being called for short OFDM) symbol.Each unit in resource lattice is called Resource Unit, by index, (k, l) is carried out to unique identification, wherein with be illustrated respectively in the sequence number of frequency-domain and time-domain.Resource Block is for describing the mapping of physical channel to Resource Unit, and Resource Block is defined as in time domain on individual continuous OFDM symbol and frequency domain individual continuous subcarrier in frequency domain, in its involved and resource lattice, downlink resource lattice as shown in Figure 1.In the down link of LTE system, in the time adopting regular circulation prefix and subcarrier spacing to be 15KHz, in a time slot, comprise 7 OFDM symbols, 12 Resource Block; In the time adopting extended cyclic prefix and subcarrier spacing to be 15KHz, in a time slot, comprise 6 OFDM symbols, 12 Resource Block, in the time that subcarrier spacing is 7.5KHz, in a time slot, comprise 3 OFDM symbols, 24 Resource Block wherein comprise multiple Resource Units on each OFDM symbol, a signal of each Resource Unit carrying.

Digital Modulation be the frequency spectrum shift of digital baseband signal to high frequency treatment, form the bandpass signal that is adapted at transmitting in channel.Basic digital modulation mode has amplitude-shift keying (being called for short ASK), frequency shift keying (being called for short FSK), phase shift keying (being called for short PSK).MPSK (Multiple Phase Shift Keying, multi-system phase shift keying) is the one in phase shift keying, is to utilize multiple outs of phase (or phase difference) of carrier wave to represent the modulation system of digital information.MPSK modulation is a kind of permanent mould modulation, and after carrying out MPSK modulation for a series of data, different modulation meets corresponding to different phase constellations, and the MPSK modulation signal on planisphere can be expressed as: wherein A is an integer of being determined by modulation system, for example, at the binary phase shift keying in LTE system (BPSK) when modulation A=1 by MPSK modulated applications, and A=1 etc. in quarternary phase-shift keying (QPSK) (QPSK) modulation; Q is the number that is subject to phase modulation, q=0, and 1 ..., M-1; M gets 2 positive integer time power.

A signal after MPSK modulation is signal after modulation through the reception signal after channel is wherein, argument θ there are differences and the impact such as Doppler frequency shift produces due to the local oscillator frequencies of time delay and transmitter and receiver.Wherein receive the e in signal R (n) ^{j θ}it is exactly phase pushing figure of the present invention.

A kind of method of calculating phase pushing figure is provided in prior art, has obtained by cell special reference, specifically comprised the steps:

Step 1, receives through time-domain signal after transmission at receiving terminal

Step 2, carries out denoising elimination to received signal;

Step 3, by the time-domain signal r receiving _ttransform to frequency domain r by FFT _f;

Step 4, by slightly synchronously obtaining in cell set and identify respectively with cell ID group

Step 5, by obtain community No. ID by No. ID of community obtain the particular location of cell special reference, thereby obtain cell special reference r _f(n);

Step 6, then by the cell special reference r of frequency domain _fand local frequency domain cell special reference (n) carry out conjugate multiplication, according to the result Y of conjugate multiplication (n)=r _f(n) s _{f this locality} ^*(n)=A ²e ^{j θ}obtain phase pushing figure

If said method in employing prior art, needs to preserve all local reference signals.Taking the tdd mode in LTE system as example, in regular circulation prefix, in a Resource Block, there are 4 cell special reference.A subframe comprises 2 time slots, again because the configuration of maximum downlink bandwidth the cell special reference number of a subframe is 110*4*2=880 so.In the uplink-downlink configuration of the LTE system of tdd mode, the rarest 2 subframes for descending, have at most 8 subframes for descending, minimumly in a radio frames so will produce 880*2=1760 local reference signal s _{f this locality}(n), produce at most 880*8=7040 local reference signal s _{f this locality}, but have and comprise multiple radio frames, the local reference signal s that produces so and store in the signal receiving (n) _{f this locality}(n) can be more, the memory consumption of this algorithm is very large, produces local reference signal s _{f this locality}(n) calculate and add the computings such as conjugate multiplication and cause computational complexity very high, computational speed is very slow.

Summary of the invention

Embodiment of the present invention technical problem to be solved is to provide a kind of computational methods and system of phase pushing figure, a kind of computational methods of timing offset value and system, a kind of computational methods of exemplary frequency deviation values and system, for simplifying the numerical procedure of phase pushing figure, timing offset value and exemplary frequency deviation values, improve computational speed.

The embodiment of the present invention provides a kind of computational methods of phase pushing figure, comprising:

Receive the time-domain signal r through transmission _t;

To described time-domain signal r _tcarry out noise elimination;

To the time-domain signal r receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

Select described frequency-region signal r _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _m, and determine described frequency-region signal r _fthe M value of the multi-system phase shift keying MPSK modulation system adopting and the value of constant A;

To described frequency-region signal r _mcarrying out M auto-correlation obtains ${\left(r_M\right)}^M=A^M{\left(e^{j\frac{2\mathrm{\&pi;}}{M}n}\right)}^M{\left(e^{\mathrm{j\&theta;}}\right)}^M=A^M{\left(e^{\mathrm{j\&theta;}}\right)}^M,$ Can obtain phase pushing figure

Accordingly, the embodiment of the present invention also provides a kind of computational methods of timing offset value, comprising:

Adopt the computational methods of phase pushing figure described above to calculate respectively the phase pushing figure e of two signals on same OFDM symbol ^{j θ 1}and e ^{j θ 2}; The alternate position spike of described two signals on frequency domain is Δ n;

By described phase pushing figure e ^{j θ 2}and e ^{j θ 2}carry out conjugate multiplication and obtain z (u):

$z\left(u\right)=e^{\mathrm{j\&theta;}1}{\left(e^{\mathrm{j\&theta;}2}\right)}^{*}=e^{j(\mathrm{\&theta;}1-\mathrm{\&theta;}2)}=e^{j\frac{2\mathrm{\&pi;d}\&times;\mathrm{\&Delta;n}}{N}};$

Calculate timing offset estimated value d according to z (u):

Wherein u represents described e ^{j θ 1}the position of respective signal in frequency domain, represent downlink bandwidth configuration, N represents the included number of subcarriers of described OFDM symbol.

Accordingly, the embodiment of the present invention also provides a kind of computational methods of exemplary frequency deviation values, comprising:

Adopt phase pushing figure described above computational methods calculate respectively the phase pushing figure e on different OFDM symbols in same time slot with two signals of identical initial phase ^{j θ 3}and e ^{j θ 4}; The alternate position spike of described two signals in time domain is m2-m1;

By described phase pushing figure e ^{j θ 3}and e ^{j θ 4}carry out conjugate multiplication and obtain z (m):

$z\left(m\right)=e^{\mathrm{j\&theta;}3}{\left(e^{\mathrm{j\&theta;}4}\right)}^{*}=e^{j(\mathrm{\&theta;}3-\mathrm{\&theta;}4)}=e^{j2\mathrm{\&pi;\&Delta;f}(m2-m1)(N+N_{\mathrm{CP}})T_s};$

Calculate frequency offset estimation value Δ f according to z (m):

$\mathrm{\&Delta;f}=\frac{1}{2\mathrm{\&pi;}(m2-m1)(N+N_{\mathrm{CP}})T_s}\arg \left(\underset{m=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(m\right)\right);$

Wherein N represents the included number of subcarriers of described OFDM symbol, N _cPrepresent circulating prefix-length, T _srepresent the sampling interval; Wherein m represents described e ^{j θ 3}the position of respective signal in frequency domain, $m=\mathrm{0,1},...,12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1.$

Accordingly, the embodiment of the present invention also provides a kind of computing system of phase pushing figure, comprising:

Receiving element, for receiving the time-domain signal r through transmission _t;

Noise is eliminated unit, for to described time-domain signal r _tcarry out noise elimination;

Converter unit, for the time-domain signal r to receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

Selected cell, for selecting described frequency-region signal r _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _m, and determine described frequency-region signal r _mthe M value of the multi-system phase shift keying MPSK modulation system adopting and the value of constant A;

The first computing unit, for to described frequency-region signal r _mcarrying out M auto-correlation obtains ${\left(r_M\right)}^M=A^M{\left(e^{j\frac{2\mathrm{\&pi;}}{M}n}\right)}^M{\left(e^{\mathrm{j\&theta;}}\right)}^M=A^M{\left(e^{\mathrm{j\&theta;}}\right)}^M,$ Thereby the phase pushing figure of obtaining

Accordingly, the embodiment of the present invention also provides a kind of computing system of timing offset value, comprising:

The computing system of phase pushing figure as above, for calculating respectively the phase pushing figure e of two signals on same OFDM symbol ^{j θ 1}and e ^{j θ 2}; The alternate position spike of described two signals on frequency domain is Δ n;

First unit that multiplies each other, for by described phase pushing figure e ^{j θ 2}and e ^{j θ 2}carry out conjugate multiplication and obtain z (u):

$z\left(u\right)=e^{\mathrm{j\&theta;}1}{\left(e^{\mathrm{j\&theta;}2}\right)}^{*}=e^{j(\mathrm{\&theta;}1-\mathrm{\&theta;}2)}=e^{j\frac{2\mathrm{\&pi;d}\&times;\mathrm{\&Delta;n}}{N}};$

The second computing unit, calculates timing offset estimated value d according to z (u): $d=\frac{N}{2\mathrm{\&pi;}\&times;\mathrm{\&Delta;n}}\arg \left(\underset{u=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(u\right)\right);$

Wherein u represents described e ^{j θ 1}the position of respective signal in frequency domain, represent downlink bandwidth configuration, N represents the included number of subcarriers of described OFDM symbol.

Accordingly, the embodiment of the present invention also provides a kind of computing system of exemplary frequency deviation values, comprising:

The computing system of phase pushing figure as above, for calculating respectively the phase pushing figure e on the different OFDM symbols of same time slot with two signals of identical initial phase ^{j θ 3}and e ^{j θ 4}; The alternate position spike of described two signals in time domain is m2-m1;

Second unit that multiplies each other, for by described phase pushing figure e ^{j θ 3}and e ^{j θ 4}carry out conjugate multiplication and obtain z (m):

$z\left(m\right)=e^{\mathrm{j\&theta;}3}{\left(e^{\mathrm{j\&theta;}4}\right)}^{*}=e^{j(\mathrm{\&theta;}3-\mathrm{\&theta;}4)}=e^{j2\mathrm{\&pi;\&Delta;f}(m2-m1)(N+N_{\mathrm{CP}})T_s};$

The 3rd computing unit, for calculate frequency offset estimation value Δ f according to z (m):

$\mathrm{\&Delta;f}=\frac{1}{2\mathrm{\&pi;}(m2-m1)(N+N_{\mathrm{CP}})T_s}\arg \left(\underset{m=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(m\right)\right);$

Wherein N represents the included number of subcarriers of described OFDM symbol, N _cPrepresent circulating prefix-length, T _srepresent the sampling interval; Wherein m represents described e ^{j θ 3}the position of respective signal in frequency domain, $m=\mathrm{0,1},...,12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1.$

The computational methods of the phase pushing figure providing in the embodiment of the present invention and system, by time-domain signal is carried out to auto-correlation, can calculate rapidly phase pushing figure, technical scheme compared to existing technology, without preserving a large amount of local reference signals, also simplified account form, improved computational speed simultaneously.

Brief description of the drawings

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, to the accompanying drawing of required use in embodiment be briefly described below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, do not paying under the prerequisite of creative work, can also obtain according to these accompanying drawings other accompanying drawing.

Fig. 1 is the structural representation of downlink resource lattice;

Figure 1A is the schematic flow sheet of the computational methods of the phase pushing figure that provides of the embodiment of the present invention one;

Fig. 2 is the schematic flow sheet of the computational methods of the phase pushing figure that provides of the embodiment of the present invention two;

Fig. 3 is the schematic flow sheet of the computational methods of the timing offset value that provides of the embodiment of the present invention three;

Fig. 4 is the schematic flow sheet of the computational methods of the exemplary frequency deviation values that provides of the embodiment of the present invention four;

Fig. 5 is the structural representation of the computing system of the phase pushing figure that provides of the embodiment of the present invention five;

Fig. 6 is the structural representation of the computing system of the phase pushing figure that provides of the embodiment of the present invention six;

Fig. 7 is the structural representation of the computing system of the timing offset value that provides of the embodiment of the present invention seven;

Fig. 8 is the structural representation of the computational methods of the exemplary frequency deviation values that provides of the embodiment of the present invention eight.

Embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiment.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtaining under creative work prerequisite, belong to the scope of protection of the invention.

A kind of computational methods and system of phase pushing figure are provided in the embodiment of the present invention, a kind of computational methods of timing offset value and system, a kind of computational methods of exemplary frequency deviation values and system, for simplifying the calculating of phase pushing figure, and then further simplify the calculating of timing offset value and exemplary frequency deviation values.Below be elaborated respectively.

Embodiment mono-:

The invention provides a kind of computational methods of phase pushing figure, as shown in Figure 1A, comprising:

101, receive the time-domain signal r through transmission _t;

102, to above-mentioned time-domain signal r _tcarry out noise elimination;

The concrete grammar that in the present embodiment, noise is eliminated can adopt method of the prior art to process, as adopt document: clear Pan in river, Chang Taihua etc. " Performance Ratio of adaptive noise elimination algorithm comparatively and emulation ". science and technology and engineering .2009,9.9 (19): in 5838-5839, disclosed method is carried out denoising elimination, certainly those skilled in the art also can select other suitable noise cancellation methods to process, and do not do concrete restriction at this;

103, to the time-domain signal r receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

104, select above-mentioned frequency-region signal r _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _m, and determine above-mentioned frequency-region signal r _mthe M value of the multi-system phase shift keying MPSK modulation system adopting;

105, to above-mentioned frequency-region signal r _mcarrying out M auto-correlation obtains ${\left(r_M\right)}^M=A^M{\left(e^{j\frac{2\mathrm{\&pi;}}{M}n}\right)}^M{\left(e^{\mathrm{j\&theta;}}\right)}^M=A^M{\left(e^{\mathrm{j\&theta;}}\right)}^M,$ Thereby obtain phase pushing figure e ^{j θ}:

Further, if employing said method calculates the phase pushing figure of multiple signals, can be to multiple phase pushing figure averagings, to improve the precision of the phase pushing figure calculating.

Embodiment bis-:

The invention provides a kind of computational methods of phase pushing figure, as shown in Figure 2, comprising:

201, receive the time-domain signal r through transmission _t;

202, to above-mentioned time-domain signal r _tcarry out noise elimination;

The concrete grammar that in the present embodiment, noise is eliminated can adopt method of the prior art to process, as adopt document: clear Pan in river, Chang Taihua etc. " Performance Ratio of adaptive noise elimination algorithm comparatively and emulation ". science and technology and engineering .2009,9.9 (19): in 5838-5839, disclosed method is carried out denoising elimination, certainly those skilled in the art also can select other suitable noise cancellation methods to process, and do not do concrete restriction at this;

203, to the time-domain signal r receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

204, judge above-mentioned frequency-region signal r _fin one of them signal be whether the signal r that adopts MPSK modulation system _m, if so, execution step 205, otherwise continue frequency-region signal r _fin next signal judge;

In the present embodiment, judge frequency-region signal r _fin any one signal through the concrete grammar of MPSK modulation can be whether: first calculate the normalization Cyclic Spectrum density of this signal and extract position and the peak value at spectrum peak, if peak value number is greater than 1, for MFSK modulation, if peak value number equals 1, calculate the standard deviation of instantaneous amplitude _aif, σ _a>0.1 is for MASK modulation, if σ _a<=0.1 is MPSK modulation, specifically can list of references: Xu Bin, and Lei Qing etc. " a kind of digital signal modulation mode recognition methods ". the communication technology, 2011,11 (44): 23-24; Certainly those skilled in the art also can select other suitable determination methods to judge, do not do concrete restriction at this;

205, select the signal r of above-mentioned employing MPSK modulation system _m, and determine above-mentioned frequency-region signal r _mthe M value of the multi-system phase shift keying MPSK modulation system adopting and the value of constant A;

In step 204, judge frequency-region signal r _mfor MPSK modulation (σ _a<=0.1) time, can be determined by a single spectral line of fixed position in spectrogram the value of M, because the value of A is determined by modulation system, be determined by M, after having obtained M value, can obtain A value; In like manner, those skilled in the art also can select other suitable methods to determine M value and A value, do not do concrete restriction at this;

206, to above-mentioned frequency-region signal r _mcarrying out M auto-correlation obtains ${\left(r_M\right)}^M=A^M{\left(e^{j\frac{2\mathrm{\&pi;}}{M}n}\right)}^M{\left(e^{\mathrm{j\&theta;}}\right)}^M=A^M{\left(e^{\mathrm{j\&theta;}}\right)}^M,$ Can obtain phase pushing figure e ^{j θ}:

Further, if employing said method calculates the phase pushing figure of multiple signals, can be to multiple phase pushing figure averagings, to improve the precision of the phase pushing figure calculating.

Embodiment tri-:

The invention provides a kind of computational methods of phase pushing figure, as shown in Figure 3, comprising:

301, receive the time-domain signal r through transmission _t;

302, to above-mentioned time-domain signal r _tcarry out noise elimination;

The concrete grammar that in the present embodiment, noise is eliminated can adopt method of the prior art to process, as adopt clear Pan in document river, Chang Taihua etc. " Performance Ratio of adaptive noise elimination algorithm comparatively and emulation ". science and technology and engineering .2009,9.9 (19): in 5838-5839, disclosed method is carried out denoising elimination, certainly those skilled in the art also can select other suitable noise cancellation methods to process, and do not do concrete restriction at this;

303, to the time-domain signal r receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

304, undertaken slightly synchronously obtaining in cell set and identifying respectively by master sync signal PSS and auxiliary synchronous signals SSS with cell ID group

305, by obtain community No. ID by No. ID of community obtain the particular location of frequency-region signal small area DRS (Dedicated Reference Signal), thereby select to adopt the cell special reference of quarternary phase-shift keying (QPSK) modulation system;

Cell special reference, owing to adopting QPSK modulation, can determine that its M value is 4, and A=1;

306, to above-mentioned cell special reference r _mcarrying out 4 auto-correlations obtains ${\left(r_M\right)}^4=A^4{\left(e^{j\frac{2\mathrm{\&pi;}}{4}n}\right)}^4{\left(e^{\mathrm{j\&theta;}}\right)}^4=A^4{\left(e^{\mathrm{j\&theta;}}\right)}^4,$ Thereby obtain phase pushing figure e ^{j θ}:

Further, if employing said method calculates the phase pushing figure of multiple signals, can be to multiple phase pushing figure averagings, to improve the precision of the phase pushing figure calculating.

Embodiment tetra-:

The present invention also provides a kind of computational methods of timing offset value, as shown in Figure 4, comprising:

401, calculate respectively the phase pushing figure e of two signals on same OFDM symbol ^{j θ 1}and e ^{j θ 2}; The alternate position spike on frequency domain of these two signals is Δ n;

In the present embodiment, can adopt the method for describing in embodiment mono-or two or three to calculate phase pushing figure e ^{j θ 1}and e ^{j θ 2}, do not repeat them here;

402, by above-mentioned phase pushing figure e ^{j θ 2}and e ^{j θ 2}carry out conjugate multiplication and obtain z (u):

$z\left(u\right)=e^{\mathrm{j\&theta;}1}{\left(e^{\mathrm{j\&theta;}2}\right)}^{*}=e^{j(\mathrm{\&theta;}1-\mathrm{\&theta;}2)}=e^{j\frac{2\mathrm{\&pi;d}\&times;\mathrm{\&Delta;n}}{N}};$

Wherein u represents e described above ^{j θ 1}the position of respective signal in frequency domain, represent downlink bandwidth configuration, N represents the included number of subcarriers of above-mentioned OFDM symbol;

403, calculate timing offset estimated value d according to z (u): wherein arg is for asking phase angle oeprator.

Embodiment five:

The present invention also provides a kind of computational methods of exemplary frequency deviation values, as shown in Figure 5, comprising:

501, calculate respectively the phase pushing figure e on different OFDM symbols in same time slot with two signals of identical initial phase ^{j θ 3}and e ^{j θ 4}; This two signals position in time domain is respectively m1 and m2, and its alternate position spike is m2-m1;

In the present embodiment, can adopt the method described in embodiment mono-or two or three to calculate phase pushing figure e ^{j θ 3}and e ^{j θ 4}, do not repeat them here;

502, by above-mentioned phase pushing figure e ^{j θ 3}and e ^{j θ 4}carry out conjugate multiplication and obtain z (m):

$z\left(m\right)=e^{\mathrm{j\&theta;}3}{\left(e^{\mathrm{j\&theta;}4}\right)}^{*}=e^{j(\mathrm{\&theta;}3-\mathrm{\&theta;}4)}=e^{j2\mathrm{\&pi;\&Delta;f}(m2-m1)(N+N_{\mathrm{CP}})T_s};$

Wherein N represents the number of subcarriers that OFDM symbol is included, N _cPrepresent circulating prefix-length, T _srepresent the sampling interval; Wherein m represents above-mentioned e ^{j θ 3}the position of respective signal in frequency domain, $m=\mathrm{0,1},...,12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1;$

503, calculate frequency offset estimation value Δ f according to z (m):

$\mathrm{\&Delta;f}=\frac{1}{2\mathrm{\&pi;}(m2-m1)(N+N_{\mathrm{CP}})T_s}\arg \left(\underset{m=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(m\right)\right),$ Wherein arg is for asking phase angle oeprator.

Embodiment six:

The present invention also provides a kind of computing system 60 of phase pushing figure, as shown in Figure 6, comprising:

Receiving element 601, for receiving the time-domain signal r through transmission _t;

Noise is eliminated unit 602, for to above-mentioned time-domain signal r _tcarry out noise elimination;

The concrete grammar that in the present embodiment, noise elimination unit 602 adopts can be: adopt method of the prior art to process, as adopt clear Pan in document river, Chang Taihua etc. " Performance Ratio of adaptive noise elimination algorithm comparatively and emulation ". science and technology and engineering .2009,9.9 (19): in 5838-5839, disclosed method is carried out denoising elimination, certainly those skilled in the art also can select other suitable noise cancellation methods to process, and do not do concrete restriction at this;

Converter unit 603, for the time-domain signal r to receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

Selected cell 604, for selecting above-mentioned frequency-region signal r _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _m, and determine above-mentioned frequency-region signal r _mthe M value of the multi-system phase shift keying MPSK modulation system adopting and the value of constant A;

The first computing unit 605, for to above-mentioned frequency-region signal r _fcarrying out M auto-correlation obtains ${\left(r_f\right)}^M=A^M{\left(e^{j\frac{2\mathrm{\&pi;}}{M}n}\right)}^M{\left(e^{\mathrm{j\&theta;}}\right)}^M,$ Thereby obtain phase pushing figure e ^{j θ}:

For instance, in a kind of execution mode, selected cell 604 can specifically comprise therein:

Judgment sub-unit, for judging above-mentioned frequency-region signal r _fin one of them whether adopt MPSK modulation system;

The first chooser unit, is judged as by above-mentioned judgment sub-unit the frequency-region signal r that adopts MPSK modulation for selecting _m;

First determines subelement, for determining above-mentioned frequency-region signal r _mthe M value of the MPSK modulation system adopting and the value of constant A;

It should be noted that, above-mentioned judgment sub-unit can adopt following determination methods: first calculate the normalization Cyclic Spectrum density of this signal and extract position and the peak value of composing peak, if peak value number is greater than 1, for MFSK modulation, if peak value number equals 1, calculate the standard deviation of instantaneous amplitude _aif, σ _a>0.1 is for MASK modulation, if σ _a<=0.1 is MPSK modulation, specifically can list of references: Xu Bin, and Lei Qing etc. " a kind of digital signal modulation mode recognition methods ". the communication technology, 2011,11 (44): 23-24; Certainly those skilled in the art also can select other suitable determination methods to judge, do not do concrete restriction at this.

It should be noted that, above-mentioned definite subelement determines that the method for M value and A can be: judging frequency-region signal r _mfor MPSK modulation (σ _a<=0.1) time, can be determined by a single spectral line of fixed position in spectrogram the value of M, because the value of A is determined by modulation system, be determined by M, after having obtained M value, can obtain A value; In like manner, those skilled in the art also can select other suitable methods to determine M value and A value, do not do concrete restriction at this.

For instance, in another embodiment, selected cell 604 also can specifically comprise:

Thick synchronous subelement, for being undertaken slightly synchronously obtaining in cell set and identifying respectively by master sync signal PSS and auxiliary synchronous signals SSS with cell ID group

The second chooser unit, for by obtain community No. ID by No. ID of community obtain frequency-region signal r _fthe particular location of small area DRS (Dedicated Reference Signal), thus selection adopts the cell special reference of quarternary phase-shift keying (QPSK) modulation system;

Second determines subelement, is that 4, A value is 1 for determining that above-mentioned cell special reference adopts the M value of MPSK modulation.

It should be noted that, cell special reference is QPSK modulation due to what adopt, and its M value is 4, corresponding, A value is 1.

Further, if adopt said system to calculate the phase pushing figure of multiple signals, this system also can comprise: averaging unit (not shown), and for to multiple phase pushing figure averagings, to improve the precision of the phase pushing figure calculating.

Embodiment seven:

The present invention also provides a kind of computing system of timing offset value, as shown in Figure 7, comprising:

As the computing system 60 of embodiment six described phase pushing figures, for calculating respectively the phase pushing figure e of two signals on same OFDM symbol ^{j θ 1}and e ^{j θ 2}; The alternate position spike on frequency domain of these two signals is Δ n;

The particular content of computing system 60 please refer to embodiment six, does not repeat them here;

First unit 701 that multiplies each other, for by above-mentioned phase pushing figure e ^{j θ 1}and e ^{j θ 2}carry out conjugate multiplication and obtain z (u):

$z\left(u\right)=e^{\mathrm{j\&theta;}1}{\left(e^{\mathrm{j\&theta;}2}\right)}^{*}=e^{j(\mathrm{\&theta;}1-\mathrm{\&theta;}2)}=e^{j\frac{2\mathrm{\&pi;d}\&times;\mathrm{\&Delta;n}}{N}};$

Wherein u represents above-mentioned e ^{j θ 1}the position of respective signal in frequency domain, represent downlink bandwidth configuration, N represents the included number of subcarriers of above-mentioned OFDM symbol;

The second computing unit 702, for calculating timing offset estimated value d according to z (u): $d=\frac{N}{2\mathrm{\&pi;}\&times;\mathrm{\&Delta;n}}\arg \left(\underset{u=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(u\right)\right),$ Wherein arg is for asking phase angle oeprator.

Embodiment eight:

The present invention also provides a kind of computing system of exemplary frequency deviation values, as shown in Figure 8, comprising:

As the computing system 60 of embodiment six described phase pushing figures, for calculating respectively the phase pushing figure e on different OFDM symbols in same time slot with two signals of identical initial phase ^{j θ 3}and e ^{j θ 4}; This two signals position in time domain is respectively m1 and m2, and its alternate position spike is m2-m1;

The particular content of computing system 60 please refer to embodiment six, does not repeat them here;

Second unit 801 that multiplies each other, for by above-mentioned phase pushing figure e ^{j θ 3}and e ^{j θ 4}carry out conjugate multiplication and obtain z (m): $z\left(m\right)=e^{\mathrm{j\&theta;}3}{\left(e^{\mathrm{j\&theta;}4}\right)}^{*}=e^{j(\mathrm{\&theta;}3-\mathrm{\&theta;}4)}=e^{j2\mathrm{\&pi;\&Delta;f}(m2-m1)(N+N_{\mathrm{CP}})T_s};$

Wherein N represents the included number of subcarriers of above-mentioned OFDM symbol, N _cPrepresent circulating prefix-length, T _srepresent the sampling interval; Wherein m represents above-mentioned e ^{j θ 3}the position of respective signal in frequency domain, $m=\mathrm{0,1},...,12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1;$

The 3rd computing unit 803, for calculate frequency offset estimation value Δ f according to z (m):

$\mathrm{\&Delta;f}=\frac{1}{2\mathrm{\&pi;}(m2-m1)(N+N_{\mathrm{CP}})T_s}\arg \left(\underset{m=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(m\right)\right),$ Wherein arg is for asking phase angle oeprator.

One of ordinary skill in the art will appreciate that all or part of step in the whole bag of tricks of above-described embodiment is can carry out the hardware that instruction is relevant by program to complete, this program can be stored in a computer-readable recording medium, storage medium can comprise: flash disk, read-only memory (Read-Only Memory, ROM), random access device (Random Access Memory, RAM), disk or CD etc.

The computational methods of the phase pushing figure above embodiment of the present invention being provided and system, the computational methods of timing offset value and system, computational methods and the system of exemplary frequency deviation values are described in detail, applied specific case herein principle of the present invention and execution mode are set forth, the explanation of above embodiment is just for helping to understand method of the present invention and core concept thereof; , for one of ordinary skill in the art, according to thought of the present invention, all will change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention meanwhile.

Claims (10)

1. computational methods for phase pushing figure, is characterized in that, comprising:

Receive the time-domain signal r through transmission _t;

To described time-domain signal r _tcarry out noise elimination;

To the time-domain signal r receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

Select described frequency-region signal r _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _m, and determine described frequency-region signal r _mthe M value of the MPSK modulation system adopting and the value of constant A;

To described frequency-region signal r _mcarrying out M auto-correlation obtains ${\left(r_M\right)}^M=A^M{\left(e^{j\frac{2\mathrm{\&pi;}}{M}n}\right)}^M{\left(e^{\mathrm{j\&theta;}}\right)}^M=A^M{\left(e^{\mathrm{j\&theta;}}\right)}^M,$ Thereby the phase pushing figure of obtaining

2. method according to claim 1, is characterized in that, the described frequency-region signal r of described selection _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _mcomprise:

Judge described frequency-region signal r _fin one of them signal whether adopt MPSK modulation system, if so, select this signal.

3. method according to claim 1, is characterized in that, the described frequency-region signal r of described selection _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _mcomprise:

Undertaken slightly synchronously obtaining in cell set and identifying respectively by master sync signal PSS and auxiliary synchronous signals SSS with cell ID group

By obtain community No. ID by No. ID of community obtain described frequency-region signal r _fthe particular location of small area DRS (Dedicated Reference Signal), thus selection adopts the cell special reference of quarternary phase-shift keying (QPSK) modulation system.

4. computational methods for timing offset value, is characterized in that, comprising:

The method of employing as described in claims 1 to 3 any one calculates respectively the phase pushing figure e of two signals on same OFDM symbol ^{j θ 1}and e ^{j θ 2}; The alternate position spike of described two signals on frequency domain is Δ n;

By described phase pushing figure e ^{j θ 2}and e ^{j θ 2}carry out conjugate multiplication and obtain z (u):

$z\left(u\right)=e^{\mathrm{j\&theta;}1}{\left(e^{\mathrm{j\&theta;}2}\right)}^{*}=e^{j(\mathrm{\&theta;}1-\mathrm{\&theta;}2)}=e^{j\frac{2\mathrm{\&pi;d}\&times;\mathrm{\&Delta;n}}{N}};$

Calculate timing offset estimated value d according to z (u): $d=\frac{N}{2\mathrm{\&pi;}\&times;\mathrm{\&Delta;n}}\arg \left(\underset{u=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(u\right)\right);$

Wherein u represents described e ^{j θ 1}the position of respective signal in frequency domain, represent downlink bandwidth configuration, N represents the included number of subcarriers of described OFDM symbol.

5. computational methods for exemplary frequency deviation values, is characterized in that, comprising:

Adopt the method as described in claims 1 to 3 any one to calculate respectively the phase pushing figure e on different OFDM symbols with two signals of identical initial phase ^{j θ 3}and e ^{j θ 4}; The alternate position spike of described two signals in time domain is m2-m1;

By described phase pushing figure e ^{j θ 3}and e ^{j θ 4}carry out conjugate multiplication and obtain z (m):

$z\left(m\right)=e^{\mathrm{j\&theta;}3}{\left(e^{\mathrm{j\&theta;}4}\right)}^{*}=e^{j(\mathrm{\&theta;}3-\mathrm{\&theta;}4)}=e^{j2\mathrm{\&pi;\&Delta;f}(m2-m1)(N+N_{\mathrm{CP}})T_s};$

Calculate frequency offset estimation value Δ f according to z (m):

$\mathrm{\&Delta;f}=\frac{1}{2\mathrm{\&pi;}(m2-m1)(N+N_{\mathrm{CP}})T_s}\arg \left(\underset{m=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(m\right)\right);$

Wherein N represents the included number of subcarriers of described OFDM symbol, N _cPrepresent circulating prefix-length, T _srepresent the sampling interval; Wherein m represents described e ^{j θ 3}the position of respective signal in frequency domain, $m=\mathrm{0,1},...,12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1.$

6. a computing system for phase pushing figure, is characterized in that, comprising:

Receiving element, for receiving the time-domain signal r through transmission _t;

Noise is eliminated unit, for to described time-domain signal r _tcarry out noise elimination;

Converter unit, for the time-domain signal r to receiving _tcarry out fast Fourier transform and obtain frequency-region signal r _f;

Selected cell, for selecting described frequency-region signal r _fat least one signal r of middle employing multi-system phase shift keying MPSK modulation system _m, and determine described frequency-region signal r _mthe M value of the multi-system phase shift keying MPSK modulation system adopting and the value of constant A;

The first computing unit, for to described frequency-region signal r _mcarrying out M auto-correlation obtains ${\left(r_M\right)}^M=A^M{\left(e^{j\frac{2\mathrm{\&pi;}}{M}n}\right)}^M{\left(e^{\mathrm{j\&theta;}}\right)}^M=A^M{\left(e^{\mathrm{j\&theta;}}\right)}^M,$ Thereby the phase pushing figure of obtaining

7. system according to claim 6, is characterized in that, described selected cell comprises:

Judgment sub-unit, for judging described frequency-region signal r _fin one of them whether adopt MPSK modulation system;

The first chooser unit, is judged as by described judgment sub-unit the frequency-region signal r that adopts MPSK modulation for selecting _m;

First determines subelement, for determining described frequency-region signal r _mthe M value of the MPSK modulation system adopting and the value of constant A.

8. system according to claim 6, is characterized in that, described selected cell comprises:

Thick synchronous subelement, for being undertaken slightly synchronously obtaining in cell set and identifying respectively by master sync signal PSS and auxiliary synchronous signals SSS with cell ID group

The second chooser unit, for by obtain community No. ID by No. ID of community obtain described frequency-region signal r _fthe particular location of small area DRS (Dedicated Reference Signal), thus selection adopts the cell special reference of quarternary phase-shift keying (QPSK) modulation system;

Second determines subelement, is that 4, A value is 1 for determining that described cell special reference adopts the M value of MPSK modulation.

9. a computing system for timing offset value, is characterized in that, comprising:

System as described in claim 6-8 any one, for calculating respectively the phase pushing figure e of two signals on same OFDM symbol ^{j θ 1}and e ^{j θ 2}; The alternate position spike of described two signals on frequency domain is Δ n;

First unit that multiplies each other, for by described phase pushing figure e ^{j θ 2}and e ^{j θ 2}carry out conjugate multiplication and obtain z (u):

$z\left(u\right)=e^{\mathrm{j\&theta;}1}{\left(e^{\mathrm{j\&theta;}2}\right)}^{*}=e^{j(\mathrm{\&theta;}1-\mathrm{\&theta;}2)}=e^{j\frac{2\mathrm{\&pi;d}\&times;\mathrm{\&Delta;n}}{N}};$

The second computing unit, calculates timing offset estimated value d according to z (u): $d=\frac{N}{2\mathrm{\&pi;}\&times;\mathrm{\&Delta;n}}\arg \left(\underset{u=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(u\right)\right);$

Wherein u represents described e ^{j θ 1}the position of respective signal in frequency domain, represent downlink bandwidth configuration, N represents the included number of subcarriers of described OFDM symbol.

10. a computing system for exemplary frequency deviation values, is characterized in that, comprising:

System as described in claim 6 to 8 any one, for calculating respectively the phase pushing figure e on the different OFDM symbols of same time slot with two signals of identical initial phase ^{j θ 3}and e ^{j θ 4}; The alternate position spike of described two signals in time domain is m2-m1;

Second unit that multiplies each other, for by described phase pushing figure e ^{j θ 3}and e ^{j θ 4}carry out conjugate multiplication and obtain z (m):

$z\left(m\right)=e^{\mathrm{j\&theta;}3}{\left(e^{\mathrm{j\&theta;}4}\right)}^{*}=e^{j(\mathrm{\&theta;}3-\mathrm{\&theta;}4)}=e^{j2\mathrm{\&pi;\&Delta;f}(m2-m1)(N+N_{\mathrm{CP}})T_s};$

The 3rd computing unit, for calculate frequency offset estimation value Δ f according to z (m):

$\mathrm{\&Delta;f}=\frac{1}{2\mathrm{\&pi;}(m2-m1)(N+N_{\mathrm{CP}})T_s}\arg \left(\underset{m=0}{\overset{12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1}{\mathrm{\&Sigma;}}}z\left(m\right)\right);$

Wherein N represents the included number of subcarriers of described OFDM symbol, N _cPrepresent circulating prefix-length, T _srepresent the sampling interval; Wherein m represents described e ^{j θ 3}the position of respective signal in frequency domain, $m=\mathrm{0,1},...,12\&times;N_{\mathrm{RB}}^{\mathrm{DL}}-1.$