Reciprocal calibration for channel estimation based on second-order statistics

Information

  • Patent Grant
  • 11329848
  • Patent Number
    11,329,848
  • Date Filed
    Thursday, June 13, 2019
    5 years ago
  • Date Issued
    Tuesday, May 10, 2022
    2 years ago
Abstract
A wireless communication method includes receiving, by a first wireless device during a training phase, reference tones using a first number of resource elements from a transmitter of a second wireless device, wherein the first wireless device comprises multiple receiving antennas, estimating, by the first wireless device, from the receiving the reference tones, a second order statistics of wireless channels between the multiple receiving antennas and the transmitter of the second wireless device, and performing channel estimation, during an operational phase subsequent to the training phase, using the second order statistics and reference tones received on a second number of resource elements, wherein the second number is less than the first number.
Description
TECHNICAL FIELD

The present document relates to wireless communication, and more particularly, precoding of wireless signals for transmission.


BACKGROUND

Due to an explosive growth in the number of wireless user devices and the amount of wireless data that these devices can generate or consume, current wireless communication networks are fast running out of bandwidth to accommodate such a high growth in data traffic and provide high quality of service to users.


Various efforts are underway in the telecommunication industry to come up with next generation of wireless technologies that can keep up with the demand on performance of wireless devices and networks. Many of those activities involve situations in which a large number of user devices may be served by a network.


SUMMARY

This document discloses techniques for channel estimation. In one example, a reciprocal calibration technique is used. In another example, a minimal set of resource elements and an estimate of second order statistics of the channel are used.


In an example aspect, a wireless communication method is disclosed. The method includes receiving, by a first wireless device during a training phase, reference tones using a first number of resource elements from a transmitter of a second wireless device, wherein the first wireless device comprises multiple receiving antennas; estimating, by the first wireless device, from the receiving the reference tones, a second order statistics of wireless channels between the multiple receiving antennas and the transmitter of the second wireless device; and performing channel estimation, during an operational phase subsequent to the training phase, using the second order statistics and reference tones received on a second number of resource elements, wherein the second number is less than the first number.


In another example aspect, a method of wireless communication is disclosed. The method includes transmitting, to multiple receive antennas of a first wireless device from a transmit antenna of a second wireless device, during a training phase, reference tones using a first number of resource elements of a wireless channel between the transmit antenna and the multiple receive antennas; receiving, at an end of the training phase, an instruction to transmit reference tones using a second number of resource elements; and transmitting, during an operational phase after the training phase, reference tones to the multiple receive antennas of the first wireless device using the second number of resource elements, wherein the second number is different from the first number and wherein the second number is based on an estimated second order statistics of the wireless channel.


In yet another example aspect, a method of wireless communication is disclosed. The method includes estimating, during a training phase, a second order statistics for a first wireless channel and a second wireless channel between a transmitter and a receiver comprising multiple antennas, wherein the second order statistics is estimated using reference tones transmitted on a first number of resource elements; predicting, during an operational phase subsequent to the training phase, an estimate of the second wireless channel based on the second order statistics and an estimate of the first wireless channel calculated using reference tones transmitted on a second number of resource elements, where the second number is less than the first number; and communicating, during the operational phase, over the second wireless channel using the estimate of the second wireless channel resulting from the predicting; wherein the first wireless channel and the second wireless channel include non-overlapping frequencies.


In yet another example aspect, a method of wireless communication is disclosed. The method includes receiving, by a first communication device, a number of subcarriers from a second communication device, each including a corresponding reference signal; calculating an inversion factor for each subcarrier based on a received value of the corresponding reference signal; and transmitting by the first communication device to the second communication device, at least some of the subcarriers by scaling a pilot signal using a corresponding inversion factor.


In yet another example aspect, a method of wireless communication is disclosed. The method includes transmitting, to a first communication device, from a second communication device, a number of subcarriers, each subcarrier including a corresponding reference signal; receiving, from the first communication device, at least some of the subcarriers carrying pilot signals scaled by inversions factors for the subcarriers; and estimating a communication channel between the second communication device to the first communication device using the inversion factors.


In yet another example aspect, a wireless communication apparatus that implements the above-described methods is disclosed.


In yet another example aspect, the methods may be embodied as processor-executable code and may be stored on a computer-readable program medium.


In yet another example aspect, a wireless communication system that operates by providing a single pilot tone for channel estimation is disclosed.


These, and other, features are described in this document.





DESCRIPTION OF THE DRAWINGS

Drawings described herein are used to provide a further understanding and constitute a part of this application. Example embodiments and illustrations thereof are used to explain the technology rather than limiting its scope.



FIG. 1 shows an example of a fixed wireless access network.



FIG. 2 shows another example of a fixed wireless access network.



FIG. 3A shows an example of a wireless channel between a transmitter and a receiver device.



FIG. 3B shows another example of a wireless channel between a transmitter and a receiver device.



FIG. 3C shows an example of a wireless channel between a single transmit antenna and multiple receive antennas.



FIG. 4 is a graph showing an example of effect of number of tones used for reference signal transmissions on average error covariance.



FIG. 5 shows an example of a mobile wireless communication network with a multi-antenna base-station and multiple user antennas.



FIG. 6 shows an example of FDD prediction setup, wherein the base-station uses the channel response of the uplink to predict the channel response of the downlink in a different frequency and on a later subframe.



FIG. 7 shows examples of bandwidth partitions.



FIG. 8 shows an example of a bandwidth partition with the same time interval for reference signal (RS) 1 and RS2.



FIG. 9 shows an example of a bandwidth partition with different time intervals for reference signal (RS) 1 and RS2.



FIG. 10 shows an example of channel prediction over the same time interval.



FIG. 11 shows an example of channel prediction over a different time interval.



FIG. 12 shows an example of a wireless channel with reciprocity.



FIG. 13 shows an example antenna configuration in which four transmit and four receive antennas are used at a network-side apparatus.



FIG. 14 shows an example antenna configuration in a user-side communications apparatus.



FIG. 15 shows a block diagram for an example implementation of reciprocity calculation.



FIG. 16 shows an example of a wireless transceiver apparatus.



FIG. 17 shows an example communication network.



FIG. 18 is a flowchart for an example wireless communication method.



FIG. 19 is a flowchart for another example wireless communication method.



FIG. 20 is a flowchart for another example wireless communication method.



FIG. 21 is a flowchart for another example wireless communication method.



FIG. 22 is a flowchart for another example wireless communication method.





DETAILED DESCRIPTION

To make the purposes, technical solutions and advantages of this disclosure more apparent, various embodiments are described in detail below with reference to the drawings. Unless otherwise noted, embodiments and features in embodiments of the present document may be combined with each other.


Section headings are used in the present document to improve readability of the description and do not in any way limit the discussion or embodiments (and/or implementations) to the respective sections only.


Channel knowledge is an important component in wireless communication, whether it is for a receiver to equalize and decode the received signal, or for a multi-antenna transmitter to generate a more efficient precoded transmission.


Typically, channel knowledge is acquired by transmitting known reference signals (pilots) and interpolating observations and results obtained by receiving them at the receiver over the entire bandwidth and time of interest. The density of the reference signals may depend on the characteristics of the channel. Higher delay spreads and higher Doppler spreads may require denser reference signals along frequency and/or time respectively, thus, occupying a considerable amount of the available capacity.


For example, consider a downlink multi-carrier transmission from a base-station to multiple users. In this case, depending on channel conditions, the reference signals should be dense enough in the time-frequency grid to allow the users to correctly estimate the channel and equalize the data (e.g., the estimation error due to interpolation of channel estimates along the time-frequency grid should be minimized). Another example is that of an uplink transmission from multi-users to a base station equipped with multi-antennas, where each user needs to allocate resource elements for its own reference signals.


In this document, we propose (i) an efficient method for estimating the channel from a very small set of resource elements, using the second-order statistics of the channel, and (ii) using the second-order statistics to reciprocally calibrate the wireless channel. The effectiveness of (i) increases with increased number of antennas at the base-station. As further described below, in some cases, it may be possible to send a single reference signal (tone) from the transmit antenna to each receive antenna at the base station, and be able to determine channel at different times and in different frequency bands based on the channel estimate obtained from the single (or, in general, reduced number of) reference tones. The efficacy of (ii) is based on deriving knowledge of the non-reciprocal parts of a channel between two nodes, and adjusting the channel response such that precoding may be used for the transmissions.


In the description, the example of a fixed wireless access (FWA) system is used only for illustrative purpose and the disclosed techniques can apply to other wireless networks, such as cellular or mobile communication networks further described in the present document.



FIG. 1 shows an example of a fixed wireless access system 100. A hub 102, that includes a transmission facility such as a cell tower, is configured to send and receive transmissions to/from multiple locations 104. For example, the locations could be user premises or business buildings. As described throughout this document, the disclosed embodiments can achieve very high cell capacity fixed wireless access, when compared to traditional fixed access technology. Some techniques disclosed herein can be embodied in implementations at the hub 102 or at transceiver apparatus located at the locations 104.



FIG. 2 shows yet another configuration of a fixed access wireless communication system 200 in which hops are used to reach users. For example, one cell tower may transmit/receive from another cell tower, which would then relay the transmissions between the principle cell tower and the users, thus extending range of the fixed wireless access system. A backhaul may connect the transmission tower 202 with an aggregation router. For example, in one configuration, a 10 Gbps fiber connection may be used to feed data between a base station at a hub and a fiber hub aggregation router. In one advantageous aspect, deployment of this technology can be achieved without having to change any network bandwidth characteristics for harder to reach areas by using the hub/home access point (AP) configuration as a launch point. Some techniques disclosed herein can be embodied in implementations at the macro tower 202 or at transceiver apparatus located at the other locations.


1. Introduction to Embodiments for Reciprocal Calibration for Channel Estimation Based on Second-Order Statistics


FIGS. 3A, 3B and 3C show simplified wireless networks to highlight certain aspects of the disclosed technology. A transmitter transmits wireless signals to a receiver in the wireless network. Some transmissions in the network, variously called as downstream or downlink transmissions, a network-side node such as a base station acts as a transmitter of wireless signals and one or more user devices act as the receiver of these wireless signals. For some other transmissions, as depicted in FIG. 3A, the direction of transmission may reverse. Such transmissions are often called uplink or upstream transmissions. For such transmissions, one or more user devices act as transmitters of the wireless signals and a network-side node such as the base station acts as the receiver of these signals (as depicted in FIG. 3A). Other type of transmissions in the network may include device-to-device transmissions, sometimes called direct or sideband transmissions. In frequency division multiplexing (FDM) networks (also called frequency division duplexing or FDD networks), the transmissions to a base station and the transmissions from the base station may occupy different frequency bands (each of which may occupy continuous or discontinuous spectrum). In time division multiplexing (TDM) networks (also called time division duplexing, or TDD, networks), the transmissions to a base station and the transmissions from the base station occupy a same frequency band, but are separated in time domain using a TDM mechanism such as time slot based transmissions. Other types of multiplexing are also possible (e.g., code division multiplexing, orthogonal time frequency space, or OTFS, multiplexing, spatial multiplexing, etc.). In general, the various multiplexing schemes can be combined with each other.


From the early days of single antenna systems in which the transmitting user device had a single transmission antenna and the receiving base station had a single receiving antenna for that signal, wireless systems have now evolved to configurations in which transmitters have multiple transmit antennas and receivers have multiple receive antennas. Such systems are often called multiple input multiple output (MIMO) systems. In a MIMO system, a wireless channel may therefore be representable as a matrix that acts on a vector representing the multiple inputs and produces the vector representing the multiple outputs. Depending on channel conditions, some entries of the channel matrix may change over time.


The wireless systems described in this document, subjected to some interference, are known to be able to achieve the channel capacity if the interference is known to the transmitter. By using a technique known as dirty paper coding (DPC), the transmitter can precode a transmission such that the signal arriving at the receiver is interference-free. In reality, a more implementation-friendly approximation of DPC known as Tomlinson-Harashima Precoding (THP) is likely to be used (and is discussed in a subsequent section of this document).


A precoded transmission is based on the knowledge of the exact channel response between the transmitting antenna(s) of a first terminal denoted by A—typically a base-station (BS)—to the receiving antenna(s) of a second terminal denoted by B—typically a piece of Consumer Premises Equipment (CPE) or a user equipment (UE). This channel response can be considered to be composed of three different parts as illustrated in FIG. 3B. First, the channel response of the transmitter in terminal A. Second, the channel response of the different reflectors. Third, the channel response of the receiver in terminal B. The transmitter channel responses may be due to the transmit chain circuitry such as power amplifiers, digital to analog converters (DACs), frequency upconverters and so on. The receiver channel response may be due to receiver-side circuitry such as low noise block (LNB), frequency downconverter, analog to digital conversion circuitry (ADC).


1.A. Second-Order Statistics

Without any loss of generality, let's consider a basic setup with a transmitter (e.g., a user device) using a single antenna and a receiver (e.g., a base station) with L antennas. If the transmitter (e.g., user device) has more than one antenna, then each antenna can be referred to as a separate transmitter.


Prior to the transmission of data, the system performs a preliminary training phase, where the second-order statistics of the channels between the transmitter and each one of the L receiver antennas are computed.



FIG. 3C shows an example of the basic setup in which the receiver comprises L receive antennas and the transmitter comprises a single transmit antenna. The signal transmitted from the transmitter is reflected (and dispersed) by a number of reflectors that give the wireless channel its characteristics. As further discussed in this document, the reflectors may be static or mobile, and have time-variant reflective characteristics.


In this initial phase, or training phase, the transmitter sends known reference signals, enabling the receiver to estimate the channel across all the bandwidth of interest and all the antennas. One example is in a multi-carrier scenario where the bandwidth of interest is represented with Nf discrete tones, at time instance k the receiver will estimate the channel response hk as a complex column vector of size N=Nf·L. This step may be repeated NSOS times and afterwards, the second-order statistics may be computed. Two different methods for this computation are proposed here as an example:


