CN106304313A - The processing method and processing device of timing slip - Google Patents
The processing method and processing device of timing slip Download PDFInfo
- Publication number
- CN106304313A CN106304313A CN201510264229.0A CN201510264229A CN106304313A CN 106304313 A CN106304313 A CN 106304313A CN 201510264229 A CN201510264229 A CN 201510264229A CN 106304313 A CN106304313 A CN 106304313A
- Authority
- CN
- China
- Prior art keywords
- sequence
- value
- metric
- noise
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W56/00—Synchronisation arrangements
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W56/00—Synchronisation arrangements
- H04W56/0035—Synchronisation arrangements detecting errors in frequency or phase
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W56/00—Synchronisation arrangements
- H04W56/0055—Synchronisation arrangements determining timing error of reception due to propagation delay
- H04W56/006—Synchronisation arrangements determining timing error of reception due to propagation delay using known positions of transmitter and receiver
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
The invention discloses the processing method and processing device of a kind of timing slip, the processing method of described timing slip comprises the following steps: obtain H1 sequence；According to described H1 retrieval metric sequence P (k), and obtain maximal metric value Power of described metric sequence P (k)；Obtain signal subset S (x1) and noise subsets NI (x2) of described metric sequence P (k) according to described maximal metric value Power, calculate characteristic measure value NIPower of described noise subsets NI (x2)；First path position value is calculated according to described maximal metric value Power, characteristic measure value NIPower, noise jamming threshold T 1 and search thresholding threshold T 2；Calculate the timing offset of described subscriber equipment according to described first path position value, the present invention is applicable to the estimation of positive and negative timing offset, and timing offset accuracy is higher.
Description
Technical field
The present invention relates to mobile communication technology field, particularly relate to processing method and the dress of a kind of timing slip
Put.
Background technology
Long Term Evolution (Long Term Evolution, LTE) and senior Long Term Evolution (LTE-Advanced,
LTE-A) in system, base station (evolved NodeB, eNB or eNodeB) in order to realize uplink synchronous,
The upward signal needing each subscriber equipment (User Equipment is called for short UE) to send can be synchronized to
Reach base station, it is therefore necessary to estimate the timing slip of each UE, then adjusted by down order notice UE
Whole launch time.The timing slip of UE is typically at Physical Random Access Channel (Physical Random
Access Channel, PRACH), Physical Uplink Shared Channel (Physical Uplink Shared
Channel, PUSCH), carry out in detection reference channel (Sounding reference signal, SRS)
Measure.
The method obtaining timing slip mainly has frequency domain and time domain two kinds, and frequency domain method is to utilize channel to estimate
Frequency domain phase place inclined when calculating, time domain approach be by by channel estimation value by certain method migration
To time domain, calculate timing slip according to the amplitude of time domain, but by determining that these methods are calculated
The precision of hour offset is the highest, and for time domain approach, is also only applicable to timing inclined, for inclined time negative
Do not process.
Foregoing is only used for auxiliary and understands technical scheme, does not represent and recognizes that foregoing is
Prior art.
Summary of the invention
Present invention is primarily targeted at the processing method and processing device that a kind of timing slip is provided, it is intended to solve
Existing timing slip is estimated to be only applicable to the technical problem that timing is inclined and accuracy is low.
For achieving the above object, the present invention provides the processing method of a kind of timing slip, described timing slip
Processing method comprise the following steps:
Obtaining H1 sequence, wherein, described H1 sequence is that the subscriber equipment that receives of base station is at predetermined channel
On the sequence that obtains with reference signal sequence conjugate multiplication of frequency-region signal；
According to described H1 retrieval metric sequence P (k), and obtain the maximum of described metric sequence P (k)
Metric Power；
According to described maximal metric value Power obtain described metric sequence P (k) signal subset S (x1) and
Noise subsets NI (x2), calculates characteristic measure value NIPower of described noise subsets NI (x2)；
According to described maximal metric value Power, characteristic measure value NIPower, noise jamming threshold T 1
And search thresholding threshold T 2 calculates first path position value；
The timing offset of described subscriber equipment is calculated according to described first path position value.
Preferably, the described letter obtaining described metric sequence P (k) according to described maximal metric value Power
Work song collection S (x1) and noise subsets NI (x2), calculate the characteristic measure value of described noise subsets NI (x2)
The step of NIPower includes:
Obtain the long w1 of left window of described metric sequence P (k), the long w2 of right window and described maximal metric value
Minimum index value MaxPidx corresponding for Power, according to the described long w1 of left window, the long w2 of right window and minimum
Index value MaxPidx obtains signal subset S (x1)；
Described noise subsets NI (x2) is obtained according to described signal subset S (x1)；
Obtain the maximal metric value of described noise subsets NI (x2) or calculate described noise subsets NI (x2)
Average degree value is as described characteristic measure value NIPower.
Preferably, described according to described maximal metric value Power, characteristic measure value NIPower, noise
Interference threshold value T1 and search thresholding threshold T 2 calculate the step of first path position value and include:
Judge that whether described maximal metric value Power and the ratio of characteristic measure value NIPower are less than described
Noise jamming threshold T 1；
The most then using described minimum index value MaxPidx as described first path position value；
If it is not, then obtain in described signal subset S (x1) more than minimum index value corresponding for Power/T2
T2idx, obtains described first path position value according to described minimum index value T2idx.
Preferably, described according to described H1 retrieval metric sequence P (k), and obtain described tolerance sequence
The step of maximal metric value Power of row P (k) including:
Increase frequency at the two ends of described H1 sequence, obtain H2 (k) sequence, wherein, described H2 (k)
The a length of M of sequence；
Described H2 (k) sequence is changed to time domain, obtains time domain h1 (k) sequence, wherein, k=0 ..., M-1；
According to described time domain h1 (k) retrieval metric sequence P (k), and obtain described metric sequence P (k)
Maximal metric value Power.
Preferably, the described step according to described time domain h1 (k) retrieval metric sequence P (k) includes:
Search sequence h2 (k) according to described time domain h1 (k) retrieval subscriber equipment；
Using the mould square of described search sequence h2 (k) as described metric sequence P (k).
Additionally, for achieving the above object, the present invention also provides for the processing means of a kind of timing slip, described
The processing means of timing slip includes:
First acquisition module, is used for obtaining H1 sequence, and wherein, described H1 sequence is that base station receives
The sequence that subscriber equipment frequency-region signal on predetermined channel obtains with reference signal sequence conjugate multiplication；
Second acquisition module, is used for according to described H1 retrieval metric sequence P (k), and obtains described
Maximal metric value Power of metric sequence P (k)；
First computing module, for obtaining described metric sequence P (k) according to described maximal metric value Power
Signal subset S (x1) and noise subsets NI (x2), calculate the characteristic measure value of described noise subsets NI (x2)
NIPower；
Second computing module, for according to described maximal metric value Power, characteristic measure value NIPower,
Noise jamming threshold T 1 and search thresholding threshold T 2 calculate first path position value；
3rd computing module, for calculating the timing slip of described subscriber equipment according to described first path position value
Amount.
Preferably, described first computing module includes:
First acquiring unit, is used for obtaining the long w1 of left window of described metric sequence P (k), the long w2 of right window
And minimum index value MaxPidx that described maximal metric value Power is corresponding, according to the described long w1 of left window,
The long w2 of right window and minimum index value MaxPidx obtain signal subset S (x1)；
Second acquisition unit, obtains described noise subsets for obtaining according to described signal subset S (x1)
NI(x2)；
Computing unit, is used for obtaining described in the maximal metric value of described noise subsets NI (x2) or calculating and makes an uproar
The average degree value of phonon collection NI (x2) is as described characteristic measure value NIPower.
Preferably, described second computing module includes:
Judging unit, for judging the ratio of described maximal metric value Power and characteristic measure value NIPower
Whether value is less than described noise jamming threshold T 1；
First processing unit, for the most then using described minimum index value MaxPidx as described first footpath
Positional value；
Second processing unit, for if it is not, then obtain in described signal subset S (x1) more than Power/T2
Corresponding minimum index value T2idx, obtains described first path position value according to described minimum index value T2idx.
Preferably, described second acquisition module includes:
Increase unit, for increasing frequency at the two ends of described H1 sequence, obtain H2 (k) sequence,
Wherein, a length of M of described H2 (k) sequence；
Converting unit, for changing described H2 (k) sequence to time domain, obtains time domain h1 (k) sequence
Row, wherein, k=0 ..., M-1；
3rd acquiring unit, is used for according to described time domain h1 (k) retrieval metric sequence P (k), and
Obtain maximal metric value Power of described metric sequence P (k).
Preferably, described 3rd acquiring unit is specifically for using according to described time domain h1 (k) retrieval
Search sequence h2 (k) of family equipment；Using the mould square of described search sequence h2 (k) as described metric sequence
P(k)。
The processing method and processing device of a kind of timing slip of the present invention, by subscriber equipment on predetermined channel
Frequency-region signal with the H1 sequence that reference signal sequence conjugate multiplication obtains to obtain corresponding metric sequence
P (k), maximal metric value Power of metric sequence P (k), noise subsets NI (x2) of metric sequence P (k)
Characteristic measure value NIPower, by maximal metric value Power and characteristic measure value NIPower and noise
Interference threshold value T1 or search thresholding threshold T 2 are compared to calculate first path position value, according to first footpath
Positional value calculates the timing offset of subscriber equipment, and the present invention can be calculated positive timing offset,
Can be calculated negative timing offset, and use the calculated timing of method institute of the present invention inclined
Shifting amount accuracy is higher.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of processing method one embodiment of timing slip of the present invention；
Fig. 2 is the refinement schematic flow sheet of step S102 in Fig. 1；
Fig. 3 is the refinement schematic flow sheet of step S103 in Fig. 1；
Fig. 4 is the refinement schematic flow sheet of step S104 in Fig. 1；
Fig. 5 is the high-level schematic functional block diagram of processing means one embodiment of timing slip of the present invention；
Fig. 6 is the refinement high-level schematic functional block diagram of the second acquisition module in Fig. 5；
Fig. 7 is the refinement high-level schematic functional block diagram of the first computing module in Fig. 5；
Fig. 8 is the refinement high-level schematic functional block diagram of the second computing module in Fig. 5.
The realization of the object of the invention, functional characteristics and advantage will in conjunction with the embodiments, do referring to the drawings further
Explanation.
Detailed description of the invention
Should be appreciated that specific embodiment described herein, only in order to explain the present invention, is not used to limit
Determine the present invention.
The present invention provides the processing method of a kind of timing slip, and with reference to Fig. 1, in one embodiment, this is fixed
The processing method of hour offset includes:
Step S101, obtains H1 sequence；
In the present embodiment, H1 sequence is the subscriber equipment that receives of base station frequency-region signal on predetermined channel
The sequence obtained with reference signal sequence conjugate multiplication.
In the present embodiment, predetermined channel is Physical Random Access Channel or Physical Uplink Shared Channel or detection
Reference channel.Reference signal sequence is preferably Zadoff-Chu sequence.
Step S102, according to described H1 retrieval metric sequence P (k), and obtains described metric sequence
Maximal metric value Power of P (k)；
In the present embodiment, the two ends of H1 sequence are processed, increase corresponding frequency at two ends, make
Must increase a length of M of the sequence after frequency, the sequence after then increasing frequency is transformed into time domain, so
Afterwards according to the configuration of system, by being transformed into the sequence of time domain to determine search sequence, in the present embodiment, logical
Cross search sequence and obtain metric sequence, such as, can obtain tolerance system by the modulus value calculating search sequence
Row, or power sequence can be selected as metric sequence: i.e. next by calculating the mould square of search sequence
Obtain tolerance series.
In the present embodiment, in maximal metric value Power obtaining metric sequence P (k), metric sequence P (k)
Maximal metric value Power may be to there being one or more index value.By searching out metric sequence
Maximal metric value Power of P (k), can carry out drawing of signal and noise by the signal that subscriber equipment sends
Point.
Step S103, obtains the signal of described metric sequence P (k) according to described maximal metric value Power
Subset S (x1) and noise subsets NI (x2), calculate the characteristic measure value of described noise subsets NI (x2)
NIPower；
In the present embodiment, the position that metric sequence P (k) moderate value maximum sign signal is the strongest, therefore,
Metric sequence can be determined according to index value corresponding at maximal metric value Power in metric sequence P (k)
Signal subset S (x1) of P (k) and noise subsets NI (x2).
In the present embodiment, or can shift obtain left window onto according to practical situation by the way of emulation long
And right window is long, according to corresponding at maximal metric value Power in left window length, right window length and metric sequence P (k)
Minimum index value determine signal subset S (x1), then except signal subset S (x1) in metric sequence P (k)
Outer sequence takes a part therein as noise subsets NI (x2).
In the present embodiment, in characteristic measure value NIPower calculating noise subsets NI (x2), Ke Yiqu
The maximal metric value of noise subsets NI (x2), or calculate the average degree value conduct of noise subsets NI (x2)
Characteristic measure value NIPower, characteristic measure value NIPower is used for characterizing noise subsets NI (x2) correspondence
The size of noise signal.
Step S104, does according to described maximal metric value Power, characteristic measure value NIPower, noise
Disturb threshold T 1 and search thresholding threshold T 2 calculates first path position value；
In the present embodiment, noise jamming threshold T 1 and search thresholding threshold T 2 are preferably by emulation
Obtain, it is also possible to determine noise jamming threshold T 1 and the search thresholding door of correspondence according to different scenes
Limit value T2, can obtain noise jamming threshold T 1 and search thresholding door when calculating first path position for the first time
Limit value T2, for follow-up use.
In the present embodiment, if the noise power of the signal of subscriber equipment transmission is relatively big, signal power is less,
Then determine first path position value Fidx according to noise jamming threshold T 1；If the letter that subscriber equipment sends
Number noise power is less and signal power is relatively big, then determine first footpath by search thresholding threshold T 2
Positional value Fidx.
Step S105, calculates the timing offset of described subscriber equipment according to described first path position value.
In the present embodiment, if first path position value Fidx < M/2, then timing offset TAest=
Fidx/M*2048；If first path position value Fidx >=M/2, then timing offset TAest=
(Fidx-M)/M*2048.Thus, it is possible to find out, the present embodiment i.e. can be calculated positive timing slip
Amount, can be calculated negative timing offset, and uses the method institute of the present embodiment calculated
Timing offset accuracy is higher.
Wherein, the unit of timing offset TAest is Ts, the 1Ts=0.001/30720 second.According to reality
The requirement of applicable cases, can use for timing offset TAest and suitable round mode, for example with
The mode rounded up rounds.
Compared with prior art, the present embodiment passes through subscriber equipment frequency-region signal on predetermined channel and ginseng
Examine the H1 sequence that signal sequence conjugate multiplication obtains and obtain metric sequence P (k) of correspondence, metric sequence
Maximal metric value Power of P (k), the characteristic measure value of noise subsets NI (x2) of metric sequence P (k)
NIPower, by maximal metric value Power and characteristic measure value NIPower and noise jamming threshold T 1
Or search thresholding threshold T 2 is compared to calculate first path position value, calculates according to first path position value and use
The timing offset of family equipment, the present embodiment can be calculated positive timing offset, can calculate
To negative timing offset, and use the calculated timing offset of method institute of the present embodiment accurate
Spend higher.
In a preferred embodiment, as in figure 2 it is shown, on the basis of the embodiment of above-mentioned Fig. 1, on
State step S102 to include:
Step S1021, increases frequency, obtains H2 (k) sequence at the two ends of described H1 sequence, wherein,
The a length of M of described H2 (k) sequence；
Step S1022, changes described H2 (k) sequence to time domain, obtains time domain h1 (k) sequence,
Wherein, k=0 ..., M-1；
Step S1023, according to described time domain h1 (k) retrieval metric sequence P (k), and obtains institute
State maximal metric value Power of metric sequence P (k).
In the present embodiment, H1 sequence two ends are processed, obtains H2 (k) sequence.Wherein, H1 sequence
It is that the subscriber equipment received base station frequency-region signal on predetermined channel is conjugated with Zadoff-Chu sequence
The result being multiplied, k=0 ..., M-1, M are the length of H2 (k) sequence.
In the present embodiment, the process to H1 sequence is mainly counted at two ends increase frequency domain, the point of increase
Value, preferably 0, it would however also be possible to employ other modes fill out value.Believe for different channels and subscriber equipment
Number sub-carrier number, the frequency domain increasing H1 sequence is counted also different: on the one hand the frequency domain of increase counts
It is the complexity considering to realize, the on the other hand resolution of timing offset to be considered.
In the present embodiment, for Physical Random Access Channel, according to LTE and LTE-A agreement, user
The sub-carrier number that the signal that equipment is sent by Physical Random Access Channel takies is 864 or 144, corresponding
Zadoff-Chu sequence length should be 839 or 139 mutually, i.e. the H1 sequence of Physical Random Access Channel is long
Degree is 839 or 139, and the length of H2 (k) sequence obtained after process preferably can take 864 or 144,
Suitable length can also be taken according to other rules.
For Physical Uplink Shared Channel and detection reference channel, need H1 sequence two ends are processed,
Make M >=N/16, wherein, N is that FFT corresponding to system bandwidth counts.Processing mode is similar to PRACH,
Being all to increase frequency at H1 sequence two ends, preferably increasing value a little is 0, or fills out other values.Increase
Count when being 0, H1=H2 (k).
In the present embodiment, H2 (k) is transformed into time domain, obtains h1 (k), k=0 ..., M-1.Frequency domain is then
The conversion in territory can have various ways, preferably inverse discrete Fourier transform (Inverse Discrete
Fourier Transform, IDFT).Configure according to system, h1 (k) obtain the search sequence of subscriber equipment
h2(k).Preferably, the code division number of users of multiplexing and code division ginseng on the frequency domain position taken according to subscriber equipment
Number determines search sequence h2 (k).If this frequency domain position is except this subscriber equipment, other users are also had to set
Standby, then the value in h1 (k) corresponding for other subscriber equipmenies is reset, obtain h2 (k).If this frequency domain position
Put only this subscriber equipment, there is no other subscriber equipmenies, then h2 (k)=h1 (k).
Wherein, h2 (k) obtain metric sequence P (k): metric sequence P (k) preferably uses power sequence,
Now P (k)=h2 (k) .* [h2 (k)] H, wherein [h2 (k)] H is the conjugation of h2 (k).Metric sequence P (k) also may be used
Calculate in other ways to adopt, such as the amplitude etc. for h2 (k).
In the present embodiment, the point of the maximum of metric sequence P (k) may have multiple, metric sequence P (k)
In the institute a little [Power, MaxPidx] of maximum, maximal metric value Power=P (MaxPidx), take wherein
Minimum index value is as minimum index value MaxPidx of the present embodiment.
In a preferred embodiment, as it is shown on figure 3, on the basis of the embodiment of above-mentioned Fig. 1, on
State step S103 to include:
Step S1031, obtain the long w1 of left window of described metric sequence P (k), the long w2 of right window and described
Big minimum index value MaxPidx the most corresponding for metric Power is long according to the described long w1 of left window, right window
W2 and minimum index value MaxPidx obtain signal subset S (x1)；
Step S1032, obtains described noise subsets NI (x2) according to described signal subset S (x1)；
Step S1033, obtains the maximal metric value of described noise subsets NI (x2) or calculates described noise
The average degree value of subset NI (x2) is as described characteristic measure value NIPower.
In the present embodiment, obtain the long w1 of left window, the long w2 of right window of metric sequence P (k), can be by imitative
Genuine mode or shift onto obtain the long w1 of left window and the long w2 of right window according to practical situation.Long according to left window
W1, the long w2 of right window and minimum index value MaxPidx obtain signal subset S (x1), wherein, x1=0 ...,
w1-1.Then sequence in addition to signal subset S (x1) takes in metric sequence P (k) part therein make
For noise subsets NI (x2), x2=0 ..., ni-1, wherein, ni is the length of NI sequence.Asked by NI (x2)
Obtain characteristic measure value NIpower.Preferably, characteristic measure value NIpower can take noise subsets NI (x2)
Maximal metric value, or calculate noise subsets NI (x2) average degree value as characteristic measure value
NIPower (this average degree value must be non-zero value).
In the present embodiment, the mode preferably calculated of signal subset S (x1) and noise subsets NI (x2) is as follows:
As MaxPidx >=w1 time:
S (x1)=P (MaxPidx-w1+1:1:MaxPidx)；If MaxPidx+w1≤M, then
NI (x2)=[P (0:1:MaxPidx-w1)；P(MaxPidx+w2:1:M-1)]；If MaxPidx+w1 > M,
Then NI (x2)=P (MaxPidx+w2-M:1:MaxPidx-w1)；
When MaxPidx < during w1:
S (x1)=[P (M-(w1-MaxPidx)+2:1:M-1)；P (0:1:MaxPidx)],
NI (x2)=P (MaxPidx+w2:1:M-(w1-MaxPidx)+1)；
Wherein (a:1:b) represents each index value between index value a to index value b.
In a preferred embodiment, as shown in Figure 4, on the basis of the embodiment of above-mentioned Fig. 3, on
State step S104 to include:
Step S1041, it is judged that described maximal metric value Power and the ratio of characteristic measure value NIPower
Whether less than described noise jamming threshold T 1；The most then enter step S1042, otherwise, enter step
Rapid S1043；
Step S1042, using described minimum index value MaxPidx as described first path position value；
Step S1043, obtains in described signal subset S (x1) more than minimum index corresponding for Power/T2
Value T2idx, obtains described first path position value according to described minimum index value T2idx.
In the present embodiment, if Power/NIpower < T1, then first path position value is Fidx=MaxPidx；
If Power/NIpower >=T1, then obtain in signal subset S (x1) corresponding more than Power/T2
Little index value T2idx, and obtain the S (T2idx) the index TP2idx corresponding to tolerance P (k), first path position
Fidx=TP2idx.
In the present embodiment, noise jamming threshold T 1 and search thresholding threshold T 2 are preferably by emulation
Obtain, it is also possible to determine noise jamming threshold T 1 and the search thresholding door of correspondence according to different scenes
Limit value T2.
The present invention also provides for the processing means of a kind of timing slip, as it is shown in figure 5, described timing slip
Processing means includes:
First acquisition module 101, is used for obtaining H1 sequence, and wherein, described H1 sequence is that base station receives
To the sequence that obtains with reference signal sequence conjugate multiplication of subscriber equipment frequency-region signal on predetermined channel
Row；
In the present embodiment, H1 sequence is the subscriber equipment that receives of base station frequency-region signal on predetermined channel
The sequence obtained with reference signal sequence conjugate multiplication.
In the present embodiment, predetermined channel is Physical Random Access Channel or Physical Uplink Shared Channel or detection
Reference channel.Reference signal sequence is preferably Zadoff-Chu sequence.
Second acquisition module 102, is used for according to described H1 retrieval metric sequence P (k), and obtains institute
State maximal metric value Power of metric sequence P (k)；
In the present embodiment, the two ends of H1 sequence are processed, increase corresponding frequency at two ends, make
Must increase a length of M of the sequence after frequency, the sequence after then increasing frequency is transformed into time domain, so
Afterwards according to the configuration of system, by being transformed into the sequence of time domain to determine search sequence, in the present embodiment, logical
Cross search sequence and obtain metric sequence, such as, can obtain tolerance system by the modulus value calculating search sequence
Row, or power sequence can be selected as metric sequence: i.e. next by calculating the mould square of search sequence
Obtain tolerance series.
In the present embodiment, in maximal metric value Power obtaining metric sequence P (k), metric sequence P (k)
Maximal metric value Power may be to there being one or more index value.By searching out metric sequence
Maximal metric value Power of P (k), can carry out drawing of signal and noise by the signal that subscriber equipment sends
Point.
First computing module 103, for obtaining described metric sequence according to described maximal metric value Power
Signal subset S (x1) of P (k) and noise subsets NI (x2), calculate the feature of described noise subsets NI (x2)
Metric NIPower；
In the present embodiment, the position that metric sequence P (k) moderate value maximum sign signal is the strongest, therefore,
Metric sequence can be determined according to index value corresponding at maximal metric value Power in metric sequence P (k)
Signal subset S (x1) of P (k) and noise subsets NI (x2).
In the present embodiment, or can shift obtain left window onto according to practical situation by the way of emulation long
And right window is long, according to corresponding at maximal metric value Power in left window length, right window length and metric sequence P (k)
Minimum index value determine signal subset S (x1), then except signal subset S (x1) in metric sequence P (k)
Outer sequence takes a part therein as noise subsets NI (x2).
In the present embodiment, in characteristic measure value NIPower calculating noise subsets NI (x2), Ke Yiqu
The maximal metric value of noise subsets NI (x2), or calculate the average degree value conduct of noise subsets NI (x2)
Characteristic measure value NIPower, characteristic measure value NIPower is used for characterizing noise subsets NI (x2) correspondence
The size of noise signal.
Second computing module 104, for according to described maximal metric value Power, characteristic measure value
NIPower, noise jamming threshold T 1 and search thresholding threshold T 2 calculate first path position value；
In the present embodiment, noise jamming threshold T 1 and search thresholding threshold T 2 are preferably by emulation
Obtain, it is also possible to determine noise jamming threshold T 1 and the search thresholding door of correspondence according to different scenes
Limit value T2.
In the present embodiment, if the noise power of the signal of subscriber equipment transmission is relatively big, signal power is less,
Then determine first path position value Fidx according to noise jamming threshold T 1；If the letter that subscriber equipment sends
Number noise power is less and signal power is relatively big, then determine first footpath by search thresholding threshold T 2
Positional value Fidx.
3rd computing module 105, for calculating the timing of described subscriber equipment according to described first path position value
Side-play amount.
In the present embodiment, if first path position value Fidx < M/2, then timing offset TAest=
Fidx/M*2048；If first path position value Fidx >=M/2, then timing offset TAest=
(Fidx-M)/M*2048.Thus, it is possible to find out, the present embodiment i.e. can be calculated positive timing slip
Amount, can be calculated negative timing offset, and uses the method institute of the present embodiment calculated
Timing offset accuracy is higher.
Wherein, the unit of timing offset TAest is Ts, the 1Ts=0.001/30720 second.According to reality
The requirement of applicable cases, can use for timing offset TAest and suitable round mode, for example with
The mode rounded up rounds.
In a preferred embodiment, as shown in Figure 6, on the basis of the embodiment of above-mentioned Fig. 5, the
Two acquisition modules 102 include:
Increase unit 1021, for increasing frequency at the two ends of described H1 sequence, obtain H2 (k) sequence
Row, wherein, a length of M of described H2 (k) sequence；
Converting unit 1022, for changing described H2 (k) sequence to time domain, obtains time domain h1 (k)
Sequence, wherein, k=0 ..., M-1；
3rd acquiring unit 1023, is used for according to described time domain h1 (k) retrieval metric sequence P (k),
And obtain maximal metric value Power of described metric sequence P (k).
In the present embodiment, H1 sequence two ends are processed, obtains H2 (k) sequence.Wherein, H1 sequence
It is that the subscriber equipment received base station frequency-region signal on predetermined channel is conjugated with Zadoff-Chu sequence
The result being multiplied, k=0 ..., M-1, M are the length of H2 (k) sequence.
In the present embodiment, the process to H1 sequence is mainly counted at two ends increase frequency domain, the point of increase
Value, preferably 0, it would however also be possible to employ other modes fill out value.Believe for different channels and subscriber equipment
Number sub-carrier number, the frequency domain increasing H1 sequence is counted also different: on the one hand the frequency domain of increase counts
It is the complexity considering to realize, the on the other hand resolution of timing offset to be considered.
In the present embodiment, for Physical Random Access Channel, according to LTE and LTE-A agreement, user
The sub-carrier number that the signal that equipment is sent by Physical Random Access Channel takies is 864 or 144, corresponding
Zadoff-Chu sequence length should be 839 or 139 mutually, i.e. the H1 sequence of Physical Random Access Channel is long
Degree is 839 or 139, and the length of H2 (k) sequence obtained after process preferably can take 864 or 144,
Suitable length can also be taken according to other rules.
For Physical Uplink Shared Channel and detection reference channel, need H1 sequence two ends are processed,
Make M >=N/16, wherein, N is that FFT corresponding to system bandwidth counts.Processing mode is similar to PRACH,
Being all to increase frequency at H1 sequence two ends, preferably increasing value a little is 0, or fills out other values.Increase
Count when being 0, H1=H2 (k).
In the present embodiment, H2 (k) is transformed into time domain, obtains h1 (k), k=0 ..., M-1.Frequency domain is then
The conversion in territory can have various ways, preferably inverse discrete Fourier transform (Inverse Discrete
Fourier Transform, IDFT).Configure according to system, h1 (k) obtain the search sequence of subscriber equipment
h2(k).Preferably, the code division number of users of multiplexing and code division ginseng on the frequency domain position taken according to subscriber equipment
Number determines search sequence h2 (k).If this frequency domain position is except this subscriber equipment, other users are also had to set
Standby, then the value in h1 (k) corresponding for other subscriber equipmenies is reset, obtain h2 (k).If this frequency domain position
Put only this subscriber equipment, there is no other subscriber equipmenies, then h2 (k)=h1 (k).
Wherein, h2 (k) obtain metric sequence P (k): metric sequence P (k) preferably uses power sequence,
Now P (k)=h2 (k) .* [h2 (k)] H, wherein [h2 (k)] H is the conjugation of h2 (k).Metric sequence P (k) also may be used
Calculate in other ways to adopt, such as the amplitude etc. for h2 (k).
In the present embodiment, the point of the maximum of metric sequence P (k) may have multiple, metric sequence P (k)
In the institute a little [Power, MaxPidx] of maximum, maximal metric value Power=P (MaxPidx), take wherein
Minimum index value is as minimum index value MaxPidx of the present embodiment.
In a preferred embodiment, as it is shown in fig. 7, on the basis of the embodiment of above-mentioned Fig. 5, institute
State the first computing module 103 to include:
First acquiring unit 1031, is used for obtaining the long w1 of left window of described metric sequence P (k), right window long
W2 and minimum index value MaxPidx corresponding to described maximal metric value Power are long according to described left window
W1, the long w2 of right window and minimum index value MaxPidx obtain signal subset S (x1)；
Second acquisition unit 1032, obtains described noise for obtaining according to described signal subset S (x1)
Collection NI (x2)；
Computing unit 1033, is used for obtaining the maximal metric value of described noise subsets NI (x2) or calculating institute
State the average degree value of noise subsets NI (x2) as described characteristic measure value NIPower.
In the present embodiment, obtain the long w1 of left window, the long w2 of right window of metric sequence P (k), can be by imitative
Genuine mode or shift onto obtain the long w1 of left window and the long w2 of right window according to practical situation.Long according to left window
W1, the long w2 of right window and minimum index value MaxPidx obtain signal subset S (x1), wherein, x1=0 ...,
w1-1.Then sequence in addition to signal subset S (x1) takes in metric sequence P (k) part therein make
For noise subsets NI (x2), x2=0 ..., ni-1, wherein, ni is the length of NI sequence.Asked by NI (x2)
Obtain characteristic measure value NIpower.Preferably, characteristic measure value NIpower can take noise subsets NI (x2)
Maximal metric value, or calculate noise subsets NI (x2) average degree value as characteristic measure value
NIPower (this average degree value must be non-zero value).
In the present embodiment, the mode preferably calculated of signal subset S (x1) and noise subsets NI (x2) is as follows:
As MaxPidx >=w1 time:
S (x1)=P (MaxPidx-w1+1:1:MaxPidx)；If MaxPidx+w1≤M, then
NI (x2)=[P (0:1:MaxPidx-w1)；P(MaxPidx+w2:1:M-1)]；If MaxPidx+w1 > M,
Then NI (x2)=P (MaxPidx+w2-M:1:MaxPidx-w1)；
When MaxPidx < during w1:
S (x1)=[P (M-(w1-MaxPidx)+2:1:M-1)；P (0:1:MaxPidx)],
NI (x2)=P (MaxPidx+w2:1:M-(w1-MaxPidx)+1)；
Wherein (a:1:b) represents each index value between index value a to index value b.
In a preferred embodiment, as shown in Figure 8, on the basis of the embodiment of above-mentioned Fig. 7, institute
State the second computing module 104 to include:
Judging unit 1041, is used for judging described maximal metric value Power and characteristic measure value NIPower
Ratio whether less than described noise jamming threshold T 1；
First processing unit 1042, for the most then using described minimum index value MaxPidx as described
First path position value；
Second processing unit 1043, for if it is not, then obtain in described signal subset S (x1) more than Power/
Minimum index value T2idx corresponding for T2, obtains described first path position according to described minimum index value T2idx
Value.
In the present embodiment, if Power/NIpower < T1, then first path position value is Fidx=MaxPidx；
If Power/NIpower >=T1, then obtain in signal subset S (x1) corresponding more than Power/T2
Little index value T2idx, and obtain the S (T2idx) the index TP2idx corresponding to tolerance P (k), first path position
Fidx=TP2idx.
In the present embodiment, noise jamming threshold T 1 and search thresholding threshold T 2 are preferably by emulation
Obtain, it is also possible to determine noise jamming threshold T 1 and the search thresholding door of correspondence according to different scenes
Limit value T2.
These are only the preferred embodiments of the present invention, not thereby limit the scope of the claims of the present invention, every
Utilize equivalent structure or equivalence flow process conversion that description of the invention and accompanying drawing content made, or directly or
Connect and be used in other relevant technical fields, be the most in like manner included in the scope of patent protection of the present invention.
Claims (10)
1. the processing method of a timing slip, it is characterised in that the processing method bag of described timing slip
Include following steps:
Obtaining H1 sequence, wherein, described H1 sequence is that the subscriber equipment that receives of base station is at predetermined channel
On the sequence that obtains with reference signal sequence conjugate multiplication of frequency-region signal；
According to described H1 retrieval metric sequence P (k), and obtain the maximum of described metric sequence P (k)
Metric Power；
According to described maximal metric value Power obtain described metric sequence P (k) signal subset S (x1) and
Noise subsets NI (x2), calculates characteristic measure value NIPower of described noise subsets NI (x2)；
According to described maximal metric value Power, characteristic measure value NIPower, noise jamming threshold T 1
And search thresholding threshold T 2 calculates first path position value；
The timing offset of described subscriber equipment is calculated according to described first path position value.
2. the processing method of timing slip as claimed in claim 1, it is characterised in that described according to institute
State maximal metric value Power and obtain signal subset S (x1) and the noise subsets of described metric sequence P (k)
NI (x2), the step of characteristic measure value NIPower calculating described noise subsets NI (x2) includes:
Obtain the long w1 of left window of described metric sequence P (k), the long w2 of right window and described maximal metric value
Minimum index value MaxPidx corresponding for Power, according to the described long w1 of left window, the long w2 of right window and minimum
Index value MaxPidx obtains signal subset S (x1)；
Described noise subsets NI (x2) is obtained according to described signal subset S (x1)；
Obtain the maximal metric value of described noise subsets NI (x2) or calculate described noise subsets NI (x2)
Average degree value is as described characteristic measure value NIPower.
3. the processing method of timing slip as claimed in claim 2, it is characterised in that described according to institute
State maximal metric value Power, characteristic measure value NIPower, noise jamming threshold T 1 and search thresholding
Threshold T 2 calculates the step of first path position value and includes:
Judge that whether described maximal metric value Power and the ratio of characteristic measure value NIPower are less than described
Noise jamming threshold T 1；
The most then using described minimum index value MaxPidx as described first path position value；
If it is not, then obtain in described signal subset S (x1) more than minimum index value corresponding for Power/T2
T2idx, obtains described first path position value according to described minimum index value T2idx.
4. the processing method of timing slip as claimed in claim 1, it is characterised in that described according to institute
State H1 retrieval metric sequence P (k), and obtain maximal metric value Power of described metric sequence P (k)
Step include:
Increase frequency at the two ends of described H1 sequence, obtain H2 (k) sequence, wherein, described H2 (k)
The a length of M of sequence；
Described H2 (k) sequence is changed to time domain, obtains time domain h1 (k) sequence, wherein, k=0 ..., M-1；
According to described time domain h1 (k) retrieval metric sequence P (k), and obtain described metric sequence P (k)
Maximal metric value Power.
5. the processing method of timing slip as claimed in claim 4, it is characterised in that described according to institute
The step stating time domain h1 (k) retrieval metric sequence P (k) includes:
Search sequence h2 (k) according to described time domain h1 (k) retrieval subscriber equipment；
Using the mould square of described search sequence h2 (k) as described metric sequence P (k).
6. the processing means of a timing slip, it is characterised in that the processing means bag of described timing slip
Include:
First acquisition module, is used for obtaining H1 sequence, and wherein, described H1 sequence is that base station receives
The sequence that subscriber equipment frequency-region signal on predetermined channel obtains with reference signal sequence conjugate multiplication；
Second acquisition module, is used for according to described H1 retrieval metric sequence P (k), and obtains described
Maximal metric value Power of metric sequence P (k)；
First computing module, for obtaining described metric sequence P (k) according to described maximal metric value Power
Signal subset S (x1) and noise subsets NI (x2), calculate the characteristic measure value of described noise subsets NI (x2)
NIPower；
Second computing module, for according to described maximal metric value Power, characteristic measure value NIPower,
Noise jamming threshold T 1 and search thresholding threshold T 2 calculate first path position value；
3rd computing module, for calculating the timing slip of described subscriber equipment according to described first path position value
Amount.
7. the processing means of timing slip as claimed in claim 6, it is characterised in that described first meter
Calculation module includes:
First acquiring unit, is used for obtaining the long w1 of left window of described metric sequence P (k), the long w2 of right window
And minimum index value MaxPidx that described maximal metric value Power is corresponding, according to the described long w1 of left window,
The long w2 of right window and minimum index value MaxPidx obtain signal subset S (x1)；
Second acquisition unit, obtains described noise subsets for obtaining according to described signal subset S (x1)
NI(x2)；
Computing unit, is used for obtaining described in the maximal metric value of described noise subsets NI (x2) or calculating and makes an uproar
The average degree value of phonon collection NI (x2) is as described characteristic measure value NIPower.
8. the processing means of timing slip as claimed in claim 7, it is characterised in that described second meter
Calculation module includes:
Judging unit, for judging the ratio of described maximal metric value Power and characteristic measure value NIPower
Whether value is less than described noise jamming threshold T 1；
First processing unit, for the most then using described minimum index value MaxPidx as described first footpath
Positional value；
Second processing unit, for if it is not, then obtain in described signal subset S (x1) more than Power/T2
Corresponding minimum index value T2idx, obtains described first path position value according to described minimum index value T2idx.
9. the processing means of timing slip as claimed in claim 6, it is characterised in that described second obtains
Delivery block includes:
Increase unit, for increasing frequency at the two ends of described H1 sequence, obtain H2 (k) sequence,
Wherein, a length of M of described H2 (k) sequence；
Converting unit, for changing described H2 (k) sequence to time domain, obtains time domain h1 (k) sequence
Row, wherein, k=0 ..., M-1；
3rd acquiring unit, is used for according to described time domain h1 (k) retrieval metric sequence P (k), and
Obtain maximal metric value Power of described metric sequence P (k).
10. the processing means of timing slip as claimed in claim 9, it is characterised in that the described 3rd
Acquiring unit is specifically for search sequence h2 (k) according to described time domain h1 (k) retrieval subscriber equipment；
Using the mould square of described search sequence h2 (k) as described metric sequence P (k).
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510264229.0A CN106304313A (en) | 2015-05-21 | 2015-05-21 | The processing method and processing device of timing slip |
| PCT/CN2016/071364 WO2016184155A1 (en) | 2015-05-21 | 2016-01-19 | Method of processing timing offset and device utilizing same |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510264229.0A CN106304313A (en) | 2015-05-21 | 2015-05-21 | The processing method and processing device of timing slip |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106304313A true CN106304313A (en) | 2017-01-04 |
Family
ID=57319372
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510264229.0A Withdrawn CN106304313A (en) | 2015-05-21 | 2015-05-21 | The processing method and processing device of timing slip |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN106304313A (en) |
| WO (1) | WO2016184155A1 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114884538B (en) * | 2022-07-11 | 2022-10-25 | 四川创智联恒科技有限公司 | Time offset estimation method and system for NR uplink control channel |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070002728A1 (en) * | 2005-06-15 | 2007-01-04 | Ntt Docomo, Inc. | Radio receiver and radio signal receiving method |
| CN101997804A (en) * | 2009-08-21 | 2011-03-30 | 大唐移动通信设备有限公司 | Method and device for estimating synchronous timing deviation |
| CN102104951A (en) * | 2009-12-17 | 2011-06-22 | 中兴通讯股份有限公司 | Method and device for timing shift estimation |
| CN102196553A (en) * | 2011-03-31 | 2011-09-21 | 重庆重邮信科通信技术有限公司 | Timing synchronization correcting method of mobile terminal |
| US20110279658A1 (en) * | 2010-05-11 | 2011-11-17 | Freescale Semiconductor, Inc | Method for synchronizing remote device |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101133582B (en) * | 2005-07-06 | 2011-02-02 | 中兴通讯股份有限公司 | Ascending synchronization establishing method based on associated transmission technology and transmission matrix thereof |
| CN101185270B (en) * | 2005-11-16 | 2010-10-13 | 中兴通讯股份有限公司 | Uplink synchronization method in mobile communication system |
| CN103298136B (en) * | 2012-02-29 | 2016-11-23 | 华为技术有限公司 | A kind of accidental access method, terminal, base station and system |
| CN104254100B (en) * | 2013-06-25 | 2017-12-15 | 普天信息技术研究院有限公司 | A kind of measuring method of uplink timing advance |
| CN104159270B (en) * | 2014-08-08 | 2018-05-01 | 大唐移动通信设备有限公司 | Control the method and apparatus of user equipment access |
-
2015
- 2015-05-21 CN CN201510264229.0A patent/CN106304313A/en not_active Withdrawn
-
2016
- 2016-01-19 WO PCT/CN2016/071364 patent/WO2016184155A1/en active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070002728A1 (en) * | 2005-06-15 | 2007-01-04 | Ntt Docomo, Inc. | Radio receiver and radio signal receiving method |
| CN101997804A (en) * | 2009-08-21 | 2011-03-30 | 大唐移动通信设备有限公司 | Method and device for estimating synchronous timing deviation |
| CN102104951A (en) * | 2009-12-17 | 2011-06-22 | 中兴通讯股份有限公司 | Method and device for timing shift estimation |
| US20110279658A1 (en) * | 2010-05-11 | 2011-11-17 | Freescale Semiconductor, Inc | Method for synchronizing remote device |
| CN102196553A (en) * | 2011-03-31 | 2011-09-21 | 重庆重邮信科通信技术有限公司 | Timing synchronization correcting method of mobile terminal |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2016184155A1 (en) | 2016-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20190239249A1 (en) | Processing of random access preamble sequences | |
| CN101295999B (en) | Detection method for accidental precursor access | |
| KR101817290B1 (en) | Method and apparatus for detecting physical random access channel preambles in a long term evolution communication system | |
| CN109076609B (en) | Method and receiver for preamble detection and time of arrival estimation | |
| US9615386B2 (en) | Random access channel false alarm control | |
| US7564909B2 (en) | Apparatus and method for detecting ranging signal in an orthogonal frequency division multiple access mobile communication system | |
| KR102248486B1 (en) | Apparatus and method for frequency offset estimation of received signal | |
| US9756657B2 (en) | Wireless network signal to interference plus noise ratio estimation for a random access channel | |
| Sanguinetti et al. | A random access algorithm for LTE systems | |
| US10834762B2 (en) | Random access preamble receiver | |
| CN112585874B (en) | Noise floor estimation for signal detection | |
| CN106792778A (en) | The method and device of uplink SRS signal power in a kind of measurement LTE system | |
| CN105072060A (en) | DTX determination method and apparatus of physical uplink control channel, and base station | |
| CN102238123B (en) | Ranging code detecting method and device | |
| US9503996B2 (en) | Method and apparatus for cell search and synchronization in mobile communication | |
| US8923376B2 (en) | Apparatus and method for estimating noise and interference power in communication system | |
| CN106304313A (en) | The processing method and processing device of timing slip | |
| US20150071258A1 (en) | Random access algorithm for lte systems | |
| US10212679B1 (en) | Method and apparatus for delay spread estimation | |
| CN102695183B (en) | Signature detection and time bias estimation | |
| CN106105126B (en) | Frame transmission, frame detection method, sending ending equipment and receiving device | |
| CN107566310A (en) | Electronic equipment and targeting signal generation, sending, receiving method and device | |
| Karthik et al. | Point Sum Average Peak Algorithm Detection of LTE Preamble | |
| Lee et al. | A Design of Low Latency Random Access Preamble Detector for LTE Uplink Receiver |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20170104 |
|
| WW01 | Invention patent application withdrawn after publication |