CN106878205A - A kind of timing offset method of estimation, device and terminal - Google Patents

Abstract

The invention discloses a kind of method of estimation of timing offset, device and terminal, be used to solve adjacent area timing and local area exist it is larger in advance or it is delayed be likely to result in effective diameter exit window in the case of, influence the calculating of effective diameter thresholding and the problem of adjacent area Timing measurement.Method is：CRS according to adjacent cell determines time domain channel estimated result, and the first effective diameter thresholding is determined according to time domain channel estimated result the first noise window power and peak power, and the first timing offset value of current area and adjacent cell is determined according to the first effective diameter thresholding；The timing offset of time domain channel estimated result is adjusted according to the first timing offset value, the second noise window power and peak power according to the time domain channel estimated result after adjustment determine the second effective diameter thresholding；The second timing offset value of current area and adjacent cell is determined according to the second effective diameter thresholding, the timing offset between current area and adjacent cell is determined according to the first timing offset value and the second timing offset value.

Description

A kind of timing offset method of estimation, device and terminal

Technical field

The present invention relates to communication technical field, more particularly to a kind of timing offset method of estimation, device and terminal.

Background technology

Inter-Cell Interference Coordination (Further Enhanced ICIC, the FeICIC) technology for further enhancing needs reply Occur under isomery scene adjacent area public reference signal (Cell-specific reference signals, CRS, it is also referred to as small The specific reference signal in area) and/or Physical Broadcast Channel (Physical Broadcast Channel, PBCH) to this cell Strong jamming, such as macro base station-femto base station (also referred to as macro-pico) deployment scenario or macro base station-base station from home are (also referred to as Macro-femto) deployment scenario, therefore, terminal (UE) needs to carry out interference elimination to CRS and/or PBCH.

Generally, interfered cell reach the time of UE and this cell reach UE time exist it is certain relatively inclined Difference, frequency deviation there is also certain relative deviation, and relative time delay deviation possible range is in [- 3,3] microsecond (us), and relative frequency deviation is inclined Difference possible range is within [- 300,300] hertz (Hz).General UE is in time-frequency convert using only a Fast Fourier Transform (FFT) (Fast Fourier Transformation, FFT) arithmetic element, therefore in order to ensure the demodulation performance of local area, UE can only base Timing windowing position is set in time domain in the reception signal of this cell, i.e., can not possibly tracking the timing of adjacent area, to carry out adjacent area real-time It is synchronous, carry out needing after channel estimation the delay spread for estimating adjacent area to enter line frequency according to local area timing receipt signal and to adjacent area Domain filters.

With Long Term Evolution (Long Term Evolution, LTE) it is descending as a example by, as shown in figure 1, to existing time delay expand Exhibition algorithm for estimating is described as follows：It is S that known descending pilot frequency sends sequence, and the signal for receiving is Y, and spatial frequency domain channel is H+n, wherein, H represents the impulse response of wireless fading channel, and n represents white Gaussian noise, then Y=(H+n) S is met, so that can To estimate that channel estimation in frequency domain is：Channel estimation in frequency domain is transformed into time domain carry out delay spread and estimate Meter, it is specific as follows：Frequency domain channel is transformed into time domain, according to Cyclic Prefix (Cyclic Prefix, CP) length or prior information Determine the position of useful signal time delay maximum diameter, take out the noise in time domain footpath without available signal power and maximum diameter, calculate Noise power P_noiseAnd peak power P_max, determine effective diameter threshold value by setting peak factor and noise factor, so that Effective diameter search space is obtained, delay spread value is determined.

Wherein, noise factor and peak factor can be configured based on experience value, and effective diameter threshold value is expressed as：Γ= (P_noise·β_noise,P_max·β_max), wherein β_noiseNoise gate is represented, is the product of noise power and noise factor, β_maxTable Show peak value threshold, be the product of peak power and peak factor.It is according to the method that CP length sets noise window：CP length is made It is the when elongatedness of time delay maximum diameter in system, calculates the position of noise window.CP is defined in LTE system under different bandwidth The T of protection interval CP_SNumber is N_CP, T_SSampling time interval is represented, FFT points are N_FFT, sub-carrier number is N_SCIf, pilot tone FFT length when frequency domain accordingly transforms to time domain uses N_SC/ 6, then the time-domain position of time delay maximum diameter can be expressed as：

Become in view of using inverse discrete Fourier transform (Inverse Discrete Fourier Transform, IDFT) Frequency-domain transform change commanders to time domain it may happen that the power leakage of signal, can retain N_FFTIn N_leakBar footpath, last effective diameter Scope be [N_FFT-N_leak:N-1,0:N_τ], noise window is [N_τ:N_FFT-N_leak]。

If adjacent area timing and local area exist it is larger in advance or it is delayed, be all likely to result in effective diameter exit window, i.e. noise window In the search space comprising effective diameter and effective diameter it is inaccurate, the calculating of effective diameter thresholding will be directly affected, and then influence adjacent Area's Timing measurement.

The content of the invention

The embodiment of the present invention provides a kind of method of estimation of timing offset, device and terminal, is used to solve in adjacent area regularly With local area exist it is larger in advance or it is delayed, in the case of being likely to result in effective diameter exit window, directly affect the meter of effective diameter thresholding Calculate, and then influence the problem of adjacent area Timing measurement.

Concrete technical scheme provided in an embodiment of the present invention is as follows：

The embodiment of the present invention provides a kind of timing offset method of estimation, including：

The public reference signal CRS of the adjacent cell according to current area carries out channel estimation and determines that time domain channel estimates knot Really, the first noise window power and peak power of the time domain channel estimated result are determined；

First effective diameter thresholding is determined according to the first noise window power and the peak power, according to the time domain Channel estimation results and the first effective diameter thresholding, determine the first timing offset value of current area and the adjacent cell；

The timing offset of the time domain channel estimated result is adjusted according to the first timing offset value, it is determined that adjusting Second noise window power of the time domain channel estimated result after whole, according to the second noise window power and the peak power Determine the second effective diameter thresholding；

According to the time domain channel estimated result after the adjustment and the second effective diameter thresholding, to tap search space Scan for determining the second timing offset value of current area and the adjacent cell, according to the first timing offset value and described Second timing offset value determines the timing offset between current area and the adjacent cell.

In possible implementation method, methods described also includes：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding and the timing offset, Scan for determining the delay spread value of current area and the adjacent cell to tap search space.

In possible implementation method, it is true that the public reference signal CRS of the adjacent cell according to current area carries out channel estimation Domain channel estimation results, including：

CRS according to the adjacent cell carries out least-squares estimation and obtains channel estimation in frequency domain result, and the frequency domain is believed Road estimated result carries out inverse discrete Fourier transform, obtains the time domain channel estimated result.