Direct method—In this method, the covariance matrix is calculated by arranging all channel responses as columns of the transfer matrix H and computing the cross correlation as:

H=[h1|h2| . . . |hNSOS]  (1)
RHH=H·H*  (2)


Average method—using this method, an average over NSOS of the variance of each channel response are calculated.










R
HH

=


1

N
SOS







k
=
1


N
SOS









h
k

·

h
k
*








(
3
)







The covariance matrix RHH is large with dimensions N×N. For example, in a 4 antenna system in which 8 tones are used for reference signal transmissions, N=32, and thus RHH is a 32×32 matrix. However, the matrix encodes only a very small amount of information originating from the delay-Doppler profile and angle of arrival/departure of the reflectors of the wireless channel. These essential parameters of the channel can be extracted from the covariance matrix using a mathematical operation such as the Principal Component Analysis (PCA). For example, these parameters may be identified by computing the dominant eigenvalues and their corresponding eigenvectors of RHH. For extracting K dominant reflectors, RHH may be approximated by

RHH≈V·D·V*  (4)


where V is an N×K eigenvectors matrix and D is the corresponding K×K eignenvalues diagonal matrix.


1.B. Channel Estimation Using a Small Set of Resource Elements

Once the second-order statistics of the channels between the transmitter and the receivers' L antennas are computed, the transmitter (e.g., a single transmit antenna of a user device) may transmit only a very small set of known reference signals, from which the receiver will be able to estimate all the channels (to all the antennas and across all the bandwidth). The general rule of thumb is that the receiver will use at least K measurements to sense the parameters of the K dominant reflectors.


More formally, let custom character={m1, m2, . . . } be a set of frequency elements (tones), where mi takes values in the range [1, 2, . . . , Nf]. Note that the size of custom character, denoted by |custom character|, depends on different parameters such as the number of antennas L, the number of reflectors K, their delay profile and their angle of arrival (AOA) to the receiver antennas. Typically, |custom character| is very small and in some cases, may even be equal one (a single tone).


The transmission of a small set of known reference signals over the tones of custom character is received over the L receiver's antennas. Let us define by M=|custom character|·L the number of received reference values, corresponding to this transmission. Using the second-order statistics, an interpolation MMSE filter, C, can be computed from these M received elements to estimate the channel over all the N elements. Let HM be the channel response over these M received elements. Then, C may be computed as

C=RHHM·(RHMHM)−1  (5)


where RHHM and RHMHM corresponds to decimated versions of the full RHH matrix (one or two of the dimensions is reduced to M).


The PCA approximation for the second-order statistics can be useful to reduce the complexity of such a computation. Let VM be the matrix V decimated along its N dimension to M. Then, C can be expressed as

C=V·D·VM*·(VM·D·VM)−1=V·(VM)−1  (6)


The last computation requires inverting VM, which has the dimensions of M×K. However, by algebraic manipulation, the inverted matrix may be even further reduced to a size of K×K

C=V·(VM*·VM)−1·VM*  (7)


Finally, the full channel response is obtained by computing

H=C·HM  (8)


1.C. Estimating the Sufficient Number of Resource Elements

The receiver may use different methods, to estimate the required number of (time-frequency) resource elements in custom character, to provide a sufficient quality of channel estimation. One method to do this, is to compute the error covariance matrix

