At least one of the present embodiments generally relates to a method and a device for an uplink transmission delay estimation between a user equipment and a base station in a wireless communication network.
In the last years, the demand for broadband high-speed, reliable and low latency services and applications, e.g. Internet Of Things devices, increased. Thus, new standards and technologies, known as 5G, have emerged to enable all these new services and applications.
For years, terrestrial and satellite networks have evolved almost independently from each other. The next generation of satellites will integrate with terrestrial 5G networks to manage connectivity to cars, vessels, airplanes and other IoT devices in remote and rural areas. 5G will encompass an evolution of today's 4G networks and the addition of a new, globally standardized radio-access technology known as “New Radio” (NR). Satellite communications will thus complement terrestrial 5G networks. To this aim, Low Earth Orbit Satellites (LEO) will play a key part in extending terrestrial 5G networks to remote areas not covered by small cell networks. In low earth orbits (LEO) satellites based mobile communication environments, the high-speed relative motion between the satellite and user equipment (UE) generates a large Doppler frequency shift. The large Doppler frequency shift is partly corrected at the UE side. However, it complicates the handling of the timing uncertainty at the base station side.
5G wireless communication system are based on orthogonal frequency-division multiplexing (OFDM). In an OFDM based system, orthogonal messages have to be sent, thus the major challenge in such a system is to maintain uplink orthogonality among users. Hence both frequency and time synchronization of the transmitted signals from the users are needed. A fine timing estimation has to be performed by the base station when the signals coming from users are detected. To this aim, it is known to use a Physical Random Access Channel (PRACH) method to obtain an initial fine time synchronization, which will allow to inform the user equipment (UE) how to compensate for the round trip delay. This compensation in the user equipment is also known as timing advance. After a successful random access procedure, in fact, the base station and the user equipment should be synchronized. In this way, the subsequent uplink signals could be correctly decoded and would not interfere with other users connected to the network.
PRACH method consists in the transmission of a set of random access preambles in order to allocate different resources to different users. A random access preamble comprises at least one cyclic prefix (CP) and Ns preamble sequences, Ns being an integer equal to or greater than 1.
Due to good zero autocorrelation and cross-correlation performance, Zadoff-Chu (ZC) sequence is generally used as the preamble sequence for timing synchronization in many systems. As an example, ZC sequences are used in current terrestrial 4G and 5G PRACH method.
At S300, a base station receives at least one random access preamble from a UE: y(t)=x(t)*h(t)+w(t), where x(t) is the transmitted random access preamble, h(t) is a transmission filter which includes equipment filters, neglected here, and a transmission delay τ and w(t) is an additive noise, mainly consisting of the thermal noise, and possibly of some interferences. The base station removes the CP(s) from the received random access preamble at S302 to obtain the Ns preamble sequences. It must be noted that this operation is approximate. Indeed, at this step the receiver of the base station has an imprecise knowledge of the time of arrival of the signal.
At S304, a N-FFT (Fast Fourier Transform) is applied on each of the Ns obtained preamble sequences to obtain Ns sequences in the frequency domain, N being the size of the FFT.
At S306, each obtained sequence in the frequency domain is multiplied by a complex conjugate of a corresponding transmitted preamble sequence. The multiplication by the complex conjugate of the transmitted preamble sequence corresponds to a convolution in the time domain. In practice, M subcarriers (M<N) are modulated for the constitution of the PRACH preamble. The other subcarriers are set to zero.
At S308, for each sequence, a N-IFFT is applied to come back in the time domain. In practice, the IFFT could be performed with a size of M, or any size between M and N, without a great impact on performance. Therefore, the result is an estimation of the transmission filter h(t) in the time domain. In a satellite context, in the absence of noise and impairments caused by the Doppler (frequency shift), this filter is roughly equivalent to a simple delay: h(t)=δ(t−τ), where δ is the Dirac pulse and τ is the transmission delay. In practice, h(t)=g(t−τ) where g(t) corresponds to all emission and reception filters, including the ones implied by the IFFT/FFT processes.
Finally, the delay can be estimated at step S310 by calculating the amplitude of each element of the estimated time response, and looking for the maximum value. The estimation in the time domain provides good results. However, the complexity is high because a N-IFFT is applied for each preamble sequence. The time latency of this method is also high which is an issue for time critical applications.
It is thus desirable to find a method for estimating a time delay that is precise and whose complexity and latency are limited.
A method for an uplink transmission delay estimation between a user equipment and a base station in a wireless communication network is disclosed. The method being implemented in the base station comprises:
the method being characterized in that it further comprises:
where l[c,s] is a phase correcting factor and Snorm[s] is a normalized version of S[s];
In an embodiment, determining the delay responsive to the determined phase comprises calculating the delay {circumflex over (τ)} as follows:
where φopt is the determined phase.
In an embodiment, determining the phase responsive to the selected candidate comprises calculating the phase φopt as follows:
φopt=−arg Γ[copt].
where Copt is the selected candidate, arg(x) an operator returning the argument of x.
In an embodiment, selecting the candidate Copt among the C candidates with a smallest associated error calculated with respect to the sums S[s] is determined as follows:
where Snorm[s] is a normalized version of S[s].
In an embodiment, the wireless communication network is a 5G network.
In an embodiment, each obtained sequence in the frequency domain is a Zadoff-Chu sequence or a M-sequence.
In an embodiment, Ns is equal to 1 or 2. In an embodiment, M is equal to 839.
A time synchronization method between a user equipment and a base station is also disclosed. The time synchronization method comprises:
A device of the base station type in communication with a user equipment in a wireless communication network is further disclosed. The device is configured to estimate an uplink transmission delay estimation from the user equipment, the device comprising:
the method being characterized in that it further comprises:
where l[c,s] is a phase correcting factor and Snorm[s] is a normalized version of S[s];
A computer program product is disclosed. It comprises program code instructions that can be loaded in a programmable device, the program code instructions causing implementation of the method according to one the previous embodiments when the program code instructions are run by the programmable device.
A storage medium, possibly non-transitory, is disclosed. It stores a computer program comprising program code instructions, the program code instructions causing implementation of the method according to one the previous embodiments when the program code instructions are read from the storage medium and run by the programmable device.
The characteristics of the invention will emerge more clearly from a reading of the following description of at least one example of embodiment, said description being produced with reference to the accompanying drawings.
The transmission time between the first user equipment 110 and the base station 100 is equal to δ1 and the transmission time between the second user equipment 120 and the base station 100 is equal to δ2. These values need to be estimated in order to synchronize the base station with these equipments and thus avoid interference with other UEs connected to the network. Indeed, the UEs have different transmission delays to the base station, and, therefore, a method to compensate for these differences is required so as to prevent that data transmitted from these different UEs to collide at the base station 100.
To this aim, a random access method may be used such as the one disclosed with respect to
At S100, the user equipment 110 transmits in an uplink a Random Access (RA) preamble to the base station 100. The RA preamble is randomly selected. The RA preamble comprises at least one cyclic prefix (CP) and at least one preamble sequence, e.g. a ZC sequence. However, the present principles are not limited to this specific RA preamble. In particular, sequences other than ZC sequences may be used such as M-sequences.
At S102, the base station estimates a delay, e.g. a round trip transmission delay. In another embodiment, the base station estimates a time difference between an average transmission delay known by the base station and the current transmission delay. From the average transmission delay and the time difference the base station can estimate the exact transmission delay. In an embodiment, the base station knows its distance from the center of the cell to which the user equipment belongs. From this distance, the base station can estimate the average transmission delay.
In another embodiment, the base station estimates a time difference between a previously estimated transmission delay and the current transmission delay. From the previously estimated transmission delay and the time difference the base station can estimate the current transmission delay. Indeed, after the UE has become synchronized, the synchronization may be lost. This can, for example, be the case if the transmission delay changes because for example the UE moves further from or closer to the base station. Therefore, a time difference may be estimated between the previous transmission delay corresponding to the former position of the UE with respect to the base station and the current one corresponding for example to the new position of the UE. The base station may also move which is the case with a LEO satellite. Knowing its displacement, the base station can estimate at least partially the effect of its displacement on the delay and thus determine a time correction to be applied. However, this time correction depends on the user equipment position and is thus not perfect.
At S104, the base station 100 transmits in a downlink a RA response to the user equipment 110. The RA response comprises the value of the transmission delay estimated by the base station. In one embodiment, the RA response indicates the successfully received preamble(s) along with the timing advance (TA) value and uplink resource allocation information to the UE. The timing advance (TA) value corresponds to the estimated transmission delay.
At S106, the UE uses the TA value to adjust its uplink timing. More precisely, the UE utilizes the received TA value to determine a suitable transmit timing.
At S108, an optional contention resolution step may apply. Indeed, the randomly selected RA preamble does not enable unique identification of the UE. Therefore, multiple UEs may try a random access with the same RA preamble sequence on the same RA channel. The contention resolution step S108 helps uniquely identify the UE that has been selected.
At S110, the user equipment 110 and the base station 100 may exchange data.
The method starts at S400. An integer s is set to the value 1. The base station 100 receives at least one random access preamble from a UE 110: y(t)=x(t)*g(t−τ)+w(t). The random access preamble comprises at least one cyclic prefix CP and Ns preamble sequences occupying M subcarriers, M being an integer and Ns being an integer equal or greater than one. In one embodiment, the received random access sequence comprises a single cyclic prefix CP followed by Ns preamble sequences. In a variant, the received random access sequence comprises Ns cyclic prefixes CPs, each CP being followed by a single preamble sequence.
The received random access sequence may comprise additional data.
At S402, the base station 100 removes the at least one CP . In one embodiment, the base station 100 removes the single cyclic prefix CP. In a variant, the base station removes the Ns cyclic prefixes. Thus, Ns preamble sequences, constituted by N samples in the time domain, occupying M subcarriers are obtained. Each one of the Ns preamble sequences is an OFDM block. Indeed, each one of the Ns preamble sequences in the time domain corresponds to a counterpart preamble sequence in the frequency domain. Said otherwise, the Ns time domain preamble sequences received by the base station were obtained, on the transmitter side, from counterpart preamble sequences in the frequency domain by applying a N-FFT (OFDM modulation). The difference N-M corresponds to N-M non-modulated subcarriers, i.e. to zero values at the input of the N-FFT.
At S404, a N-FFT, i.e. a Fast Fourier Transform of size N, is applied on each of the Ns obtained preamble sequences to obtain Ns sequences in the frequency domain, e.g. ZC sequences or M-sequences. Each sequence in the frequency domain is identified by an index p varying from 1 to Ns. N is equal or proportional to a nominal value N′ defined by the standard corresponding to the sent sequence (4G, 5G . . . ). At S406, in the frequency domain, each obtained sequence in the frequency domain, identified by an index p, is multiplied by a complex conjugate of a corresponding transmitted preamble sequence. Thus, for each obtained sequence in the frequency domain of index p, M values Zkp are obtained where k is an index identifying the subcarriers, k thus varying from 1 to M.
In the frequency domain, with a delay τ, after multiplication by the conjugate and omitting the noise and interference due to Doppler, Zkp is obtained for the subcarrier k: Zkp=e−j2πkτ/N. The phase of the obtained sequence Zkp in the frequency domain is thus representative of the transmission delay which is to be estimated.
At S408, a sum S[s] is calculated as follows from the values Zkp and its conjugate(Zkp)*:
where M is the number of subcarriers occupied by the sequence (839 in 4G and 5G), the index s refers to a specific carrier separation, known as “the step” between the subcarriers involved in the sum calculation and the variable ηs corresponds to the noise plus interference for example caused by Doppler. ρs is equal to the amplitude of S[s]. Said otherwise, Snorm[s]=, where Snorm[s] is a normalized version of S[s]. Each sum S [s] after normalization is thus an estimate of
.
Thus, each sum corresponds to an average of a differential demodulation in the frequency domain.
At S410, the value of s is compared with smax. If s<smax, the method continues at S412 and s is increased by 1. Otherwise, the method continues at S414. Thus at S414, K sums S[s] are obtained with sϵ[1; smax]. From each sum S[s], a phase {circumflex over (φ)}s may be estimated, wherein the phase {circumflex over (φ)}s is a multiple of the delay τ. However, each phase {circumflex over (φ)}s presents a phase ambiguity of
which limits the delay range if only an estimation from S[s] is performed with s>1. Indeed, if
represents the real phase to be estimated, the estimation {circumflex over (φ)}s corresponds to
{circumflex over (φ)}s=sψ+noise (Eq. 2)
In the general case, for a given step s, an estimation {circumflex over (φ)}s corresponds to a set of s values:
Thus, from one estimation {circumflex over (φ)}s obtained from a single sum S[s], it is not possible to know the real phase ψ in the case where s is different of 1. There is thus an ambiguity.
To resolve these ambiguities, the different estimations of the phase are combined in step S414 in order to get a unique estimate of the phase and thus of the delay.
At S414, for a given set of steps s, C candidate sets are calculated from the sums S[s]. A given candidate set is composed of L (which is equal to the number of steps) estimations of ψ, each such estimations corresponding to a particular value of
l=l[c,s] is a phase correcting factor that depends on the step s and on the candidate c. l[c, s] is an integer. This phase correcting factor is used to cope with the inherent phase ambiguity
For example, if s∈{1,2,3,4}, i.e. smax=4, C=6 candidate sets are calculated as indicated in the table 1.
For example, if c=4, the L=4 estimations of {circumflex over (ψ)}[4, s] correspond to
For each candidate set, the L estimations are combined into one. It could be possible to directly combine these phases {circumflex over (ψ)}[c, s]. However, it would have the disadvantage that the corresponding noise will no longer be Gaussian. In particular, this would make more problematic to select the best combination with a least square criterion.
Therefore, the base station 100 combines the
instead of the
In the above calculation, the normalized version of S[s] may be used. In this latter case,
The normalized version Snorm[s] of S[s] is obtained by dividing S[s] by its amplitude in order to obtain a complex number of amplitude 1. Instead of directly adding the different Γ[c, s] in order to get a simple estimate Γ[c] for each candidate set c, a sum weighted by W[s] is used. Thus, Γ[c] for each candidate set c is calculated as follows:
Therefore, for each candidate (c=1 to 6 in the previous example), Γ[c] corresponds to a linear combination of the
W[s] is a weighting factor that depends on s and increases with s. The weights W[s] may take the following values according to various embodiments: W[s]=s2, W[s]=s2, W[s]=1.
At the end of S414, C estimated values Γ[c] are obtained. The base station 100 has to select the most reliable estimation among the C estimations Γ[c]. To this aim, a least square criterion may be used, taking benefit that the noise/interference is still Gaussian or almost Gaussian.
According to equations (1) and (2), S[s] is an estimate of e−jsψ. Besides, each Γ[c] is an estimate of e−jψ. Therefore, for each candidate c, the base station calculates the following square error:
Error[c]=Σs|Snorm[s]−Γs[c]|2 where Γs[c] means Γ[c] raised to the power of s.
At S416, the base station 100 selects a best candidate Copt as the one minimizing Error[c]. Minimizing the square error is equivalent to maximizing the reliability of the estimation under the Gaussian assumption.
At S418, a phase φopt is determined from the selected best candidate Copt as follows: φopt=−arg Γ[copt].
At S420, the base station 100 estimates the corresponding delay {circumflex over (τ)} from the determined phase φopt as follows:
The method ends at S422.
Each combination Γ[c] corresponds to a different way to resolve the different phase ambiguities. The reliability, i.e. Error[c], of each combination Γ[c] is evaluated with respect to the observations, i.e. the original sums S[s]. The combination with the best reliability, i.e. the lowest Error[c], is selected in order to provide the final delay estimation. By combining the different estimations, the phase ambiguity is thus resolved and thus a unique estimate of the phase and therefore of the delay is obtained. The method may be used in a 5G satellite system and works with very large Doppler. Thus, the method makes it possible to estimate the delay despite the difficulty caused by Doppler and by the low SNR inherent to satellite transmission.
The method starts at step S500. An integer r is set to the value 0.
At S502, the Least Common Multiple (LCM) β of all the steps s is determined. In the case where the steps are equal to the followings values {1, 2, 3, 4}, β=12. The phase range [0; 2π] is divided into β intervals of size
At S506, a phase correcting factor corresponding to a step s in calculated as follows:
for the interval
For the calculation of l[s] the division by β is an integer division.
At S508, r is compared with β−1. If r<β−1, the method continues at S510 and r is increased by one. Otherwise the method continues at S512.
At step S512, the redundant sets are eliminated. At the end each set l[c,s] is unique.
In the above example (for steps sϵ{1, 2, 3, 4} and β=12), the sets of correcting factors are defined as follows before elimination of redundant sets:
There are 6 redundant sets underlined in the above table 2. After elimination of the redundant sets, 6 candidate sets remains as defined in the table 1.
In one embodiment, this method is performed off-line. The phase correcting factors l[s;c] determined off-line may thus be stored in memory for further use by a receiver. Assuming the step s varies from 1 to smax, the number of candidate sets is provided in Table 2.
In a specific embodiment, the receiver of the base station chooses the smax value, and thus retrieves the stored phase correcting factors l[s;c] from memory. In a variant, the receiver may compute the phase correcting factors l[s;c] using the above method.
The device 100 comprises, connected by a communication bus 110: a processor or CPU (acronym of “Central Processing Unit”) 101; a random access memory RAM 102; a read only memory ROM 103; a storage unit or a storage medium reader such as an SD (acronym of “Secure Digital”) card reader 104 or a HDD (acronym of “Hard Disk Drive”); and at least one set of communication interfaces COM 105 enabling the device 100 to transmit and receive data.
The processor 101 is capable of executing instructions loaded into the RAM 102 from the ROM 103, from an external memory (such as an SD card), from a storage medium (such as the HDD), or from a communication network. When the device 100 is powered up, the processor 101 is capable of reading instructions from the RAM 102 and executing them. These instructions form a computer program causing the implementation, by the processor 101, of the method described in relation to
The methods described in relation to
Number | Date | Country | Kind |
---|---|---|---|
20185334.8 | Jul 2020 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2021/017315 | 4/23/2021 | WO |