In possible implementation method, the first noise window power and its peak work of the time domain channel estimated result are determined Rate, including：

Calculate the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, it is expressed as：Wherein,Represent the of port pThe time domain channel of row pilot tone estimates knot Really, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the first noise window power is：Wherein, N_τ Represent the time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the peak power is：

In possible implementation method, the first effective diameter is determined according to the first noise window power and the peak power Thresholding, including：

The product for calculating the first noise window power and the first coefficient obtains the first product, and calculates its peak work Rate obtains the second product with the product of the second coefficient；

Select the maximum in first product and second product, however, it is determined that the maximum is more than or equal to institute Peak power is stated, determines that the first effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak value Power, determines that the first effective diameter thresholding is the maximum.

In possible implementation method, according to the time domain channel estimated result and the first effective diameter thresholding, it is determined that First timing offset value of current area and the adjacent cell, including：

Determine that current area is with the first timing offset value of the adjacent cell： Wherein, P represents total port number, and R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.

In possible implementation method, the timing according to the first timing offset value to the time domain channel estimated result is inclined Difference is adjusted, including：

By when the corresponding time-domain taps of the first timing offset are moved to zero described in the time domain channel estimated result Domain tap, is successively translated each time-domain taps in the time domain channel estimated result, and the time domain channel after being adjusted is estimated Counting result is：Wherein, k=0 ..., N_FFT- 1, N_FFTTable Show that FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel of row pilot tone estimates knot Really.

In possible implementation method, the second effective diameter is determined according to the second noise window power and the peak power Thresholding, including：

The product for calculating the second noise window power and the 3rd coefficient obtains the 3rd product, and calculates its peak work Rate obtains the 4th product with the product of the 4th coefficient；

Select the maximum in the 3rd product and the 4th product, however, it is determined that the maximum is more than or equal to institute Peak power is stated, determines that the second effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak value Power, determines that the second effective diameter thresholding is the maximum.