RE=V·D·V*−C·VM·D·V*=(V−C·VMD·V*  (9)


and use it to compare the mean energy along the diagonal to a threshold










d

R
E


=

diag


(

R
E

)






(
10
)








R
E

_



=
Δ





1
N






i
=
1

N











d

R
E




(
i
)




2



<
TH





(
11
)







This criterion may be used to determine whether the number of elements in custom character is sufficient. Afterwards, the receiver may instruct the transmitter to adjust the number of required reference signals to transmit.


To illustrate this concept, an extreme case of channels with many reflectors and small angle-of-arrivals was examined with different number of tones allocated for reference signals. White Gaussian noise was added to the received signal, creating a signal-to-noise (SNR) ratio of 60 dB. The graph in FIG. 4, shows for a varying number of tones with reference signals (1 to 5), the computed values of the average error covariance and the average estimation SNR (inverse of mean-square-error compared to the true channel). The computed average error covariance, RE, can be roughly compared to a threshold of 235 dB, afterwards the noise contribution of the estimated channel is smaller than that of the AWGN. Therefore, for this example, it would be recommended to use 5 tones with reference signals.


1.D. Embodiments for Prediction in FDD/TDD Networks

The methods described in the previous sections may also be applied directly for efficient channel prediction. Channel prediction may include a training phase, where reference signals are transmitted over two different time-frequency resources sets and the second-order statistics is computed for the corresponding channels. These sets may include two non-overlapping frequency bands, two different time instances, or a combination of both. For example, in FDD networks, downlink and uplink transmissions use different frequency bands, and channel characteristics over one of downlink/uplink channel could be predicted based on observed channel estimates over the other (uplink or downlink) channel.


Let hk,1 represent a column vector of the channel response over Nf1 frequency tones and L antennas, across a first frequency band and at time instance Tk (first set) and hk,2 as a column vector of the channel response over Nf2 frequency tones and L antennas, across a second frequency band and at time instance Tk+Δ (second set). A matrix H can be formulated as the concatenation of these column vectors from multiple training sessions and compute the covariance matrix RHH.









H
=

[





h

1
,
1







h

2
,
1













h

2
,
1







h

2
,
2

















h


N
SOS

,
1







h


N
SOS

,
2





]











(
12
)









R
HH
=H·H*  (13)


The procedure described in previous sections, may be fully applied here as well. PCA may be applied to RHH for K dominant reflectors. In the operational phase, the transmitter may send a small set of reference signals, m, over part of the first time-frequency set and use the computed interpolation filter

C=RHHM·(RHMHM)−1  (14)
H=C·HM  (15)


to estimate the channel across the entire second time-frequency set.


For example, after a training session, multiple user devices are transmitting their small set of resources to the base station, where it is used for predicting future channels in a different band (FDD) or in the same band (TDD). The base station uses the predicted channels to compute a precoded transmission back to the users.


This method may also be applied for non-direct prediction techniques. For example, as a first step, the second-order statistics of the channel may be computed from full-band reference signals. Then, small sets of reference signals at both time-frequency sets are transmitted. The received signals in each set, consisting of M elements may be further compressed to K elements by computing

HK=(VM)−1·HM  (16)


Then, the covariance matrices of the two sets and the cross-covariance matrix may be computed along with an interpolation filter from one set to the other

R(1)=HK(1)·(HK(1))*  (17)
R(2,1)=HK(2)·(HK(1))*  (18)
C(2,1)=R(2,1)·(R(1))−1  (19)


Finally, the channel in the second set is predicted from the channel in the first set

HK(2)=C(2,1)·HK(1)  (20)

and computed over all the N elements

H=V·HK(2)  (21)


1.E. Channel Responses and Channel Reciprocity

There are two main differences between the channel responses at terminals A and B and the channel response of the wireless channel reflectors:


1. The channel response of the wireless channel reflectors in a time-division duplex (TDD) system is reciprocal whereas the channel response of the terminals is not.


2. The channel response of the wireless channel reflectors may change rapidly (e.g., in 1-10 milliseconds, depending on the Doppler of the reflectors and terminals), but the channel response of the terminals changes slowly, mostly with temperature.


There are several methods for obtaining the complete channel response from terminal A to B described in the literature. For example, an explicit method would be to send known reference signals from terminal A to B and have terminal B transmit back the values of the received reference signals to terminal A. This is often referred to as explicit feedback. However, each value must be represented with multiple bits, and in a system where terminal A has many antennas, there are many user terminals and significant Doppler effects causing the propagation channel to change rapidly, the amount of information that needs to be transmitted can severely reduce the overall system efficiency. In the extreme case with high levels of Doppler, it is simply not possible to feedback all the required Channel State Information (CSI) quickly enough, resulting in stale CSI and suboptimal precoding.


Instead, a TDD system can use an approach known as “reciprocity calibration” to obtain the relationship between the non-reciprocal parts of the channel response in both transmission directions: the AB (from A to B) and the BA (from B to A). Terminal B first transmits known reference signals that allow terminal A to compute the AB channel response. Using knowledge of the non-reciprocal relationship, terminal A can adjust the BA channel response to make it suitable for precoding a transmission back to terminal B.


More formally, for a multi-carrier TDD system that uses multi-carrier modulation, where the channel can be described as a complex value in the frequency domain for a specific sub-carrier (tone), the three components of the AB channel response can be denoted as HATX, HCH and HBRX. Similarly, the three components of the BA channel response are HBTX, HCH and HARX. The overall downlink (AB) channel response is

HAB=HATX·HCH·HBRX  (22)


and the overall uplink (BA) channel response is

HBA=HBTX·HCH·HARX  (23)


From HAB and HBA, the reciprocity calibration factor can be written as









α
=



H
A
TX

·

H
B
TX




H
B
TX

·

H
A
RX







(
24
)







Therefore, if HBA is known at terminal A, it can compute HAB=αHBA. The question that remains is how to obtain α. Note that for the multi-carrier system, the above Equations (22) to (24) will provide reciprocity calibration values and channel responses on a per sub-carrier basis for sub-carriers on which reference signals are transmitted.


Different methods exist within the literature for computing the reciprocity calibration factor. The most straight forward of these is to utilize explicit feedback as described above, but only feed back HAB when α is re-calculated. Since the transmitter and receiver channel responses change relatively slowly, the rate of feedback is typically in the order of minutes and thus represents negligible overhead for a modest number of terminals and antennas. However, when the number of antennas in terminal A and the number of CPEs (terminal B) is large, as can be the case in a massive multiple-input multiple-output (MIMO) system with many subscribers, the feedback overhead can consume a considerable portion of the system capacity.


Another approach is to have terminal A transmit reference signals between its own antennas and calculate calibration factors for only HATX and HARX. That is, obtain:










α
A

=


H
A
TX


H
A
RX






(
25
)








which results in

{tilde over (H)}ABA·HBA=HATX·HBTX·HCH  (26)


Terminal A will then precode one reference symbol using {tilde over (H)}AB that terminal B can use to remove its HBTX and HBRX contributions from all subsequent precoded transmissions. This technique may be called relative calibration. Whilst this approach entirely removes the need for feedback of HBA, the need for terminal A to transmit to itself during a calibration procedure and then to CPEs that could be located many hundreds of meters or even kilometers away can create dynamic range challenges. It is typically desirable to use the same hardware gain stages in the transmit chain when calibrating as those used for transmission, since having to switch gain stages between calibration and transmission can change the nature of HATX and HARX.


This document describes, amongst other embodiments and approaches, a new approach for computing the reciprocity calibration factor that avoids the dynamic range concern of relative calibration whilst maintaining high levels of efficiency when scaling to a larger number of antennas and terminals. As described herein, the reference signals transmitted for calibration and at the same power level as typical signal transmissions, and hence are better suited to capture and calibrate the distortions introduced by transmit/receive circuitry.


1.F. Reciprocity Calibration Via Receiver-Side Inversion

Let Terminal A transmit known reference signals over a subset of multi-carrier tones and P be a specific reference signal at one of these tones. For example, Terminal A may use every Mth subcarrier for reference signal transmission, where M is an integer. For example, M may be 8 or 16 in practical systems. Terminal B receives

YB=HAB·P+W  (27)

where W is additive white Gaussian noise with zero mean and variance N0. Note that the above equation is a scalar equation because the equation represents the received signal at a single subcarrier. For each subcarrier on which a reference is transmitted, there will be one such equation. Terminal B estimates HAB from YB and inverts it. To avoid singularities and cope with a large dynamic range, regularized zero forcing may be used to compute the inversion:











H
~

AB

-
1


=



H
AB
*




H
AB
*

·

H
AB


+

N
0





H
AB

-
1







(
28
)







Terminal B then transmits {tilde over (H)}AB−1 back to terminal A over the same tone. This transmission should quickly follow the first one—especially in the presence of Doppler—to ensure HCH remains relatively constant. Terminal A then receives

YB=HBA·{tilde over (H)}AB−1+W  (29)


Ignoring the noise term, which may be averaged out over multiple transmissions, it can be seen that YB is the inverse of the reciprocity calibration factor:











Y
B





H
B
TX

·

H
CH

·

H
A
RX




H
A
TX

·

H
CH

·

H
B
RX




=




H
B
TX

·

H
A
RX




H
A
TX

·

H
B
RX



=

α

-
1







(
30
)







Since these are scalar values, the inversion processing is for both HAB and YB is straightforward. Here, the inverse reciprocity calibration factor represents a ratio of circuitry channel from Terminal B to Terminal A, and a circuitry channel from Terminal A to Terminal B.


In multi-carrier systems, the above-described procedure may be repeated over multiple tones and the result interpolated to yield the full set of calibration factors over the bandwidth of interest. This full set may be obtained, for example, by averaging or interpolating the calibration factors are the subcarriers at which reference signals were transmitted. Since the Tx and Rx contributions of both terminal A and B will be relatively flat across frequency, it should be possible to use a sparse subgrid of tones with the appropriate interpolation to obtain an accurate level of calibration.


The results of the channel estimation as above may be combined with channel estimation of the HCH channel to obtain an estimate of the overall channel HAB and HBA.


2. Wireless Precoded Communication for FDD Systems

A wireless system, with a multi-antenna base-station and multiple user antennas, is shown in FIG. 5. Each transmission from a user antenna to one of the base-station antennas (or vice versa), experiences a different channel response (assuming the antennas are physically separated enough). For efficient communication, the base-station improves the users' received Signal-to-Interference-Noise-Ratio (SINR) by means of precoding. However, to precode, the base-station needs to have an accurate estimation of the downlink channels to the users during the transmission time.


In FDD systems, two directions of transmissions may use two different (typically non-overlapping) frequency bands. These transmissions may include, for example, transmissions from a network-side node such as a base station or an access point to multiple user devices, often called the downlink direction, and transmissions from the multiple user devices to the network-side node, often called the upstream or uplink direction. The various embodiments described in the present document, and the attached appendices, perform signal processing to improve communication by, for example, estimating second order statistics of a wireless channel and using the estimate to perform pre-coding in an FDD system.


In the following subsections, an efficient system and a method for predicting the downlink channel for precoding in an FDD system is described.


2.A. Second-Order Statistics Training

For simplicity, this section focuses on a single user antenna and the L base-station antennas. This can be easily extended to any number of users. The setup of the system is shown in FIG. 6. The base-station predicts from the uplink channel response, the downlink channel response in a different frequency band and Nlatency subframes later.


To achieve this, the system preforms a preliminary training phase, consisting of multiple sessions, where in each session i=1, 2, . . . , Ntraining, the following steps are taken:

    • At subframe n, the user equipment transmits reference signals (RS) in the uplink. The base-station receives them and estimate the uplink channel HUL(i) over the L base-station antennas.
    • At subframe n+Nlatency, the base-station transmits reference signals in the downlink from all its antennas. The user equipment receives it and sends it back as uplink data in a later subframe. The base-station computes the downlink channel estimation for it, HDL(i). In a different implementation, it is possible that the UE will compute the channel estimation and send it to the base-station as uplink data.
    • The base-station computes the second-order statistics

      RUL(i)=HUL(i)·(HUL(i))H
      RDL,UL(i)=HDL(i)·(HUL(i))H
      RDL(i)=HDL(i)·(HDL(i))H


Herein, (⋅)H is the Hermitian operator. For the case that the channel has non-zero-mean, both the mean and the covariance matrix should be determined.


When the training sessions are completed, the base-station averages out the second order statistics:







R
UL

=


1

N
training







i
=
1


N
training








R
UL

(
i
)











R

DL
,
UL


=


1

N
training







i
=
1


N
training








R

DL
,
UL


(
i
)











R
DL

=


1

N
training







i
=
1


N
training








R
DL

(
i
)








Then, it computes the prediction filter and the covariance of the estimation error:

Cprediction=RDL,UL·(RUL)−1
RE=RDL−Cprediction·(RDL,UL)H


The inversion of RUL may be approximated using principal component analysis techniques. We compute {λ}, the K most dominant eigenvalues of RUL, arranged in a diagonal matrix D=diag(λ1, λ2, . . . , λK) and their corresponding eigenvectors matrix V. Typically, K will be in the order of the number of reflectors along the wireless path. The covariance matrix can then be approximated by RUL≈V·D·(V)H and the inverse as RUL≈V·D−1·(V)H.


Note, that there is a limited number of training sessions and that they may be done at a very low rate (such as one every second) and therefore will not overload the system too much.


To accommodate for possible future changes in the channel response, the second-order statistics may be updated later, after the training phase is completed. It may be recomputed from scratch by initializing again new Ntraining sessions, or by gradually updating the existing statistics.


The interval at which the training step is to be repeated depends on the stationarity time of the channel, i.e., the time during which the second-order statistics stay approximately constant. This time can be chosen either to be a system-determined constant, or can be adapted to the environment. Either the base-station or the users can detect changes in the second-order statistics of the channel and initiate a new training phase. In another embodiment, the base-station may use the frequency of retransmission requests from the users to detect changes in the channel, and restart the process of computing the second-order statistics of the channel.


2.B. Scheduling a Downlink Precoded Transmission

For each subframe with a precoded downlink transmission, the base-station should schedule all the users of that transmission to send uplink reference signals Nlatency subframes before. The base-station will estimate the uplink channel responses and use it to predict the desired downlink channel responses

HDL=Cprediction·HUL


Then, the downlink channel response HDL and the prediction error covariance RE will be used for the computation of the precoder.


3. Efficient Channel Estimation Using Second-Order Statistics

In some embodiments, channel knowledge is typically acquired by transmitting known reference signals (pilots) and interpolating them at the receiver over the entire bandwidth and time of interest. Typically, the density of the pilots depends on characteristics of the channel. Higher delay spreads require more dense pilots along frequency and higher Doppler spreads require more dense pilots along time. However, the pilots are typically required to cover the entire bandwidth of interest and, in some cases, also the entire time interval of interest.


Herein, a method based on the computation of the second-order statistics of the channel, where after a training phase, the channel can be estimated over a large bandwidth from reference signals in a much smaller bandwidth is proposed. Even more, the channel can also be predicted over a future time interval.


3.A. Second-Order Statistics Training for Channel Estimation

Similar to the embodiments described in Section 2.A, a system preforms a preliminary training phase, consisting of multiple sessions, where in each session i=1, 2, . . . , Ntraining, the following steps are taken:

    • The transmitter sends reference signals to the receiver. We partition the entire bandwidth of interest into two parts BW1 and BW2, as shown in FIG. 7, where typically the size of BW1 will be smaller or equal to BW2. Note, that these two parts do not have to from a continuous bandwidth. The transmitter may send reference signals at both parts at the same time interval (FIG. 8) or at different time intervals (FIG. 9).
    • The receiver receives the reference signals and estimates the channel over their associated bandwidth, resulting in channel responses H1(i) and H2(i).
    • The receiver computes the second-order statistics of these two parts:

      R1(i)=H1(i)·(H1(i))H
      R2,1(i)=H2(i)·(H1(i))H
      R2(i)=H2(i)·(H2(i))H


Herein, (⋅)H is the Hermitian operator. For the case that the channel has non-zero-mean, both the mean and the covariance matrix should be determined, as is further discussed in the attached appendices. When the training sessions are completed, the base-station averages out the second order statistics:







R
1

=


1

N
training







i
=
1


N
training








R
1

(
i
)











R

2
,
1


=


1

N
training







i
=
1


N
training








R

2
,
1


(
i
)











R
2

=


1

N
training







i
=
1


N
training








R
2

(
i
)








Then, the receiver computes a prediction filter:

Cprediction=R2,1·(R1)−1


and if needed, also the covariance of the estimation error

RE=R2−Cprediction·(R2,1)H


The inversion of R1 may be approximated using principal component analysis techniques. We compute {λ}, the K most dominant eigenvalues of R1, arranged in a diagonal matrix D=diag(λ1, λ2, . . . , λK) and their corresponding eigenvectors matrix V. Typically, K will be in the order of the number of reflectors along the wireless path. The covariance matrix can then be approximated by R1≈V·D·(V)H and the inverse as R1−1≈V·D−1·(V)H.


Note, that there is a limited number of training sessions and that they may be done at a very low rate (such as one every second) and therefore will not overload the system too much.


To accommodate for possible future changes in the channel response, the second-order statistics may be updated later, after the training phase is completed. It may be recomputed from scratch by initializing again new Ntraining sessions, or by gradually updating the existing statistics.


The interval at which the training step is to be repeated depends on the stationarity time of the channel, i.e., the time during which the second-order statistics stay approximately constant. This time can be chosen either to be a system-determined constant, or can be adapted to the environment. Either the base-station or the users can detect changes in the second-order statistics of the channel and initiate a new training phase. In another embodiment, the base-station may use the frequency of retransmission requests from the users to detect changes in the channel, and restart the process of computing the second-order statistics of the channel.


3.B. Embodiments for Efficient Channel Estimation

After the training phase is completed, the transmitter may only send reference signals corresponding to BW1. The receiver, estimated the channel response H1 and use it to compute (and predict) and channel response H2 over BW2 using the prediction filter:

H2=Cprediction·H1



FIGS. 10 and 11 describe example two prediction scenarios (same time interval and future time interval, respectively).


4. Reciprocal Calibration for Reverse Channel Estimation

In the operation of the wireless systems described in this document, n−1 signals, intended for n−1 individual UEs, will act as interference for the target UE. A transmit pre-coder cancels the interference generated at the target UE by the n−1 un-intended signals meant for other UEs. To build a pre-coder, down link channel state information (CSI) is used.


In an extrinsic beamforming technique, CSI is fed back from the UE to BS through a feedback up-link channel. However, considerable amount of data BW is used for this, thus affecting the overall system throughput efficiency.


For Time Division Duplex (TDD) systems, the physical channel in the air (sometimes called the radio channel) is reciprocal within the channel coherence time. i.e., the up-link (UE to BS) and down-link (BS to UE) are identical (in SISO (transpose in MIMO). However, when the transceiver front-end (FE) hardware is also taken into account, channel reciprocity no longer holds. This is due to the non-symmetric characteristics of the RF hardware. It includes PA non-linearity, RF chain crosstalk, phase noise, LNA non-linearity and noise figure, carrier and clock drifts etc.


In some embodiments, a calibration mechanism can be designed to calibrate for the nonreciprocal components of the wireless link such that embodiments can estimate the down link by observing the up-link with the help of these calibration coefficients. If this is feasible, no CSI feedback is necessary (as in the case of extrinsic beam forming), thus improving the overall system throughput efficiency. The associated beamforming is also sometimes called intrinsic beamforming.



FIG. 12 shows an example block diagram of a communication channel with reciprocity. The composite wireless channel from A to B may be represented as: ĤA,B=CRX,B·HA,B·CTX,A. For the reciprocal channel, it may be assumed that HAB=λHBAT for a complex scalar (λ).


In the case of a non-reciprocal channel, with analog and RF components, Non-reciprocal analog and RF components: CTX, A, CRX, A, CRX, B, CTX, B, ideally for simplicity, it is beneficial if each matrix is a diagonal matrix. Such an embodiment may also use a design that minimizes the coupling between Tx and Rx paths. Similarly, the composite channel from B to A is given by ĤB,A=CRX, A·HB,A·CTX,B.


If all the C matrices can be estimated a priori, the BS to UE channel can be estimated from the UE to BS channel. In such a case, feeding back channel state information for transmit beamforming may not be needed, thereby making the upstream bandwidth available for data instead of having to transmit channel state information. Estimation of the C matrices may also improve system efficiency. In some embodiments disclosed herein, the reciprocity calibration may be performed by calibrating Tx and Rx of the BS and UE side during a startup or a pre-designated time. The diagonal matrices CTX, A, CRX, A, CRX, B, CTX, B may be estimated. These matrices may be re-estimated and updated periodically. The rate of change of the C matrices will typically be slow and may be related to factors such as the operating temperature of the electronics used for Tx and Rx.


4.A. Notation

In the description herein, ha1a2 denotes the channel from transmitter (TX) a1 to receiver (RX) a2. This notation is different from the conventional MIMO channel notation. In the conventional methods, this will be denoted as ha2a1)·Also, conjugate of a complex quantity is represented with a*, e.g., conj(h)=h*.


4.B. Downlink Channel Estimation from Uplink Channel and Calibration Coefficients

While the disclosed techniques are more generally applicable, for the ease of explanation, the following assumptions are made:


[1] cross talk between the TX-TX, RX-RX and TX-RX RF chains is negligible


[2] Antenna mutual coupling in negligible


[3] TX and RX at BS are working with clocks generated from the same PLL so that carrier and symbol clocks are synchronized in wireless transmissions among transceivers at A or B (and not necessarily between A and B). Here, “A” and “B” represent communication devices at two ends of a communication medium. For example, A may refer to a base station (or user equipment) and correspondingly B may refer to a user equipment (or base station). As another example, A may refer to a hub and B may refer to a remote station. Without loss of generality, the channel from A to B may be called the downlink (DL) channel and the channel from B to A may be called the uplink (UL) channel. See, e.g., FIG. 13 and FIG. 14.


[4] same assumptions as above for UE.


Inventors' measurement on some existing equipment has verified that a) the coupling between different RF paths is typically of the order of −30 dB. A careful design of the RF front end can ensure even lesser levels of cross talk. b) The isolation between the cross polarizations of the antenna is of the order of 15 to 20 dB. This means that if a signal of xdB power is sent on the vertical polarization of a cross polarized antenna, an image with (x−15) dB power will appear on the horizontal polarization. This isolation cannot be improved much even under improved antenna design. So, for the below calibration mechanism to work properly, embodiments should either use i) antenna with single polarization is used or ii) if dual polarized antenna are used, take care that simultaneous transmission on both the polarizations is never happening.


However, these assumptions can be relaxed, as described herein, and modifications in the below described calibration algorithm will be presented as well. If dual polarized antenna is the design choice, modifications to the disclosed algorithm as described herein could be used in some embodiments.


Some embodiments of a calibration algorithm are described herein for a 4×4 MIMO system. This is to keep the description simple and easy to comprehend. The same mechanism can be generalized to systems with any number of BS and UE antenna.


With reference to FIG. 13 and FIG. 14, a four antenna system 1300 at A and a four antenna system 1400 at B are depicted, respectively. For example, the configuration 1300 at A may represent a base station and the configuration 1400 at B may represent a UE (or vice versa). Let halal denote the channel from TX a1 to RX a2. It is constituted by 2 components a). The reciprocal radio channel from antenna a1 to antenna a2 and b) the non-reciprocal components at TX a1 and RX a2. Non-reciprocal components are captured in two memory-less complex scalars denoted by ta1 and ra2 corresponding to the TX and RX respectively. In this patent document, modifications when delay is involved will be demonstrated as well.


Thus, ĥa1a2 can be written as:

ĥa1a2=ta1·ha1a2·ra2  (31)


Similarly, channel between TX at a2 to RX at a1, ĥa2a1, can be written as

ĥa2a1=ta2·ha2a1·ra1  (32)


Taking the ratio between ĥa1a2 and ĥa2a1, and noting that ha1 a2 and ha2a1 are identical, the following can be written:













c

a





1

a





2


