High resolution timing advance estimation based on PRACH and sparse IFFT algorithm for LTE PRACH
11102044 · 2021-08-24
Assignee
Inventors
- Hemanth Palally (Bangalore, IN)
- Saikat Senapati (Bangalore, IN)
- Somasekhar Pemmasani (Bangalore, IN)
- Ramesh Annavajjala (Nashua, NH, US)
Cpc classification
H04W64/00
ELECTRICITY
G01S5/0036
PHYSICS
H04W56/0045
ELECTRICITY
G01S13/765
PHYSICS
International classification
Abstract
Systems, methods and computer software are disclosed for providing high resolution timing advance estimation based on Physical Random Access Channel (PRACH). An example method includes receiving a preamble signal r(n); performing signal conditioning on r(n); down sampling the signal and performing antialiasing filtering to provide a y(n) signal; correlating y(n) with a reference preamble with a reference preamble sequence c(n) to provide correlation output Ryc; using a peak value P of the correlation output Ryc to detect a preamble ID and a timing advance; constructing a sequence s(n) by segmenting r_centered(n) for length L around an index P*24; performing time domain interpolation of c(n) around index P to obtain a sequence c_interpolated(n); performing time domain interpolation between sequences s(n) and c_interpolated(n); detecting a peak position Q of the correlation; and deriving TA as P*24−L/2+q in terms of Ts.
Claims
1. A method for enhancing UE position accuracy using Physical Random Access Channel (PRACH), the method comprising: receiving a preamble signal r(n) having a predetermined sampling frequency and a predetermined length; performing signal conditioning on r(n) to provide a r_centered(n) signal; down sampling the r_centered(n) signal by a factor of 24 and performing antialiasing filtering to provide a y(n) signal; correlating the down sampled version of the received preamble y(n) with a reference preamble sequence c(n) using a Fast Fourier transform (FFT) method to provide correlation output Ryc; using a peak value P of the correlation output Ryc to detect a preamble ID and a timing advance at a resolution of 16Ts, wherein a Ts is 1/30720000 seconds; constructing a sequence s(n) by segmenting r_centered(n) for length L around an index P*24, wherein L is a number of Lags; performing time domain interpolation of c(n) around index P to obtain a sequence c_interpolated(n); performing time domain interpolation between sequences s(n) and c_interpolated(n); detecting a peak position Q of the correlation; and deriving TA as P*24 L/2+q in terms of Ts, wherein TA is a Timing Advance and q is an offset value.
2. The method of claim 1 wherein the preamble signal predetermined sampling frequency comprises 30.72 Million samples per second (sps).
3. The method of claim 1 wherein the preamble signal predetermined length comprises 27744 samples in length.
4. The method of claim 1 further comprising conditioning the preamble signal after the receiving a preamble signal.
5. The method of claim 4 wherein conditioning the preamble includes removing a cyclic prefix and extracting only a preamble sequence.
6. The method of claim 4 further comprising positioning of the signal in the frequency domain by a factor n_PRB_Offset, wherein the signal is shifted in frequency by the n_PRB_Offset so that the signal is centered about Direct Center (DC).
7. The method of claim 1 further comprising correlating the sequence y(n) with the reference preamble sequence c(n) using the FFT method Ryc=ifft{Y(k).Math.C*(k)}, where Y(k) and C(k) are 1024-point FFT of y(n) and c(n) respectively.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
DETAILED DESCRIPTION
(19) We are looking for Positioning based on Cell-ID (CID) which uses geographical knowledge of a UE's serving cell. To improve the accuracy, measurements made by the UE and/or the eNodeB can be utilized in addition.
(20) Basic CID positioning estimates the location of a UE using only the coordinates of its serving eNodeB. Typically, basic CID positioning provides only coarse estimation of the UE location, with accuracy of roughly the same order as the cell radius.
(21) Enhanced CID positioning uses additional information beyond the identity of the eNodeB that is serving the UE. The distance of a UE from its serving eNodeB or cell can be estimated from the Round Trip Time (RTT). Two measurements are defined in LTE Release 9 by which an eNodeB can indicate the RTT to the E-SMLC, namely ‘Timing Advance Type 1’ and ‘Timing Advance Type 2’. A diagram of type 1 and type 2 timing advance are shown in
(22) Type 1 and type 2 Timing Advance definition is mentioned in 3GPP spec 36.214 section 5.2.4, which is hereby incorporated by reference.
(23) TABLE-US-00001 Definition Type1: Timing advance (T.sub.ADV) type 1 is defined as the time difference T.sub.ADV = (eNB Rx − Tx time difference) + (UE Rx − Tx time difference), where the eNB Rx − Tx time difference corresponds to the same UE that reports the UE Rx − Tx time difference. Type2: Timing advance (T.sub.ADV) type 2 is defined as the time difference T.sub.ADV = (eNB Rx − Tx time difference), where the eNB Rx − Tx time difference corresponds to a received uplink radio frame containing PRACH from the respective UE.
(24) Physical Random Access Channel (PRACH) is used to achieve uplink time synchronization for a UE which either has not yet acquired, or has lost, its uplink synchronization. A successful PRACH attempt should allow subsequent UE transmissions to be inserted among the scheduled synchronized transmissions of other UEs. This sets the required timing estimation accuracy which must be achievable from the PRACH decoding at eNB
(25) PRACH Transmitter
(26) Most convenient implementation of PRACH transmitter in the UE is to generate the preamble using smallest possible IFFT and shifting the preamble to the required frequency location through time-domain up sampling and filtering (commonly known as hybrid frequency/time domain PRACH generation shown in the
(27) A PRACH Preamble structure is shown in
(28) PRACH Formats
(29) TABLE-US-00002 Preamble format T.sub.CP T.sub.SEQ 0 3168 .Math. T.sub.s .sup. 24576 .Math. T.sub.s 1 21024 .Math. T.sub.s .sup. 24576 .Math. T.sub.s 2 6240 .Math. T.sub.s 2 .Math. 24576 .Math. T.sub.s 3 21024 .Math. T.sub.s 2 .Math. 24576 .Math. T.sub.s 4* 448 .Math. T.sub.s .sup. 4096 .Math. T.sub.s
(30) As described above Type 2 Measurement report of LPPa protocol requires 2Ts/8Ts resolution timing advance to better locate UE positioning. Physical layer is responsible for decoding PRACH Preamble and its corresponding Timing Advance (TA) and report to higher layers. Traditional measurement report gets 16Ts timing resolution from Physical layer (PHY) which boils down to approx. 80 meters range of UE positioning. As part of LPPa requirement, PHY is required to report Timing Advance with higher accuracy of 2Ts. This 2Ts resolution will give UE positioning accuracy up to 10 meters.
Solution to Problem
(31) The current implementation of the PRACH detector estimates the timing advance with the accuracy of 16Ts. The proposed method works in tandem with the existing implementation to estimate the TA with higher resolution of 2Ts. A key insight is that the timing offset can be determined from the correlation of the received PRACH preamble with the reference received preamble. In other words, the transmitted preamble is known to be an interpolated or upsampled version of the original short preamble sequence.
(32) The processing steps required in the method are as follows, and are shown in
(33) The current implementation of the PRACH detector estimates the timing advance with the accuracy of 16Ts. The proposed method works in tandem with the existing implementation to estimate the TA with higher resolution of 2Ts.
(34) This method removes the need for large IFFT as we are performing correlation and interpolation in time domain.
(35) The processing steps required in the method are as follows
(36) The received preamble signal r(n) has the sampling frequency of 30.72 Msps and is 27744 samples long. The following signal conditioning operations are performed on r(n)
(37) At 301, cyclic prefix is removed, and only Preamble sequence is extracted. The resulting sequence will be 24576 samples long
(38) The positioning of the signal in the frequency domain is determined by the factor n_PRB_Offset. The signal is shifted in frequency by n_PRB_Offset so that it is centered about DC.
(39) Anti-alias filtering in time domain to remove out of band noise. Let the resulting signal be denoted as r_centred(n)
(40) At 302, the frequency centered signal is down sampled by a factor of 24 after performing antialiasing filtering. Let this signal be called Y(n). At 303, a reference preamble may be used.
(41) At 304, the down sampled version of the received preamble is correlated with the reference preamble sequence c(n) using the FFT method as below
Ryc=ifft{Y(k).Math.C*(k)}
(42) Where Y(k) and C(k) are 1024-point FFT of y(n) and c(n) respectively
(43) At 305, the peak position P of the correlation output Ryc is used to detect the preamble ID and the Timing advance at a resolution of 16Ts
(44) At 306, we construct a sequence s(n) by segmenting r_centred(n) for length L centered around the index P*24.
(45) At 307, we perform time domain interpolation of the reference preamble c(n) around the index P to get a sequence c_interpolated(n). At 308, a table for interpolation coefficients may be used
(46) As we know that the transmitted preamble is constructed by interpolating only the 839 samples preamble sequence, the interpolation coefficients can be derived beforehand and stored in a LUT. Here we are exploiting the fact that the preamble at sampling rate of Fs=30.72 Msps is highly correlated in nature.
(47) Perform time domain interpolation between the sequences s(n) and c_interpolated(n).
(48) At 310, detect the peak position Q of the correlation from 309.
(49) At 311, TA can be derived as P*24−L/2+Q in terms of Ts.
(50) The insight that the timing offset can be determined from the correlation of the received PRACH preamble with the reference received preamble also means that the equivalent correlation could be performed in any other RAT with an equivalent PRACH. The PRACH preamble structure used in 5G NR is exactly the same as that used in 4G. So all high resolution timing advance (HRTA) related methods included herein can be used for 5G.
(51) The standard LTE conformance requires the timing advance to be done with an Accuracy of 16Ts. We are estimating it at a higher resolution of 2Ts. This is useful for public safety networks using LTE or 5G technology base stations.
(52) In this method we have to store the sequence c centered(n) in the memory till we find the peak P. So, we are making a tradeoff between memory usage and computation complexity.
(53) This method is designed to be performed at the base station.
(54) In various embodiments, this method can be designed to be performed at the base station or could be performed anywhere where sufficient processing power is available, including at a baseband unit colocated or remotely located relative to the base station. No modification is required to the UE.
(55) It is noted that the PRACH preamble structure used in 5G NR (New Radio) is exactly the same as that used in 4G, so the methods disclosed herein could be used by one having skill in the art at a 5G base station or in a 5G network. It is noted that different PRACH preamble structures could be standardized and used, and, when used in combination with the presently disclosed methods and new RF technologies, could provide UE positioning with greater resolution than 2 Ts.
(56) The timing advance values derived using the presently described methods could be used, in some embodiments, for UE location/positioning; responding to UE paging requests; reporting UE location to the core network; plotting UE location on a map; using UE location for compliance with governmental requirements for location services and/or emergency services; optimizing beamforming for MIMO applications; optimizing UE location for purposes of assessing handover and inter-cell interference coordination (ICIC), etc. The timing advance values derived using the presently described methods could, in the case of a UE that is a smart watch or smart fob or smart tag, be used to locate the UE within a small radius.
(57)
(58) The block diagram of a typical LTE PRACH receiver with High resolution timing advance estimation capability is as shown in
(59) First, we explain the properties of the received PRACH signal that we use in the proposed method and in the second section we provide the steps involved in the proposed method
(60) Property1: PRACH signal is inherently Windowed in Frequency Domain
(61) It can be seen from
(62) However, DFT of the transmitted preamble (
W(k)=1, −419≤k≤419
(63) W(k)=0, otherwise
(64) The IFFT of W(k).fwdarw.w(n) is a sinc function defined as below
(65)
(66) M=839, N=16384
(67)
(68) Due to the ideal autocorrelation property of the basic preamble its circular correlation (
(69) This is because windowing in frequency domain with a rectangular pulse equates to convolution in time domain with sinc pulse
(70) Property2: Magnitudes of the peaks of circular correlation of a preamble and its up sampled version are proportional by a factor determined by the up-sampling factor
(71) Let P.sub.mag be the magnitude of the peak P from 1024-point circular correlation obtained at module (1) of
(72) Then
(73)
(74) Algorithm
(75) The steps involved in the proposed method are as follows
(76) Define a search window L within which the peak Q is located, with the center of the search window S=16P
(77) Define two points A and B that are equidistant on either side of S
(78) Use Grouped FFT method to compute the transformed outputs at points A and B
(79) Choose between point A and B, whichever has greater magnitude
(80) Use the magnitude and index from step 4 and the equations Eq1 and Eq2 to define a narrow search zone Z
(81) Compute transformed outputs with in the Z
(82) Find the position of the max from above step
(83) Case 1: When peak is in the center of search window, shown in
(84) Case 2: When peak is in the right side of search window, shown in
(85) Case 3: when peak is in the left side of search window, shown in
(86) Step-Wise Peak Search
(87) Under the Rayleigh and Rician channel conditions, the central lobe will not be symmetric about the peak. Under such conditions, stepwise peak search can be used to find the peak. This method is also computationally efficient, and is shown in
(88) Stepwise Peak Search Steps
(89) Let S be the center of the search window. S=16P, where P is the peak position obtained from 1024 point correlation.
(90) Choose a Point A (value A is empirically chosen) that is at the left of the S such a way that we definitely capture a point in the left half of the main lobe in spite of the movement of the peak anywhere in the search window.
(91) From Point A move to the right in steps of A and compute transformed output at each of those steps.
(92) Stop the movement to the right when the current transformed output is lesser than the previous. This indicates that we have just passed the peak.
(93) Move backwards in smaller steps or step of one till we find the peak.
(94) Step wise peak search reduces the number of transformed outputs to be computed to less than one-third for each peak.
(95) Example of Complexity Reduction
(96) In a typical real time PRACH processing we need to detect multiple peaks as shown in
(97) Let L=60, Num Peaks to be detected=16, Z=4
(98) For each peak we need only 6 transformed outputs to be computed that include 4 points in the narrower search zone Z and two for points A and B.
(99) So totally only 6*16=96 transformed outputs need to be computed
(100)
(101)
(102) Processor 1702 and baseband processor 1706 are in communication with one another. Processor 1702 may perform routing functions, and may determine if/when a switch in network configuration is needed. Baseband processor 1706 may generate and receive radio signals for both radio transceivers 1712 and 1714, based on instructions from processor 1702. The steps described herein could be performed at the baseband processor 1706, or, in some embodiments, at processor 1702. In some embodiments, processors 1702 and 1706 may be on the same physical logic board. In other embodiments, they may be on separate logic boards.
(103) Processor 1702 may identify the appropriate network configuration, and may perform routing of packets from one network interface to another accordingly. Processor 1702 may use memory 1704, in particular to store a routing table to be used for routing packets. Baseband processor 1706 may perform operations to generate the radio frequency signals for transmission or retransmission by both transceivers 1710 and 1712. Baseband processor 1706 may also perform operations to decode signals received by transceivers 1712 and 1714. Baseband processor 1706 may use memory 1708 to perform these tasks.
(104) The first radio transceiver 1712 may be a radio transceiver capable of providing LTE eNodeB functionality, and may be capable of higher power and multi-channel OFDMA. The second radio transceiver 1714 may be a radio transceiver capable of providing LTE UE functionality. Both transceivers 1712 and 1714 may be capable of receiving and transmitting on one or more LTE bands. In some embodiments, either or both of transceivers 1712 and 1714 may be capable of providing both LTE eNodeB and LTE UE functionality. Transceiver 1712 may be coupled to processor 1702 via a Peripheral Component Interconnect-Express (PCI-E) bus, and/or via a daughtercard. As transceiver 1714 is for providing LTE UE functionality, in effect emulating a user equipment, it may be connected via the same or different PCI-E bus, or by a USB bus, and may also be coupled to SIM card 1718. First transceiver 1712 may be coupled to first radio frequency (RF) chain (filter, amplifier, antenna) 1722, and second transceiver 1714 may be coupled to second RF chain (filter, amplifier, antenna) 1724.
(105) SIM card 1718 may provide information required for authenticating the simulated UE to the evolved packet core (EPC). When no access to an operator EPC is available, a local EPC may be used, or another local EPC on the network may be used. This information may be stored within the SIM card, and may include one or more of an international mobile equipment identity (IMEI), international mobile subscriber identity (IMSI), or other parameter needed to identify a UE. Special parameters may also be stored in the SIM card or provided by the processor during processing to identify to a target eNodeB that device 1700 is not an ordinary UE but instead is a special UE for providing backhaul to device 1700.
(106) Wired backhaul or wireless backhaul may be used. Wired backhaul may be an Ethernet-based backhaul (including Gigabit Ethernet), or a fiber-optic backhaul connection, or a cable-based backhaul connection, in some embodiments. Additionally, wireless backhaul may be provided in addition to wireless transceivers 1712 and 1714, which may be Wi-Fi 1702.11a/b/g/n/ac/ad/ah, Bluetooth, ZigBee, microwave (including line-of-sight microwave), or another wireless backhaul connection. Any of the wired and wireless connections described herein may be used flexibly for either access (providing a network connection to UEs) or backhaul (providing a mesh link or providing a link to a gateway or core network), according to identified network conditions and needs, and may be under the control of processor 1702 for reconfiguration.
(107) A GPS module 1730 may also be included, and may be in communication with a GPS antenna 1732 for providing GPS coordinates, as described herein. When mounted in a vehicle, the GPS antenna may be located on the exterior of the vehicle pointing upward, for receiving signals from overhead without being blocked by the bulk of the vehicle or the skin of the vehicle. Automatic neighbor relations (ANR) module 1732 may also be present and may run on processor 1702 or on another processor, or may be located within another device, according to the methods and procedures described herein.
(108) Other elements and/or modules may also be included, such as a home eNodeB, a local gateway (LGW), a self-organizing network (SON) module, or another module. Additional radio amplifiers, radio transceivers and/or wired network connections may also be included.
(109) The system may include 5G equipment. 5G networks are digital cellular networks, in which the service area covered by providers is divided into a collection of small geographical areas called cells. Analog signals representing sounds and images are digitized in the phone, converted by an analog to digital converter and transmitted as a stream of bits. All the 5G wireless devices in a cell communicate by radio waves with a local antenna array and low power automated transceiver (transmitter and receiver) in the cell, over frequency channels assigned by the transceiver from a common pool of frequencies, which are reused in geographically separated cells. The local antennas are connected with the telephone network and the Internet by a high bandwidth optical fiber or wireless backhaul connection.
(110) 5G uses millimeter waves which have shorter range than microwaves, therefore the cells are limited to smaller size. Millimeter wave antennas are smaller than the large antennas used in previous cellular networks. They are only a few inches (several centimeters) long. Another technique used for increasing the data rate is massive MIMO (multiple-input multiple-output). Each cell will have multiple antennas communicating with the wireless device, received by multiple antennas in the device, thus multiple bitstreams of data will be transmitted simultaneously, in parallel. In a technique called beamforming the base station computer will continuously calculate the best route for radio waves to reach each wireless device, and will organize multiple antennas to work together as phased arrays to create beams of millimeter waves to reach the device.
(111) The foregoing discussion discloses and describes merely exemplary embodiments of the present invention. In some embodiments, software that, when executed, causes a device to perform the methods described herein may be stored on a computer-readable medium such as a computer memory storage device, a hard disk, a flash drive, an optical disc, or the like. As will be understood by those skilled in the art, the present invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. Various components in the devices described herein may be added, removed, split across different devices, combined onto a single device, or substituted with those having the same or similar functionality.
(112) Although the present disclosure has been described and illustrated in the foregoing example embodiments, it is understood that the present disclosure has been made only by way of example, and that numerous changes in the details of implementation of the disclosure may be made without departing from the spirit and scope of the disclosure, which is limited only by the claims which follow. Various components in the devices described herein may be added, removed, or substituted with those having the same or similar functionality. Various steps as described in the figures and specification may be added or removed from the processes described herein, and the steps described may be performed in an alternative order, consistent with the spirit of the invention. Features of one embodiment may be used in another embodiment. Other embodiments are within the following claims.