In possible implementation method, according to the time domain channel estimated result after the adjustment, the second effective diameter thresholding And the timing offset, scan for determining the delay spread value of current area and the adjacent cell to tap search space, Including：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding, it is determined that first path position is：I=-N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFTFFT points are represented, p represents port numbers, The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath Bar number, (τ '₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that tail footpath position is： Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,Represent port The call number of the pilot tone column of p,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ’₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that the tail footpath position subtracts the result obtained by the first path position, by the result be defined as current area with The delay spread value of the adjacent cell.

A kind of timing offset estimation unit is additionally provided in the embodiment of the present invention, including：

First processing module, for the public reference signal CRS of the adjacent cell according to current area, to carry out channel estimation true Domain channel estimation results, determine the first noise window power and peak power of the time domain channel estimated result；

Second processing module, for determining the first effective diameter according to the first noise window power and the peak power Thresholding, according to the time domain channel estimated result and the first effective diameter thresholding, determines current area with the adjacent cell The first timing offset value；

3rd processing module, it is inclined for the timing according to the first timing offset value to the time domain channel estimated result Difference is adjusted, it is determined that the second noise window power of the time domain channel estimated result after adjustment, according to the second noise window work( Rate and the peak power determine the second effective diameter thresholding；

Fourth processing module, for according to the time domain channel estimated result after the adjustment and second effective diameter door Limit, scans for determining the second timing offset value of current area and the adjacent cell to tap search space, according to described the One timing deviation and the second timing offset value determine the timing offset between current area and the adjacent cell.

In possible implementation method, also including the 5th processing module, it is used for：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding and the timing offset, Scan for determining the delay spread value of current area and the adjacent cell to tap search space.

In possible implementation method, the first processing module specifically for：

CRS according to the adjacent cell carries out least-squares estimation and obtains channel estimation in frequency domain result, and the frequency domain is believed Road estimated result carries out inverse discrete Fourier transform, obtains the time domain channel estimated result.

In possible implementation method, the first processing module specifically for：

Calculate the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, it is expressed as：Wherein,Represent the of port pThe time domain channel of row pilot tone estimates knot Really, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the first noise window power is：Wherein, N_τ Represent the time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the peak power is：

In possible implementation method, the Second processing module specifically for：

The product for calculating the first noise window power and the first coefficient obtains the first product, and calculates its peak work Rate obtains the second product with the product of the second coefficient；

Select the maximum in first product and second product, however, it is determined that the maximum is more than or equal to institute Peak power is stated, determines that the first effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak value Power, determines that the first effective diameter thresholding is the maximum.

In possible implementation method, the Second processing module specifically for：

Determine that current area is with the first timing offset value of the adjacent cell： Wherein, P represents total port number, and R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.

In possible implementation method, the 3rd processing module specifically for：

By when the corresponding time-domain taps of the first timing offset are moved to zero described in the time domain channel estimated result Domain tap, is successively translated each time-domain taps in the time domain channel estimated result, and the time domain channel after being adjusted is estimated Counting result is：Wherein, k=0 ..., N_FFT- 1, N_FFTTable Show that FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel of row pilot tone estimates knot Really.

In possible implementation method, the 3rd processing module specifically for：

The product for calculating the second noise window power and the 3rd coefficient obtains the 3rd product, and calculates its peak work Rate obtains the 4th product with the product of the 4th coefficient；

Select the maximum in the 3rd product and the 4th product, however, it is determined that the maximum is more than or equal to institute Peak power is stated, determines that the second effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak value Power, determines that the second effective diameter thresholding is the maximum.

In possible implementation method, the 5th processing module specifically for：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding, it is determined that first path position is：I=-N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFTFFT points are represented, p represents port numbers,The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFT's The bar number in footpath, (τ '₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing partially Difference；

Determine that tail footpath position is： Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,Represent end The call number of the pilot tone column of mouth p,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ’₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that the tail footpath position subtracts the result obtained by the first path position, by the result be defined as current area with The delay spread value of the adjacent cell.

The embodiment of the present invention additionally provides a kind of terminal, and the terminal mainly includes processor and memory, wherein, memory In preserve default program, processor is used to read the program preserved in memory, and procedure below is performed according to the program：

The public reference signal CRS of the adjacent cell according to current area carries out channel estimation and determines that time domain channel estimates knot Really, the first noise window power and peak power of the time domain channel estimated result are determined；

First effective diameter thresholding is determined according to the first noise window power and the peak power, according to the time domain Channel estimation results and the first effective diameter thresholding, determine the first timing offset value of current area and the adjacent cell；

The timing offset of the time domain channel estimated result is adjusted according to the first timing offset value, it is determined that adjusting Second noise window power of the time domain channel estimated result after whole, according to the second noise window power and the peak power Determine the second effective diameter thresholding；

According to the time domain channel estimated result after the adjustment and the second effective diameter thresholding, to tap search space Scan for determining the second timing offset value of current area and the adjacent cell, according to the first timing offset value and described Second timing offset value determines the timing offset between current area and the adjacent cell.

In possible implementation method, processor is according to the time domain channel estimated result after the adjustment, described second effective Footpath thresholding and the timing offset, scan for determining that current area is expanded with the time delay of the adjacent cell to tap search space Exhibition value.

In possible implementation method, processor carries out least-squares estimation and obtains frequency domain letter according to the CRS of the adjacent cell Road estimated result, inverse discrete Fourier transform is carried out to the channel estimation in frequency domain result, is obtained the time domain channel and is estimated knot Really.

In possible implementation method, processor determines the first noise window power and peak of the time domain channel estimated result It is worth power, detailed process is：

Calculate the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, it is expressed as：Its In,Represent the of port pThe time domain channel of row pilot tone is estimated As a result, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the first noise window power is：Wherein, N_τ Represent the time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the peak power is：

In possible implementation method, processor calculates the first noise window power and the product of the first coefficient obtains first Product, and the calculating peak power obtains the second product with the product of the second coefficient；

Select the maximum in first product and second product, however, it is determined that the maximum is more than or equal to the peak Value power, determines that the first effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, determines that first has Effect footpath thresholding is the maximum.

In possible implementation method, processor determines that current area is with the first timing offset value of the adjacent cell：Wherein,P tables Show total port number, R represents reception antenna sum, wherein Γ^{(r, p)}Represent the first effective diameter thresholding.

In possible implementation method, processor is according to the first timing offset value to the time domain channel estimated result Timing offset is adjusted, and is with process：

By when the corresponding time-domain taps of the first timing offset are moved to zero described in the time domain channel estimated result Domain tap, is successively translated each time-domain taps in the time domain channel estimated result, and the time domain channel after being adjusted is estimated Counting result is：, wherein, k=0 ..., N_FFT- 1, N_FFTTable Show that FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel of row pilot tone estimates knot Really.

In possible implementation method, processor determines that the process of the second effective diameter thresholding is：

The product for calculating the second noise window power and the 3rd coefficient obtains the 3rd product, and calculates its peak work Rate obtains the 4th product with the product of the 4th coefficient；The maximum in the 3rd product and the 4th product is selected, if really The fixed maximum is more than or equal to the peak power, the second effective diameter thresholding is determined for the peak power, if really The fixed maximum is less than the peak power, determines that the second effective diameter thresholding is the maximum.

In possible implementation method, processor is according to the time domain channel estimated result after the adjustment, described second effective Footpath thresholding, it is determined that first path position is： I=-N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFTRepresent FFT points, p Represent port numbers,The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent The N of reservation_FFTFootpath bar number, (τ '₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent institute State the second timing offset；

Determine that tail footpath position is： Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,Represent end The call number of the pilot tone column of mouth p,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ’₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that the tail footpath position subtracts the result obtained by the first path position, by the result be defined as current area with The delay spread value of the adjacent cell.

Based on above-mentioned technical proposal, in the embodiment of the present invention, there is the situation of time-delay deviation in adjacent cell and current area Under, the CRS according to adjacent cell carries out channel estimation and obtains time domain channel estimated result, determines the of the time domain channel estimated result One noise window power and peak power, after determining the first effective diameter thresholding according to the first noise window power and peak power, adopt Thick synchronization gain is carried out with the first effective diameter thresholding carries out the first timing offset value of thick synchronization, using the first timing offset value pair The timing offset of time domain channel estimated result carries out thick synchronous adjustment, it is determined that time domain channel estimated result after adjustment second is made an uproar Sound window power, the second effective diameter thresholding is determined according to the second noise window power and peak power, using the second effective diameter thresholding Scan for determining the second timing offset value of precise synchronization to tap search space, deposited with local area such that it is able to the timing in adjacent area Shift to an earlier date or delayed larger, in the case of being likely to result in effective diameter exit window, time domain channel is estimated by coarse timing synchronization carry out Synchronous adjustment so that effective diameter obtains effective effective diameter thresholding without falling into noise window, and then obtain and relatively accurate determine When estimation of deviation measurement and delay spread measurement.

Brief description of the drawings

Fig. 1 is the implementation process schematic diagram of existing delay spread algorithm for estimating；

Fig. 2 is the method flow schematic diagram of timing offset estimation in the embodiment of the present invention；

Fig. 3 is the structural representation of timing offset estimation unit in the embodiment of the present invention；

Fig. 4 is terminal structure schematic diagram in the embodiment of the present invention.

Specific embodiment

In order that the object, technical solutions and advantages of the present invention are clearer, below in conjunction with accompanying drawing the present invention is made into One step ground is described in detail, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole implementation Example.Based on the embodiment in the present invention, what those of ordinary skill in the art were obtained under the premise of creative work is not made All other embodiment, belongs to the scope of protection of the invention.

In the embodiment of the present invention, as shown in Fig. 2 the method detailed flow that timing offset is estimated is as follows：

Step 201：The CRS of the adjacent cell according to current area carries out channel estimation and determines time domain channel estimated result, really The first noise window power and peak power of domain channel estimation results.

Specifically, the CRS of the adjacent cell according to current area carries out channel estimation and determines time domain channel estimated result, specifically For：

CRS according to adjacent cell carries out least square (Least Square, LS) and estimates to obtain channel estimation in frequency domain result, Channel estimation in frequency domain result is expressed as formula (1)：

Wherein, r represents that reception antenna number is indexed, and p represents that antenna end slogan is indexed,The pilot tone columns of port p is represented, Assuming that reception antenna sum is R, total port number is P；

To channel estimation in frequency domain resultInverse discrete Fourier transform is carried out, time domain channel estimated result is obtainedIt is expressed as formula (2)：

Wherein,Represent the of port pRow pilot tone Time domain channel estimated result, N_FFTFFT points are represented, r represents that reception antenna number is indexed.

Specifically, it is determined that the detailed process of the first noise window power of time domain channel estimated result and peak power is such as Under：

Calculate time domain channel estimated result in each port each column pilot tone time domain channel estimated result mould square, It is expressed as formula (3)：

Formula (3)

Wherein,Represent the of port pRow pilot tone when Domain channel estimation results, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According in time domain channel estimated result the mould of the time domain channel estimated result of each column pilot tone of each port square, Determine the first noise window power for shown in formula (4)：

Wherein, N_τRepresent the time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According in time domain channel estimated result the mould of the time domain channel estimated result of each column pilot tone of each port square, Determine peak power for shown in formula (5)：

Step 202：First effective diameter thresholding is determined according to the first noise window power and peak power, according to time domain channel Estimated result and the first effective diameter thresholding, determine the first timing offset value of current area and adjacent cell.

Specifically, it is determined that the process of the first effective diameter thresholding is：Calculate the product of the first noise window power and the first coefficient The first product is obtained, and calculating peak power obtains the second product with the product of the second coefficient；Select the first product and second Maximum in product, however, it is determined that the maximum is more than or equal to the peak power, determines that the first effective diameter thresholding is the peak value Power, however, it is determined that the maximum is less than the peak power, determines that the first effective diameter thresholding is the maximum.

Wherein, the first coefficient and the second coefficient can be determined by emulating.

First effective diameter thresholding formula (6) is expressed as：

Wherein,The first noise window power is represented,Represent peak power, β_noiseRepresent the first coefficient, β_maxTable Show the second coefficient.If met according to the first effective diameter thresholding that formula (6) determinesThen have to first again Effect footpath thresholding is entered as：

Specifically, according to time domain channel estimated result and the first effective diameter thresholding, determine current area with adjacent cell First timing offset value, detailed process is：

Determine current area with the first timing offset value of adjacent cell for shown in formula (7)：

Wherein,I=-N_leak,…, N_τ, P represents total port number, and R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.Wherein,'s Physical significance is in i=-N_leak,…,N_τHunting zone in find causeMore than or equal to Γ^(r,p)Minimum rope Draw.

Step 203：The timing offset of time domain channel estimated result is adjusted according to the first timing offset value, it is determined that adjusting Second noise window power of the time domain channel estimated result after whole, second is determined according to the second noise window power and peak power Effective diameter thresholding.

Wherein, the span of noise window is to preset.

Specifically, the timing offset of time domain channel estimated result is adjusted according to the first timing offset value, specially：

The corresponding time-domain taps of first timing offset are moved to zero time-domain taps and obtain the estimation of the second time domain channel Result is shown in formula (8)：

Wherein, k=0 ..., N_FFT- 1, N_FFTRepresent FFT points, τ₁The first timing offset value is represented,Represent port The of pThe time domain channel estimated result of row pilot tone.

Respectively according to formula (4) and formula (5), according to the time domain channel estimated result after adjustment, to every reception antenna Each port is in noise window [N_τ:N_FFT-N_leak] in the range of carry out noise power calculation and obtain the second noise window powerIf Putting the 3rd coefficient isSetting the 4th coefficient isWherein, the 3rd coefficient and the 4th coefficient can be determined by emulating.

Specifically, it is determined that the process of the second effective diameter thresholding is：Calculate the product of the second noise window power and the 3rd coefficient The 3rd product is obtained, and calculating peak power obtains the 4th product with the product of the 4th coefficient；

Select the maximum in the 3rd product and the 4th product, however, it is determined that the maximum is more than or equal to peak power, really Fixed second effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, determines the second effective diameter thresholding It is the maximum.

Second effective diameter thresholding formula (9) is expressed as：

If met according to the second effective diameter thresholding that formula (9) determinesThen again to the second effective diameter Thresholding is entered as：

Step 204：According to the time domain channel estimated result after adjustment and the second effective diameter thresholding, to tap search space Scan for determining the second timing offset value of current area and adjacent cell, according to the first timing offset value and the second timing offset Value determines the timing offset between current area and adjacent cell.

Wherein, the span of tap search space is to preset, specifically, the span of tap search space with The span of noise window is complementary relationship, it is assumed that search space total length is N_FFTIf, part point setting in search space It is spatial noise, then another part is tap search space.

Specifically, it is determined that to tap search space [N after the second effective diameter thresholding_FFT-N_leak:N_FFT-1,0:N_τ] searched Rope, obtains the second timing offset value τ 1 ' of current area and adjacent cell, is expressed as with formula (10) and formula (11)：

Wherein, i=-N_leak..., N_τ, P represents total port number, and R represents reception antenna sum, whereinRepresent the Two effective diameter thresholdings.Wherein, the peak power of the time domain channel estimated result after adjustment is estimated with the time domain channel before adjustment The peak power of result is identical.

Specifically, calculate the first timing offset value and the second timing offset value and, by gained and value be defined as currently Timing offset between cell and adjacent cell.

In implementation, can also determining according to the time domain channel estimated result after adjustment, the second effective diameter thresholding and determination When deviation, scan for determining the delay spread value of current area and adjacent cell to tap search space.

Specifically, according to the time domain channel estimated result after adjustment, the second effective diameter thresholding, it is determined that first path position is formula (12) shown in：

Wherein, i=-N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFT FFT points are represented, p represents port numbers,The call number of the pilot tone column of port p is represented,Represent the second effective diameter Thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ '₁+τ₁) represent timing offset, τ₁Represent the first timing offset, τ '₁Table Show the second timing offset；

Determine tail footpath position for shown in formula (13)：

Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFT's The bar number in footpath, (τ '₁+τ₁) represent timing offset, τ₁Represent the first timing offset, τ '₁Represent the second timing offset；

Determine that tail footpath position subtracts the result obtained by first path position, by the result be defined as current area and adjacent cell when Prolong expanding value.

Wherein, in formula (12) and formula (13) Wherein,Represent the of port pAfter the adjustment of row pilot tone when Domain channel estimation results, r represents that reception antenna number is indexed, N_FFTRepresent FFT points.

The calculating of timing offset and delay spread below by way of two specific embodiments to being provided in the embodiment of the present invention Process is illustrated.

Specific embodiment one, under LTE FEICIC macro-pico deployment scenarios, has one with a width of 10MHz macrocells A, has Pico cells a B, UE with a width of 10MHz to be in the edge of Pico cells, and capable and experienced by the CRS of macrocell Disturb, A cells reach UE and are advanced by 3us with the time delay of B cells arrival UE, and specific implementation step is as follows：

Step 1, the LS estimations of CRS are carried out to adjacent area A, obtain channel estimation in frequency domain result

Wherein,R represents that reception antenna number is indexed, and p represents antenna end slogan rope Draw,Represent the pilot tone columns of port p, it is assumed that reception antenna sum is R, total port number is P.

Step 2, to channel estimation in frequency domain resultCarry out IDFT conversion and obtain time domain channel estimated result It is expressed asWherein, Represent the of port pThe time domain channel estimated result of row pilot tone, N_FFTFFT points are represented, r represents that reception antenna number is indexed, Wherein, N_FFT=256.

Step 3, according toCalculate the first noise window general powerAnd peak powerWherein,

Represent the of port pThe time domain channel of row pilot tone Estimated result, r represents that reception antenna number is indexed, N_FFTRepresent FFT points.N_τRepresent the time-domain position of time delay maximum diameter, N_leakTable Show the N of reservation_FFTFootpath bar number.

Wherein,

Step 4, carries out coarse timing synchronization, obtains the first timing offset value τ₁。

Specifically, the first effective diameter thresholding is setWherein,Table Show the first noise window power,Represent the first peak power, β_noiseRepresent the first coefficient, β_maxRepresent the second coefficient.Wherein, β_maxArranges value is more than β_noiseSo that there is the first effect footpath thresholding of a relatively high, first footpath power is maximum under General System, effective diameter door The limit relative precision that can ensure timing estimation higher.β_maxCan be determined by emulating.If determine according to formula first has Effect footpath thresholding meetsThe first effective diameter thresholding is entered as again then：

First timing offset value τ₁It is expressed as：Wherein,P tables Show total port number, R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.

Step 5, according to the first timing offset value τ₁Calculate the second effective diameter thresholding

Specifically, according to τ₁AdjustmentPhysical meaning is that coarse timing synchronization is estimated into corresponding time-domain taps movement Position to the 0th time-domain taps obtainsFinally give the time domain channel estimated result after adjustmentSpecifically It is as follows：Wherein, k=0 ..., N_FFT- 1, N_FFTRepresent FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel estimated result of row pilot tone.

SetWithTo every reception antenna, each portIn noise window [N_τ:N_FFT-N_leak] model Being calculated for noise power is carried out in enclosingAccording to every reception antenna, each portIn peak powerAnd noise powerIt is determined thatWherein,WithCan With by emulating determination.If met according to the second effective diameter thresholding that formula determinesThen have to second again Effect footpath thresholding is entered as：

Step 6, estimates time domain channel to be timed synchro measure and delay spread measurement.

Specifically, to tap search space [N_FFT-N_leak:N_FFT-1,0:N_τ] scan for, obtain Timing Synchronization measured value Or delay spread measured value, wherein, can be respectively provided with for timing offset measured value and delay spread measured value differentWithTo meet different demands.

Specifically, timing offset measured value is expressed as τ '₁+τ₁, τ '₁Acquisition process be： Wherein, P represents total port number, and R represents reception antenna sum, whereinRepresent the second effective diameter thresholding.Wherein,

Specifically, the acquisition process of delay spread measured value is：

It is determined that first path position is

Wherein, i =-N_leak,…,N_τ, N_τRepresent the time-domain position of time delay maximum diameter；

Determine that tail footpath position is

Wherein, I=-N_leak,…,N_τ；

Determine that tail footpath position subtracts the result obtained by first path position, the result is current area and expands with the time delay of adjacent cell Exhibition value.

Specific embodiment two, under LTE FEICIC macro-pico deployment scenarios, has one with a width of 10MHz macrocells A, has Pico cells a B, UE with a width of 10MHz to be in the edge of Pico cells, and capable and experienced by the CRS of macrocell Disturb, A cells reach time delay delayed 3us, specific implementation process and the step 1 in specific embodiment one that UE and B cells reach UE The description of~step 6 is identical, and here is omitted.

A kind of timing offset estimation unit is additionally provided based on same inventive concept, in the embodiment of the present invention, the device Specific implementation can be found in the description of above method embodiment, repeats part and repeats no more, as shown in figure 3, the device mainly includes：

First processing module 301, the public reference signal CRS for the adjacent cell according to current area carries out channel estimation Determine time domain channel estimated result, determine the first noise window power and peak power of the time domain channel estimated result；

Second processing module 302, for determining that first has according to the first noise window power and the peak power Effect footpath thresholding, according to the time domain channel estimated result and the first effective diameter thresholding, determines current area with the neighbour First timing offset value of cell；

3rd processing module 303, for being determined the time domain channel estimated result according to the first timing offset value When deviation be adjusted, it is determined that adjustment after time domain channel estimated result the second noise window power, according to second noise Window power and the peak power determine the second effective diameter thresholding；

Fourth processing module 304, for according to the time domain channel estimated result after the adjustment and described second effective Footpath thresholding, scans for determining the second timing offset value of current area and the adjacent cell, according to institute to tap search space State the timing offset that the first timing offset value and the second timing offset value are determined between current area and the adjacent cell.

In implementation, also including the 5th processing module 305, it is used for：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding and the timing offset, Scan for determining the delay spread value of current area and the adjacent cell to tap search space.

In implementation, the first processing module specifically for：

CRS according to the adjacent cell carries out least-squares estimation and obtains channel estimation in frequency domain result, and the frequency domain is believed Road estimated result carries out inverse discrete Fourier transform, obtains the time domain channel estimated result.

In implementation, the first processing module specifically for：

Calculate the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, it is expressed as：Its In,Represent the of port pThe time domain channel of row pilot tone is estimated As a result, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the first noise window power is：Wherein, N_τ Represent the time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the peak power is：

In implementation, the Second processing module specifically for：

The product for calculating the first noise window power and the first coefficient obtains the first product, and calculates its peak work Rate obtains the second product with the product of the second coefficient；

Select the maximum in first product and second product, however, it is determined that the maximum is more than or equal to the peak Value power, determines that the first effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, determines that first has Effect footpath thresholding is the maximum.

In implementation, the Second processing module specifically for：

Determine that current area is with the first timing offset value of the adjacent cell： Wherein, P represents total port number, and R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.

In implementation, the 3rd processing module specifically for：

By when the corresponding time-domain taps of the first timing offset are moved to zero described in the time domain channel estimated result Domain tap, is successively translated each time-domain taps in the time domain channel estimated result, and the time domain channel after being adjusted is estimated Counting result is：Wherein, k=0 ..., N_FFT- 1, N_FFTTable Show that FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel of row pilot tone estimates knot Really.

In implementation, the 3rd processing module specifically for：

The product for calculating the second noise window power and the 3rd coefficient obtains the 3rd product, and calculates its peak work Rate obtains the 4th product with the product of the 4th coefficient；

Select the maximum in the 3rd product and the 4th product, however, it is determined that the maximum is more than or equal to institute Peak power is stated, determines that the second effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak value Power, determines that the second effective diameter thresholding is the maximum.

In implementation, the 5th processing module specifically for：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding, it is determined that first path position is：I=-N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFTFFT points are represented, p represents port numbers, The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath Bar number, (τ '₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that tail footpath position is： Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,Represent end The call number of the pilot tone column of mouth p,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ’₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that the tail footpath position subtracts the result obtained by the first path position, by the result be defined as current area with The delay spread value of the adjacent cell.

A kind of terminal is additionally provided based on same inventive concept, in the embodiment of the present invention, the specific implementation of the terminal can be joined See the specific descriptions of above method embodiment, repeat part and repeat no more, as shown in figure 4, the terminal mainly includes processor 401 With memory 402, wherein, preserve default program in memory 402, processor 401 is used to reading and is preserved in memory 402 Program, according to the program perform procedure below：

The public reference signal CRS of the adjacent cell according to current area carries out channel estimation and determines that time domain channel estimates knot Really, the first noise window power and peak power of the time domain channel estimated result are determined；

First effective diameter thresholding is determined according to the first noise window power and the peak power, according to the time domain Channel estimation results and the first effective diameter thresholding, determine the first timing offset value of current area and the adjacent cell；

The timing offset of the time domain channel estimated result is adjusted according to the first timing offset value, it is determined that adjusting Second noise window power of the time domain channel estimated result after whole, according to the second noise window power and the peak power Determine the second effective diameter thresholding；

According to the time domain channel estimated result after the adjustment and the second effective diameter thresholding, to tap search space Scan for determining the second timing offset value of current area and the adjacent cell, according to the first timing offset value and described Second timing offset value determines the timing offset between current area and the adjacent cell.

In implementation, processor according to the time domain channel estimated result after the adjustment, the second effective diameter thresholding and The timing offset, scans for determining the delay spread value of current area and the adjacent cell to tap search space.

In implementation, processor carries out least-squares estimation and obtains channel estimation in frequency domain result according to the CRS of the adjacent cell, Inverse discrete Fourier transform is carried out to the channel estimation in frequency domain result, the time domain channel estimated result is obtained.

In implementation, processor determines the first noise window power and peak power of the time domain channel estimated result, tool Body process is：

Calculate the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, it is expressed as：Its In,Represent the of port pThe time domain channel of row pilot tone is estimated As a result, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the first noise window power is：Wherein, N_τ Represent the time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According to the mould of the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result Square, determine that the peak power is：

In implementation, processor calculates the first noise window power and the product of the first coefficient obtains the first product, and Calculate the peak power and obtain the second product with the product of the second coefficient；

Select the maximum in first product and second product, however, it is determined that the maximum is more than or equal to the peak Value power, determines that the first effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, determines that first has Effect footpath thresholding is the maximum.

In implementation, processor determines that current area is with the first timing offset value of the adjacent cell： Wherein, P represents total port number, and R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.

In implementation, processor enters according to the first timing offset value to the timing offset of the time domain channel estimated result Row adjustment, be with process：

By when the corresponding time-domain taps of the first timing offset are moved to zero described in the time domain channel estimated result Domain tap, is successively translated each time-domain taps in the time domain channel estimated result, and the time domain channel after being adjusted is estimated Counting result is：Wherein, k=0 ..., N_FFT- 1, N_FFTTable Show that FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel of row pilot tone estimates knot Really.

In implementation, processor determines that the process of the second effective diameter thresholding is：

The product for calculating the second noise window power and the 3rd coefficient obtains the 3rd product, and calculates its peak work Rate obtains the 4th product with the product of the 4th coefficient；The maximum in the 3rd product and the 4th product is selected, if really The fixed maximum is more than or equal to the peak power, the second effective diameter thresholding is determined for the peak power, if really The fixed maximum is less than the peak power, determines that the second effective diameter thresholding is the maximum.

In implementation, processor according to the time domain channel estimated result after the adjustment, the second effective diameter thresholding, it is determined that First path position is：I=- N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFTFFT points are represented, p is represented Port numbers,The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent and retain N_FFTFootpath bar number, (τ '₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent described Two timing offset；

Determine that tail footpath position is： Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,Represent end The call number of the pilot tone column of mouth p,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ’₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that the tail footpath position subtracts the result obtained by the first path position, by the result be defined as current area with The delay spread value of the adjacent cell.

Wherein, connected by bus between processor and memory, bus architecture can include any number of interconnection The various circuits of the memory that bus and bridge, the one or more processors for specifically being represented by processor and memory are represented are linked Together.Bus architecture can also be by various other circuits of such as ancillary equipment, voltage-stablizer and management circuit or the like Link together, these are all it is known in the art, therefore, being no longer described further to it herein.EBI is carried For interface.Processor is responsible for bus architecture and common treatment, and memory can store processor and perform operation when institute The data for using.

Based on above-mentioned technical proposal, in the embodiment of the present invention, there is the situation of time-delay deviation in adjacent cell and current area Under, the CRS according to adjacent cell carries out channel estimation and obtains time domain channel estimated result, determines the of the time domain channel estimated result One noise window power and peak power, after determining the first effective diameter thresholding according to the first noise window power and peak power, adopt Thick synchronization gain is carried out with the first effective diameter thresholding carries out the first timing offset value of thick synchronization, using the first timing offset value pair The timing offset of time domain channel estimated result carries out thick synchronous adjustment, it is determined that time domain channel estimated result after adjustment second is made an uproar Sound window power, the second effective diameter thresholding is determined according to the second noise window power and peak power, using the second effective diameter thresholding Scan for determining the second timing offset value of precise synchronization to tap search space, deposited with local area such that it is able to the timing in adjacent area Shift to an earlier date or delayed larger, in the case of being likely to result in effective diameter exit window, time domain channel is estimated by coarse timing synchronization carry out Synchronous adjustment so that effective diameter obtains effective effective diameter thresholding without falling into noise window, and then obtain and relatively accurate determine When estimation of deviation measurement and delay spread measurement.

It should be understood by those skilled in the art that, embodiments of the invention can be provided as method, system or computer program Product.Therefore, the present invention can be using the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware Apply the form of example.And, the present invention can be used and wherein include the computer of computer usable program code at one or more The shape of the computer program product implemented in usable storage medium (including but not limited to magnetic disk storage and optical memory etc.) Formula.

The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product Figure and/or block diagram are described.It should be understood that every first-class during flow chart and/or block diagram can be realized by computer program instructions The combination of flow and/or square frame in journey and/or square frame and flow chart and/or block diagram.These computer programs can be provided The processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce A raw machine so that produced for reality by the instruction of computer or the computing device of other programmable data processing devices The device of the function of being specified in present one flow of flow chart or multiple one square frame of flow and/or block diagram or multiple square frames.

These computer program instructions may be alternatively stored in can guide computer or other programmable data processing devices with spy In determining the computer-readable memory that mode works so that instruction of the storage in the computer-readable memory is produced and include finger Make the manufacture of device, the command device realize in one flow of flow chart or multiple one square frame of flow and/or block diagram or The function of being specified in multiple square frames.

These computer program instructions can be also loaded into computer or other programmable data processing devices so that in meter Series of operation steps is performed on calculation machine or other programmable devices to produce computer implemented treatment, so as in computer or The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one The step of function of being specified in individual square frame or multiple square frames.

Obviously, those skilled in the art can carry out various changes and modification without deviating from essence of the invention to the present invention God and scope.So, if these modifications of the invention and modification belong to the scope of the claims in the present invention and its equivalent technologies Within, then the present invention is also intended to comprising these changes and modification.

Claims (18)

1. a kind of timing offset method of estimation, it is characterised in that including：

The public reference signal CRS of the adjacent cell according to current area carries out channel estimation and determines time domain channel estimated result, really The first noise window power and peak power of the fixed time domain channel estimated result；

First effective diameter thresholding is determined according to the first noise window power and the peak power, according to the time domain channel Estimated result and the first effective diameter thresholding, determine the first timing offset value of current area and the adjacent cell；

The timing offset of the time domain channel estimated result is adjusted according to the first timing offset value, it is determined that after adjustment Time domain channel estimated result the second noise window power, determined according to the second noise window power and the peak power Second effective diameter thresholding；

According to the time domain channel estimated result after the adjustment and the second effective diameter thresholding, tap search space is carried out Search determines the second timing offset value of current area and the adjacent cell, according to the first timing offset value and described second Timing offset value determines the timing offset between current area and the adjacent cell.

2. the method for claim 1, it is characterised in that methods described also includes：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding and the timing offset, to taking out Head search space scans for determining the delay spread value of current area and the adjacent cell.

3. method as claimed in claim 1 or 2, it is characterised in that the public reference signal of the adjacent cell according to current area CRS carries out channel estimation and determines time domain channel estimated result, including：

CRS according to the adjacent cell carries out least-squares estimation and obtains channel estimation in frequency domain result, and the frequency domain channel is estimated Meter result carries out inverse discrete Fourier transform, obtains the time domain channel estimated result.

4. method as claimed in claim 1 or 2, it is characterised in that determine the first noise of the time domain channel estimated result Window power and peak power, including：

Calculate the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result mould square, It is expressed as： ${\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}=\[|h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(0\right){|}^2,|h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(1\right){|}^2,...,|h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}(N_{FFT}-1){|}^2\],$ Wherein,Represent the of port pThe time domain channel of row pilot tone estimates knot Really, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According in the time domain channel estimated result mould of the time domain channel estimated result of each column pilot tone of each port square, Determine that the first noise window power is： $P_{noise}^{(r,p)}=\underset{j=0}{\overset{L_{RS}^{\left(p\right)}-1}{\Σ}}\underset{k=-N_{\mathrm{\τ}}}{\overset{N_{FFT}-N_{leak}}{\Σ}}{\mathrm{pow}}_{cir,j}^{(r,p)}\left(\right(k+N_{FFT}\left){\mathrm{\%N}}_{FFT}\right),$ Wherein, N_τRepresent The time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According in the time domain channel estimated result mould of the time domain channel estimated result of each column pilot tone of each port square, Determine that the peak power is：

5. method as claimed in claim 4, it is characterised in that according to the first noise window power and the peak power Determine the first effective diameter thresholding, including：

The product for calculating the first noise window power and the first coefficient obtains the first product, and calculate the peak power with The product of the second coefficient obtains the second product；

Select the maximum in first product and second product, however, it is determined that the maximum is more than or equal to the peak Value power, determines that the first effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, Determine that the first effective diameter thresholding is the maximum.

6. method as claimed in claim 5, it is characterised in that had according to the time domain channel estimated result and described first Effect footpath thresholding, determines the first timing offset value of current area and the adjacent cell, including：

Determine that current area is with the first timing offset value of the adjacent cell： ${\mathrm{\τ}}_1=\underset{r=0}{\overset{R}{\Σ}}\underset{p=0}{\overset{P}{\Σ}}\underset{l_{RS}^{\left(p\right)}=0}{\overset{L_{RS}^{\left(p\right)}}{\Σ}}{\mathrm{\τ}}_{l_{RS}^{\left(p\right)}}^{(r,p)}/(L_{RS}^{\left(p\right)}\·P\·R),$ Wherein, ${\mathrm{\τ}}_{l_{RS}^{\left(p\right)}}^{(r,p)}=\arg \underset{i}{min}\{{\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(i+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right)\≥{\mathrm{\Γ}}^{(r,p)}\},i=-N_{leak},...,N_{\mathrm{\τ}},$ P represents total port number, and R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.

7. method as claimed in claim 1 or 2, it is characterised in that believed the time domain according to the first timing offset value The timing offset of road estimated result is adjusted, including：

The corresponding time-domain taps of first timing offset described in the time domain channel estimated result are moved into zero time domain to take out Head, is successively translated each time-domain taps in the time domain channel estimated result, and the time domain channel after being adjusted estimates knot It is really： ${\tilde{h}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(k\right)=h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(k+{\mathrm{\τ}}_1+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right),$ Wherein, k=0 ..., N_FFT- 1, N_FFTRepresent FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel estimated result of row pilot tone.

8. method as claimed in claim 7, it is characterised in that according to the second noise window power and the peak power Determine the second effective diameter thresholding, including：

The product for calculating the second noise window power and the 3rd coefficient obtains the 3rd product, and calculate the peak power with The product of the 4th coefficient obtains the 4th product；

Select the maximum in the 3rd product and the 4th product, however, it is determined that the maximum is more than or equal to the peak Value power, determines that the second effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, Determine that the second effective diameter thresholding is the maximum.

9. method as claimed in claim 2, it is characterised in that according to the time domain channel estimated result after the adjustment, described Second effective diameter thresholding and the timing offset, scan for determining current area and the adjacent cell to tap search space Delay spread value, including：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding, it is determined that first path position is： ${\mathrm{\τ}}_{l_{RS}^{\left(p\right)}}^{(r,p)}=\arg \underset{i}{min}\{{\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(i+\mathrm{\τ}{1}^{,}+\mathrm{\τ}1+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right)\≥{\widetilde{\mathrm{\Γ}}}^{(r,p)}\},$ I=-N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFTFFT points are represented, p represents port numbers, The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath Bar number, (τ '₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that tail footpath position is： $\mathrm{\τ}{2}_{l_{RS}^{\left(p\right)}}^{(r,p)}=\arg \underset{i}{\max }\{{\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(i+{{\mathrm{\τ}}^{,}}_1+{\mathrm{\τ}}_1+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right)\≥{\widetilde{\mathrm{\Γ}}}^{(r,p)}\},$ Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,Represent port The call number of the pilot tone column of p,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ’₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that the tail footpath position subtracts the result obtained by the first path position, by the result be defined as current area with it is described The delay spread value of adjacent cell.

10. a kind of timing offset estimation unit, it is characterised in that including：

First processing module, when the public reference signal CRS for the adjacent cell according to current area carries out channel estimation and determines Domain channel estimation results, determine the first noise window power and peak power of the time domain channel estimated result；

Second processing module, for determining the first effective diameter door according to the first noise window power and the peak power Limit, according to the time domain channel estimated result and the first effective diameter thresholding, determines current area with the adjacent cell First timing offset value；

3rd processing module, for being entered to the timing offset of the time domain channel estimated result according to the first timing offset value Row adjustment, it is determined that adjustment after time domain channel estimated result the second noise window power, according to the second noise window power with And the peak power determines the second effective diameter thresholding；

Fourth processing module, for according to the time domain channel estimated result after the adjustment and the second effective diameter thresholding, Scan for determining the second timing offset value of current area and the adjacent cell to tap search space, it is fixed according to described first When deviation and the second timing offset value timing offset that determines between current area and the adjacent cell.

11. devices as claimed in claim 10, it is characterised in that also including the 5th processing module, be used for：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding and the timing offset, to taking out Head search space scans for determining the delay spread value of current area and the adjacent cell.

12. device as described in claim 10 or 11, it is characterised in that the first processing module specifically for：

CRS according to the adjacent cell carries out least-squares estimation and obtains channel estimation in frequency domain result, and the frequency domain channel is estimated Meter result carries out inverse discrete Fourier transform, obtains the time domain channel estimated result.

13. device as described in claim 10 or 11, it is characterised in that the first processing module specifically for：

Calculate the time domain channel estimated result of each column pilot tone of each port in the time domain channel estimated result mould square, It is expressed as： ${\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}=\[|h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(0\right){|}^2,|h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(1\right){|}^2,...,|h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}(N_{FFT}-1){|}^2\],$ Wherein,Represent the of port pThe time domain channel of row pilot tone estimates knot Really, r represents that reception antenna number is indexed, N_FFTRepresent FFT points；

According in the time domain channel estimated result mould of the time domain channel estimated result of each column pilot tone of each port square, Determine that the first noise window power is： $P_{noise}^{(r,p)}=\underset{j=0}{\overset{L_{RS}^{\left(p\right)}-1}{\Σ}}\underset{k=-N_{\mathrm{\τ}}}{\overset{N_{FFT}-N_{leak}}{\Σ}}{\mathrm{pow}}_{cir,j}^{(r,p)}\left(\right(k+N_{FFT}\left){\mathrm{\%N}}_{FFT}\right),$ Wherein, N_τRepresent The time-domain position of time delay maximum diameter, N_leakRepresent the N for retaining_FFTFootpath bar number；

According in the time domain channel estimated result mould of the time domain channel estimated result of each column pilot tone of each port square, Determine that the peak power is：

14. devices as claimed in claim 13, it is characterised in that the Second processing module specifically for：

The product for calculating the first noise window power and the first coefficient obtains the first product, and calculate the peak power with The product of the second coefficient obtains the second product；

Select the maximum in first product and second product, however, it is determined that the maximum is more than or equal to the peak Value power, determines that the first effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, Determine that the first effective diameter thresholding is the maximum.

15. devices as claimed in claim 14, it is characterised in that the Second processing module specifically for：

Determine that current area is with the first timing offset value of the adjacent cell： ${\mathrm{\τ}}_1=\underset{r=0}{\overset{R}{\Σ}}\underset{p=0}{\overset{P}{\Σ}}\underset{l_{RS}^{\left(p\right)}=0}{\overset{L_{RS}^{\left(p\right)}}{\Σ}}{\mathrm{\τ}}_{l_{RS}^{\left(p\right)}}^{(r,p)}/(L_{RS}^{\left(p\right)}\·P\·R),$ Wherein, ${\mathrm{\τ}}_{l_{RS}^{\left(p\right)}}^{(r,p)}=\arg \underset{i}{min}\{{\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(i+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right)\≥{\mathrm{\Γ}}^{(r,p)}\},i=-N_{leak},...,N_{\mathrm{\τ}},$ P represents total port number, and R represents reception antenna sum, wherein Γ^(r,p)Represent the first effective diameter thresholding.

16. device as described in claim 10 or 11, it is characterised in that the 3rd processing module specifically for：

The corresponding time-domain taps of first timing offset described in the time domain channel estimated result are moved into zero time domain to take out Head, is successively translated each time-domain taps in the time domain channel estimated result, and the time domain channel after being adjusted estimates knot It is really： ${\tilde{h}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(k\right)=h_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(k+{\mathrm{\τ}}_1+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right),$ Wherein, k=0 ..., N_FFT- 1, N_FFTRepresent FFT counts, τ₁The first timing offset value is represented,Represent the of port pThe time domain channel estimated result of row pilot tone.

17. devices as claimed in claim 16, it is characterised in that the 3rd processing module specifically for：

The product for calculating the second noise window power and the 3rd coefficient obtains the 3rd product, and calculate the peak power with The product of the 4th coefficient obtains the 4th product；

Select the maximum in the 3rd product and the 4th product, however, it is determined that the maximum is more than or equal to the peak Value power, determines that the second effective diameter thresholding is the peak power, however, it is determined that the maximum is less than the peak power, Determine that the second effective diameter thresholding is the maximum.

18. devices as claimed in claim 11, it is characterised in that the 5th processing module specifically for：

According to the time domain channel estimated result after the adjustment, the second effective diameter thresholding, it is determined that first path position is： ${\mathrm{\τ}}_{l_{RS}^{\left(p\right)}}^{(r,p)}=\arg \underset{i}{min}\{{\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(i+\mathrm{\τ}{1}^{,}+\mathrm{\τ}1+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right)\≥{\widetilde{\mathrm{\Γ}}}^{(r,p)}\},$ I=-N_leak,…,N_τ, N_τThe time-domain position of time delay maximum diameter is represented, r represents that reception antenna number is indexed, N_FFTFFT points are represented, p represents port numbers, The call number of the pilot tone column of port p is represented,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath Bar number, (τ '₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that tail footpath position is： $\mathrm{\τ}{2}_{l_{RS}^{\left(p\right)}}^{(r,p)}=\arg \underset{i}{\max }\{{\mathrm{pow}}_{cir,l_{RS}^{\left(p\right)}}^{(r,p)}\left(\right(i+{{\mathrm{\τ}}^{,}}_1+{\mathrm{\τ}}_1+N_{FFT}\left)\mathrm{mod}{N}_{FFT}\right)\≥{\widetilde{\mathrm{\Γ}}}^{(r,p)}\},$ Wherein, i=-N_leak,…,N_τ, r represent reception antenna number index, N_FFTFFT points are represented, p represents port numbers,Represent end The call number of the pilot tone column of mouth p,Represent the second effective diameter thresholding, N_leakRepresent the N for retaining_FFTFootpath bar number, (τ’₁+τ₁) represent the timing offset, τ₁Represent first timing offset, τ '₁Represent second timing offset；

Determine that the tail footpath position subtracts the result obtained by the first path position, by the result be defined as current area with it is described The delay spread value of adjacent cell.