=



t

a





1


·

r

a





2





r

a





1


·

t

a





2










=


1

c

a





2

a





1



.








(
33
)







The coefficient ca1a2 is referred to as the calibration coefficient from a1 to a2. Similarly, the calibration coefficient from a2 to a3, ca2a3, can be written as:













c

a





2

a





3


=



t

a





2


·

r

a





3





r

a





2


·

t

a





3










=




t

a





1


·

r

a





3





r

a





1


·

t

a





3




·



r

a





1


·

t

a





2





t

a





1


·

r

a





2











=


c

a





1

a





3



c

a





1

a





2










(
34
)







From Eq. (34), it can be seen that the calibration coefficient between antenna a2 (TX) and antenna a3 (RX) can be written as the ratio of the calibration coefficients of the reference antenna a1 to that of a3 and a2. Similar relations can be derived for TXs and RXs at B, depicted in configuration 1400 of FIG. 14. Equations for calibration coefficients that involve antennas from both sides A and B can be derived as below.













c

a





2

b





2


=



t

a





2


·

r

b





2





t

b





2


·

r

a





2










=




t

a





1


·

r

b





1





t

b





1


·

r

a





1




·



r

a





1


·

t

a





2





t

a





1


·

r

a





2




·



t

b





1


·

r

b





2





t

b





2


·

r

b





1











=



c

a





1

b





1



c

a





1

a





2



·

c

b





1

b





2










(
35
)







From Eq. (35), it should be clear that any wireless channel between A and B can be written in terms of a) a calibration coefficient with respect to a reference antenna at A, b) a calibration coefficient with respect to a reference antenna at B, and c) a calibration coefficient between the same reference antennas at A and B.


Similarly, the downstream channels (BS-UE) can be represented in terms of the upstream channels (UE-BS) and calibration coefficients.

ĥa1b1=ca1b1·ĥb1a1
ĥa1b2=ca1b1·cb1b2·ĥb2a1 (ca1b2·ĥb2a1)
ĥa1b3=ca1b1·cb1b3·ĥb3a1 (ca1b3·ĥb3a1)
ĥa1b4=ca1b1·cb1b4·ĥb4a1  (36)


Similarly,












h
^


a





2

b





1


=



c

a





1

b





1



c

a





1

a





2



·


h
^


b





1

a





2












h
^


a





2

b





2


=



c

a





1

b





1



c

a





1

a





2



·

c

b





1

b





2


·


h
^


b





2

a





2












h
^


a





2

b





3


=



c

a





1

b





1



c

a





1

a





2



·

c

b





1

b





3


·


h
^


b





3

a





2












h
^


a





2

b





4


=



c

a





1

b





1



c

a





1

a





2



·

c

b





1

b





4


·



h
^


b





4

a





2


.







(
37
)







Furthermore












h
^


a





3

b





1


=



c

a





1

b





1



c

a





1

a





3



·


h
^


b





1

a





3












h
^


a





3

b





2


=



c

a





1

b





1



c

a





1

a





3



·

c

b





1

b





2


·


h
^


b





2

a





3












h
^


a





3

b





3


=



c

a





1

b





1



c

a





1

a





3



·

c

b





1

b





3


·


h
^


b





3

a





3













h
^


a





3

b





4


=



c

a





1

b





1



c

a





1

a





3



·

c

b





1

b





4


·


h
^


b





4

a





3




,





(
38
)







And finally












h
^


a





4

b





1


=



c

a





1

b





1



c

a





1

a





4



·


h
^


b





1

a





4












h
^


a





4

b





2


=



c

a





1

b





1



c

a





1

a





4



·

c

b





1

b





2


·


h
^


b





2

a





4












h
^


a





4

b





3


=



c

a





1

b





1



c

a





1

a





4



·

c

b





1

b





3


·


h
^


b





3

a





4












h
^


a





4

b





4


=



c

a





1

b





1



c

a





1

a





4



·

c

b





1

b





4


·


h
^


b





4

a





4








(
39
)







Using the results from Eq. (36) to Eq. (39), and denoting ca1b1 as ζ, a complex constant, the downlink MIMO channel can be expressed in terms of the uplink MIMO channel using the following equation:










[





h
^


a





1

b





1






h
^


a





2

b





1






h
^


a





3

b





1






h
^


a





4

b





1








h
^


a





1

b





2






h
^


a





2

b





2






h
^


a





3

b





2






h
^


a





4

b





2








h
^


a





1

b





3






h
^


a





2

b





3






h
^


a





3

b





3






h
^


a





4

b





3








h
^


a





1

b





4






h
^


a





2

b





4






h
^


a





3

b





4






h
^


a





4

b





4





]

=

ζ
·

[





h
^


a





1

b





1






1

c

a





1

a





2






h
^


b





1

a





2







1

c

a





1

a





2






h
^


b





1

a





3







1

c

a





1

a





2






h
^


b





1

a





4









c

b





1

b





2





h
^


b





2

a





1








c

a





1

b





1



c

a





1

a





2






h
^


b





2

a





2








c

a





1

b





1



c

a





1

a





2






h
^


b





2

a





3








c

a





1

b





1



c

a





1

a





2






h
^


b





2

a





4









c

b





1

b





3





h
^


b





3

a





1








c

a





1

b





3



c

a





1

a





2






h
^


b





3

a





2








c

a





1

b





3



c

a





1

a





2






h
^


b





3

a





3








c

a





1

b





3



c

a





1

a





2






h
^


b





3

a





4









c

b





1

b





4





h
^


b





4

a





1








c

a





1

b





4



c

a





1

a





2






h
^


b





4

a





2








c

a





1

b





4



c

a





1

a





2






h
^


b





4

a





3








c

a





1

b





4



c

a





1

a





2






h
^


b





4

a





4






]






(

40

a

)







The right hand side of Eq. (40a) can be further decomposed as:









=

ζ
·

[



1


0


0


0




0



c

b





1

b





2




0


0




0


0



c

b





1

b





3




0




0


0


0



c

b





1

b





4





]

·




[





h
^


a





1

b





1






h
^


a





1

b





2






h
^


a





1

b





3






h
^


a





1

b





4








h
^


a





2

b





1






h
^


a





2

b





2






h
^


a





2

b





3






h
^


a





2

b





4








h
^


a





3

b





1






h
^


a





3

b





2






h
^


a





3

b





3






h
^


a





3

b





4








h
^


a





4

b





1






h
^


a





4

b





2






h
^


a





4

b





3






h
^


a





4

b





4





]

·

[



1


0


0


0




0



1

c

a





1

a





2





0


0




0


0



1

c

a





1

a





3





0




0


0


0



1

c

a





1

a





3






]








(

40

b

)







Eq. (40b) is of the form ζ·KB·HU·KA. Elements in the calibration coefficient matrix KA and KB are obtained by calibrations performed at A (BS) and B (UE) and later by transferring UE coefficients to the BS. Note that calibration coefficient estimation at the BS may involve transmission and reception of calibration signals among BS antennas (a.k.a. local calibration). Similarly, estimation of calibration coefficients at UE can be performed using local transmission and reception of calibration signals among UE antennas.


In some embodiments, the TX and RX timing at a device (BS or UE) may be operated from the same PLL. This eliminates the carrier and/or clock offset impairments that is often associated with the detector at B for RF transmissions from A and vice versa. This is because A and B will be deriving all their internal clock frequencies, in general, from 2 different PLLs, one at A and another one at B. If these impairments are a part of the calibration coefficients, e.g., manifest themselves as a time varying phase rotation, then the coefficients will vary more frequently due to time varying carrier or clock errors in addition to its own time variability. Since KA and KB are obtained from measurements exclusively at BS or UE (and not using transmissions from BS to UE or vice versa), they vary relatively slowly.


Local calibrations, thus, are generally stable and do not change much over a period of several minutes (e.g., 30 minutes).


Eq. (40b) further reveals that the reverse MIMO channel can be mathematically modelled as the composition of a) a complex scalar, b) calibration coefficients at B, c) MIMO uplink channel transfer function, and d) calibration coefficients at A.


A pre-coder can be built at A (or at B), by acquiring calibration coefficients of the receiver side. The pre-coder implementation could be based on any of several pre-coders available in the literature. Some examples include: an MMSE pre-coder, a regularized MMSE precoder, a zero forcing precoder, a Tomlinson-Harashima pre-coder, and so on. As will be appreciated by one of skill in the art, either linear pre-coders (first three examples above) or non-linear pre-coders (the last example above) may be used.


To illustrate this point, an example embodiment of a reciprocity based zero-forcing (ZF) pre-coder using the data above (ref. Eq. 40b) is described below. The configuration of this pre-coder is depicted in FIG. 15. Note that, the same set of data can be used to design any type of pre-coder.


A ZF pre-coder may have the following form.









W
=



K
A

-
1


·

H
U

-
1


·

K
B

-
1




E


{





K
A

·

H
U



·
B




F

}







(

40

c

)







Here, ∥⋅∥F is the Frobenious norm of a given matrix. Note that ζ, the calibration coefficient between A and B does not appear in the pre-coder. The effect of ζ will be counteracted by the equalizer (single tap or MIMO) at the UEs.


Some embodiments use the fact that KB and KA are typically slow-varying, so that their inverse can be implemented in software (instead of implementing in hardware). In some embodiments, W2 (HU−1 or inverse of HU) may be fast varying and implemented in hardware circuits.


In some embodiments, W2 may be obtained as a by-product of the receiver Equalizer at the BS. For example BS equalizer often implements a variant of the HU−1


Therefore, using the techniques disclosed herein, some embodiments may model the downlink channel as a composition of the uplink channel and slow varying local coefficients. This enables to build a variety of different pre-coders with minimal-feedback overhead, and certain linear precoders, use receiver equalizer computations performed at the base station.


Estimation of ζ or ca1b1 involves calibrating across BS and UE. For the reasons discussed above, these coefficients can be frequently varying and the estimation and feedback of these coefficients could consume a lot of bandwidth. Advantageously, a transmit-side pre-coder to cancel the multi-user interference can be designed without the knowledge of ζ. It can be designed from the upstream channel measurements KA and KB, as described in the present document.


Several estimation methods to determine the calibration coefficients are described in the present document. When the number of antenna is relatively small, the method described in Section 4.C may be used. When a large number of antennas are involved, the method described in Section 4.D may be used.


4.C. Estimation of Calibration Coefficients (Cx1x2)—Method 1—Iterative Algorithm

When the number of antenna is relatively small, such as 4 as in this example, an iterative algorithm can be used to compute the calibration coefficients. This, however, will entail a large amount of calibrations in a massive MIMO scenario and may become impractical. In a massive MIMO scenario, estimation described in sec. III may be used.


Iteration 1:

ca1a3=ca1a2˜ca2a3  (41)
ca1a4=ca1a2·ca2a4  (42)


Make an initial estimate of ca1a2, ca1a3, ca1a3, as described herein. Also estimate ca2a3 and ca2a4. Using the latter, make an alternative estimate of ca1a3 and ca1a4 (eq. 41 and 42). From these two independent estimates of ca1a3 and ca1a4, make a refined estimate of ca1a3 and ca1a4 (e.g., the arithmetic mean). A refined estimate of ca1a2 can be obtained from Eq. (41) or (42)

ca1a4=ca1a3·ca3a4  (43)


A refined estimate of ca1a4 can be obtained (as described above), from a measurement of ca3a4 using Eq. (43). Using the updated ca1a4, estimates of ca1a2 and ca1a3 can be updated.


More iterations, in the above manner, can be performed to further refine the above coefficients.


The number of calibrations to be performed in this method is 4C2=6. Number of calibrations will grow exponentially in N, the number of antenna; hence not practical for large array of antenna.


4.D. Estimation of Calibration Coefficients (Cx1x2)—Method 2—Large Number of Antenna

When the number of antenna is high, e.g., 64, Method 1 is not practically feasible. This section describes a simple method for the same.


The following equations have been previously seen derived.

ĥa1a2=ca1a2ĥa2a1  (44)
ya2a2a1sa1+na2  (45)
ya1a2a1sa2+na1  (46)


where sa1, sa2 are known training symbols, ya1 and ya2 are the received signal and n is additive noise. It is assumed that the noise is white in the frequency band of interest. The LS estimates of








h
^


a





1

a





2


=





y

a





2


·

s

a





1

*






s

a





1




2







and







h
^


a





2

a





1



=




y

a





1


·

s

a





2

*






s

a





2




2


.







Note that the estimation error is







CN


(

0
,


σ
n
2





s
b



2



)


.





Now







c

a





1

a





2


=




h
^


a





1

a





2




h
^


a





2

a





1



.






Similar expressions hold good for other calibration coefficients.


For a 64 antenna system, this method would use a total of 64 calibration signal transmissions.


4.E. Estimation of Calibration Coefficients (Cx1x2)—Method 3—Using Total Least Squares (TLS)

What is described below is an algorithm to estimate the calibration coefficients using the method of total least squares.


Refer to Method 2. Take k (say 4) such LS estimates and form the following matrix

Ha1a2=[ĥa1a2,1a1a2,2 ⋅ ⋅ ⋅ ĥa1a2,4]  (47)
Ha2a1=[ĥa2a1,1a2a1,2 ⋅ ⋅ ⋅ ĥa2a1,4]  (48)


Let D=[Ha2a1, Ha1a2]. Note that D is a 4×2 matrix in this example.








c
^


a





1

a





2


=



argmin


Δ






H

a





1

a





2



,

Δ






H

a





2

a





1










Δ






H

a





2

a





1





F
2


+




Δ






H

a





1

a





2





F
2







subject to Ha1a2+ΔHa1a2=(Ha2a1+ΔHa2a1)·ca1a2.


Solution for the above is obtained as below.


Let D=UΣVH where Σ=diag (σ1, σ2), σ12, where σs are singular values of D; U and V are the left and right eigen-vectors respectively. Let









V
=


[




v
11




v
12






v
21




v
22




]

.





(
49
)







To have a solution for this Total Least Squares problem, υ22≠0. Only if σ1≠σ2, the solution will be unique.


The solution is given by,








c
^



a





1

a





2

-
opt


=



-

v
12



v
22


.





In some embodiments, total least squares error optimization criteria may be used. Total least squares assumes that both h.(..) and c.(..) are subject to estimation errors. This is more realistic compared to assuming that h.(..) is impervious to estimation errors.


Note that three different methods for estimating the calibration coefficients are described above. During the implementation/prototyping phase, effectiveness of each method can be evaluated separately and the one with the best merit need to be selected for the final implementation. Alternatively, the decision may be made based on the number of antennas for which the calibration is performed. For example, method 1 may be suitable for up to 8 receive and/or transmit antennas, while methods 2 or 3 may be used for higher number of antennas.


4.F. Effect of Delays in the TX/RX Path

Let delays associated with ta1, ta2, ra1 and ra2 be τta2, τta2, τra1, τra2. Therefore,













c

a





1

a





2


=






t

a





1


·

r

a





2









r

a





1


·

t

a





2






·



e


-
j






2

π






fr

a





2




·

e


-
j






2

π






fr

a





1







e


-
j






2

π






fr

a





1




·

e


-
j






2

π






fr

a





2













=




c

a





1

a





2




·

e


-
j






2

π





f





δ





r










(
50
)







where δτ=(τta1ra2)−(τta2ra1), the delay difference between the a1 and a2. It is clear that, the effect of δτ is to cause a linear phase on the calibration coefficients across the frequency.


Calibration coefficient as can be seen is a function of δτ. The delays in the digital sections of TX and RX will not affect the calibration coefficients, as these will cancel out in the expression above (Note that digital delays are identical for all transmitters. This is true for receivers as well). It is clear that, a delay effect, as described above, will be manifested if the synchronizing clocks at BS or UE are routed using long traces.


It has been assumed that ta1 etc. is memory-less. The case were ta1 for example is say a 2 tap filter is not considered in this formulation.


4.G. Calibration Coefficients when there is Mutual Coupling Between TX/TX, TX/RX Paths

Non-negligible coupling between different TX/RX path can exist due to a) Imperfect hardware design or more importantly b) if cross-polarized antenna is used in the design. Modifications to the calibration algorithm under these conditions is described below. It is assumed that there is no coupling between the TX and RX paths.


Referring to FIG. 13 and FIG. 14, assume that the receive paths have mutual coupling between them. For e.g., ra1a2 is the mutual coupling between receive paths r1 and r2. Similarly ra1a3 refers to the coupling between the receive paths r1 and r3. Similar explanation applies to mutual coupling between the transmit paths. Note that ra1a1 refers to ra1 of the Section 4.B.


Assume that antenna a1 transmits a calibration sequence, sa1. It has been received by receive antenna a2, a3 and a4. Similar to the formulation in Sections 4.B. and 4.C., the following expressions can be derived.










[




y

a





2








y

a





3








y

a





4






]

=




[




r

a





2

a





2





r

a





2

a





3





r

a





2

a





4







r

a





3

a





2





r

a





3

a





3





r

a





3

a





4







r

a





4

a





2





r

a





4

a





3





r

a





4

a





4





]



[




h

a





1

a





2







h

a





1

a





3







h

a





1

a





4





]





t

a





1


·

s

a





1




+

[




n

a





2







n

a





3







n

a





4





]






(
51
)















[




y

a





1

n






y

a





3

n






y

a





4

n




]

=




[




r

a





1

a





1





r

a





1

a





3





r

a





1

a





4







r

a





3

a





1





r

a





3

a





3





r

a





3

a





4







r

a





4

a





1





r

a





4

a





3





r

a





4

a





4





]



[




h

a





2

a





1







h

a





2

a





3







h

a





2

a





4





]





t

a





2


·

s

a





2




+

[




n

a





1







n

a





3







n

a





4





]






(
52
)







[




y

a





1

m






y

a





2

m






y

a





4

m




]

=




[




r

a





1

a





1





r

a





1

a





2





r

a





1

a





4







r

a





2

a





1





r

a





2

a2





r

a





2

a





4







r

a





4

a





1





r

a





4

a2





r

a





4

a





4





]



[




h

a





3

a





1







h

a





3

a





2







h

a





3

a





4





]





t

a





3


·

s

a





3




+

[




n

a





1







n

a





2







n

a





4





]






(
53
)







[




y

a





1

mn






y

a





2

mn






y

a





3

mn




]

=




[




r

a





1

a





1





r

a





1

a





2





r

a





1

a





3







r

a





2

a





1





r

a





2

a





2





r

a





2

a





3







r

a





3

a





1





r

a





3

a





2





r

a





3

a





2





]



[




h

a





4

a





1







h

a





4

a





2







h

a





4

a





3





]





t

a





4


·

s

a





4




+

[




n

a





1







n

a





2







n

a





3





]






(
54
)







Proceeding as in Section 4.D., LS estimates of ĥa1a2, ĥa2a3 and ĥa1a4 can be obtained from Eq. (51). The term ĥa1a2 is the non-reciprocal (baseband) channel from antenna a1 to antenna a2, whereas ha1a2 is the reciprocal radio channel from antenna a1 to antenna a2g. Similarly, LS estimates ĥa2a1{circumflex over ( )}, ĥa3a1{circumflex over ( )} and ĥa1a1{circumflex over ( )} can be obtained from equations 52, 53 and 54.










[





h
^


a





1

a





2








h
^


a





1

a





3








h
^


a





1

a





4





]

=


t

a





1


·


[




r

a





2

a





2





r

a





2

a





3





r

a





2

a





4







r

a





3

a





2





r

a





3

a





3





r

a





3

a





4







r

a





4

a





2





r

a





4

a





3





r

a





4

a





4





]



[




h

a





1

a





2







h

a





1

a





3







h

a





1

a





4





]







(
55
)







[





h
^


a





2

a





1








h
^


a





3

a





1








h
^


a





4

a





1





]

=


[





r

a





1

a





1


·

t

a





2







r

a





1

a





3


·

t

a





2







r

a





1

a





4


·

t

a





2









r

a





1

a





1


·

t

a





3







r

a





1

a





2


·

t

a





3







r

a





1

a





4


·

t

a





3









r

a





1

a





1


·

t

a





4







r

a





1

a





2


·

t

a





4







r

a





1

a





3


·

t

a





4






]



[




h

a





2

a





1







h

a





3

a





1







h

a





4

a





1





]






(
56
)







From Eq. 55 and 56, we can write the following










[





h
^


a





2

a





1








h
^


a





3

a





1








h
^


a





4

a





1





]

=




1

t

a





1



·



[





r

a





1

a





1


·

t

a





2







r

a





1

a





3


·

t

a





2







r

a





1

a





4


·

t

a





2









r

a





1

a





1


·

t

a





3







r

a





1

a





2


·

t

a





3







r

a





1

a





4


·

t

a





3









r

a





1

a





1


·

t

a





4







r

a





1

a





2


·

t

a





4







r

a





1

a





3


·

t

a





4






]



[




r

a





2

a





2





r


a





2

a3










r

a





2

a





4







r

a





3

a





2





r

a





3

a





3





r

a





3

a





4







r

a





4

a





2





r

a





4

a





3





r

a





4

a





4





]



-
1





c






[





h
^


a





1

a





2








h
^


a





1

a





3








h
^


a





1

a





4





]







(
57
)







where C is the calibration coefficient to be estimated. It is interesting to note that when the mutual coupling coefficients are set to 0, the above equation gives raise to eq. 33. Similarly ĥa2a3 ĥa2a4 and ĥa3a4 can be expressed in terms of the product of a calibration matrix and ĥa3a2, ĥa4a2 and ĥa4a3.


These calibration coefficients can be computed using Total Least Squares method (method 3).


Note that there is no mutual coupling between sides A and B. This enables embodiments to write downstream link exactly as in Eq. 40c. However, in this case, KA and KB will no longer be diagonal matrices.


5. Embodiments and Method for the Disclosed Technology


FIG. 16 is a block diagram representation of a wireless hardware platform 1600 which may be used to implement the various methods described in the present document. The hardware platform 1600 may be incorporated within a base station or a user device. The hardware platform 1600 includes a processor 1602, a memory 1604 and a transceiver circuitry 1606. The processor 1602 may execute instructions, e. g., by reading from the memory 1604, and control the operation of the transceiver circuitry 1606 and the hardware platform 1600 to perform the methods described herein. In some embodiments, the memory, processor and transceiver may be included in a single chip solution.



FIG. 17 shows an example of a wireless communication system 1700 in which a transmitter device 1702 transmits signals to a receiver 1704. The signals may undergo various wireless channels and multipaths, as depicted. Some reflectors such as buildings and trees may be static, while others such as cars, may be moving scatterers. The transmitter device 1702 may be, for example, a mobile phone, a tablet, a computer, or another Internet of Things (IoT) device such as a smartwatch, a camera, and so on. The receiver device 1704 may be a network device such as the base station. The signals transmitted from the base station to the transmitter 1702 may experience similar channel degradations produced by static or moving scatterers. The techniques described in the present document may be implemented by the devices in the wireless communication system 1700.



FIG. 18 is a flowchart for an example method 1800 for wireless communication. The method 1800 may be implemented by a base station in a wireless network. The method 1800 includes receiving (1802), by a first wireless device (e.g., the base station) during a training phase, reference tones using a first number of resource elements from a transmitter of a second wireless device, wherein the first wireless device comprises multiple receiving antennas, estimating (1804), by the first wireless device, from the receiving the reference tones, a second order statistics of wireless channels between the multiple receiving antennas and the transmitter of the second wireless device, and performing (1806) channel estimation, during an operational phase subsequent to the training phase, using the second order statistics and reference tones received on a second number of resource elements, wherein the second number is less than the first number.


As previously described, the training phase may be based on a pre-defined number of repetitions of the reference signal transmission, or may be dependent on passage of a time, such as 10 to 20 seconds, during which the training is performed. The estimation operation 1804 may be performed by solving for the various equations described herein.


Upon the base station deciding that the training phase is over, then the base station may begin operating in the operational phase. The transition from the training phase to the operational phase may be communicated to user devices via an over-the-air message transmission. Alternative, the user devices and the base station may keep track of number of reference signal transmissions performed, and enter operational phase after a pre-defined threshold is exceeded.


As described herein, during operational phase, a reduced number of resource elements may be used for reference signal transmissions. In some embodiments, this number may be selected to be a fixed number that is known to both the transmitter and the receiver. In general, this number is less (or far less) than the first number of reference tones used during the training phase. As previously discussed, preferably sufficient reference tones are transmitted and received so that the essential parameters of the channel can be extracted.


The second order statistics may be measured in one of two different techniques as described here. In one technique, called the averaging method, the channel matrix is formulated as an average of autocovariance of an estimated channel response over all the receiving antennas of the first wireless device, at different time instances. In another technique, called, the direct method, includes estimating covariance of a channel matrix using a direct method in which the channel matrix is formulated to comprise a number of columns, wherein each column comprises an estimated channel response over all the receiving antennas of the first wireless device at a time instance. Further details are described with reference to equations (1) to (3).


During the operational phase, the method 700 may include performing channel estimation using an interpolation filter that interpolates from a decimated version of the channel estimate onto the entire channel estimate. Some example embodiments are provided with respect to Equations (5) to (9).


The method 1800 may further include estimating an error covariance matrix representative of an error in the estimated channel response over all of the receiving antennas computed by interpolating channel estimates; and selectively revising, for future use, the second number of resource elements of reference tones using a measure of the error covariance matrix. The measure of error covariance may be the comparison between a sum of squares of diagonal entries with a threshold, as described with respect to Equation (10) and (11).



FIG. 19 is a flowchart for an example method 1900 of wireless communication. The method 1900 may be implemented by a user device (user equipment) in a wireless system. The method includes transmitting (1902), to multiple receive antennas of a first wireless device from a transmit antenna of a second wireless device, during a training phase, reference tones using a first number of resource elements of a wireless channel between the transmit antenna and the multiple receive antennas, receiving (1904), at an end of the training phase, an instruction to transmit reference tones using a second number of resource elements, and transmitting (1906), during an operational phase after the training phase, reference tones to the multiple receive antennas of the first wireless device using the second number of resource elements, wherein the second number is different from the first number and wherein the second number is based on an estimated second order statistics of the wireless channel.


In some embodiments, the method 1900 further includes receiving, during the operational phase, another instruction to transmit reference tones using a third number of resource elements, wherein the third number is different from the first number and the second number, and transmitting, after receiving the another instruction, reference tones to the first wireless device using the third number of resource elements. For example, the base station may decide that the reduced number of reference tones may have to be adjusted upwards (if seeing too many errors, see equations (9) and (10)) or adjusted lower when the number of reflectors in the channel is less than what was observed in the training phase. As previously described, the lower number of reference signal overhead can improve the efficiency of a wireless channel by making greater bandwidth for data traffic.



FIG. 20 discloses another method 2000 of wireless communication. The method 2000 includes estimating (2002), during a training phase, a second order statistics for a first wireless channel and a second wireless channel between a transmitter and a receiver comprising multiple antennas, wherein the second order statistics is estimated using reference tones transmitted on a first number of resource elements, predicting (2004), during an operational phase subsequent to the training phase, an estimate of the second wireless channel based on the second order statistics and an estimate of the first wireless channel calculated using reference tones transmitted on a second number of resource elements, where the second number is less than the first number, and communicating (2006), during the operational phase, over the second wireless channel using the estimate of the second wireless channel resulting from the predicting, wherein the first wireless channel and the second wireless channel include non-overlapping frequencies.


In some embodiments, the second order statistics may include a cross-covariance between an estimate of the first wireless channel and an estimate of the second wireless channel. The estimates of the first and second wireless channels may be obtained during different time periods (e.g., training phases that are offset in time with respect to each other). Sections 2 and 3 have provided additional details used by the method 2000.


In the methods 1800, 1900 and 2000, the resource elements may represent time slots or subcarriers. The reference tone transmission during the training phase may occur on non-contiguous frequencies.


It will be appreciated that technique for reducing the amount of transmission resources used for reference signal transmissions are disclosed. Using the disclosed techniques, e.g., determination of second order statistics, calculation of an interpolation filter, and interpolating based on the interpolation filter, can be used to reduce the number of pilot transmissions after an initial training phase in which an estimate of the channel and its second order statistics are obtained. The reduced number of transmissions is sufficient as long as these transmissions provide sufficient information via diversity in angle of arrival, number of transmissions and resources used for transmissions such that the receiver is able to discern the wireless channel characteristics, e.g., dominant reflectors in the channel. This reduced set of number may be fixed a priori (e.g., based on a general knowledge of the expected channel under which the wireless communication is to operate) or may be changed from time to time based on the disclosed error measurements. As an illustrative example, if a wireless channel is characterized by 8 dominant reflectors, it may be sufficient to send a single reference tone to 8 receive antennas of a base station from a transmit antenna (e.g., a user device) to estimate and predict the entire channel both in the uplink and in the downlink directions, and at future transmission times.



FIG. 21 is a flowchart for an example method 600 of wireless communication. At 2102, a first communication device receives a number of subcarriers from a second communication device. Each subcarrier includes a corresponding reference signal. For example, the first communication device may be Terminal B in FIG. 3B. In some embodiments, the first communication device may be a user equipment (UE) or a wireless terminal that is incorporated into a phone, camera, laptop, tablet, Internet of Things (IoT) device or another device capable of received and transmitting wireless signals. In some embodiments, the second communication device may be the above-discussed Terminal A. In some embodiments, the second communication device may be a base station, an access point, Hub 102 (in FIG. 1), macro tower 202 (in FIG. 2), or another network-side device that is providing communication connectivity to multiple devices such as the first communication device.


In some embodiments, the receiving operation 2102 may occur in a single time slot of a TDD communication system. Alternatively, the reference signal transmissions over the subcarriers may be received over multiple time slots. Similarly, in some implementations, multiple reference signal transmissions on a same subcarrier may be received over multiple timeslots and the results may be averaged to reduce effect of noise on the various calculations described herein.


At 2102, the first communication device may calculate an inversion factor for each subcarrier based on the received value of the reference signal received on respective subcarrier. As described in the present document, e.g., equation (28), a zero forcing technique to ensure numerical stability of the calculation may be used. The regularized zero forcing technique, for example, eliminates singularity caused by zero or low valued denominator of Equation (28).


At 2106, the first communication device may transmit pilot signals back to the second communication device. The pilot signals are scaled by the inversion factor calculated for the subcarrier on which the pilot is being transmitted.


As described in the present document, one of the advantageous aspect of the method 2100 is that each inversion factor is a scalar. Inversion factors may be rational numbers or, in general, complex numbers representing both a scale and a phase shift for the subcarrier.


In some embodiments, the pilot signal and the reference signals may both be unit signals. For example, the pilot and reference signals may simply multiply a pre-determined subcarrier level by a factor of 1. Alternatively, other scaling factors may be used. Similarly, the reference signal may also simply multiply a nominal subcarrier level by a scale factor of 1. Other scale factors may be used for multiplying a nominal subcarrier value both to the first communication device, and from the subcarrier device. These scale factors may, for example, be selected such that the product of “to” and “from” scale factors for nominal subcarrier signal level scaling is 1.


As previously described, averaging may be performed by the first communication device, by repeating inversion factor calculations for subcarriers over multiple receptions.



FIG. 22 is a flowchart of an example method 2200 of wireless communication. The method 2200 may transmit reference signals and precoded signals to a first communication device from a second communication device. The first and the second communication devices may be embodied, for example, as described with respect to method 2100.


At 2202, the second communication device transmits reference signals using a number of subcarriers. The reference signals may simply scale by unity a nominal signal level of the subcarrier. Other scale factors may be used, as described with respect to method 2100. The reference signal transmissions may be performed in a single time slot. Alternatively, multiple reference signal transmissions may be performed over multiple time slots. The reference signal transmissions may be repeated as often as needed, e.g., based on the effectiveness of the subsequent pre-coding to achieve or maintain low error communication.


At 2204, the second communication device receives inversion factors in the form or scaled pilot signals from the first communication device, where the subcarriers that were used for reference signal transmissions are scaled by the corresponding inversion factor calculated by the first communication device.


At 2206, the second communication device may estimate the communication channel from the second communication device to the first communication device by generating values on all the subcarriers based on the inversion factors that were received on the subset of the subcarriers. For example, simple linear interpolation may be used for reduced computational complexity.


After an estimate of the communication channel based on the inversion factors is obtained, the second communication device may then computer a full channel response that also includes channel response of the reflectors. In some embodiments, second order statistics of the effect of the reflectors may be used for obtaining the full channel response.


In some embodiments, the first communication device is Terminal B described in the present document and depicted in FIG. 3B. Typically, multiple Terminal Bs may be present in a wireless communication system, with a single Terminal A, and possibly a secondary base station, will provide wireless access to multiple Terminal B's. In some embodiments, the second communication device is Terminal A described in the present document and depicted in FIG. 3B. With respect to the methods 600 and 700, the sparse set of subcarriers over which the reference signals (and subsequently pilot signals) may be transmitted may include every Mth subcarrier, such as every 8th or every 16th subcarrier. In some embodiments, a non-uniform group of subcarriers may be used.


In some embodiments, a wireless communication device (e.g., 1600) may include a processor and a transceiver circuitry such that the processor controls the operation of inversion factor calculations and the transceiver circuitry is used to transmit or receive the reference signals and the pilot signals.


It will be appreciated by one of skill in the art that the described receiver-side inversion is a novel approach for computing reciprocity calibration factors that has low feedback overhead and more favorable dynamic range requirements compared to existing solutions. For a given tone in the frequency domain, only one symbol transmission in each direction is required to obtain the reciprocity calibration factor, which greatly reduces overhead compared to feeding back channels as data samples. As such, it is well suited to enabling efficient bidirectional transmission in TDD wireless systems operating with a large number of both BS antennas and CPEs.


The following listing of examples provide embodiments that can addressed the technical problems described in the present document, among other problems.


1. A wireless communication method, comprising: receiving, by a first wireless device during a training phase, reference tones using a first number of resource elements from a transmitter of a second wireless device, wherein the first wireless device comprises multiple receiving antennas; estimating, by the first wireless device, from the receiving the reference tones, a second order statistics of wireless channels between the multiple receiving antennas and the transmitter of the second wireless device; and performing channel estimation, during an operational phase subsequent to the training phase, using the second order statistics and reference tones received on a second number of resource elements, wherein the second number is less than the first number.


2. The method of example 1, wherein the estimating the second order statistics includes estimating covariance of a channel matrix using a direct method in which the channel matrix is formulated to comprise a number of columns, wherein each column comprises an estimated channel response over all the receiving antennas of the first wireless device at a time instance.


3. The method of example 1, wherein the estimating the second order statistics includes estimating covariance of a channel matrix using an averaging method in which the channel matrix is formulated as an average of autocovariance of an estimated channel response over all the receiving antennas of the first wireless device, at different time instances.


4. The method of example 1, wherein the performing the channel estimation includes: performing the channel estimation by interpolating channel estimates from the reference tones received on the second number of resource elements using an interpolation filter.


5. The method of example 1, further comprising: determining, using the second order statistics, an interpolation filter.


6. The method of example 5, wherein the determining the interpolation filter includes using a decimated version of a covariance matrix calculated during the training phase.


7. The method of example 6, wherein the determining the interpolation filter is performed by: determining an eigenvector factorization of the covariance matrix, wherein the eigenvector factorization is represented as: RHH≈V·D·V*, wherein RHH is the covariance matrix, V is an eigenvectors matrix and D is a diagonal matrix of eigenvalues, and * is a conjugate transpose operation; and estimating the interpolation filter as function of the eigenvectors matrix and a decimated version of the eigenvectors matrix.


8. The method of example 7, wherein the function is represented as C=V·(VM)−1 where C represents the interpolation filter, VM is the decimated version of V.


9. The method of example 5, further including: estimating an error covariance matrix representative of an error in the estimated channel response over all of the receiving antennas computed by interpolating channel estimates; and selectively revising, for future use, the second number of resource elements of reference tones using a measure of the error covariance matrix.


10. The method of example 9, wherein the measure of the error covariance matrix includes a comparison between a sum of squares of diagonal entries of the error covariance matrix and a threshold.


11. A method of wireless communication, comprising: transmitting, to multiple receive antennas of a first wireless device from a transmit antenna of a second wireless device, during a training phase, reference tones using a first number of resource elements of a wireless channel between the transmit antenna and the multiple receive antennas; receiving, at an end of the training phase, an instruction to transmit reference tones using a second number of resource elements; and transmitting, during an operational phase after the training phase, reference tones to the multiple receive antennas of the first wireless device using the second number of resource elements, wherein the second number is different from the first number and wherein the second number is based on an estimated second order statistics of the wireless channel.


12. The method of example 11, further including: receiving, during the operational phase, another instruction to transmit reference tones using a third number of resource elements, wherein the third number is different from the first number and the second number; and transmitting, after receiving the another instruction, reference tones to the first wireless device using the third number of resource elements.


13. The method of any of examples 1 to 12, wherein the resource elements comprise subcarriers.


14. The method of any of examples 1 to 12, wherein the resource elements comprise time slots.


15. The method of any of examples 1 to 14, wherein the second number of resource elements is equal to 1.


16. A wireless communication method, implementable by a wireless communication apparatus, comprising: estimating, during a training phase, a second order statistics for a first wireless channel and a second wireless channel between a transmitter and a receiver comprising multiple antennas, wherein the second order statistics is estimated using reference tones transmitted on a first number of resource elements; predicting, during an operational phase subsequent to the training phase, an estimate of the second wireless channel based on the second order statistics and an estimate of the first wireless channel calculated using reference tones transmitted on a second number of resource elements, where the second number is less than the first number; and communicating, during the operational phase, over the second wireless channel using the estimate of the second wireless channel resulting from the predicting; wherein the first wireless channel and the second wireless channel include non-overlapping frequencies.


17. The method of example 16, wherein the second order statistics comprises a cross-covariance between an estimate of the first wireless channel and an estimate of the second wireless channel.


18. The method of example 17, wherein the estimate of the first wireless channel and the estimate of the second wireless channel are obtained at different times.


19. The method of example 16 or 17, wherein the method is implemented by a base station in a frequency division duplexing (FDD) network.


20. The method of example 16 or 17, wherein the method is implemented by a base station in a time division duplexing (TDD) network.


21. The method of any of examples 1 to 20, wherein the reference tones during the training phase are transmitted over a non-contiguous spectrum.


22. A method of wireless communication, comprising: receiving, by a first communication device, a number of subcarriers from a second communication device, each including a corresponding reference signal; calculating an inversion factor for each subcarrier based on a received value of the corresponding reference signal; and transmitting by the first communication device to the second communication device, at least some of the subcarriers by scaling a pilot signal using a corresponding inversion factor.


23. The method of example 22, wherein the subcarriers are received in a single time slot.


24. The method of example 22, wherein the subcarriers are received over multiple time slots.


25. The method of any of examples 22 to 24, wherein the calculating the inversion factor includes using a regularized zero forcing technique, thereby avoiding singularities in calculations.


26. The method of example 22, wherein the inversion factor is a complex number.


27. The method of any of examples 22 to 26, wherein the scaling the pilot signal includes multiplying the pilot signal by the inversion factor.


28. The method of example 22, wherein the first communication device is a user terminal in a wireless network and the second communication device is a network device in the wireless network.


29. The method of example 25, wherein the calculating the inversion factor includes evaluating:









H
~

AB

-
1


=


H
AB
*




H
AB
*

·

H
AB


+

N
0




,





wherein {tilde over (H)}AB−1 represents the inversion factor, HAB* represents complex conjugate of the received value, and N0 represents noise variance in received reference signals.


30. The method of example 29, wherein the calculating is repeated and averaged over multiple received reference signal transmissions for each subcarrier.


31. The method of any of examples 22 to 29, wherein the reference signal and the pilot signal are inverse functions of each other.


32. A method of wireless communication, comprising: transmitting, to a first communication device, from a second communication device, a number of subcarriers, each subcarrier including a corresponding reference signal; receiving, from the first communication device, at least some of the subcarriers carrying pilot signals scaled by inversions factors for the subcarriers; and estimating a communication channel between the second communication device to the first communication device using the inversion factors.


33. The method of example 32, wherein the estimating the communication channel includes interpolating inversion factors at intermediate subcarriers for which no inversion factors were received from the first communication device.


34. The method of example 32, wherein the corresponding reference signals transmitted on each subcarriers are identical.


35. The method of example 32, further including: performing a subsequent transmission from the second communication device to the first communication device by pre-coding using a result of the estimating the communication channel.


36. The method of example 35, wherein the pre-coding includes Tomlison-Harashima precoding.


37. The method of example 32, wherein the subcarriers on which reference signals are transmitted include every Mth subcarrier of the communication channel, where M is an integer greater than 1.


38. The method of example 32, wherein the estimating the communication channel further includes estimating contributions of reflectors to the communication channel using second order statistics.


39. The method of example 32, wherein the inversion factors are complex scalar numbers.


40. A wireless communication apparatus comprising: a processor; a non-transitory memory; and a wireless transceiver, wherein a method recited in one or more of examples 1 to 39 is stored in the non-transitory memory, and wherein the processor is configured to perform the method using the wireless transceiver for transmitting or receiving signals.


41. The wireless communication apparatus of example 22, wherein the wireless communication apparatus is a base station of a multi-user multi-input multi-output (MU-MIMO) wireless system.


The disclosed and other embodiments, modules and the functional operations described in this document can be implemented in digital electronic circuitry, or in computer software, firmware, or hardware, including the structures disclosed in this document and their structural equivalents, or in combinations of one or more of them. The disclosed and other embodiments can be implemented as one or more computer program products, i.e., one or more modules of computer program instructions encoded on a computer readable medium for execution by, or to control the operation of, data processing apparatus. The computer readable medium can be a machine-readable storage device, a machine-readable storage substrate, a memory device, a composition of matter effecting a machine-readable propagated signal, or a combination of one or more them. The term “data processing apparatus” encompasses all apparatus, devices, and machines for processing data, including by way of example a programmable processor, a computer, or multiple processors or computers. The apparatus can include, in addition to hardware, code that creates an execution environment for the computer program in question, e.g., code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of one or more of them. A propagated signal is an artificially generated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal, that is generated to encode information for transmission to suitable receiver apparatus.


A computer program (also known as a program, software, software application, script, or code) can be written in any form of programming language, including compiled or interpreted languages, and it can be deployed in any form, including as a standalone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program does not necessarily correspond to a file in a file system. A program can be stored in a portion of a file that holds other programs or data (e.g., one or more scripts stored in a markup language document), in a single file dedicated to the program in question, or in multiple coordinated files (e.g., files that store one or more modules, sub programs, or portions of code). A computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network.


The processes and logic flows described in this document can be performed by one or more programmable processors executing one or more computer programs to perform functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application specific integrated circuit).


Processors suitable for the execution of a computer program include, by way of example, both general and special purpose microprocessors, and any one or more processors of any kind of digital computer. Generally, a processor will receive instructions and data from a read-only memory or a random access memory or both. The essential elements of a computer are a processor for performing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto optical disks, or optical disks. However, a computer need not have such devices. Computer readable media suitable for storing computer program instructions and data include all forms of non-volatile memory, media and memory devices, including by way of example semiconductor memory devices, e.g., EPROM, EEPROM, and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks; magneto optical disks; and CD ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.


While this patent document contains many specifics, these should not be construed as limitations on the scope of an invention that is claimed or of what may be claimed, but rather as descriptions of features specific to particular embodiments. Certain features that are described in this document in the context of separate embodiments can also be implemented in combination in a single embodiment. Conversely, various features that are described in the context of a single embodiment can also be implemented in multiple embodiments separately or in any suitable sub-combination. Moreover, although features may be described above as acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a sub-combination or a variation of a sub-combination. Similarly, while operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results.


Only a few examples and implementations are disclosed. Variations, modifications, and enhancements to the described examples and implementations and other implementations can be made based on what is disclosed.

Claims
  • 1. A wireless communication method, comprising: receiving, by a first wireless device during a training phase, reference tones using a first number of resource elements from a transmitter of a second wireless device, wherein the first wireless device comprises multiple receiving antennas;estimating, by the first wireless device, from the receiving the reference tones, a second order statistics of wireless channels between the multiple receiving antennas and the transmitter of the second wireless device; andperforming channel estimation, during an operational phase subsequent to the training phase, using the second order statistics and reference tones received on a second number of resource elements, wherein the second number is less than the first number,wherein performing the channel estimation comprises: interpolating channel estimates from the reference tones received on the second number of resource elements using an interpolation filter, andwherein the estimating the second order statistics includes estimating covariance of a channel matrix using (a) a direct method that formulates the channel matrix to comprise a number of columns, wherein each column comprises an estimated channel response over all the receiving antennas of the first wireless device at a time instance, or (b) an averaging method that formulates the channel matrix as an average of autocovariance of an estimated channel response over all the receiving antennas of the first wireless device at different time instances.
  • 2. The method of claim 1, further comprising: determining, using the second order statistics, the interpolation filter.
  • 3. The method of claim 2, further including: estimating an error covariance matrix representative of an error in the estimated channel response over all of the receiving antennas computed by interpolating the channel estimates; andselectively revising, for future use, the second number of resource elements of reference tones using a measure of the error covariance matrix.
  • 4. The method of claim 3, wherein the measure of the error covariance matrix includes a comparison between a sum of squares of diagonal entries of the error covariance matrix and a threshold.
  • 5. A wireless communication method, comprising: receiving, by a first wireless device during a training phase, reference tones using a first number of resource elements from a transmitter of a second wireless device, wherein the first wireless device comprises multiple receiving antennas;estimating, by the first wireless device, from the receiving the reference tones, a second order statistics of wireless channels between the multiple receiving antennas and the transmitter of the second wireless device;performing channel estimation, during an operational phase subsequent to the training phase, using the second order statistics and reference tones received on a second number of resource elements, wherein the second number is less than the first number,wherein performing the channel estimation comprises: interpolating channel estimates from the reference tones received on the second number of resource elements using an interpolation filter; anddetermining, using the second order statistics, the interpolation filter, wherein the determining the interpolation filter includes using a decimated version of a covariance matrix calculated during the training phase.
  • 6. The method of claim 5, wherein the determining the interpolation filter is performed by: determining an eigenvector factorization of the covariance matrix, wherein the eigenvector factorization is represented as: RHH≈V·D·V* wherein RHH is the covariance matrix, V is an eigenvectors matrix and D is a diagonal matrix of eigenvalues, and * is a conjugate transpose operation; andestimating the interpolation filter as function of the eigenvectors matrix and a decimated version of the eigenvectors matrix.
  • 7. The method of claim 6, wherein the function is represented as C=V·(VM)−1, where C represents the interpolation filter, VM is the decimated version of V.
  • 8. A method of wireless communication, comprising: transmitting, to multiple receive antennas of a first wireless device from a transmit antenna of a second wireless device, during a training phase, reference tones using a first number of resource elements of a wireless channel between the transmit antenna and the multiple receive antennas;receiving, at an end of the training phase, an instruction to transmit reference tones using a second number of resource elements;transmitting, during an operational phase after the training phase, reference tones to the multiple receive antennas of the first wireless device using the second number of resource elements, wherein the second number is different from the first number, wherein the second number is based on an estimated second order statistics of the wireless channel, wherein the first wireless device is configured to perform channel estimation during the operational phase, and wherein performing the channel estimation comprises interpolating channel estimates from the reference tones received on the second number of resource elements using an interpolation filter;receiving, during the operational phase, another instruction to transmit reference tones using a third number of resource elements, wherein the third number is different from the first number and the second number; andtransmitting, after receiving the another instruction, reference tones to the first wireless device using the third number of resource elements.
  • 9. The method of claim 8, wherein the resource elements comprise subcarriers.
  • 10. The method of claim 8, wherein the resource elements comprise time slots.
  • 11. The method of claim 8, wherein the second number is less than the first number.
  • 12. The method of claim 8, wherein the second number is a pre-defined number.
  • 13. The method of claim 8, wherein the training phase may be based on a pre-defined number of repetitions of reference signal transmissions or on passage of a pre-defined time.
  • 14. The method of claim 8, wherein the transition from the training phase to the operational phase is transmitted from the first wireless device to the second wireless device.
  • 15. A first wireless communication device, comprising: a transceiver configured to receive, during a training phase, reference tones using a first number of resource elements from a transmitter of a second wireless device, wherein the first wireless device comprises multiple receiving antennas; anda processor configured to:estimate, from the receiving the reference tones, a second order statistics of wireless channels between the multiple receiving antennas and the transmitter of the second wireless device; andperform channel estimation, during an operational phase subsequent to the training phase, using the second order statistics and reference tones received on a second number of resource elements, wherein the second number is less than the first number,wherein the channel estimation is performed by: interpolating channel estimates from the reference tones received on the second number of resource elements using an interpolation filter, andwherein the estimating the second order statistics includes estimating covariance of a channel matrix using (a) a direct method that formulates the channel matrix to comprise a number of columns, wherein each column comprises an estimated channel response over all the receiving antennas of the first wireless device at a time instance, or (b) an averaging method that formulates the channel matrix as an average of autocovariance of an estimated channel response over all the receiving antennas of the first wireless device at different time instances.
  • 16. The first wireless communication device of claim 15, wherein the processor is further configured to: determine, using the second order statistics, the interpolation filter.
  • 17. The first wireless communication device of claim 15, wherein the processor is further configured to: estimate an error covariance matrix representative of an error in the estimated channel response over all of the receiving antennas computed by interpolating the channel estimates; andselectively revise, for future use, the second number of resource elements of reference tones using a measure of the error covariance matrix.
  • 18. The first wireless communication device of claim 17, wherein the measure of the error covariance matrix includes a comparison between a sum of squares of diagonal entries of the error covariance matrix and a threshold.
Parent Case Info

This patent document is a 371 of International Application No. PCT/US2019/037095, filed Jun. 13, 2019, which claims priority to and benefits of U.S. Provisional Application No. 62/684,594 entitled “SECOND ORDER STATISTICS FOR EFFICIENT CHANNEL ESTIMATION,” filed on Jun. 13, 2018 and U.S. Provisional Application No. 62/726,822 entitled “RECIPROCITY CALIBRATION IN WIRELESS COMMUNICATIONS,” filed on Sep. 4, 2018. The entire contents of the aforementioned patent applications are incorporated by reference as part of the disclosure of this patent document.

PCT Information
Filing Document Filing Date Country Kind
PCT/US2019/037095 6/13/2019 WO 00
Publishing Document Publishing Date Country Kind
WO2019/241589 12/19/2019 WO A
US Referenced Citations (251)
Number Name Date Kind
4754493 Coates Jun 1988 A
5083135 Nagy et al. Jan 1992 A
5182642 Gersdorff et al. Jan 1993 A
5623511 Bar-David et al. Apr 1997 A
5831977 Dent Nov 1998 A
5872542 Simons et al. Feb 1999 A
5956624 Hunsinger et al. Sep 1999 A
6212246 Hendrickson Apr 2001 B1
6289063 Duxbury Sep 2001 B1
6356555 Rakib et al. Mar 2002 B1
6388621 Lynch May 2002 B1
6426983 Rakib et al. Jul 2002 B1
6608864 Strait Aug 2003 B1
6631168 Izumi Oct 2003 B2
6704366 Combes et al. Mar 2004 B1
6956814 Campanella Oct 2005 B1
6973144 Zhu Dec 2005 B1
7327812 Auer Feb 2008 B2
7392018 Ebert et al. Jun 2008 B1
7773685 Tirkkonen et al. Aug 2010 B2
7864877 Hottinen Jan 2011 B2
8229017 Lee et al. Jul 2012 B1
8259845 Dent Sep 2012 B2
8401131 Fety et al. Mar 2013 B2
8547988 Hadani et al. Oct 2013 B2
8619892 Vetter et al. Dec 2013 B2
8879378 Rakib et al. Nov 2014 B2
8892048 Turner Nov 2014 B1
8976851 Hadani et al. Mar 2015 B2
9031141 Hadani et al. May 2015 B2
9071285 Hadani et al. Jun 2015 B2
9071286 Hadani et al. Jun 2015 B2
9083483 Rakib et al. Jul 2015 B1
9083595 Rakib et al. Jul 2015 B2
9130638 Hadani et al. Sep 2015 B2
9282528 Hashimoto Mar 2016 B2
9294315 Hadani et al. Mar 2016 B2
9444514 Hadani et al. Sep 2016 B2
9548840 Hadani et al. Jan 2017 B2
9553984 Krause et al. Jan 2017 B2
9590779 Hadani et al. Mar 2017 B2
9634719 Rakib et al. Apr 2017 B2
9660851 Hadani et al. May 2017 B2
9667307 Hadani et al. May 2017 B1
9668148 Hadani et al. May 2017 B2
9712354 Hadani et al. Jul 2017 B2
9722741 Rakib et al. Aug 2017 B1
9729281 Hadani et al. Aug 2017 B2
10667148 Hadani et al. May 2020 B1
10681568 Hadani et al. Jun 2020 B1
10693581 Rakib et al. Jun 2020 B2
10693692 Hadani et al. Jun 2020 B2
10716095 Rakib et al. Jul 2020 B2
10749651 Hebron et al. Aug 2020 B2
10826728 Tsatsanis et al. Nov 2020 B2
10855425 Kons et al. Dec 2020 B2
10873418 Kons et al. Dec 2020 B2
20010031022 Petrus et al. Oct 2001 A1
20010033614 Hudson Oct 2001 A1
20010046205 Easton et al. Nov 2001 A1
20020001308 Heuer Jan 2002 A1
20020034191 Shattil Mar 2002 A1
20020181388 Jain et al. Dec 2002 A1
20020181390 Mody et al. Dec 2002 A1
20020181607 Izumi Dec 2002 A1
20030073464 Giannakis et al. Apr 2003 A1
20030185295 Yousef Oct 2003 A1
20030235147 Walton et al. Dec 2003 A1
20040044715 Aldroubi et al. Mar 2004 A1
20040174812 Murakami et al. Sep 2004 A1
20040189581 Sako et al. Sep 2004 A1
20040218523 Varshney et al. Nov 2004 A1
20050157778 Trachewsky et al. Jul 2005 A1
20050157820 Wongwirawat et al. Jul 2005 A1
20050180517 Abe Aug 2005 A1
20050207334 Hadad Sep 2005 A1
20050251844 Martone et al. Nov 2005 A1
20060008021 Bonnet Jan 2006 A1
20060039270 Strohmer et al. Feb 2006 A1
20060189280 Goldberg Aug 2006 A1
20070014272 Palanki et al. Jan 2007 A1
20070038691 Candes et al. Feb 2007 A1
20070078661 Sriram et al. Apr 2007 A1
20070104283 Han et al. May 2007 A1
20070110131 Guess et al. May 2007 A1
20070177695 Salem Aug 2007 A1
20070211952 Faber et al. Sep 2007 A1
20070237181 Cho et al. Oct 2007 A1
20070253465 Muharemovic et al. Nov 2007 A1
20070253504 Hasegawa Nov 2007 A1
20080043857 Dias et al. Feb 2008 A1
20080117999 Kadous et al. May 2008 A1
20080186843 Ma et al. Aug 2008 A1
20080187062 Pan et al. Aug 2008 A1
20080232504 Ma et al. Sep 2008 A1
20080310383 Kowalski Dec 2008 A1
20090080403 Hamdi Mar 2009 A1
20090092259 Jot et al. Apr 2009 A1
20090103593 Bergamo Apr 2009 A1
20090122854 Zhu et al. May 2009 A1
20090161804 Chrabieh et al. Jun 2009 A1
20090204627 Hadani Aug 2009 A1
20090222226 Baraniuk et al. Sep 2009 A1
20090303961 Popovic et al. Dec 2009 A1
20100001901 Baraniuk et al. Jan 2010 A1
20100002816 Mody Jan 2010 A1
20100008432 Kim et al. Jan 2010 A1
20100027608 Priotti Feb 2010 A1
20100111138 Hosur et al. May 2010 A1
20100142476 Jiang et al. Jun 2010 A1
20100017884 Nassiri-Toussi et al. Jul 2010 A1
20100187914 Rada et al. Jul 2010 A1
20100238787 Guey Sep 2010 A1
20100277308 Potkonjak Nov 2010 A1
20100303136 Ashikhmin et al. Dec 2010 A1
20100322349 Lee et al. Dec 2010 A1
20110007789 Garmany Jan 2011 A1
20110110532 Svendsen May 2011 A1
20110116489 Grandhi May 2011 A1
20110116516 Hwang et al. May 2011 A1
20110126071 Han et al. May 2011 A1
20110131463 Gunnam Jun 2011 A1
20110195670 Dakshinamurthy et al. Aug 2011 A1
20110216808 Tong et al. Sep 2011 A1
20110286502 Adachi et al. Nov 2011 A1
20110287778 Levin et al. Nov 2011 A1
20110292971 Hadani et al. Dec 2011 A1
20110293030 Rakib et al. Dec 2011 A1
20110299379 Sesia et al. Dec 2011 A1
20110305267 Riu et al. Dec 2011 A1
20120021769 Lindoff et al. Jan 2012 A1
20120051457 Ma et al. Mar 2012 A1
20120140716 Baldemair et al. Jun 2012 A1
20120170684 Yim et al. Jul 2012 A1
20120201322 Rakib et al. Aug 2012 A1
20120213098 Sun Aug 2012 A1
20120235795 Liao et al. Sep 2012 A1
20120269201 Atungsiri et al. Oct 2012 A1
20120272117 Stadelmeier et al. Oct 2012 A1
20120320994 Loghin et al. Dec 2012 A1
20130003802 Doron Jan 2013 A1
20130021977 Yang et al. Jan 2013 A1
20130058390 Haas et al. Mar 2013 A1
20130077579 Cho et al. Mar 2013 A1
20130083661 Gupta et al. Apr 2013 A1
20130121497 Smaragdis et al. May 2013 A1
20130230010 Kim et al. Sep 2013 A1
20130260787 Hashimoto Oct 2013 A1
20130279627 Wu et al. Oct 2013 A1
20130315133 Wang et al. Nov 2013 A1
20140143639 Loghin et al. May 2014 A1
20140161154 Hadani et al. Jun 2014 A1
20140169385 Hadani et al. Jun 2014 A1
20140169406 Hadani et al. Jun 2014 A1
20140169433 Hadani et al. Jun 2014 A1
20140169436 Hadani et al. Jun 2014 A1
20140169437 Hadani et al. Jun 2014 A1
20140169441 Hadani et al. Jun 2014 A1
20140247803 Arambepola et al. Sep 2014 A1
20140348252 Siohan et al. Nov 2014 A1
20140364128 Lee et al. Dec 2014 A1
20150117395 Hadani et al. Apr 2015 A1
20150304130 Logothets et al. Oct 2015 A1
20150326273 Rakib et al. Nov 2015 A1
20150327085 Hadani et al. Nov 2015 A1
20150382231 Jabbar et al. Dec 2015 A1
20160043835 Hadani et al. Feb 2016 A1
20160135132 Donepudi et al. May 2016 A1
20160182269 Hadani et al. Jun 2016 A1
20160191217 Hadani et al. Jun 2016 A1
20160191280 Hadani et al. Jun 2016 A1
20160254889 Shattil Sep 2016 A1
20160277225 Frenne et al. Sep 2016 A1
20160309345 Tehrani et al. Oct 2016 A1
20160380743 Rakib Dec 2016 A1
20160381576 Hadani et al. Dec 2016 A1
20170012749 Rakib et al. Jan 2017 A1
20170012810 Rakib et al. Jan 2017 A1
20170019297 Rakib Jan 2017 A1
20170033899 Rakib et al. Feb 2017 A1
20170040711 Rakib et al. Feb 2017 A1
20170078054 Hadani et al. Mar 2017 A1
20170099122 Hadani et al. Apr 2017 A1
20170099607 Hadani et al. Apr 2017 A1
20170149594 Rakib et al. May 2017 A1
20170149595 Rakib et al. May 2017 A1
20170201354 Hadani et al. Jul 2017 A1
20170207817 Hadani et al. Jul 2017 A1
20170222700 Hadani et al. Aug 2017 A1
20170230215 Rakib et al. Aug 2017 A1
20170244524 Hadani et al. Aug 2017 A1
20170272177 Seyama Sep 2017 A1
20170288913 Rakib et al. Oct 2017 A1
20170289961 Rakib et al. Oct 2017 A1
20170302353 Rahman Oct 2017 A1
20180076871 Rahman Mar 2018 A1
20180109284 Hadani et al. Apr 2018 A1
20180167165 Kons et al. Jun 2018 A1
20180205481 Rakib et al. Jul 2018 A1
20180227159 Rakib et al. Aug 2018 A1
20180242170 Hadani et al. Aug 2018 A1
20180262306 Hadani et al. Sep 2018 A1
20180288809 Delfeld et al. Oct 2018 A1
20180302802 Fanfelle Oct 2018 A1
20190013914 Xu Jan 2019 A1
20190036577 Delfeld et al. Jan 2019 A1
20190036741 Hadani et al. Jan 2019 A1
20190044682 Hebron et al. Feb 2019 A1
20190075551 Hadani et al. Mar 2019 A1
20190081836 Hadani et al. Mar 2019 A1
20190159177 Rakib et al. May 2019 A1
20190173617 Kons et al. Jun 2019 A1
20190173630 Kons et al. Jun 2019 A1
20190215109 Hadani et al. Jul 2019 A1
20190238189 Delfeld et al. Aug 2019 A1
20190327054 Kons et al. Oct 2019 A1
20190342126 Hadani et al. Nov 2019 A1
20190342136 Hadani et al. Nov 2019 A1
20190379422 Hadani et al. Dec 2019 A1
20200045562 Hadani et al. Feb 2020 A1
20200119868 Rakib et al. Apr 2020 A1
20200137774 Molisch et al. Apr 2020 A1
20200145273 Rakib et al. May 2020 A1
20200153107 Rakib May 2020 A1
20200186397 Tsatsanis et al. Jun 2020 A1
20200204309 Namboodiri Jun 2020 A1
20200204410 Hadani Jun 2020 A1
20200228170 Delfeld et al. Jul 2020 A1
20200259604 Hadani et al. Aug 2020 A1
20200259692 Hadani et al. Aug 2020 A1
20200259697 Delfeld Aug 2020 A1
20200280138 Fanfelle et al. Sep 2020 A1
20200287672 Namboodiri et al. Sep 2020 A1
20200288333 Rakib et al. Sep 2020 A1
20200305010 Hadani et al. Sep 2020 A1
20200313695 Namboodiri et al. Oct 2020 A1
20200313949 Hadani Oct 2020 A1
20200322185 Kons et al. Oct 2020 A1
20200322202 Hadani et al. Oct 2020 A1
20200351836 Rakib et al. Nov 2020 A1
20200367252 Hebron et al. Nov 2020 A1
20200389268 Sathyanarayan et al. Dec 2020 A1
20200403829 Namboodiri et al. Dec 2020 A1
20210028877 Rakib et al. Jan 2021 A1
20210036823 Hebron et al. Feb 2021 A1
20210058114 Molisch et al. Feb 2021 A1
20210077625 Kons et al. Mar 2021 A1
20210105155 Kons et al. Apr 2021 A1
20210126750 Kons et al. Apr 2021 A1
20210135905 Kons et al. May 2021 A1
20210194732 Kons et al. Jun 2021 A1
Foreign Referenced Citations (56)
Number Date Country
2011150315 Dec 2011 WO
2013148546 Oct 2013 WO
2014004585 Jan 2014 WO
2016014596 Jan 2016 WO
2016014598 Jan 2016 WO
2016176642 Nov 2016 WO
2016183230 Nov 2016 WO
2016183240 Nov 2016 WO
2016209848 Dec 2016 WO
2017003952 Jan 2017 WO
2017011455 Jan 2017 WO
2017011478 Jan 2017 WO
2017044501 Mar 2017 WO
2017087706 May 2017 WO
2017100666 Jun 2017 WO
2017147439 Aug 2017 WO
2017165697 Sep 2017 WO
2017173160 Oct 2017 WO
2017173389 Oct 2017 WO
2017201467 Nov 2017 WO
2018031367 Feb 2018 WO
2018031938 Feb 2018 WO
2018031952 Feb 2018 WO
2018032016 Feb 2018 WO
2018064327 Apr 2018 WO
2018064587 Apr 2018 WO
2018064605 Apr 2018 WO
2018086696 May 2018 WO
2018106731 Jun 2018 WO
2018129554 Jul 2018 WO
2018140837 Aug 2018 WO
2018191309 Oct 2018 WO
2018195548 Oct 2018 WO
2018200567 Nov 2018 WO
2018200577 Nov 2018 WO
2019014332 Jan 2019 WO
2019032142 Feb 2019 WO
2019032605 Feb 2019 WO
2019036492 Feb 2019 WO
2019051093 Mar 2019 WO
2019051427 Mar 2019 WO
2019055861 Mar 2019 WO
2019068053 Apr 2019 WO
2019060596 May 2019 WO
2019089986 May 2019 WO
2019113046 Jun 2019 WO
2019157230 Aug 2019 WO
2019173775 Sep 2019 WO
2019241436 Dec 2019 WO
2019241589 Dec 2019 WO
2020142520 Jul 2020 WO
2020206304 Oct 2020 WO
2020227619 Nov 2020 WO
2020247768 Dec 2020 WO
2021026212 Feb 2021 WO
2021062354 Apr 2021 WO
Non-Patent Literature Citations (12)
Entry
International Search Report and Written Opinion for PCT/US19/037095, dated Sep. 2, 2019, 13 pages.
Banelli, P. et al., “Modulation Formats and Waveforms for 5G Networks: Who Will Be the Heir of OFDM?,” IEEE Signal Processing Magazine, vol. 81, pp. 80-93, Nov. 2014.
El Hattachi, R. et al., “NGMN 5G Initiative White Paper,” NGMN Alliance, Feb. 17, 2015. [Online]. Available: https://www.ngmn.org/uploads/media/NGMN_5G_White_Paper_V1_0.pdf, 125 pages.
Rusek, F. et al., “Scaling Up MIMO, Opportunities and Challenges with Very Large Arrays,” IEEE Signal Processing Magazine, pp. 40-60 (2013).
Vodafone, “Cellular Internet of Things: Architectural Aspects,” RP-150869, 3GPP RAN#68, Malmo, Sweden (Jun. 9, 2015), 19 pages.
Gurevich, S. et al. “Group Representation Design of Digital Signals and Sequences,” S.W. Golomb et al. (eds.), SETA 2008, LNCS 5203, pp. 153-166, Springer-Verlag Berlin Heidelberg (2008).
“AT&T Annual Report 2014,” Opening Our Network [Online]. Retrieved from the Internet: Sep. 22, 2016. <URL: http://www.att.com/Investor/ATT_Annual/2014/att_introduces_new_concepts_for_telecom_network.html>, 5 pages.
Catt, “UL ACK/NACK transmission methods for LTE-A,” 3GPP TSG RAN WG1 Meeting #60bis, R1-102453, Beijing, China, Apr. 12-16, 2010, 8 pages.
Toskala, A. et al., “Physical Layer,” Chapter 5 In: “LTE for UMTS: OFDMA and SC-FDMA Based Radio Access,” Holma, H. et al. (eds.), John Wiley & Sons, Ltd., United Kingdom, 2009, pp. 83-135.
Mecklenbrauker, W., “A Tutorial on Non-Parametric Bilinear Time-Frequency Signal Representations,” In: Time and Frequency Representation of Signals and Systems, Longo, G. et al. (eds.), Springer-Verlag Wien, vol. 309, pp. 11-68 (1989).
Nehorai, A. et al., “MURI: Adaptive waveform design for full spectral dominance (2005-2010),” AFOSR FA9550-05-1-0443, Final Report, [online], Mar. 11, 2011 Retrieved on May 11, 2013, Retrieved from the Internet <URL: http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADA565420>, 103 pages.
Extended European Search Report dated Jun. 25, 2021 for EP Application No. 19819211.4, 10 pages.
Related Publications (1)
Number Date Country
20210135905 A1 May 2021 US
Provisional Applications (2)
Number Date Country
62726822 Sep 2018 US
62684594 Jun 2018 US