# Patent

## Title
Time-reversal wireless systems having asymmetric architecture

## Abstract

In an asymmetric time-reversal wireless system, a base station includes an input circuit configured to, during a hand-shake period, receive a channel response signal derived from a probe signal sent from a first terminal device to the apparatus through multiple wireless propagation paths, and during an uplink transmission period, receive combined signals that include a signal from the first terminal device and a signal from a second terminal device. The base station includes a data processor configured to calculate a signature waveform for the first terminal device based on the channel response, and determine the signal sent from the first terminal device during the uplink transmission period based on the combined signals and the signature waveform for the first terminal device.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application is the U.S. National Stage Application of International Application No. PCT/US14/51148, filed on Aug. 14, 2014, which claims priority to U.S. patent application Ser. No. 13/969,271, filed on Aug. 16, 2013, and U.S. patent application Ser. No. 13/969,320, filed on Aug. 16, 2013, now U.S. Pat. No. 9,559,874. The entire disclosures of the above applications are incorporated herein by reference.
TECHNICAL FIELD
This disclosure generally relates to time-reversal wireless systems having asymmetric architecture.
BACKGROUND
In some time-reversal communication systems, when a transceiver A intends to transmit information to a transceiver B, transceiver B first sends a waveform that propagates through a scattering and multi-path environment, and the signals are received by transceiver A. Transceiver A transmits time-reversed signals back through the same channel to transceiver B. Based on channel reciprocity, a time-reversal communication system leverages the multi-path channel as a matched filter, i.e., treats the environment as a facilitating matched filter computing machine, and focuses the wave at the receiver in both space and time domains.
SUMMARY
In general, in one aspect, a method for time-reversal wireless communication comprises at a transceiver, receiving a channel response signal derived from a probe signal sent from a first device to the transceiver through multiple wireless propagation paths; calculating a signature waveform for the first device based on the channel response; receiving combined signals that include a signal from the first device and a signal from a second device; and determining the signal from the first device based on the combined signals and the signature waveform for the first device.
Implementations of the method may include one or more of the following features. Determining the signal from the first device comprises performing a convolution operation, or an operation equivalent to the convolution operation, between the combined signals and the signature waveform for the first device.
Determining the signal from the first device can include filtering the combined signals using the signature waveform as a filtering parameter.
The signature waveform for the first device can be a time-reversed waveform of the channel response signal derived from the probe signal sent from the first device.
In some examples, the channel response signal can be represented as h(t) and the signature waveform for the first device can be represented as h*(T−t), t representing time, and T representing a delay spread of the multiple propagation paths. In some examples, the channel response signal can be represented as h[k] and the signature waveform for the first device can be represented as h*[L−k], k being an integer, and L representing a delay spread of the multiple propagation paths.
The method can include transmitting, from the transceiver to the first receiver and a second receiver, a combined downlink signal derived from a combination of a first downlink waveform and a second downlink waveform, the first downlink waveform including data intended for the first device, the second downlink waveform including data intended for the second device.
The method can include generating the first downlink waveform based on the data intended for the first device and the signature waveform for the first device.
Generating the first downlink waveform can include performing a convolution of the data intended for the first device and the signature waveform for the first device.
The combined downlink signal can be configured to enable each of the first device and the second device to receive multipath signals that can be used determine a data signal intended for the device, the first and second devices receiving the downlink signal through different propagation paths and determining different data signals.
The transceiver can be part of a base station.
Receiving combined signals can include receiving combined signals that include signals from three or more devices communicating with the transceiver, each of the devices being associated with a corresponding signature waveform.
In some examples, the first device does not calculate the signature waveform. The first device does not perform any signal processing that uses the signature waveform.
In general, in another aspect, an apparatus for time-reversal wireless communication comprises an input circuit configured to, during a hand-shake period, receive a channel response signal derived from a probe signal sent from a first device to the apparatus through multiple wireless propagation paths, and during an uplink transmission period, receive combined signals that include a signal from the first device and a signal from a second device; and a data processor configured to calculate a signature waveform for the first device based on the channel response, and determine the signal sent from the first device during the uplink transmission period based on the combined signals and the signature waveform for the first device.
Implementations of the apparatus may include one or more of the following features. The apparatus can include a base station that includes the input circuit and the data processor.
The data processor can determine the signal sent from the first device by performing a convolution operation, or an operation equivalent to the convolution operation, between the combined signals and the signature waveform for the first device.
The data processor can determine the signal sent from the first device by filtering the combined signals using the signature waveform as a filtering parameter.
The data processor can determine the signature waveform for the first device by generating a time-reversed waveform of the channel response signal derived from the probe signal sent from the first device.
In some examples, the channel response signal can be represented as h(t) and the signature waveform for the first device can be represented as h*(T−t), t representing time, and T representing a delay spread of the multiple propagation paths. In some examples, the channel response signal can be represented as h[k] and the signature waveform for the first device can be represented as h*[L−k], k being an integer, and L representing a delay spread of the multiple propagation paths.
The data processor can be configured to determine a combined downlink signal derived from a combination of a first downlink waveform and a second downlink waveform, the first downlink waveform including data intended for the first device, the second downlink waveform including data intended for the second device.
The data processor can be configured to determine the first downlink waveform based on the data intended for the first device and the signature waveform for the first device.
The data processor can determine the first downlink waveform by performing a convolution, or an operation equivalent to the convolution, of the data intended for the first device and the signature waveform for the first device.
The combined downlink signal can be configured to enable each of the first device and the second device to receive multipath signals that can be used determine a data signal intended for the device, the first and second devices receiving the downlink signal through different propagation paths and determining different data signals.
The input circuit can be configured to receive combined signals that include signals from three or more devices that communicate with the apparatus, and the data processor is configured to calculate a signature waveform for each of the devices.
In some examples, the first device does not calculate the signature waveform. The first device does not perform any signal processing that uses the signature waveform.
In general, in another aspect, a method of communication between a transceiver and multiple devices using wireless multi-path signals is provided. The method comprises performing a hand-shake process between the transceiver and the devices, in which for each device, the transceiver receives a channel response signal derived from a probe signal sent from the device to the transceiver through multiple propagation paths; calculating a signature waveform for each of the devices based on the channel response derived from the probe signal sent from the device; receiving a combined uplink signal that includes uplink signals from two or more devices; and determining the uplink signal from each one of the devices based on the combined uplink signal and the signature waveform for the device.
Implementations of the method may include one or more of the following features. The transceiver can be part of a base station.
In some examples, the devices do not calculate any signature waveform. The devices do not perform any signal processing that uses the signature waveform.
In general, in another aspect, a method of asymmetric time-reversal communication between a plurality of devices using wireless multi-path signals is provided. The method comprises performing a first hand-shake process between the first and second devices, in which the first device receives a channel response signal derived from a probe signal sent from the second device to the first device through multiple wireless propagation paths; performing a second hand-shake process between the first and third devices, in which the first device receives a channel response signal derived from a probe signal sent from the third device to the first device through multiple wireless propagation paths; at the first device, calculating a signature waveform for the second device based on the channel response derived from the probe signal sent from the second device; at the first device, calculating a signature waveform for the third device based on the channel response derived from the probe signal sent from the third device; at the first device, receiving a combined uplink signal that includes uplink signals from the second and third devices; at the first device, determining the uplink signal from the second device based on the combined uplink signal and the signature waveform for the second device; and at the first device, determining the uplink signal from the third device based on the combined uplink signal and the signature waveform for the third device.
Implementations of the method may include one or more of the following features. The second and third devices do not calculate any signature waveform.
The method can include transmitting, from the first device to the second and third devices, a combined downlink signal derived from a combination of a second downlink waveform and a third downlink waveform, the second downlink waveform including data intended for the second device, the third downlink waveform including data intended for the third device.
The second downlink waveform can be embedded with the signature waveform corresponding to the second device.
The method can include at the second device, receiving a signal intended for the second device without performing any signal processing that uses the signature waveform corresponding to the second device.
The method can include at the first device, generating the second downlink waveform based on the data intended for the second device and the signature waveform for the second device.
Generating the second downlink waveform can include performing a convolution of the data intended for the second device and the signature waveform for the second device.
The combined downlink signal can be configured to enable each of the second and third devices to receive multipath signals that can be used determine a data signal intended for the device, the second and third devices receiving the combined downlink signal through different propagation paths and determining different data signals.
Determining the uplink signal from the second device can include performing a convolution operation, or an operation equivalent to the convolution operation, between the combined uplink signal and the signature waveform for the second device.
Determining the uplink signal from the second device can include filtering the combined uplink signal using the signature waveform corresponding to the second device as a filtering parameter.
The signature waveform for the second device can be a time-reversed waveform of the channel response signal derived from the probe signal sent from the second device.
In some examples, the channel response signal derived from the probe signal sent from the second device can be represented as h(t) and the signature waveform for the second device can be represented as h*(T−t), t representing time, and T representing a delay spread of the multiple propagation paths. In some examples, the channel response signal derived from the probe signal sent from the second device can be represented as h[k] and the signature waveform for the second device can be represented as h*[L−k], k being an integer, and L representing a delay spread of the multiple propagation paths.
The first device can include a base station.
Receiving a combined uplink signal can include receiving a combined uplink signal that includes signals from the second and third devices, and one or more additional devices communicating wirelessly with the first device, each of the devices being associated with a corresponding signature waveform.
The details of one or more implementations of an asymmetric time-reversal wireless system are set forth in the accompanying drawings and the description below. Other features, aspects, and advantages will become apparent from the description, the drawings, and the claims.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a diagram showing an exemplary environment for operating a time-reversal system.
FIG. 2A is a graph of an exemplary recorded channel response waveform.
FIG. 2B is a graph of an exemplary time-reversed waveform generated by reversing the waveform of FIG. 2A with respect to time.
FIG. 3 is a diagram showing an exemplary environment for operating a time-reversal system having multiple receivers.
FIG. 4 is a diagram of an exemplary multi-user time reversal communication system.
FIG. 5 is a diagram showing a hand-shaking process.
FIG. 6 is a diagram showing a down-link transmission process.
FIG. 7 is a diagram showing an up-link transmission process.
FIG. 8 is a graph showing the bit error rate (BER) performance of the time reversal system for downlink transmission in an outdoor environment.
FIG. 9 is a graph showing the bit error rate performance of the time reversal system for uplink transmission in an outdoor environment.
FIG. 10 is a graph showing the bit error rate performance of the time reversal system for downlink transmission in an indoor environment.
FIG. 11 is a graph showing the bit error rate performance of the time reversal system for uplink transmission in an indoor environment.
FIGS. 12-14 are flow diagrams of exemplary processes for asymmetric time-reversal wireless communication.
FIG. 15 is a diagram of an exemplary multi-user time reversal communication system.
FIGS. 16A and 16D are diagrams showing an exemplary 2D interference cancelation and decision module.
FIG. 17 is a diagram showing a multi-stage 2D interference cancelation and decision module.
FIGS. 18A and 18B are graphs showing examples of channel impulse responses.
FIGS. 19A and 19B are graphs showing examples of convolutions of channel impulse responses and time-reversal signature waveforms.
FIGS. 20 to 23 are graphs showing the bit error rate performances of the time reversal system using 2D parallel interference cancellation under various operation conditions.
FIGS. 24 to 26 are flow diagrams of exemplary processes for asymmetric time-reversal wireless communication with interference cancellation.
FIG. 27 is a diagram of an exemplary multi-user time reversal communication system that uses multiple communication protocols.
FIG. 28A is a diagram of an exemplary transmitter based on Wi-Fi protocol.
FIG. 28B is a diagram of an exemplary receiver based on Wi-Fi protocol.
FIG. 29A is a diagram of an exemplary transmitter at a base station for downlink transmission based on 4G LTE protocol.
FIG. 29B is a diagram of an exemplary receiver at a terminal device for downlink transmission based on 4G LTE protocol.
FIG. 30A is a diagram of an exemplary transmitter at the terminal device for uplink transmission based on 4G LTE protocol.
FIG. 30B is a diagram of an exemplary receiver at the base station for uplink transmission based on 4G LTE protocol.
FIG. 31A is a diagram of an exemplary downlink architecture of an asymmetrical time-reversal communication system.
FIG. 31B is a diagram of an exemplary uplink architecture of the asymmetrical time-reversal communication system.
FIG. 32 is a diagram of an exemplary symmetrical time-reversal communication system.
FIG. 33 is a diagram of exemplary system architectures of a base station and terminal devices that use both symmetrical and asymmetrical time-reversal communication.
FIG. 34 is a diagram of an implementation of a base station and terminal devices that use both symmetrical and asymmetrical time-reversal communication.
FIG. 35 is a diagram of an exemplary a mobile device having communication modules using a first communication protocol having a substantially symmetric architecture and a second communication having an asymmetric architecture.
FIG. 36 is a diagram of an exemplary group of robots that communicate using a first communication protocol that has a substantially symmetrical architecture and a second communication protocol that has an asymmetrical architecture.
FIG. 37 is a diagram of an exemplary vehicle having a communication unit that communicates using a first communication protocol that has a substantially symmetrical architecture and a second communication protocol that has an asymmetrical architecture.



Like reference symbols in the various drawings indicate like elements.
DETAILED DESCRIPTION
This disclosure describes a novel asymmetric multi-user time reversal wireless communication system architecture for both downlink and uplink communication between a base station (BS) and multiple terminal devices (TDs). In some implementations, the complexities of the base station and the terminal devices are asymmetric, such that the base station performs most of the signal processing as both a transmitter (for the downlink) and receiver (for the uplink), allowing the use of low complexity terminal devices.
For example, an intelligent house may include one or more computers that communicate with several sensors (e.g., temperature, humidity, light, and motion sensors), meters (e.g., electricity and water meters), appliances (e.g., refrigerator, oven, washing machine), electronic devices (e.g., television, digital video recorder, audio/video system, telephone, digital photo album, intelligent lamp, security system), climate control systems (e.g., fans, thermostats for air conditioning and heating, motorized shades), power generators (e.g., backyard wind turbine, solar panel, geothermal energy system).
To reduce the cost of communication, it may be preferable that the devices communicating with the computer be simple without the need to perform complex signal processing. The computer may function as a base station or be coupled to a base station that performs most of the signal processing when transmitting signals to the devices and receiving signals sent from the devices.
For example, an intelligent factory may include one or more computers that communicate wirelessly with robots working in assembly lines, vehicles that move supplies and assembled products, climate control systems, security systems, inventory control systems, and power systems. For example, a laboratory may include one or more computers that communicate wirelessly with instruments that monitor parameters when conducting experiments.
In the examples above, the computer can communicate with the devices using time-reversal division multiple access technology that uses the environment to provide spatial filtering, allowing a large number of devices to communicate with the computers simultaneously. Compared to using previous wireless communication technologies, such as Wi-Fi or Bluetooth, time-reversal division multiple access has the advantage that the additional cost for enabling each device to communicate with the computer is small because the device itself does not need to perform complicated signal processing. Most of the signal processing is performed at the computer. The power consumption by each device for enabling wireless communication may also be much smaller compared to previous wireless technologies.
Because the signals are transmitted through multiple propagation paths having various propagation lengths, there may be a large delay spread, and it may be possible for the computer to receive the tail end of the signal of an earlier symbol and the front end of a later symbol at the same time, resulting in inter-symbol interference. During an uplink cycle, multiple devices may transmit signals to the computer at the same time, so the computer may receive a combined signal having components from several devices. The time-reversal system may use signature waveforms associated with each device to identify the signal components for each device. After using the signature waveforms as a filtering mechanism to identify individual signal components, there may still be interference among signals from different devices. To increase the accuracy of communication, the computer may perform signal processing to cancel inter-symbol interference and/or inter-user interference.
Compared to previous multi-carrier techniques used to reduce inter-symbol interference, such as orthogonal frequency-division multiplexing (OFDM) or orthogonal frequency-division multiple access (OFDMA), the time-reversal divisional multiple access technology leverages the unique temporal and spatial focusing effects of the time reversal phenomenon to provide cost-effective single-carrier broadband multi-user communication.
The time-reversal division multiple access scheme uses the multi-path channel profile associated with each user's location as a location-specific signature for the user. Each path of the multi-path channel can be treated as a virtual antenna in the time-reversal division multiple access system, which collectively results in very high-resolution spatial focusing with potentially “pin-point” accuracy. Meanwhile, the temporal focusing effect can effectively suppress inter-symbol interference which significantly simplifies the terminal user's complexity and gives rise to higher-speed data transmission.
The computer may function as a base station or be coupled to a base station that performs most of the signal processing when transmitting signals to the devices and receiving signals sent from the devices. A 2D parallel interference cancellation scheme may allow multiple devices to communicate with the computer at the same time with reduced interference. For a specified bit error rate and data transmission rate, the 2D parallel interference cancellation scheme can allow more devices to communicate with the computer simultaneously. For a specified bit error rate and number of devices, the 2D parallel interference cancellation scheme can allow the devices to communicate with the computer at higher data transmission rates.
Overview of Time-Reversal System
The following provides an overview of a time-reversal multiuser downlink system. Referring to FIG. 1, a time-reversal system can be used in an environment having structures or objects that may cause one or more reflections of wireless signals. For example, a venue 102 may have a first room 104 and a second room 106. When a first device 108 in the first room 104 transmits a signal to a second device 110 in the second room 106, the signal can propagate in several directions and reach the second device 110 by traveling through several propagation paths, e.g., 112, 114, and 116. The signal traveling through multiple propagation paths may be referred to as a multipath signal. As the signal travels through the propagation paths, the signal may become distorted. The multipath signal received by the second device 110 can be quite different from the signal transmitted by the first device 108.
For example, referring to FIG. 2A, when the first device 108 sends a pulse signal, the signal received by the second device 110 may have a waveform 120. The waveform 120 may be referred to as the channel impulse response signal.
Referring to FIG. 2B, a time-reversed waveform 130 can be generated by reversing the waveform 120 with respect to time. If the second device 110 sends a signal having the waveform 130, the signal will propagate in various directions, including through propagation paths 112, 114, and 116 (in reverse direction relative to the propagation direction of the impulse signal), and reach the first device 108. The multipath signal received at the first device 108 may form an impulse signal that is similar to the impulse signal previously sent from the first device 108 to the second device 110.
The waveforms 120 and 130 shown in FIGS. 2A and 2B are merely examples. The waveforms in time-reversal systems can vary depending on, e.g., the environment and the information or data being transmitted. In addition, the initial signal sent from the first device 108 can be any sort of signal, of which an impulse is just one example. The initial signal can be any waveform.
When the second device 110 intends to transmit a data stream to the first device 108, the second device 110 may use a normalized time-reversed conjugate of the signal received from device 108 as a basic transmission waveform. The second device 110 may encode the data stream on the basic waveform and transmit the signal through the wireless channel. The signal received at the device 108 may be described as the convolution of the transmitted signal and the channel impulse response, plus additive white Gaussian noise. Because the transmitted signal has been designed based on a time reversed version of the channel impulse response, the first device 108 may only need to perform a simple adjustment to the received signal and down-sample it to recover the data stream transmitted by the second device 110.
In some examples a transmitting device or base station or access point may send signals to two or more receiving devices at the same time. The transmitted signals travel through multiple propagation paths to each receiver. Because the receivers are positioned at different locations, the multipath signals travel through different propagation paths to reach the receivers. By carefully constructing the waveform of the signal sent from the transmitter, it is possible to allow each receiver to receive data intended for the receiver with sufficiently high quality.
Referring to FIG. 3, the first device 108 may communicate with the second device 110 and a third device 140. The second device 110 may send a probe signal that travels through propagation paths 112, 114, and 116 to the first device 108. The probe signal can be, e.g., a pulse signal, a signal that has a predetermined waveform, or a signal that includes symbols. The first device 108 may record the received waveform representing the channel response for a first multipath channel. The third device 140 may send a probe signal that travels through propagation paths 142, 144, and 146 to the first device 108. The first device 108 may record the received waveform representing the channel response for the second multipath channel.
The first device 108 may construct a downlink signal based on a first time-reversed multipath channel response, a second time-reversed multipath channel response, a first data stream intended for the second device receiver 110, and a second data stream intended for the third device receiver 140. The first device 108 may transmit the downlink signal so that a first portion of the downlink signal travels through propagation paths 112, 114, and 116 to reach the second device 110. A second portion of the downlink signal may travel through propagation paths 142, 144, and 146 to reach the third device 140. The first signal portion received at the second device 110 may form a first data stream. The second signal portion received at the third device 140 may form a second data stream.
In the example of FIG. 3, the device operating as a transmitter may use either an omnidirectional antenna or a directional antenna for broadcasting the downlink signal, as long as the downlink signal reaches each of the receivers through multiple propagation paths. In some systems, the transmitting antenna may use beam-forming techniques to launch a signal that undergoes a large number of reflections. That is, the transmitter may launch a wireless signal that preferentially experiences multiple transmission paths.
In some examples, a multiple input multiple output (MIMO) system may be used in which the device operating as a transmitter has multiple transmit antennas, and each of the devices operating as a receiver has multiple receive antennas. In some systems, the MIMO transmitting antennas may launch signals that undergo a large number of reflections. That is, the transmitter may launch a wireless signal that preferentially experiences multiple transmission paths.
System Architecture
Referring to FIG. 4, an exemplary multi-user time reversal communication system 150 includes a base station (BS) 152 and multiple terminal devices (e.g., 154a, 154b, 154c, collectively 154). Each of the terminal devices 154 is associated with a multi-path wireless channel (e.g., 156a, 156b, 156c, collectively 156) between itself and the base station 152. Each multi-path wireless channel 156 in the figure represents two or more multiple signal propagation paths between the corresponding terminal and the base station. In some implementations, all the devices (including the base station 152 and the terminals 154) can operate at the same frequency band. For the downlink, the base station 152 can send multiple messages (either independent or non-independent) simultaneously to multiple selected terminals 154. For the uplink, multiple terminals 154 can send their own messages to the base station 152 simultaneously.
Hand-Shake Process
The following describes an exemplary hand-shaking process between the base station 152 and the terminals 154.
Referring to FIG. 5, in some implementations, after each terminal device 154 registers with the base station 152, a hand-shaking process may occur periodically in a round-robin manner among the registered terminal devices 154 to maintain an up-to-date record of the channel responses associated with the terminal devices 154 at the base station 152.
During the hand-shaking process, at a given time, only one selected terminal, e.g., 154a, may send a probe signal 160, and may perform radio frequency modulation 162 to produce a radio frequency signal that propagates through the associated multi-path channel 156a to the base station 152. The probe signal 160 can be, e.g., a pulse signal, a signal that has a predetermined waveform, or a signal that includes symbols to enable the base station to perform synchronization and/or other functions. In the examples below, the probe signal is a pulse signal. However, other types of probe signals can be used in the system 150. The base station 152 may receive the transmitted radio frequency signal, perform radio frequency demodulation 164 to recover the unique channel response (waveform) 166 and record it. When the selected terminal, e.g., 154a, is performing the hand-shaking process with the base station 152, the other terminals, e.g., 154b and 154c, remain silent. Upon receiving the channel response 166 associated with the terminal 154a, the base station 152 calculates a signature waveform g 168 for the terminal 154a based on its channel response received at the base station 152.
For example, one of the basic signature waveforms for a terminal associated with a channel response h(t) can be h*(T−t), which is a time-reversed waveform of the channel response, with T representing the delay spread of the multi-path channel. If h(t) is a complex value, h*(T−t) is the time reverse and conjugation of the channel response h(t). In the discrete time domain, the channel response can be represented as h[k], and the corresponding signature waveform can be h*[L−k], in which L represents the delay spread of the multi-path channel.
The delay spread can be a measure of the multipath richness of a communications channel. In general, it can be interpreted as the difference between the time of arrival of the earliest significant multipath component (typically the line-of-sight component) and the time of arrival of the latest multipath component. For the discrete time model, by the definition of delay spread, when the delay spread is L, the total channel length is L+1.
The system 150 can use any type of waveform design that is designed for a variety of performance metrics, such as designed for reducing or minimizing inter-symbol interference and/or inter-user interference.
After the hand-shaking process, the calculated signature waveforms can be used for data transmissions during both the downlink and the uplink processes.
Downlink Transmission
The following describes an exemplary downlink transmission process.
Referring to FIG. 6, assume that there are N active terminal devices 154 receiving data from the base station 152 concurrently. In some implementations, at the base station 152, there may be N independent information sources 170, each generating information bits to a corresponding terminal device 154. The information bits from the sources may be first encoded by channel-coding using encoders 172. The channel-coding process is optional and may be omitted in some implementations. The encoded bit-stream for any given terminal iε{1, 2, . . . N} may be represented by a sequence of symbols {Xi[k]}.
A parameter referred to as the rate back-off factor D may be used to match the symbol rate (signal bandwidth) with the higher system's sampling rate (channel bandwidth). For any terminal's symbol sequence {Xi[k]}, ∀iε{1, 2, . . . N}, the rate matching process at the base station 152 may be performed by up-sampling the sequence of symbols by the factor D. The rate matching processes may be performed at rate matching modules 174. The up-sampled sequence of {Xi[k]}, ∀iε{1, 2, . . . N} can be expressed as




  
    
      
        
          X
          i
          
            [
            D
            ]
          
        
        =
        
          {
          
            
              
                
                  
                    
                      X
                      i
                    
                    ⁡
                    
                      [
                      
                        k
                        ⁢
                        
                          /
                        
                        ⁢
                        D
                      
                      ]
                    
                  
                  ,
                
              
              
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      k
                      ⁢
                      
                        
                      
                      ⁢
                      mod
                      ⁢
                      
                        
                      
                      ⁢
                      D
                    
                    =
                    0
                  
                  ,
                
              
            
            
              
                
                  0
                  ,
                
              
              
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    k
                    ⁢
                    
                      
                    
                    ⁢
                    mod
                    ⁢
                    
                      
                    
                    ⁢
                    D
                  
                  ≠
                  0.
                
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          1
        
        )
      
    
  




At the waveform embedding modules 176, the up-sampled sequences may be embedded with the calculated (based on the channel responses obtained in the hand-shaking process) signature waveform gi, ∀iε{1, 2, . . . N} by taking the convolution of the up-sampled symbol sequence and the corresponding signature waveform. As described above, one possible signature waveform is the time-reversed (and conjugated if the channel response is of complex values) channel response, i.e.,




  
    
      
        
          
            g
            i
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            1
            
              
                
                  ∑
                  
                    l
                    =
                    0
                  
                  L
                
                ⁢
                
                  
                    
                    
                      
                        h
                        i
                      
                      ⁡
                      
                        [
                        l
                        ]
                      
                    
                    
                  
                  2
                
              
            
          
          ⁢
          
            
              
                h
                i
                *
              
              ⁡
              
                [
                
                  L
                  -
                  k
                
                ]
              
            
            .
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          2
        
        )
      
    
  




The output of the signature embedding block 176 for the i-th terminal is the convolution of {Xi[D][k]} and {gi[k]}, i.e.




  
    
      
        
          
            
              (
              
                
                  X
                  i
                  
                    [
                    D
                    ]
                  
                
                *
                
                  g
                  i
                
              
              )
            
            ⁡
            
              [
              k
              ]
            
          
          =
          
            
              ∑
              
                l
                =
                0
              
              L
            
            ⁢
            
              
            
            ⁢
            
              
                
                  g
                  i
                
                ⁡
                
                  [
                  l
                  ]
                
              
              ⁢
              
                
                  X
                  i
                  
                    [
                    D
                    ]
                  
                
                ⁡
                
                  [
                  
                    k
                    -
                    l
                  
                  ]
                
              
            
          
        
        ,
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          3
        
        )
      
    
  




where (*) is the linear convolution.

After the signature embedding, each of the streams may be multiplied by a scalar power gain β to implement a power control mechanism at the power control module 178. All the streams (embedded with signature waveforms) may be multiplexed together by a summation module 180 to obtain a combined signal S[k]:




  
    
      
        
          S
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            ∑
            
              i
              =
              0
            
            N
          
          ⁢
          
            
          
          ⁢
          
            
              
                
                  β
                  i
                
                ⁡
                
                  (
                  
                    
                      X
                      i
                      
                        [
                        D
                        ]
                      
                    
                    *
                    
                      g
                      i
                    
                  
                  )
                
              
              ⁡
              
                [
                k
                ]
              
            
            .
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          4
        
        )
      
    
  




The combined signal S[k] may be transmitted by the RF module 182 and may propagate through different multi-path channels {hi[k]} 156, for all iε{1, 2, . . . N}.

The signal received at the i-th terminal 154 is the multi-path channel output plus additive noise ñi.




  
    
      
        
          
            Y
            i
            
              [
              D
              ]
            
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            
              
                (
                
                  S
                  *
                  
                    h
                    i
                  
                
                )
              
              ⁡
              
                [
                k
                ]
              
            
            +
            
              
                
                  n
                  ~
                
                i
              
              ⁡
              
                [
                k
                ]
              
            
          
          =
          
            
              
                ∑
                
                  j
                  =
                  1
                
                N
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    β
                    j
                  
                  ⁡
                  
                    (
                    
                      
                        X
                        j
                        
                          [
                          D
                          ]
                        
                      
                      *
                      
                        g
                        j
                      
                      *
                      
                        h
                        i
                      
                    
                    )
                  
                
                ⁡
                
                  [
                  k
                  ]
                
              
            
            +
            
              
                
                  
                    n
                    ~
                  
                  i
                
                ⁡
                
                  [
                  k
                  ]
                
              
              .
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          5
        
        )
      
    
  




Due to the temporal focusing effect, the signal energy can be concentrated in a single time sample or close to a single time sample. Then, the i-th terminal 154 may not need to apply equalization to the received signal to recover the signal. The i-th terminal 154 may down-sample the received signal with the same factor D (to recover the signal's rate from the system's sampling rate to the symbol rate), producing Yi[k], represented as follows:




  
    
      
        
          
            
              Y
              i
            
            ⁡
            
              [
              k
              ]
            
          
          =
          
            
              
                ∑
                
                  j
                  =
                  1
                
                N
              
              ⁢
              
                
              
              ⁢
              
                
                  ∑
                  
                    l
                    =
                    
                      -
                      
                        ⌊
                        
                          L
                          ⁢
                          
                            /
                          
                          ⁢
                          D
                        
                        ⌋
                      
                    
                  
                  
                    ⌊
                    
                      L
                      ⁢
                      
                        /
                      
                      ⁢
                      D
                    
                    ⌋
                  
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      
                        β
                        j
                      
                      ⁡
                      
                        (
                        
                          
                            g
                            j
                          
                          *
                          
                            h
                            i
                          
                        
                        )
                      
                    
                    ⁡
                    
                      [
                      
                        L
                        +
                        Dl
                      
                      ]
                    
                  
                  ⁢
                  
                    
                      X
                      j
                    
                    ⁡
                    
                      [
                      
                        k
                        -
                        l
                      
                      ]
                    
                  
                
              
            
            +
            
              
                n
                i
              
              ⁡
              
                [
                k
                ]
              
            
          
        
        ,
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          6
        
        )
      
    
  




where Yi[k] is used to estimate Xi[k], and ni[k]=ñi[Dk].

Upon obtaining Yi[k] at the i-th terminal (∀iε{1, 2, . . . N}) 154, an estimation {circumflex over (X)}i[k] about Xi[k] can be made based on the value of Yi[k]. Then, the estimated/reconstructed symbol sequence {{circumflex over (X)}i[k]} can be fed into the decoder (which corresponds to the encoder) to recover the information bits at each individual terminal device.
Uplink Transmission
The following describes an exemplary uplink transmission process.
Referring to FIG. 7, for the uplink, assume that there are N active terminal devices 154 transmitting their own data to the base station 152 concurrently. Each terminal device 154 may have its own information source 190 which may be independent of others.
At each terminal device 154, the information bits may be first encoded by the channel coding block 192. Then, the coded symbols {Xi[k]} (∀iε{1, 2, . . . N}) may be up-sampled by the rate back-off factor D for rate matching at the rate matching module 194.
The up-sampled symbol sequence for the i-th terminal can be expressed as




  
    
      
        
          X
          i
          
            [
            D
            ]
          
        
        =
        
          {
          
            
              
                
                  
                    
                      X
                      i
                    
                    ⁡
                    
                      [
                      
                        k
                        ⁢
                        
                          /
                        
                        ⁢
                        D
                      
                      ]
                    
                  
                  ,
                
              
              
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      k
                      ⁢
                      
                        
                      
                      ⁢
                      mod
                      ⁢
                      
                        
                      
                      ⁢
                      D
                    
                    =
                    0
                  
                  ,
                
              
            
            
              
                
                  0
                  ,
                
              
              
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    k
                    ⁢
                    
                      
                    
                    ⁢
                    mod
                    ⁢
                    
                      
                    
                    ⁢
                    D
                  
                  ≠
                  0.
                
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          7
        
        )
      
    
  




The scaling factors {ai} may be used to implement power control at the power control module 196. After multiplying with scaling factor, the sequence of ai{Xi[D][k]}, for all iε{1, 2, . . . N}, may be transmitted by the RF module 198.
When the sequence ai{Xi[D][k]} propagates through its wireless channel {hi[k]} 156, the convolution between ai{Xi[D][k]} and the channel impulse response {hi[k]} is automatically taken as the channel output for the i-th terminal. Since each terminal's multi-path channel response is unique, the automatic convolution between the ai{Xi[D][k]} and {hi[k]} is mathematically similar to the base station's signature embedding process in the downlink.
All of the channel outputs for the N terminal devices 154 may be mixed together in the air plus the additive noise at the base station 152. As a result, the mixed signal received at the base station 152 can be written as




  
    
      
        
          Z
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            
              ∑
              
                i
                =
                0
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                
                  a
                  i
                
                ⁡
                
                  (
                  
                    
                      X
                      i
                      
                        [
                        D
                        ]
                      
                    
                    *
                    
                      g
                      i
                    
                  
                  )
                
              
              ⁡
              
                [
                k
                ]
              
            
          
          +
          
            
              
                n
                ~
              
              ⁡
              
                [
                k
                ]
              
            
            .
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          8
        
        )
      
    
  




Upon receiving the mixed signal as shown above, the base station 152 may pass this mixed signal through an information extraction module that can include a filter bank 200 of the calculated signature waveforms gi, ∀iε{1, 2, . . . N}, each of which may perform the convolution between its input signal Z[k] and the corresponding terminal device's signature waveform {gi[k]}. Such a convolution using the signature waveform may extract the useful signal component and suppress the signals of other terminal devices. As the output of the i-th filter, the convolution can be represented as




  
    
      
        
          
            Y
            i
            
              [
              D
              ]
            
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            
              
                ∑
                
                  j
                  =
                  1
                
                N
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    a
                    j
                  
                  ⁡
                  
                    (
                    
                      
                        X
                        j
                        
                          [
                          D
                          ]
                        
                      
                      *
                      
                        h
                        j
                      
                      *
                      
                        g
                        i
                      
                    
                    )
                  
                
                ⁡
                
                  [
                  k
                  ]
                
              
            
            +
            
              
                (
                
                  
                    g
                    i
                  
                  *
                  
                    n
                    ~
                  
                
                )
              
              ⁡
              
                [
                k
                ]
              
            
          
          =
          
            
              
                ∑
                
                  j
                  =
                  1
                
                N
              
              ⁢
              
                
              
              ⁢
              
                
                  ∑
                  
                    l
                    =
                    0
                  
                  
                    2
                    ⁢
                    L
                  
                
                ⁢
                
                  
                    
                      
                        a
                        j
                      
                      ⁡
                      
                        (
                        
                          
                            h
                            j
                          
                          *
                          
                            g
                            i
                          
                        
                        )
                      
                    
                    ⁡
                    
                      [
                      l
                      ]
                    
                  
                  ⁢
                  
                    
                      X
                      j
                      
                        [
                        D
                        ]
                      
                    
                    ⁡
                    
                      [
                      
                        k
                        -
                        l
                      
                      ]
                    
                  
                
              
            
            +
            
              
                
                  (
                  
                    
                      g
                      i
                      *
                    
                    ⁢
                    
                      n
                      ~
                    
                  
                  )
                
                ⁡
                
                  [
                  k
                  ]
                
              
              .
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          9
        
        )
      
    
  




If the time reversal waveform is used for {gi[k]}, in the above equation, the highest gain for the i-th terminal's symbol may be achieved at the temporal focusing time l=L, with




  
    
      
        
          
            (
            
              
                h
                i
                *
              
              ⁢
              
                g
                i
              
            
            )
          
          ⁡
          
            [
            L
            ]
          
        
        =
        
          
            
              
                ∑
                
                  l
                  =
                  0
                
                L
              
              ⁢
              
                
              
              ⁢
              
                
                  
                  
                    
                      h
                      i
                    
                    ⁡
                    
                      [
                      l
                      ]
                    
                  
                  
                
                2
              
            
          
          .
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          10
        
        )
      
    
  




The rate matching may be performed by down-sampling (with the same factor D) the filter bank's output signal at the rate matching module 202 to recover the original symbol rate of each terminal 154.
After the rate matching, the down-sampled symbol sequence {Yi[k]} can be obtained as




  
    
      
        
          
            Y
            i
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            
              ∑
              
                j
                =
                1
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                ∑
                
                  l
                  =
                  
                    -
                    
                      ⌊
                      
                        L
                        ⁢
                        
                          /
                        
                        ⁢
                        D
                      
                      ⌋
                    
                  
                
                
                  ⌊
                  
                    L
                    ⁢
                    
                      /
                    
                    ⁢
                    D
                  
                  ⌋
                
              
              ⁢
              
                
                  
                    
                      a
                      j
                    
                    ⁡
                    
                      (
                      
                        
                          g
                          i
                        
                        *
                        
                          h
                          j
                        
                      
                      )
                    
                  
                  ⁡
                  
                    [
                    
                      L
                      +
                      Dl
                    
                    ]
                  
                
                ⁢
                
                  
                    X
                    j
                  
                  ⁡
                  
                    [
                    
                      k
                      -
                      l
                    
                    ]
                  
                
              
            
          
          +
          
            
              n
              i
            
            ⁡
            
              [
              k
              ]
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          11
        
        )
      
    
  




where the colored noise ni[k]=Σl=0Lgi[l]ñi[Dk−l].

Based on Yi[k] (∀ε{1, 2, . . . N}), an estimation {circumflex over (X)}i[k] about Xi[k] can be made. Then, the estimated/reconstructed symbol sequence {{circumflex over (X)}i[k]} can be fed into the corresponding decoder 204 to recover the information bits from each terminal device at the base station 152.
Simulation Results
The following describes results of simulation of the system 150 under various conditions.
In the simulation, the IEEE 802.15.4a (broadband) channel model was used to evaluate the effectiveness of the proposed system. In the examples shown here, the typical delay spread of the outdoor channels are between 200 ns and 400 ns, and the typical delay spread of the indoor channels are between 70 ns and 120 ns. Without loss of generality, and for the exemplary embodiments described here, the following system parameters were used to in the simulation:










IEEE 802.15.4a Outdoor Channel Model


















System sampling rate
1 GHz



Rate back-off factor: D
32, 64



Eb/N0 (energy per bit to noise power
0 dB~20 dB



spectral density ratio)



Number of concurrent terminals:
2, 4



N



Coding Scheme
None (Uncoded)



Power Control Scheme
Equal-Power Allocation




















IEEE 802.15.4a Indoor Channel Model


















System sampling rate
1 GHz



Rate back-off factor: D
16, 32



Eb/N0
0 dB~20 dB



Number of concurrent terminals:
2, 4



N



Coding Scheme
None (Uncoded)



Power Control Scheme
Equal-Power Allocation










Referring to FIG. 8, a graph 210 shows the predicted bit error rate (BER) performance of the system 150 for downlink transmission in an outdoor environment. A curve 212 represents the predicted bit error rate when the rate back-off factor D is equal to 64, and the number of users N is equal to 2. A curve 214 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 2. A curve 216 represents the predicted bit error rate when the rate back-off factor D is equal to 64, and the number of users N is equal to 4. A curve 218 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 4.
Referring to FIG. 9, a graph 220 shows the predicted bit error rate performance of the system 150 for uplink transmission in an outdoor environment. A curve 222 represents the predicted bit error rate when the rate back-off factor D is equal to 64, and the number of users N is equal to 2. A curve 224 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 2. A curve 226 represents the predicted bit error rate when the rate back-off factor D is equal to 64, and the number of users N is equal to 4. A curve 228 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 4.
Referring to FIG. 10, a graph 230 shows the predicted bit error rate performance of the system 150 for downlink transmission in an indoor environment. A curve 232 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 2. A curve 234 represents the predicted bit error rate when the rate back-off factor D is equal to 16, and the number of users N is equal to 2. A curve 236 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 4. A curve 238 represents the predicted bit error rate when the rate back-off factor D is equal to 16, and the number of users N is equal to 4.
Referring to FIG. 11, a graph 240 shows the predicted bit error rate performance of the system 150 for uplink transmission in an indoor environment. A curve 242 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 2. A curve 244 represents the predicted bit error rate when the rate back-off factor D is equal to 16, and the number of users N is equal to 2. A curve 246 represents the predicted bit error rate when the rate back-off factor D is equal to 32, and the number of users N is equal to 4. A curve 248 represents the predicted bit error rate when the rate back-off factor D is equal to 16, and the number of users N is equal to 4.
Referring to FIG. 12, an exemplary process 250 for asymmetric time-reversal wireless communication is provided. For example, the process 250 can be implemented using the base station 152 shown in FIGS. 4-7. The process 250 may include, at a transceiver, receiving a channel response signal derived from a probe signal sent from a first device to the transceiver through multiple wireless propagation paths (252). For example, the transceiver can be part of the base station 152.
The process 250 may include calculating a signature waveform for the first device based on the channel response (254). For example, the first device can be the first terminal device 154a. 
The process 250 may include receiving combined signals that include a signal from the first device and a signal from a second device (256). For example, the second device can be the second terminal device 154b. 
The process 250 may include determining the signal from the first device based on the combined signals and the signature waveform for the first device (258). For example, the base station 152 can receive the combined signals from the first and second terminal devices 154a and 154b and apply the signature waveform for the first terminal device 154a to the combined signals to determine the signal from the first terminal device 154a. For example, a convolution operation, or an operation equivalent or similar to the convolution operation, can be performed between the combined signals and the signature waveform for the first device. For example, the signature waveform for the first device can be a time-reversed waveform of the channel response signal derived from the probe signal sent from the first device.
Referring to FIG. 13, an exemplary process 270 for communication between a transceiver and multiple devices using wireless multi-path signals is provided. For example, the transceiver can be the base station 152, and the multiple devices can be the terminal devices 154 in FIGS. 4-7.
The process 270 may include performing a hand-shake process between the transceiver and the devices, in which for each device, the transceiver receives a channel response signal derived from a probe signal sent from the device to the transceiver through multiple propagation paths (272). For example, after each device registers with the transceiver, a hand-shaking process can occur periodically in a round-robin manner among the registered devices to maintain an up-to-date record of the channel responses associated with the devices at the transceiver. In embodiments, the hand-shaking process may occur on a predetermined schedule, or a calculated schedule, or in response to a measured parameter of the receivers or transmitters or system performance, and the like. During the hand-shaking process, at a given time, only one selected device may perform radio frequency modulation to produce a radio frequency signal that propagates through the associated multi-path channel to the transceiver. The transceiver receives the transmitted radio frequency signal and may perform radio frequency demodulation to recover the unique channel response (waveform) and may record it. In an exemplary embodiment, when the selected device is performing the hand-shaking process with the transceiver, the other devices may remain silent.
The process 270 may include calculating a signature waveform for each of the devices based on the channel response derived from the probe signal sent from the device (274). For example, a signature waveform for a device associated with a channel response h(t) can be h*(T−t), which is a time-reversed waveform of the channel response, with T representing the delay spread of the multi-path channel. In discrete time domain, the channel response can be represented as h[k], and the corresponding signature waveform can be h*[L−k], in which L represents the delay spread of the multi-path channel.
The process 270 may include receiving a combined uplink signal that includes uplink signals from two or more devices (276), and determining the uplink signal from each one of the devices based on the combined uplink signal and the signature waveform for the device (278). For example, the base station 152 can receive a combined uplink signal that includes uplink signals from the first and second terminal devices 154a and 154b. The base station 152 can apply the signature waveform for the first terminal device 154a to the combined signals to determine the signal from the first terminal device 154a, and apply the signature waveform for the second terminal device 154b to the combined signals to determine the signal from the second terminal device 154b. For example, a convolution operation, or an operation equivalent or similar to the convolution operation, can be performed between the combined signals and the signature waveform for the first or second device.
Referring to FIG. 14, an exemplary process 290 for asymmetric time-reversal communication between a plurality of devices using wireless multi-path signals is provided. For example, the devices can include the base station 152 and the terminal devices 154 of FIGS. 4-7.
The process 290 may include performing a first hand-shake process between the first and second devices, in which the first device may receive a channel response signal derived from a probe signal sent from the second device to the first device through multiple wireless propagation paths (292). For example, the first device can be the base station 152, and the second device can be the first terminal device 154a. 
The process 290 may include performing a second hand-shake process between the first and third devices, in which the first device receives a channel response signal derived from a probe signal sent from the third device to the first device through multiple wireless propagation paths (294). For example, the third device can be the second terminal device 154b. 
The process 290 may include at the first device, calculating a signature waveform for the second device based on the channel response derived from the probe signal sent from the second device (296); and at the first device, calculating a signature waveform for the third device based on the channel response derived from the probe signal sent from the third device (298). For example, the signature waveform for the second or third terminal devices 154a or 154b can be a time-reversed waveform of the channel response signal derived from the probe signal sent from the second or third terminal devices 154a or 154b, respectively. A signature waveform for the second or third terminal devices associated with a channel response h(t) can be h*(T−t), which is a time-reversed waveform of the channel response, with T representing the delay spread of the multi-path channel. If h(t) is a complex value, h*(T−t) is the time reverse and conjugation of the channel response h(t). In discrete time domain, the channel response can be represented as h[k], and the corresponding signature waveform can be h*[L−k], in which L represents the delay spread of the multi-path channel. Note that each terminal device may have a unique, location dependent, signature waveform.
The process 290 may include at the first device, receiving a combined uplink signal that includes uplink signals from the second and third devices (300). For example, the base station 152 can receive a combined uplink signal that includes uplink signals from the second and third terminal devices 154a, 154b. 
The process 290 may include at the first device, determining the uplink signal from the second device based on the combined uplink signal and the signature waveform for the second device (302); and at the first device, determining the uplink signal from the third device based on the combined uplink signal and the signature waveform for the third device (304). For example, the base station 152 can apply the signature waveform for the second terminal device 154a to the combined signals to determine the signal from the second terminal device 154a, and apply the signature waveform for the third terminal device 154b to the combined signals to determine the signal from the third terminal device 154b. For example, a convolution operation, or an operation equivalent or similar to the convolution operation, can be performed between the combined signals and the signature waveform for the second and/or third terminal device 154a and/or 154b. 
Uplink System with Parallel Interference Cancellation
The novel time-reversal division multiple access based multi-user uplink architecture described above can also utilize a 2-dimensional (2D) parallel interference cancellation scheme to enhance the system performance. In an exemplary time-reversal division multiple access uplink architecture the processing power at a base station (BS) that has already been made available for the downlink, can be reused, thus potentially reducing the cost for the uplink modules at the terminal devices. The 2D parallel interference cancellation scheme uses the tentative decisions of detected symbols to effectively cancel the interference in both the time dimension (inter-symbol interference or ISI) and the user dimension (inter-user interference (IUI)), which may significantly improve the bit-error-rate performance and achieve a high signal-to-noise-ratio (SNR). To further improve the BER performance, a multi-stage processing can be implemented by cascading multiple stages of the proposed 2D interference cancellation, with a total delay that increases linearly with the number of stages, but independent of the number of users. Simulation results are provided for up to 3 stages of interference cancellation and compared with the basic time-reversal division multiple access system without interference cancelation.
In this description, depending on context, the term “user” may refer to a device. For example, in a system that has multiple devices communicating with a base station, the term “multi-user uplink” refers to the uplink by multiple devices, and the term “inter-user interference” refers to the interference among various devices.
As described above, a time-reversal division multiple access architecture may be described as having two parts, the downlink part and the uplink part. In a time-reversal division multiple access downlink system, a base station (BS) may transmit multiple data streams to several users simultaneously, in which each user may be associated with a unique multi-path profile in a rich-scattering environment. The time-reversal division multiple access downlink scheme may exploit the spatial degrees of freedom of the environment and focus the useful signal power mostly at the intended locations. Time reversal mirrors (TRMs) at the base station may use a time-reversed copy of the channel response of each user's channel as the user's signature waveform, and then embed these signatures into the corresponding data streams. The transmitted signal from the base station in the time-reversal division multiple access downlink may be a mixed signal that includes the data intended to be sent to several users (including different data intended for different users). When the combined signal propagates to a certain user through the corresponding multipath channel, the encoded waveform may be temporally and spatially focused at the intended user's location.
Within the time-reversal division multiple access framework, more sophisticated signature waveforms can be derived based on the multi-path channel responses to further improve the performance of the time-reversal division multiple access downlink system.
To further enhance the system performance of the time-reversal division multiple access uplink, a 2-dimensional (2D) parallel interference cancellation scheme may be provided by exploring the signal structure of both the inter-symbol interference (ISI) (time dimension) and the inter-user interference (IUI) (user dimension). Compared with existing multi-user detection (MUD) techniques for code division multiple access (CDMA) systems for example, a TRDMA system may experience more severe inter-user interference because the multi-path signature waveforms are not naturally orthogonal. Also, the time-reversal division multiple access system allows overlap between the transmitted signature waveforms to boost system throughput, which may cause inter-symbol interference, depending on the extent of overlapping.
Here we disclose a multi-user uplink scheme based on time-reversal division multiple access, and a 2D parallel interference cancellation technique. The 2D parallel interference cancellation scheme uses the tentative decisions of detected symbols to effectively reduce or cancel both the inter-symbol interference and inter-user interference at the base station. To further improve the bit error rate performance, multi-stage processing can be performed by cascading multiple stages of cancellation, with a total delay that increases linearly with the number of stages, but independent of the number of users.
Referring again to FIG. 15, an exemplary broadband multi-user time reversal communication system 1550 includes a base station 1552 and multiple users or terminal devices (e.g., 1554a, 1554b, 154c, collectively 1554). In some implementations, all the devices (including the base station 1552 and the terminals 1554) can operate at the same frequency band, and the system operates in multi-path environments. For the downlink, the base station 1552 can send multiple messages (either independent or non-independent) simultaneously to multiple selected terminals 1554. For the uplink, multiple terminals 1554 can send their own messages to the base station 1552 simultaneously.
In an exemplary embodiment, N users 1554 may simultaneously transmit independent bit streams to the base station 1552. In this example, binary phase-shift keying (BPSK) modulation is considered, and the polarity of the modulated symbols {Xi[k]ε{−1, +1}} carries the binary information for user i. Other types of modulation formats including ASK, FSK, QAM, and any known variations and combinations of those modulation formats can also be used.
For any given user i in the uplink network, the channel hi 1556 between the base station 1552 and user i is a multi-path channel that can be characterized by a unique discrete-time channel impulse response (CIR)




  
    
      
        
          
            h
            i
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            ∑
            
              l
              =
              0
            
            
              L
              -
              1
            
          
          ⁢
          
            
          
          ⁢
          
            
              h
              
                i
                ,
                l
              
            
            ⁢
            
              δ
              ⁡
              
                [
                
                  k
                  -
                  l
                
                ]
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          12
        
        )
      
    
  




where hi,lεR is the l-th tap of the channel impulse response with length L, and δ[•] is the Dirac delta function. We assume that the channels are quasi-static and reciprocal, in which information about the channels can be acquired at the base station through a channel probing phase.

During the channel probing phase, each user 1554 may take turns sending a probe signal or waveform to the base station 1552. The probe signal can be, e.g., a pulse signal, a signal that has a predetermined waveform, or a signal that includes symbols to enable the base station to perform synchronization and/or other functions. In the examples below, the exemplary probe signal is a pulse signal. However, other types of probe signals can be used in the system 1550. When the probe signals are pulse signals, the channel impulse response {hi[k]} of each user's link can be recorded by a time-reversal mirror 1558 at the base station 1552. Upon recording the channel impulse response, the time-reversal mirror 1558 may reverse the recorded waveform in the time domain and normalize it as the unique signature waveform of user i. As described above, the time-reversed waveform of user i may be used in the uplink data transmission phase to extract the desired signal from a combination of the multiple access signals that are mixed in the air.
After the channel probing phase, the users 1554 can start to transmit the statistically independent messages {X1[k], X2[k], . . . , XN[k]} to the base station 1552 through the multi-path channels 1556. Also as described above, each device 1554 may have a rate matching module 1560 that introduces a rate back-off factor D to match the symbol rate (signal bandwidth) with the higher system sampling rate (channel bandwidth). The scaling factors ai, for iε{1, 2, . . . , N} may be used to implement the transmit power control, whose values are assumed to be instructed by the base station 1552 through a feedback/control channel. After multiplying with scaling factor, the sequence of aiXi[D][k] for all iε1{1, 2, . . . , N}, may be transmitted through the corresponding multi-path channel {hi[k]}.
When the sequence {aiXi[D][k]} propagates through its wireless channel {hi[k]} 1556, the convolution between {aiXi[D][k]} and the channel impulse response {hi[k]} may be automatically taken as the channel output for user i. All of the channel outputs for the N users are mixed together in the air
Upon receiving the mixed signal, the base station 1552 may pass this mixed signal through a bank of N time-reversal mirrors 1558, each of which may perform the convolution between its input signal {S[k]} and the user's signature waveform {gi[k]}. Such a convolution using the signature waveform can extract the useful signal component and suppress the signals of the other users. As the output of the i-th time-reversal mirror 1558, the convolution of {S[k]} and the signature of user i {gi[k]} can be represented as




  
    
      
        
          
            
              Y
              i
              
                [
                D
                ]
              
            
            ⁡
            
              [
              k
              ]
            
          
          =
          
            
              
                
                  ∑
                  
                    j
                    =
                    1
                  
                  N
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      a
                      j
                    
                    ⁡
                    
                      (
                      
                        
                          g
                          i
                        
                        *
                        
                          h
                          j
                        
                        *
                        
                          X
                          j
                          
                            [
                            D
                            ]
                          
                        
                      
                      )
                    
                  
                  ⁡
                  
                    [
                    k
                    ]
                  
                
              
              +
              
                
                  (
                  
                    
                      g
                      i
                    
                    *
                    
                      n
                      ~
                    
                  
                  )
                
                ⁡
                
                  [
                  k
                  ]
                
              
            
            =
            
              
                
                  ∑
                  
                    j
                    =
                    1
                  
                  N
                
                ⁢
                
                  
                
                ⁢
                
                  
                    ∑
                    
                      l
                      =
                      0
                    
                    
                      
                        2
                        ⁢
                        
                          
                        
                        ⁢
                        L
                      
                      -
                      2
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      
                        
                          a
                          j
                        
                        ⁡
                        
                          (
                          
                            
                              g
                              i
                            
                            *
                            
                              h
                              j
                            
                          
                          )
                        
                      
                      ⁡
                      
                        [
                        l
                        ]
                      
                    
                    ⁢
                    
                      
                        X
                        j
                        
                          [
                          D
                          ]
                        
                      
                      ⁡
                      
                        [
                        
                          k
                          -
                          l
                        
                        ]
                      
                    
                  
                
              
              +
              
                
                  (
                  
                    
                      g
                      i
                    
                    *
                    
                      n
                      ~
                    
                  
                  )
                
                ⁡
                
                  [
                  k
                  ]
                
              
            
          
        
        ,
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          13
        
        )
      
    
  




in which the highest gain for user i's symbol is achieved at the temporal focusing time l=L−1, with





  
    
      
        
          
            (
            
              
                g
                i
              
              *
              
                h
                i
              
            
            )
          
          ⁡
          
            [
            
              L
              -
              1
            
            ]
          
        
        =
        
          
            
              ∑
              
                l
                =
                0
              
              
                L
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                
                
                  
                    h
                    i
                  
                  ⁡
                  
                    [
                    l
                    ]
                  
                
                
              
              2
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          14
        
        )
      
    
  




Comparing Equation 13 and the received signal at the terminal users in the downlink, the same mathematical structure can be found by switching the roles of the signature waveforms {gi} s and the channel impulse responses his in the convolution (and ignoring the scaling factor ai and noise term). Therefore, mathematically, a virtual spatial focusing effect as observed in the downlink can be seen in the user's signature domain of the uplink scheme. Unlike the physical spatial focusing effect observed in the downlink in which the useful signal power is concentrated at different physical locations, in the uplink, the signal power concentration in the users' signature waveform space is achieved mathematically at the base station 1552. Such a virtual spatial focusing effect is used to separate the useful signal from the signals of other users. The base station 1552 includes a rate matching module 1566 to perform rate matching by down-sampling (with the same factor D) the time-reversal mirror's output signal to recover the original symbol rate of the modulated symbols of each user.
After the rate matching, the down-sampled time-reversal mirror output Yi[k] can be represented as




  
    
      
        
          
            Y
            i
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            
              ∑
              
                j
                =
                1
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                ∑
                
                  l
                  =
                  
                    -
                    
                      ⌊
                      
                        
                          L
                          -
                          1
                        
                        D
                      
                      ⌋
                    
                  
                
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    
                      a
                      j
                    
                    ⁡
                    
                      (
                      
                        
                          g
                          i
                        
                        *
                        
                          h
                          j
                        
                      
                      )
                    
                  
                  ⁡
                  
                    [
                    
                      L
                      -
                      1
                      +
                      Dl
                    
                    ]
                  
                
                ⁢
                
                  
                    X
                    j
                  
                  ⁡
                  
                    [
                    
                      k
                      -
                      l
                    
                    ]
                  
                
              
            
          
          +
          
            
              n
              i
            
            ⁡
            
              [
              k
              ]
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          15
        
        )
      
    
  




where the colored noise ni[k]=Σl=0L−1gl[l]ñ[Dk−l]=gi{tilde over (n)}[k] with gi=[gi[0], gi[1], . . . , gi[L−1]] and {tilde over (n)}[k]=[ñ[k], ñ[k−1], . . . , ñ[k−L+1]]T.

A more rigorous representation of Yi[k] is Yi[k]=Yi[D][L−1+Dk], which aligns the highest temporal focusing gain (gi*hi)[L−1] in Yi[k] with the transmitted symbol Xi[k] in time for ease of simple notation.
Decomposing the signal shown in Equation 15, we have the following components as




  
    
      
        
          
            
              
                
                  Y
                  i
                
                ⁡
                
                  [
                  k
                  ]
                
              
              =
              
              ⁢
              
                
                  
                    
                      
                        a
                        i
                      
                      ⁡
                      
                        (
                        
                          
                            g
                            i
                          
                          *
                          
                            h
                            i
                          
                        
                        )
                      
                    
                    ⁡
                    
                      [
                      
                        L
                        -
                        1
                      
                      ]
                    
                  
                  ⁢
                  
                    
                      X
                      i
                    
                    ⁡
                    
                      [
                      k
                      ]
                    
                  
                
                +
              
            
          
          
            
              
              ⁢
              
                (
                Signal
                )
              
            
          
        
        
          
            
              
              ⁢
              
                
                  
                    a
                    i
                  
                  ⁢
                  
                    
                      ∑
                      
                        
                          l
                          =
                          
                            -
                            
                              ⌊
                              
                                
                                  L
                                  -
                                  1
                                
                                D
                              
                              ⌋
                            
                          
                        
                        
                          l
                          ≠
                          0
                        
                      
                      
                        ⌊
                        
                          
                            L
                            -
                            1
                          
                          D
                        
                        ⌋
                      
                    
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        
                          (
                          
                            
                              g
                              i
                            
                            *
                            
                              h
                              i
                            
                          
                          )
                        
                        ⁡
                        
                          [
                          
                            L
                            -
                            1
                            +
                            Dl
                          
                          ]
                        
                      
                      ⁢
                      
                        
                          X
                          i
                        
                        ⁡
                        
                          [
                          
                            k
                            -
                            l
                          
                          ]
                        
                      
                    
                  
                
                +
              
            
          
          
            
              
              ⁢
              
                (
                ISI
                )
              
            
          
        
        
          
            
              
              ⁢
              
                
                  
                    
                      ∑
                      N
                    
                    
                      
                        j
                        =
                        1
                      
                      
                        j
                        ≠
                        i
                      
                    
                  
                  ⁢
                  
                    
                      a
                      j
                    
                    ⁢
                    
                      
                        
                          ∑
                          
                            ⌊
                            
                              
                                L
                                -
                                1
                              
                              D
                            
                            ⌋
                          
                        
                        
                          l
                          =
                          
                            -
                            
                              ⌊
                              
                                
                                  L
                                  -
                                  1
                                
                                D
                              
                              ⌋
                            
                          
                        
                      
                      ⁢
                      
                        
                          
                            (
                            
                              
                                g
                                i
                              
                              *
                              
                                h
                                j
                              
                            
                            )
                          
                          ⁡
                          
                            [
                            
                              L
                              -
                              1
                              +
                              Dl
                            
                            ]
                          
                        
                        ⁢
                        
                          
                            X
                            j
                          
                          ⁡
                          
                            [
                            
                              k
                              -
                              l
                            
                            ]
                          
                        
                      
                    
                  
                
                +
              
            
          
          
            
              
              ⁢
              
                (
                IUI
                )
              
            
          
        
        
          
            
              
              ⁢
              
                
                  n
                  i
                
                ⁡
                
                  [
                  k
                  ]
                
              
            
          
          
            
              
              ⁢
              
                (
                Noise
                )
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          16
        
        )
      
    
  




The basic time-reversal division multiple access uplink scheme uses the signal Yi[k] in Equation 15 to estimate the transmitted symbol Xi[k]. A simple decision rule can be derived to implement the maximum-likelihood estimation (MLE) for the binary phase-shift keying symbols Xi[k]ε{+1, −1} for i=1, 2, . . . , N. By the central limit theorem, we model the total interference term (including the inter-symbol interference and inter-user interference) as a Gaussian random variable with zero mean (because the interference term is a linear combination of the zero-mean binary symmetric symbols Xi[k]ε{+1,−1}) and variance




  
    
      
        
          σ
          
            I
            ,
            i
          
          2
        
        =
        
          
            
              
                
                
                  a
                  i
                
                
              
              2
            
            ⁢
            
              
                ∑
                
                  
                    l
                    =
                    
                      -
                      
                        ⌊
                        
                          
                            L
                            -
                            1
                          
                          D
                        
                        ⌋
                      
                    
                  
                  
                    
                      l
                      1
                    
                    ⁢
                    0
                  
                
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                  
                    
                      (
                      
                        
                          g
                          i
                          *
                        
                        ⁢
                        
                          h
                          i
                        
                      
                      )
                    
                    ⁡
                    
                      [
                      
                        L
                        -
                        1
                        +
                        Dl
                      
                      ]
                    
                  
                  
                
                2
              
            
          
          +
          
            
              ∑
              
                
                  j
                  =
                  1
                
                
                  
                    j
                    1
                  
                  ⁢
                  i
                
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                
                  
                  
                    a
                    j
                  
                  
                
                2
              
              ⁢
              
                
                  ∑
                  
                    l
                    =
                    
                      -
                      
                        ⌊
                        
                          
                            L
                            -
                            1
                          
                          D
                        
                        ⌋
                      
                    
                  
                  
                    ⌊
                    
                      
                        L
                        -
                        1
                      
                      D
                    
                    ⌋
                  
                
                ⁢
                
                  
                    
                    
                      
                        (
                        
                          
                            g
                            i
                            *
                          
                          ⁢
                          
                            h
                            j
                          
                        
                        )
                      
                      ⁡
                      
                        [
                        
                          L
                          -
                          1
                          +
                          Dl
                        
                        ]
                      
                    
                    
                  
                  2
                
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          17
        
        )
      
    
  




The likelihood ratio can be derived as




  
    
      
        
          Λ
          ⁡
          
            (
            
              
                Y
                i
              
              ⁡
              
                [
                k
                ]
              
            
            )
          
        
        =
        
          
            
              L
              ⁡
              
                (
                
                  
                    
                      X
                      i
                    
                    ⁡
                    
                      [
                      k
                      ]
                    
                  
                  =
                  
                    1
                    ❘
                    
                      
                        Y
                        i
                      
                      ⁡
                      
                        [
                        k
                        ]
                      
                    
                  
                
                )
              
            
            
              L
              ⁡
              
                (
                
                  
                    
                      X
                      i
                    
                    ⁡
                    
                      [
                      k
                      ]
                    
                  
                  =
                  
                    
                      -
                      1
                    
                    ❘
                    
                      
                        Y
                        i
                      
                      ⁡
                      
                        [
                        k
                        ]
                      
                    
                  
                
                )
              
            
          
          =
          
            
              
                
                  f
                  i
                
                ⁡
                
                  [
                  
                    
                      
                        
                          Y
                          i
                        
                        ⁡
                        
                          [
                          k
                          ]
                        
                      
                      ❘
                      
                        
                          X
                          i
                        
                        ⁡
                        
                          [
                          k
                          ]
                        
                      
                    
                    =
                    1
                  
                  ]
                
              
              
                
                  f
                  i
                
                ⁡
                
                  [
                  
                    
                      
                        
                          Y
                          i
                        
                        ⁡
                        
                          [
                          k
                          ]
                        
                      
                      ❘
                      
                        
                          X
                          i
                        
                        ⁡
                        
                          [
                          k
                          ]
                        
                      
                    
                    =
                    
                      -
                      1
                    
                  
                  ]
                
              
            
            =
            
              exp
              ⁡
              
                (
                
                  
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          
                            a
                            i
                          
                          ⁡
                          
                            (
                            
                              
                                g
                                i
                                *
                              
                              ⁢
                              
                                h
                                i
                              
                            
                            )
                          
                        
                        ⁡
                        
                          [
                          
                            L
                            -
                            1
                          
                          ]
                        
                      
                    
                    
                      (
                      
                        
                          σ
                          
                            I
                            ,
                            i
                          
                          2
                        
                        +
                        
                          σ
                          
                            N
                            ,
                            i
                          
                          2
                        
                      
                      )
                    
                  
                  ⁢
                  
                    
                      Y
                      i
                    
                    ⁡
                    
                      [
                      k
                      ]
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          18
        
        )
      
    
  




where ƒi(y|x) is the conditional probability density function (pdf) of Yi[k] given that Xi[k]=x is transmitted. The simple form of the likelihood ratio shown in Equation 18 leads to a simple decision rule for the maximum-likelihood estimation, specifically,





  
    
      
        
          
            
              X
              ^
            
            i
            
              (
              0
              )
            
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            sgn
            (
            
              
                Y
                i
              
              ⁡
              
                [
                k
                ]
              
            
            )
          
          =
          
            {
            
              
                
                  
                    
                      +
                      1
                    
                    ,
                  
                
                
                  
                    
                      
                        if
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                            Y
                            i
                          
                          ⁡
                          
                            [
                            k
                            ]
                          
                        
                      
                      ≥
                      0
                    
                    ,
                  
                
              
              
                
                  
                    
                      -
                      1
                    
                    ,
                  
                
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      
                        
                          Y
                          i
                        
                        ⁡
                        
                          [
                          k
                          ]
                        
                      
                    
                    <
                    0.
                  
                
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          19
        
        )
      
    
  




In Equation 19, the superscript “(0)” of {circumflex over (X)}i(0)[k] indicates the initial stage of the interference cancellation procedure. Such a notation is useful in the following discussion of multi-stage interference cancellation schemes.
The error probability of the estimator shown in Equation 19 can be calculated based on the Gaussian approximations of the interference as follows




  
    
      
        
          
            P
            err
            
              (
              0
              )
            
          
          ⁡
          
            (
            i
            )
          
        
        =
        
          
            Q
            ⁡
            
              (
              
                
                  
                    
                      
                        
                          
                            
                              
                                
                                
                                  a
                                  i
                                
                                
                              
                              2
                            
                            ⁢
                            
                              ∑
                              
                                l
                                =
                                0
                              
                              
                                L
                                -
                                1
                              
                            
                          
                          
                        
                        ⁢
                        
                          
                            h
                            i
                          
                          ⁡
                          
                            [
                            l
                            ]
                          
                        
                      
                      
                    
                    2
                  
                  
                    
                      σ
                      
                        I
                        ,
                        i
                      
                      2
                    
                    +
                    
                      σ
                      N
                      2
                    
                  
                
              
              )
            
          
          =
          
            Q
            ⁡
            
              (
              
                
                  SINR
                  i
                  
                    (
                    0
                    )
                  
                
              
              )
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          20
        
        )
      
    
  




Where SINRi(0) is the signal-to-interference-plus-noise ratio (SINR) for user i at the initial stage. From Equation 20, one can see that the error probability decreases with the achieved signal-to-interference-plus-noise ratio, i.e., the quality of the signal before the final decision.

The base station 1552 includes a 2D interference cancellation and decision module 1568 that cancels the inter-symbol interference and the inter-user interference. The following describes the 2D interference cancellation scheme that uses the estimated symbols to effectively cancel both the inter-symbol interference and inter-user interference, and significantly improves the performance of the uplink communication.
2-Dimensional Parallel Interference Cancelation
The 2D parallel interference cancellation scheme used by the 2D interference cancelation and decision module 1568 can have one or more stages. A single stage 2D parallel interference cancellation scheme is described first. The time-reversal division multiple access system can be an interference-limited system, especially in the high signal-to-noise ratio (SNR) regime. The interference terms shown in Equation 16 have structures that can be exploited to improve the bit error rate performances. Because the channel impulse responses have been obtained at the base station 1552 during the channel probing phase, the interference terms in Equation 16 can be reconstructed if the relevant transmitted symbols are known. Note that the coefficients {ai} can be determined by the base station 1552 and sent to the users 1554. In the interference cancellation scheme described here, the estimated symbols from the previous stage can be used to approximate the interference terms in the current stage. The basic TRDMA system (without interference cancelation) may be considered as the initial (0-th) stage of the interference cancellation. The interference cancellation of the TRDMA system includes two parts belonging to two different dimensions; the inter-symbol interference due to the multi-path effect of broadband channels, which is in the time domain, and the inter-user interference caused by the simultaneous transmission of multiple users, which is in the user's signature domain. The 2D parallel interference cancellation scheme for the TRDMA uplink system targets the interference in both dimensions by exploiting the structure of interference in both dimensions.
Tentative Decision Vector
Due to the unique structure of time-reversal waveforms, each received symbol may include interference caused by the symbols transmitted before and after the current symbol. According to Equation 16, in order to ideally cancel the interference for user i's symbol Xi[k], one has to know all the other users' transmitted symbols from time




  (
  
    k
    -
    
      ⌊
      
        
          L
          -
          1
        
        D
      
      ⌋
    
  
  )




to





  (
  
    k
    +
    
      ⌊
      
        
          L
          -
          1
        
        D
      
      ⌋
    
  
  )




for the inter-user interference; and user i's own transmitted symbols from time





  (
  
    k
    -
    
      ⌊
      
        
          L
          -
          1
        
        D
      
      ⌋
    
  
  )




to (k−1) and from (k+1) to





  (
  
    k
    +
    
      ⌊
      
        
          L
          -
          1
        
        D
      
      ⌋
    
  
  )




for the inter-symbol interference.

Tentative decisions can be made in an attempt to estimate these symbols. To simplify the notation, let's define the vector {circumflex over (X)}j(0)[k] for all jε{1, 2, . . . , N}, as the stage-0 tentative decision vector for user j, where




  
    
      
        (
        12
        )
      
    
    
      
        
      
    
  
  
    
      
        
          
            
              
                
                  
                    
                      X
                      ^
                    
                    _
                  
                  j
                  
                    (
                    0
                    )
                  
                
                ⁡
                
                  [
                  k
                  ]
                
              
              ⁢
              
                =
                Δ
              
              ⁢
              
              ⁢
              
                [
                
                  
                    
                      
                        X
                        ^
                      
                      j
                      
                        (
                        0
                        )
                      
                    
                    ⁡
                    
                      [
                      
                        k
                        +
                        
                          ⌊
                          
                            
                              L
                              -
                              1
                            
                            D
                          
                          ⌋
                        
                      
                      ]
                    
                  
                  ,
                  
                    
                      
                        X
                        ^
                      
                      j
                      
                        (
                        0
                        )
                      
                    
                    ⁡
                    
                      [
                      
                        k
                        +
                        
                          ⌊
                          
                            
                              L
                              -
                              1
                            
                            D
                          
                          ⌋
                        
                        -
                        1
                      
                      ]
                    
                  
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                
              
            
          
        
        
          
            
              
              ⁢
              
                
                  
                    
                      X
                      ^
                    
                    j
                    
                      (
                      0
                      )
                    
                  
                  ⁡
                  
                    [
                    
                      k
                      -
                      
                        ⌊
                        
                          
                            L
                            -
                            1
                          
                          D
                        
                        ⌋
                      
                    
                    ]
                  
                
                T
              
            
          
        
        
          
            
              =
              
              ⁢
              
                [
                
                  
                    sgn
                    (
                    
                      
                        Y
                        j
                      
                      ⁡
                      
                        [
                        
                          k
                          +
                          
                            ⌊
                            
                              
                                L
                                -
                                1
                              
                              D
                            
                            ⌋
                          
                        
                        ]
                      
                    
                    )
                  
                  ,
                  
                    sgn
                    ⁡
                    
                      (
                      
                        
                          Y
                          j
                        
                        ⁡
                        
                          [
                          
                            k
                            +
                            
                              ⌊
                              
                                
                                  L
                                  -
                                  1
                                
                                D
                              
                              ⌋
                            
                            -
                            1
                          
                          ]
                        
                      
                      )
                    
                  
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                
              
            
          
        
        
          
            
              
                
                  
                  ⁢
                  
                    sgn
                    ⁡
                    
                      (
                      
                        
                          Y
                          j
                        
                        ⁡
                        
                          [
                          
                            k
                            -
                            
                              ⌊
                              
                                
                                  L
                                  -
                                  1
                                
                                D
                              
                              ⌋
                            
                          
                          ]
                        
                      
                      )
                    
                  
                  ]
                
                T
              
              .
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          21
        
        )
      
    
  




In an exemplary embodiment, where the stage-0 tentative decisions for user j depend only on this user's own time-reversal mirror 1558 output, then the tentative decision vectors can be obtained in parallel for all the users in the 2D parallel interference cancellation scheme.
Approximated Interference Reconstruction
Based on the tentative decisions of the transmitted symbols, the approximated interference terms in Equation 16 can be reconstructed by looking at the structure of the inter-symbol interference and the inter-user interference. In doing so, we first define the row vector Ui,j for ∀i,jε{1, 2, . . . , N} as




  
    
      
        
          U
          
            i
            ,
            j
          
        
        ⁢
        
          =
          Δ
        
        ⁢
        
          
            [
            
              
            
            ⁢
            
              
                
                  (
                  
                    
                      g
                      i
                      *
                    
                    ⁢
                    
                      h
                      j
                    
                  
                  )
                
                [
                
                  
                
                ⁢
                
                  L
                  -
                  1
                  -
                  
                    D
                    ⁢
                    
                      ⌊
                      
                        
                          L
                          -
                          1
                        
                        D
                      
                      ⌋
                    
                  
                
                ]
              
              ,
              
                
              
              ⁢
              
                
                  (
                  
                    
                      g
                      i
                      *
                    
                    ⁢
                    
                      h
                      j
                    
                  
                  )
                
                [
                
                  
                
                ⁢
                
                  L
                  -
                  1
                  -
                  
                    D
                    ⁡
                    
                      (
                      
                        
                          ⌊
                          
                            
                              L
                              -
                              1
                            
                            D
                          
                          ]
                        
                        -
                        1
                      
                      )
                    
                  
                
                ]
              
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  (
                  
                    
                      g
                      i
                      *
                    
                    ⁢
                    
                      h
                      j
                    
                  
                  )
                
                ⁡
                
                  [
                  
                    L
                    -
                    1
                    +
                    
                      D
                      ⁢
                      
                        ⌊
                        
                          
                            L
                            -
                            1
                          
                          D
                        
                        ⌋
                      
                    
                  
                  ]
                
              
            
            ]
          
          .
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          22
        
        )
      
    
  




so that the inter-symbol interference canceller vector and the inter-user interference canceller vector can be represented as described below.

Considering the inter-symbol interference to Xi[k] as a linear combination of user i's own symbols, we define the inter-symbol interference canceller vector Vi,i for user i to be

Vi,iUi,iD(1,0,1),  (Equ. 23)

where D(z) is a diagonal matrix whose diagonal elements are listed by z and 1{1, 1, . . . 1} of length





  ⌊
  
    
      L
      -
      1
    
    D
  
  ⌋




As a result, the approximated inter-symbol interference term for user i's symbol Xi[k] can be written in a compact form as the product of the inter-symbol interference canceller vector Vi,i and the tentative decision vector {circumflex over (X)}i(0)[k] as shown below





  
    
      
        ISI
        =
        
          
            
              a
              i
            
            ⁢
            
              
                ∑
                
                  
                    
                      
                        1
                        =
                        
                          -
                          
                            ⌊
                            
                              
                                L
                                -
                                1
                              
                              D
                            
                            ⌋
                          
                        
                      
                    
                  
                  
                    
                      
                        l
                        ≠
                        0
                      
                    
                  
                
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    (
                    
                      
                        g
                        i
                      
                      *
                      
                        h
                        i
                      
                    
                    )
                  
                  ⁡
                  
                    [
                    
                      L
                      -
                      1
                      +
                      Dl
                    
                    ]
                  
                
                ⁢
                
                  
                    
                      X
                      ^
                    
                    j
                  
                  ⁡
                  
                    [
                    
                      k
                      -
                      l
                    
                    ]
                  
                
              
            
          
          =
          
            
              a
              i
            
            ⁢
            
              V
              
                i
                ,
                i
              
            
            ⁢
            
              
                
                  
                    X
                    ^
                  
                  i
                  
                    (
                    0
                    )
                  
                
                ⁡
                
                  [
                  k
                  ]
                
              
              .
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          24
        
        )
      
    
  




We define the inter-user interference canceller vector for the inter-user interference caused by user j to user i as Ui,j, so that the estimated inter-user interference term to be canceled for user i's symbol Xi[k] can be obtained as




  
    
      
        IUI
        =
        
          
            
              ∑
              
                
                  
                    
                      j
                      =
                      1
                    
                  
                
                
                  
                    
                      j
                      ≠
                      i
                    
                  
                
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                a
                j
              
              ⁢
              
                
                  ∑
                  
                    l
                    =
                    
                      -
                      
                        ⌊
                        
                          
                            L
                            -
                            1
                          
                          D
                        
                        ⌋
                      
                    
                  
                  
                    ⌊
                    
                      
                        L
                        -
                        1
                      
                      D
                    
                    ⌋
                  
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      (
                      
                        
                          g
                          i
                        
                        *
                        
                          h
                          j
                        
                      
                      )
                    
                    ⁡
                    
                      [
                      
                        L
                        -
                        1
                        +
                        Dl
                      
                      ]
                    
                  
                  ⁢
                  
                    
                      
                        X
                        ^
                      
                      j
                    
                    ⁡
                    
                      [
                      
                        k
                        -
                        l
                      
                      ]
                    
                  
                
              
            
          
          =
          
            
              ∑
              
                
                  
                    
                      j
                      =
                      1
                    
                  
                
                
                  
                    
                      j
                      ≠
                      i
                    
                  
                
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                a
                j
              
              ⁢
              
                U
                
                  i
                  ,
                  j
                
              
              ⁢
              
                
                  
                    
                      
                        X
                        ^
                      
                      _
                    
                    j
                    
                      (
                      0
                      )
                    
                  
                  ⁡
                  
                    [
                    k
                    ]
                  
                
                .
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          25
        
        )
      
    
  




Single-Stage 2D Interference Cancellation
Referring to FIGS. 16A to 16C, in an exemplary embodiment, a 2D interference cancelation and decision module 1670 may use a single-stage 2D parallel interference cancellation scheme. The 2D interference cancelation and decision module 1670 may include N interference cancelation and decision module units (e.g., 1672a, 1672b, 1672c, collectively 1672) each capable of at least partially canceling the interference associated with symbols sent from one of the users 1554. In this example, the interference cancelation and decision module unit 1672a may cancel the inter-symbol interference and the inter-user interference from the symbols sent from the user 1554a, the unit 1672b may cancel the interference from the symbols sent from the user 1554b, the unit 1672c may cancel the interference from the symbols sent from the user 1554c, and so forth.
FIG. 16D shows an enlarged diagram of the exemplary interference cancelation and decision module unit 1672a. The input signal Y1[k] 1674 associated with user 1 (1554a) may be buffered by a delay chain 1676 of length




  2
  ⁢
  
    
      ⌊
      
        
          L
          -
          1
        
        D
      
      ⌋
    
    .
  




Tentative hard decision units 1678 may make tentative hard decisions in parallel for the signal associated with user 1, and the tentative hard decisions may be concatenated at a concatenation unit 1680 to obtain a vector {circumflex over (X)}1(0) 1682 using the decision rule of Equation 19. The vector {circumflex over (X)}1(0) represents the tentative hard decisions for user 1.
The matrix product of vector {circumflex over (X)}1(0) 1682 and the inter-symbol interference canceller vector V11, scaled by a power coefficient a1, can be calculated to obtain the inter-symbol interference ISI1,1 (1684a) associated with the symbols sent from user 1.
The matrix product of a vector {circumflex over (X)}2(0) 1682b and the inter-user interference canceller vector U1,2, scaled by a power coefficient a2, can be calculated to obtain the inter-user interference IUI1,2 (1684b) representing the interference of user 2 (1554b) to user 1.
In a similar manner, the matrix product of a vector {circumflex over (X)}N(0) and the inter-user interference canceller vector U1,N, scaled by a power coefficient aN, can be calculated to obtain the inter-user interference IUI1,N representing the interference of user N to user 1.
An interference cancellation module 1686 may subtract the inter-symbol interference ISI1,1 (1684a) and inter-user interferences IUI1,2 (1684b) to IUI1,N (1684c) from the time-reversal mirror output Y1[k] to obtain a soft bit representing a refined estimate of the symbol sent from user 1, represented as {tilde over (X)}1(1)[k]. The soft bit may pass through a hard decision unit 1690, which may output a hard decision of the symbol sent from user 1, represented as {circumflex over (X)}1(1)[k]. The superscript (1) indicates that this is the hard decision for stage 1 interference cancelation.
In general, at each interference cancelation and decision module unit 1672i (i=1 to N), the input signal Yi[k] associated with the i-th user can be buffered by a delay chain of length




  2
  ⁢
  
    
      ⌊
      
        
          L
          -
          1
        
        D
      
      ⌋
    
    .
  




Tentative hard decision units can make tentative hard decisions in parallel for the signal associated with user i, and the tentative hard decisions can be concatenated at a concatenation unit to obtain a vector {circumflex over (X)}i(0)[k] for ∀iε{1, 2, . . . , N} using the decision rule of Equation 19. In embodiments, the vector {circumflex over (X)}i(0) represents the tentative hard decisions for user i.

The matrix product of vector {circumflex over (X)}i(0) and the inter-symbol interference canceller vector Vii, scaled by a power coefficient ai, can be calculated to obtain the inter-symbol interference ISIi associated with the symbols sent from the user i.
The matrix product of each of vectors {circumflex over (X)}j(0) (j≠i) and the inter-user interference canceller vector Ui,j, scaled by a power coefficient aj, can be calculated to obtain the inter-user interference IUIi,j representing the interference of user j to user i.
An interference cancellation module 1686 may subtract the inter-symbol interference IUIi and inter-user interferences IUIi,j (for all j not equal to i) from the time-reversal mirror output Yi[k] to obtain a soft bit representing a refined estimate of the symbol sent from user i, represented as {tilde over (X)}i(1)[k]. The soft bit may pass through a hard decision unit, which may output a refined hard decision of the symbol sent from user i, represented as {circumflex over (X)}i(1)[k]. The superscript (1) indicates that this is the hard decision for stage 1 interference cancelation.
Mathematical derivations of the final hard decision {circumflex over (X)}i(1)[k] are provided below. As described above, upon obtaining the tentative decision vectors {circumflex over (X)}i(0)[l] for all iε{1, 2, . . . , N}, the interference terms may be reconstructed and then subtracted from the signal Yi[k] with the inter-symbol interference and inter-user interference canceller vectors.
Similar to the definition in Equation 21, denoting




  
    
      
        
          X
          _
        
        j
      
      ⁡
      
        [
        k
        ]
      
    
    ⁢
    
      =
      Δ
    
    ⁢
    
      
        [
        
          
            
              X
              j
            
            ⁢
            
              ⌊
              
                k
                +
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
              ⌉
            
          
          ,
          
            
              X
              j
            
            ⁡
            
              [
              
                k
                +
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
                -
                1
              
              ]
            
          
          ,
          ...
          ⁢
          
            
          
          ,
          
            
              X
              j
            
            ⁡
            
              [
              
                k
                -
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
              ]
            
          
        
        ]
      
      T
    
  
  ,




we can rewrite Equation 16 in a more compact form as





  
    
      
        
          
            Y
            i
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            
              a
              i
            
            ⁢
            
              
                
                  ∑
                  
                    l
                    =
                    0
                  
                  
                    L
                    -
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                    
                      
                        h
                        i
                      
                      ⁡
                      
                        [
                        l
                        ]
                      
                    
                    
                  
                  2
                
              
            
            ⁢
            
              
                X
                i
              
              ⁡
              
                [
                k
                ]
              
            
          
          +
          
            
              a
              i
            
            ⁢
            
              V
              
                i
                ,
                i
              
            
            ⁢
            
              
                
                  X
                  _
                
                i
              
              ⁡
              
                [
                k
                ]
              
            
          
          +
          
            
              ∑
              
                
                  
                    
                      j
                      =
                      1
                    
                  
                
                
                  
                    
                      j
                      ≠
                      i
                    
                  
                
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                a
                j
              
              ⁢
              
                U
                
                  i
                  ,
                  j
                
              
              ⁢
              
                
                  
                    X
                    _
                  
                  j
                
                ⁡
                
                  [
                  k
                  ]
                
              
            
          
          +
          
            
              
                n
                i
              
              ⁡
              
                [
                k
                ]
              
            
            .
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          26
        
        )
      
    
  




After the interference cancellation, the resulting soft-bit {tilde over (X)}i(1)[k] can be written as





  
    
      
        
          
            
              X
              ~
            
            i
            
              (
              1
              )
            
          
          ⁡
          
            [
            k
            ]
          
        
        =
        
          
            
              Y
              i
            
            ⁡
            
              [
              k
              ]
            
          
          -
          
            
              a
              i
            
            ⁢
            
              V
              
                i
                ,
                i
              
            
            ⁢
            
              
                
                  
                    X
                    _
                  
                  ^
                
                i
                
                  (
                  0
                  )
                
              
              ⁡
              
                [
                k
                ]
              
            
          
          -
          
            
              ∑
              
                
                  
                    
                      j
                      =
                      1
                    
                  
                
                
                  
                    
                      j
                      ≠
                      i
                    
                  
                
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                a
                j
              
              ⁢
              
                U
                
                  i
                  ,
                  j
                
              
              ⁢
              
                
                  
                    
                      X
                      _
                    
                    ^
                  
                  j
                  
                    (
                    0
                    )
                  
                
                ⁡
                
                  [
                  k
                  ]
                
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          27
        
        )
      
    
  




A hard decision {circumflex over (X)}i(1)[k]=sgn({tilde over (X)}i(1)[l]) can be made based on {tilde over (X)}i(1)[k] to achieve a more refined estimation for the transmitted symbol.

The Multi-Stage Iterative Scheme
To further improve the detection performance, multi-stage processing can be performed by cascading multiple stages of the 2D parallel interference cancellation. In this part, we generalize the basic single-stage scheme into an iterative multi-stage scheme, in which the signal quality may be refined after each stage of interference cancellation. This may be useful when, for example, the bit error rate performance of the initial stage is not sufficiently high.
Referring to FIG. 17, a 2D interference cancelation and decision module 400 may use an M-stage interference cancellation scheme. As the interference cancellation procedures proceed in each stage, the same 2D parallel interference cancellation may be performed in an attempt to remove both the inter-symbol interference and the inter-user interference with the updated estimates of the transmitted symbols. Each stage may take the soft-bits from the previous stage as input, based on which tentative decisions are made to estimate the interference and generate the soft-bits of the current stage as the output.
More specifically, in the context of multiple stages, the m-th stage's operation can be described as follows:
Delay and Buffering: The soft-bits from the previous stage (the (m−1)-th stage) are delayed and buffered to form a soft-bit vector for each user, such that for all iε{1, 2, . . . , N}.




  
    
      
        
          
            
              
                X
                ~
              
              _
            
            i
            
              (
              
                m
                -
                1
              
              )
            
          
          ⁡
          
            [
            
              k
              -
              
                m
                ⁢
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
            
            ]
          
        
        ⁢
        
          =
          Δ
        
        ⁢
        
          [
          
            
              
                
                  X
                  ~
                
                i
                
                  (
                  
                    m
                    -
                    1
                  
                  )
                
              
              ⁡
              
                [
                
                  k
                  -
                  
                    
                      (
                      
                        m
                        -
                        1
                      
                      )
                    
                    ⁢
                    
                      ⌊
                      
                        
                          L
                          -
                          1
                        
                        D
                      
                      ⌋
                    
                  
                
                ]
              
            
            ,
            
              
                
                  X
                  ~
                
                i
                
                  (
                  
                    m
                    -
                    1
                  
                  )
                
              
              ⁡
              
                [
                
                  k
                  -
                  
                    
                      (
                      
                        m
                        -
                        1
                      
                      )
                    
                    ⁢
                    
                      ⌊
                      
                        
                          L
                          -
                          1
                        
                        D
                      
                      ⌋
                    
                  
                  -
                  1
                
                ]
              
            
            ⁢
            
              
            
            ,
            
               
              
                ...
                ⁢
                
                  
                
                ,
                
                   
                  
                    
                  
                  ⁢
                  
                     
                    
                      
                        
                          X
                          ~
                        
                        i
                        
                          (
                          
                            m
                            -
                            1
                          
                          )
                        
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                         
                        
                           
                          
                            [
                            
                              k
                              -
                              
                                
                                  (
                                  
                                    m
                                    +
                                    1
                                  
                                  )
                                
                                ⁢
                                
                                  ⌊
                                  
                                    
                                      L
                                      -
                                      1
                                    
                                    D
                                  
                                  ⌋
                                
                              
                            
                            ]
                          
                          ]
                        
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          28
        
        )
      
    
  




Tentative Decision: The tentative decisions are made based on the soft-bits from the previous stage (the (m−1)-th stage) in an attempt to estimate the transmitted symbols, i.e., for all iε{1, 2, . . . , N},




  
    
      
        
          
            
              
                X
                ^
              
              _
            
            i
            
              (
              
                m
                -
                1
              
              )
            
          
          ⁡
          
            [
            
              k
              -
              
                m
                ⁢
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
            
            ]
          
        
        =
        
          sgn
          ⁡
          
            (
            
              
                
                  
                    X
                    ^
                  
                  _
                
                i
                
                  (
                  
                    m
                    -
                    1
                  
                  )
                
              
              ⁢
              
                
                
                  k
                  -
                  
                    m
                    ⁢
                    
                      ⌊
                      
                        
                          L
                          -
                          1
                        
                        D
                      
                      ⌋
                    
                  
                
                
              
            
            )
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          29
        
        )
      
    
  




where sgn(•) is the sign function applied element-wise to the operand when the operand is a vector.

Interference Cancellation: The soft-bits of the current stage (the m-th stage) are generated by subtracting the estimated interference terms from the original signals {Yi}. The soft-bit of user i generated by the m-th stage is given by




  
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          30
        
        )
      
    
    
      
        
      
    
  
  
    
      
        
          
            
              X
              ~
            
            i
            
              (
              m
              )
            
          
          [
          
            k
            -
            
              m
              ⁢
              
                ⌊
                
                  
                    L
                    -
                    1
                  
                  D
                
                ⌋
              
            
          
          ⌋
        
        =
        
          
            
              Y
              i
            
            [
            
              k
              -
              
                m
                ⁢
                
                  ⌊
                  
                    
                      L
                      -
                      1
                    
                    D
                  
                  ⌋
                
              
            
            ⌉
          
          -
          
            
              a
              i
            
            ⁢
            
              V
              
                i
                ,
                i
              
            
            ⁢
            
              
                
                  
                    X
                    _
                  
                  ^
                
                i
                
                  (
                  
                    m
                    -
                    1
                  
                  )
                
              
              ⁡
              
                [
                
                  k
                  -
                  
                    m
                    ⁢
                    
                      ⌊
                      
                        
                          L
                          -
                          1
                        
                        D
                      
                      ⌋
                    
                  
                
                ]
              
            
          
          -
          
            
              ∑
              
                
                  
                    
                      j
                      =
                      1
                    
                  
                
                
                  
                    
                      j
                      ≠
                      i
                    
                  
                
              
              N
            
            ⁢
            
              
            
            ⁢
            
              
                a
                j
              
              ⁢
              
                U
                
                  i
                  ,
                  j
                
              
              ⁢
              
                
                  
                    
                      X
                      _
                    
                    ^
                  
                  j
                  
                    (
                    
                      m
                      -
                      1
                    
                    )
                  
                
                [
                
                  k
                  -
                  
                    m
                    ⁢
                    
                      ⌊
                      
                        
                          L
                          -
                          1
                        
                        D
                      
                      ⌋
                    
                  
                
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




In particular, such an M-stage scheme is initialized by setting the soft-bits of the initial stage (the 0-th stage) as

{tilde over (X)}1(0)[k]=Y1[k],{tilde over (X)}2(0)[k]=Y2[l] . . . ,{tilde over (X)}N(0)[k]=YN[k].  (Equ. 31)

After the signal flows through all the M stages, the final decision is made for each user based on the output of the stage-M interference cancellation,




  
    
      
        
          
            
              
                X
                ^
              
              i
              
                (
                M
                )
              
            
            ⁡
            
              [
              
                k
                -
                
                  M
                  ⁢
                  
                    ⌊
                    
                      
                        L
                        -
                        1
                      
                      D
                    
                    ⌋
                  
                
              
              ]
            
          
          =
          
            sgn
            ⁡
            
              (
              
                
                  
                    X
                    ~
                  
                  i
                  
                    (
                    M
                    )
                  
                
                ⁡
                
                  [
                  
                    k
                    -
                    
                      M
                      ⁢
                      
                        ⌊
                        
                          
                            L
                            -
                            1
                          
                          D
                        
                        ⌋
                      
                    
                  
                  ]
                
              
              )
            
          
        
        ,
        
          

        
        ⁢
        
          
            for
            ⁢
            
              
            
            ⁢
            all
            ⁢
            
              
            
            ⁢
            i
          
          ∈
          
            {
            
              1
              ,
              2
              ,
              ...
              ⁢
              
                
              
              ,
              N
            
            }
          
        
      
    
    
      
        (
        
          Equ
          .
          
            
          
          ⁢
          32
        
        )
      
    
  




Comparing the timing of the estimated symbols after the final decision




  
    
      X
      ^
    
    i
    
      (
      M
      )
    
  
  ⁡
  
    [
    
      k
      -
      
        M
        ⁢
        
          ⌊
          
            
              L
              -
              1
            
            D
          
          ⌋
        
      
    
    ]
  




and the original time-reversal mirror output signal Yi[k], one can see that each stage incurs a delay of





  
    ⌊
    
      
        L
        -
        1
      
      D
    
    ⌋
  
  .




Therefore, the total detection delay grows linearly with the number of stages, but is independent of the number of users. The complexity of the 2D parallel interference cancellation increases linearly with the number of users (or devices), as opposed to previous joint detection schemes whose complexity grows exponentially with the number of users.

Simulation Results
The following describes simulation results on the bit error rate performance of the 2D parallel interference cancellation scheme. To study the interference cancellation scheme in a realistic setting, we used the practical IEEE 802.15.4a outdoor non-line-of-sight (NLOS) channel model to evaluate the bit error rate performance of the time-reversal system shown in FIGS. 15-17. In such a channel model, each channel tap is a real number.
Referring to FIG. 18A, a graph 420 shows an example of a channel impulse response for user 1 under the IEEE 802.15.4a outdoor non-line-of-sight channel model. Referring to FIG. 18B, a graph 422 shows an example of a channel impulse response for user 2 under the IEEE 802.15.4a outdoor non-line-of-sight channel model.
Referring to FIG. 19A, a graph 424 shows the convolution of the channel impulse response for user 1 with the time-reversal signature waveform for user 1.
Referring to FIG. 19B, a graph 426 shows the convolution of the channel impulse response for user 2 with the time-reversal signature waveform for user 1.
The channels used for the simulations for the examples shown in FIGS. 18A and 18B are randomly generated according to the channel model specified in IEEE 802.15.4a, with the system sampling period Ts=1 ns and the channel length truncated (because the amplitude of the remaining paths after 300 ns is typically small enough to be neglected) at LTs=300 ns (i.e., L=300).
In FIG. 19A, the convolution between user 1's channel impulse response h1 and its matched time-reversal signature waveform g1 exhibits a prominent central peak 428 at (h1*g1)[L−1], demonstrating the temporal focusing effect of the time-reversal technique.
In FIG. 19B, the amplitude of the convolution between the time-reversal signature waveform g1 and the mismatched channel impulse response h2 is significantly smaller than the central peak (h1*g1)[L−1] shown in FIG. 19A, demonstrating the virtual spatial focusing effect in the user's signature domain.
The following provides numerical evaluations on the predicted bit error rate performance of the 2D interference cancellation scheme. In the simulations, the channel impulse responses for various users are randomly and independently generated using the IEEE 8.2.15.4a channel model, with Ts=1 ns and L=300. Without loss of generality, the channel impulse response of each user is normalized so that




  
    
      
        ∑
        
          k
          =
          0
        
        
          L
          -
          1
        
      
      ⁢
      
        
      
      ⁢
      
        
          
          
            
              h
              i
            
            ⁡
            
              [
              k
              ]
            
          
          
        
        2
      
    
    =
    1
  
  ,
  
    ∀
    
      i
      ∈
      
        {
        
          1
          ,
          2
          ,
          ...
          ⁢
          
            
          
          ,
          N
        
        }
      
    
  
  ,




and we assume that all the power control coefficients ai=1, ∀iε{1, 2, . . . , N}. This provides equal power allocation among the users.

A large number of independent trials of channel realizations were conducted and averaged to characterize the average performance of the interference cancellation scheme under this channel model.
Bit Error Rate Versus Eb/N0 
The following describes the predicted bit error rate performance versus Eb/N0 (the energy-per-bit to noise-power-spectral-density ratio) with various combinations of rate back-off factor D and the total number of users N. The energy-per-bit Eb is normalized to 1 by the assumption that each BPSK symbol Xi[k]ε{−1,+1} has a unit power. Accordingly, the power of the received additive white Gaussian noise (AWGN) ñ[k] at the base station is given by




  
    σ
    N
    2
  
  =
  
    
      
        N
        0
      
      2
    
    =
    
      
        
          (
          
            
              2
              ⁢
              
                E
                b
              
            
            
              N
              0
            
          
          )
        
        
          -
          1
        
      
      .
    
  




Referring to FIG. 20, a graph 430 shows the predicted bit error rate performances of the 2D parallel interference cancellation scheme (up to 3 stages), compared with the basic time-reversal division multiple access system without interference cancellation. In the example of FIG. 20, the simulation considers that case where there are N=5 end-users accessing the base station at the same time with a rate back-off factor D=16 (about 5.3% of the channel length L=300).
Referring to FIG. 21, a graph 440 shows the predicted bit error rate performances of the 2D parallel interference cancellation scheme (up to 3 stages), compared with the basic time-reversal division multiple access system without interference cancellation. In the example of FIG. 21, the simulation considers that case where there N=10 end-users accessing the base station at the same time with a rate back-off factor D=32 (about 10.7% of the channel length).
FIGS. 20 and 21 show that significant bit error rate performance gain can be achieved by the 2D parallel interference cancellation scheme, compared with the baseline time-reversal division multiple access system without interference cancellation. Additional gain can be achieved by cascading more stages of the interference cancellation scheme, at the price of increased decoding delay that grows linearly with the total number of stages. The largest gain-per-stage is obtained by the first stage, and the marginal gain diminishes for the following stages. There is a tradeoff between the system performance and the number of stages, and different tradeoffs can be considered for different applications. The theoretical approximation derived above matches reasonably well with the simulation results, considering its complicated correlation between tentative decisions and the Gaussian assumptions that were made in the approximation. A comparison of FIGS. 20 and 21 shows that the accuracy of the approximation improves as the number of users increases, which agrees with the central limit theorem.
High Signal-to-Noise Ratio Regime
Referring to FIG. 22, a graph 450 shows how the predicted bit error rate performance changes with the number of users N with a high signal-to-noise ratio. In this example, the rate back-off factor D=32. The graph 450 shows that the predicted bit error rate increases with the number of users N due to the increased inter-user interference. For a given bit error rate level (e.g., 10−3), the 2D interference cancellation may enable more users to transmit simultaneously and therefore may increase the system capacity. The benefit of using the interference cancellation may diminish when the bit error rate of the baseline system (considered as the initial stage) is above a certain threshold. This is because the interference cancellation scheme at the receiver eventually relies on the tentative decisions of the detected symbols to cancel the interference, whose effectiveness depends on the quality of those tentative decisions.
Referring to FIG. 23, a graph 460 shows how the predicted bit error rate performance is affected by the rate back-off factor D with high signal-to-noise ratio, assuming that there are 5 users accessing the base station simultaneously. Because both the inter-symbol interference and the inter-user interference are reduced with a larger rate back-off factor D, the predicted bit error rate decreases as D increases. The graph 460 shows that by using the 2D parallel interference cancellation scheme, the time-reversal division multiple access system can use a smaller D to achieve the same bit error rate, which translates to higher throughput for each user. When D is small (e.g., for D=6 (only 2% of the channel length)), there can be more erroneous tentative decisions (i.e., high bit error rate (greater than 10−9 for the initial stage), which in turn enhances the interference power, rather than canceling the interference. However, many time-reversal division multiple access applications will operate at a much lower bit error rate level, so the interference cancellation scheme is not expected to contribute to additional interference. Note that a lower bit error rate level can be achieved by increasing the rate back-off factor D and/or reducing the number of concurrent users N, thereby increasing the accuracy of the tentative decisions at the initial stage.
Referring to FIG. 24, an exemplary process 470 for time-reversal wireless communication with interference cancellation is provided. For example, the process 470 can be implemented by the base station 1552 of FIG. 15. The process 470 may include at a transceiver, receiving a combined signal that includes a first signal from a first device and a second signal from a second device, the first signal sent from the first device to the transceiver through first multiple wireless propagation paths, the second signal sent from the second device to the transceiver through second multiple wireless propagation paths, each of the first and second signals including a plurality of symbols (472). For example, the transceiver can be part of the base station 1552, the first device can be the first terminal device 1554a, and the second device can be the second terminal device 1554b. 
The process 470 may include generating a first set of estimated symbols based on the combined signal and a first signature waveform associated with the first device, the first set of estimated symbols representing a rough estimate of the symbols sent from the first device (474). For example, the base station 1552 can apply a signature waveform associated with the first terminal device 1554a to a combined signal to generate a first set of estimated symbols, e.g., according to Equation 13.
The process 470 may include generating a second set of estimated symbols based on the combined signal and a second signature waveform associated with the second device, the second set of estimated symbols representing a rough estimate of the symbols sent from the second device (476). For example, the base station 1552 can apply a signature waveform associated with the second terminal device 1554b to the combined signal to generate a second set of estimated symbols, e.g., according to Equation 13.
The process 470 may include determining interference based on at least one of the first set of estimated symbols or the second set of estimated symbols (478). For example, the base station 1552 can determine the interference, which can include inter-symbol interference and/or inter-user interference. For the first set of estimated symbols, the inter-symbol interference can be determined based on the first set of estimated symbols, and the inter-user interference can be determined based on the second set of estimated symbols. For the second set of estimated symbols, the inter-symbol interference can be determined based on the second set of estimated symbols, and the inter-user interference can be determined based on the first set of estimated symbols. For example, the interference can be determined according to Equations 24 and 25.
The process 470 may include generating a refined first set of estimated symbols by canceling the interference from the first set of estimated symbols, the refined first set of estimated symbols representing a refined estimate of the symbols sent from the first device (480). For example, the base station 1552 can generate the refined first set of estimated symbols by subtracting the inter-symbol interference and/or the inter-user interference from the first set of estimated symbols, e.g., according to Equation 27.
Referring to FIG. 25, an exemplary process 490 for time-reversal wireless communication with interference cancellation is provided. For example, the process 490 can be implemented by the base station 1552 of FIG. 15. The process 490 may include at a transceiver, receiving a combined signal that includes a first signal from a first device and a second signal from a second device, the first signal sent from the first device to the transceiver through first multiple wireless propagation paths, the second signal sent from the second device to the transceiver through second multiple wireless propagation paths (492). For example, the transceiver can be part of the base station 1552, the first device can be the first terminal device 1554a, and the second device can be the second terminal device 1554b. 
The process 490 may include estimating the first signal sent from the first device to generate a first estimated signal based on the combined signal and a first signature waveform associated with the first device (494). For example, the base station 1552 can apply a signature waveform associated with the first terminal device 1554a to a combined signal to generate a first estimated signal.
The process 490 may include estimating the second signal sent from the second device to generate a second estimated signal based on the combined signal and a second signature waveform associated with the second device (496). For example, the base station 1552 can apply a signature waveform associated with the second terminal device 1554b to the combined signal to generate a second estimated signal.
The process 490 may include determining interference based on at least one of the first estimated signal or the second estimated signal (498). For example, the base station 1552 can determine the interference, which can include interference between signals from the same device and/or interference between signals from different devices. For the first estimated signal, the interference between signals from the same device can be determined based on the first estimated signal, and the interference between signals from different devices can be determined based on the second estimated signal. For the second estimated signal, the interference between signals from the same device can be determined based on the second estimated signal, and the interference between signals from different devices can be determined based on the first estimated signal.
The process 490 may include determining the signal sent from the first device by subtracting the interference from the first estimated signals (500). For example, the base station 1552 can determine the signal sent from the first device by subtracting the interference between signals from the same device and/or the interference between signals from different devices from the first estimated signal.
Referring to FIG. 26, an exemplary process 510 for time-reversal wireless communication with interference cancellation is provided. For example, the process 510 can be implemented by the base station 1552 of FIG. 15. The process 510 may include at a transceiver, receiving a combined signal that includes signals from a plurality of devices, each device sending a signal to the transceiver through multiple wireless propagation paths (512). For example, the transceiver can be part of the base station 1552, and the plurality of devices can include the first terminal device 1554a and the second terminal device 1554b. 
The process 510 may include estimating signals sent from each of the devices based on the combined signal and a signature waveform associated with the device (514). For example, for each of the terminal devices 1554, the base station 1552 can apply a signature waveform associated with the terminal device 1554 to a combined signal to generate estimated signals for the terminal device 1554.
The process 510 may include determining an interference associated with the estimated signal from each of the devices based on the estimated signals from the devices (516). For example, the base station 1552 can determine the interference, which can include interference between signals from the same device and/or interference between signals from different devices. For each terminal device 1554, the interference between signals from the same device can be determined based on the estimated signal for the terminal device 1554, and the interference between signals from different devices can be determined based on the estimated signals for the other devices.
The process 510 may include determining the signal sent from each of the devices by subtracting the interference associated with the estimated signal associated with the device from the estimated signal associated with the device (518). For example, for each of the terminal devices 1554, the base station 1552 can determine the signal sent from the device by subtracting the interference between signals from the same device and/or the interference between signals from different devices from the estimated signal for the terminal device.
A multi-user time-reversal division multiple access uplink architecture and a 2D parallel interference cancellation scheme are provided to enhance the system performance. The TRDMA uplink architecture described above may keep the cost of the communication components of the end-users at a low level, and may reuse the processing power at the base station that has already been made available for the downlink. The 2D parallel interference cancellation scheme may utilize the tentative decisions of detected symbols to effectively reduce or cancel the interference in both the time dimension (ISI) and the user dimension (IUI). To further improve the bit error rate performance, a multi-stage scheme may be provided by cascading multiple stages of the 2D interference cancellation, with a total delay that increases linearly with the number of stages, but independent of the number of users.
In this disclosure, an asymmetric system architecture for the time reversal wireless communication system 150 has been described. The system 100 can have one or more of the following features and advantages:
The system 150 may use a single-carrier wireless broadband communication technology using the time reversal structure and other modified waveforms, which effectively suppress the inter-symbol interference.
The system 150 may leverage the time-reversal division multiple access concept to explore the spatial degrees of freedom, and thus enable multiple simultaneous transmissions on the same frequency band.
The system architecture described above can accommodate all types of waveform design (as a natural generalization of the time-reversed channel response), and arbitrary channel coding scheme (including the uncoded scheme).
When the time-reversal waveform is used as the signature waveform, the downlink scheme described above can form the spatial focusing effect of time reversal, focusing the signal power only at the locations of the intended terminal devices. Such a physical spatial focusing effect can reduce the co-channel interference and electromagnetic pollution to the environment, and may enhance the privacy and security of the transmitted message.
The asymmetric complexity-distribution between the base station (e.g., 152) (which has the most complexity) and the terminal devices (e.g., 154) (which can be of lower complexity) may make the terminal devices more affordable, which is a desirable feature for systems having a large number of terminal devices that may be limited by size, power, and/or computational capability.
The architecture design described above does not require channel equalization at the terminal devices to enable two-way communication.
In some implementations, the base station 152 can be part of a mobile or stationary device. For example, the base station 152 can be implemented as part of a sensor module, a controller, a mobile phone, a laptop computer, or an electronic appliance that communicates wirelessly with multiple other devices. For example, a mobile phone or a laptop computer may communicate simultaneously with a television, a printer, a thermometer, a radio, a refrigerator, a lighting control system, and other devices using the techniques described above. The base station may be an access point or a router and may be associated with a wireless hot spot, or wireless cell, or picocell, or nanocell and the like.
Systems Using Multiple Communication Protocols
In some implementations, the base station and the terminal devices can communicate with one another other using multiple communication protocols, including communication protocols that use time-reversal techniques and communication protocols that do not use time-reversal techniques. For example, a mobile phone may have a cellular module configured to communicate with a wireless telephone base station using cellular communication protocols, such as Global System for Mobile Communications (GSM), Code Division Multiple Access (CDMA), Long Term Evolution (LTE), or LTE Advanced protocols. The mobile phone may have a Wi-Fi module configured to communicate with a network router using Wi-Fi protocols (e.g., IEEE 802.11, 802.11a, 802.11b, 802.11g, 802.11n, 802.11ac, 802.11ad). The mobile phone may have a Bluetooth module configured to communicate with Bluetooth enabled devices. The mobile phone may have a time-reversal module configured to communicate with sensors and other devices using the time-reversal communication protocol described above.
Having the time-reversal module in combination with other communication modules may enhance the functionality of the mobile phone. For example, the mobile phone may route calls through the wireless telephone base station when the mobile phone is in an outdoor location, and switch to routing calls through time-reversal base stations when the mobile phone is in an indoor venue where multipath interference is high. For example, if a network router has both Wi-Fi and time-reversal modules, the mobile phone may communicate with the network router using Wi-Fi protocols when executing an application that requires a higher data rate (e.g., when sending image and video data), and switch to using the time-reversal protocol when executing an application that requires a lower data rate (e.g., when sending text or numerical data).
Referring to FIG. 27, in some exemplary implementations, a communication system 520 may include a base station 522 and multiple terminal devices (e.g., 524a, 524b, collectively 524). The base station 522 may include a first communication module 526 that uses the time-reversal communication techniques described above, and a second communication module 528 that uses the Wi-Fi, Bluetooth, or cellular communication protocol.
In the example of FIG. 27, the first terminal device 524a may have a first communication module 530a that uses time-reversal communication techniques, and a second communication module 532a that uses the Wi-Fi, Bluetooth, or cellular communication protocol. The second terminal device 524b may have a first communication module 530b that uses time-reversal communication techniques and a second communication module 532b that uses the Wi-Fi, Bluetooth, or cellular communication protocol. In the following description, depending on what communication protocol is used, the communication modules 528 and 532 may be referred to as a Wi-Fi communication module, a Bluetooth communication module, or a cellular communication module. In some examples, each of the base station 522 and terminal devices 524 may have two or more modules that comply with two or more of Wi-Fi, Bluetooth, cellular, and other communication protocols.
The time-reversal communication module 526 at the base station 522 may communicate with the time-reversal communication modules 530a, 530b at the first and second terminal devices 524a, 524b, respectively. The Wi-Fi (or Bluetooth, cellular) communication module 528 at the base station 522 may communicates with the Wi-Fi (or Bluetooth, cellular) communication modules 532a, 532b at the first and second terminal devices 524a, 524b, respectively.
In some implementations, the base station 522 and the terminal device 524a may perform a handshake or channel probing process using the Wi-Fi, Bluetooth, or cellular protocol in order to obtain channel information, then use time-reversal techniques to transmit data. In this example, the time-reversal communication module 526 can continuously transmit downlink data packets or receive uplink data packets to or from the time-reversal communication module 530a without the need to periodically perform hand-shake processes to obtain channel information. Another advantage of the system 520 is that existing Wi-Fi and/or cellular communication hardware can be used for channel probing for time-reversal communications.
There are multiple channels in Wi-Fi and cellular protocols, in which each channel can be, e.g., 5 MHz for cellular protocol and 20/40 MHz for the Wi-Fi protocol. The channel information for each channel can be obtained according to the respective cellular or Wi-Fi protocol. By combining the channel information of multiple channels, we can obtain the channel information for a wider band and then generate the corresponding signature waveform for time-reversal communications.
In some implementations, both communication modules 526, 528 may be used to transmit data. For example, the base station 522 may initially use the Wi-Fi (or Bluetooth, cellular) communication module 528 to establish a Wi-Fi (or Bluetooth, cellular) communication link with the terminal device 524a, then based on certain criteria, switch to using the time-reversal communication module 526 to establish a communication link with the terminal device 524a. 
The criteria for switching may be based on one or more factors such as interference level of the communication channel, battery level of the terminal device, location of the terminal device, and the type of data being transmitted. For example, the Wi-Fi communication modules 528, 532 may be used when the interference level is below a threshold, and the time-reversal communication modules 526, 530 may be used when the interference level is above the threshold. The Wi-Fi communication modules 528, 532 may be used when the application requires a higher data rate (e.g., when transmitting image or video data), and the time-reversal communication modules 526, 530 may be used when the application requires a lower data rate (e.g., when transmitting text or numerical data). The Wi-Fi communication modules 528, 532 may be used when the terminal device 524 has sufficient battery power, and the time-reversal communication modules 526, 530 may be used when the battery level of the terminal device 524 is low. The time-reversal communication module 530 may consume less power than the Wi-Fi communication module 532, so switching to using the time-reversal communication module 530 may reduce power consumption and allow the terminal device 524 to operate for a longer period of time using battery power.
The base station 522 and the terminal device 524 may perform a handoff procedure to terminate the Wi-Fi (or Bluetooth, cellular) communication session between the modules 526, 530, and initiate a time-reversal communication session between the modules 528, 532. Prior to handoff, the base station 522 may send a command signal to the terminal device 524 to initiate the handoff procedure. The time-reversal communication modules 530 and 526 may perform a hand-shake process to obtain channel and timing information. The switching from using the Wi-Fi (or Bluetooth, cellular) protocol to using the time-reversal communication protocol, or vice versa, can be initiated by either the base station 522 or the terminal device 524.
In some implementations, prior to switching from the Wi-Fi protocol to the time-reversal communication protocol, the base station 522 and the terminal device 524 may test the signal quality of the link between the modules 526 and 530, and only switch to using the modules 526 and 530 if the time-reversal communication link meets predetermined quality criteria.
In some implementations, after switching to using the communication modules 526, 530, the base station 522 and the terminal device 524 may switch back to using the communication modules 528, 532 if certain criteria are met. For example, if switching from Wi-Fi to time-reversal is due to interference the base station 522 and the terminal device 524 may periodically test the communication link between the Wi-Fi communication modules 528 and 532, and if the signal quality satisfies predetermined criteria, the base station 522 and the terminal device 524 perform another handoff procedure to switch to using the Wi-Fi communication modules 528 and 532.
In some implementations, the base station 522 may communicate with a first group of terminal devices 524 using the Wi-Fi (or Bluetooth, cellular) protocol, and communicate with a second group of terminal devices 524 using the time-reversal communication protocol.
In FIG. 27, the communication modules 526 and 528 are shown as two separate blocks. In some examples, the modules 526 and 528 can have overlapping components, such as sharing one or more common data processors, controllers, and/or memory components.
The Wi-Fi, Bluetooth, and cellular communication protocols each has a substantially symmetric architecture, whereas the time-reversal communication protocol may have an asymmetric architecture. An asymmetric architecture of time-reversal communication can be seen from the examples shown in FIGS. 6 and 7, in which the complex signal processing (e.g., performing convolution computations) is performed at the base station for both uplink and downlink.
Referring to FIGS. 28A and 28B, for example, a device that is configured to communicate according to the Wi-Fi protocol may include both a transmitter 540 as shown in FIG. 28A, and a receiver 550 as shown in FIG. 28B. For example, the Wi-Fi communication module 528 of the base station 522 in FIG. 27 may include both the transmitter 540 and the receiver 550. Likewise, the Wi-Fi communication module 532 of the terminal device 524 may include both the transmitter 540 and the receiver 550. For the downlink, the transmitter 540 at the base station 522 may send data to the receiver 550 at the terminal device 524. For the uplink, the transmitter 540 at the terminal device 524 may send data to the receiver 550 at the base station 522.
The complexity of signal processing at the transmitter 540 may be comparable to that of the signal processing at the receiver 550, except that the operations are substantially reversed. For example, at the transmitter 540, channel coding (e.g., convolution code) may be applied to binary input data, and mapped according to a specific modulation (such as QPSK, 16QAM, or 64QAM). Pilot symbols may be inserted into the modulated data, and the inserted data may undergo serial-to-parallel (S/P) conversion, forming a vector of symbols. An inverse Fourier transform (IFFT) can be performed on the vector of symbols to generate a time-domain vector. After parallel-to-serial conversion (P/S), cyclic prefixes can be added to reduce interference. The output signal can then be modulated by a radio frequency carrier and transmitted through an antenna.
At the receiver 550, the opposite or reverse set of operations can be performed. A radio frequency signal may be received at the antenna and processed (e.g., including A/D conversion and low-pass filtering), and cyclic prefixes may be removed from the received signal. The signal may go through a serial-to-parallel (S/P) conversion to form a vector. A fast Fourier transform (FFT) may be performed on the vector, generating a frequency-domain vector. After parallel-to-serial conversion (P/S), the pilot may be deleted. The signal may then go through demodulation and channel decoding to estimate the transmitted binary data.
FIGS. 28A and 28B show that the signal processing at the transmitter 540 and the receiver 550 may be substantially symmetrical, thus the Wi-Fi protocol may be said to have a substantially symmetrical architecture. Moreover, the signal processing performed by the transmitter 540 at the base station 522 during downlink may be substantially the same as the signal processing performed by the transmitter 540 at the terminal device 524 during uplink. Likewise, the signal processing performed by the receiver 550 at the terminal device 524 during downlink may be substantially the same as the signal processing performed by the receiver 550 at the base station 522 during uplink.
Note that various modifications can be made to the transmitter 540 and the receiver 550. The transmitter 540 (or receiver 550) at the base station 522 may be implemented differently from the transmitter 540 (or receiver 550) at the terminal device 524. For example, the transmitter 540 at the base station 522 may use a data processor that is more powerful than the data processor used by the transmitter 540 at the terminal device 524. When we say that the Wi-Fi protocol has a substantially symmetrical architecture, we mean that the signal processing specified by the Wi-Fi protocol for the base station is comparable to that for the terminal device, even though the base station may have circuitry (including, e.g., data processor, memory, active/passive electronic components), software, and/or firmware different from those of the terminal device.
FIGS. 29A and 29B show the substantially symmetric architecture of a 4G LTE downlink. FIGS. 30A and 30B show the substantially symmetric architecture of a 4G LTE uplink. For example, an LTE communication module at the base station may include a transmitter 560 as shown in FIG. 29A and a receiver 590 as shown in FIG. 30B. A LTE communication module at the terminal device may include a transmitter 580 as shown in FIG. 30A and a receiver 570 as shown in FIG. 29B. In the downlink, the transmitter 560 at the base station may send data to the receiver 570 at the terminal device. In the uplink, the transmitter 580 at the terminal device may send data to the receiver 590 at the base station.
Referring to FIG. 29A, in the 4G LTE downlink, at the transmitter 560 at the base station, a sequence of bits to be transmitted to a terminal device may undergo a serial-to-parallel (S/P) conversion, forming a vector of symbols. The vector of symbols may be re-organized through a sub-carrier mapping. An inverse Fourier transform (IFFT) can be performed on the re-organized vector to generate a time-domain vector. After parallel-to-serial conversion (P/S), cyclic prefixes can be added to reduce interference. The output signal can then be modulated by a radio frequency carrier and transmitted through an antenna.
Referring to FIG. 29B, at the receiver at the terminal device, the opposite set of operations may be performed. A radio frequency signal may be received at an antenna and processed (e.g., including A/D conversion and low-pass filtering), and cyclic prefixes may be removed from the received signal. The signal may go through a serial-to-parallel (S/P) conversion to form a vector. A fast Fourier transform (FFT) can be applied to the vector, generating a frequency-domain vector. Then, sub-carrier demapping/equalization can be applied to extract the information for the terminal device. After a parallel-to-serial conversion (P/S), the transmitted data can be obtained at the terminal device.
Referring to FIGS. 30A and 30B, in the 4G LTE uplink, the signal processing procedures may be similar to the downlink except that at the transmitter 580 of the terminal device, a fast Fourier transform (FFT) may be applied before the subcarrier mapping, and at the receiver 590 of the base station, an inverse fast Fourier transform (IFFT) may be performed after the subcarrier demapping.
FIGS. 29A and 29B show that in the 4G LTE downlink, the signal processing at the transmitter 560 and the receiver 570 may be substantially symmetrical. FIGS. 30A and 30B show that the signal processing at the transmitter 580 and the receiver 590 may be substantially symmetrical. For both downlink and uplink, the complexity of signal processing at the terminal device may be comparable to that of the signal processing at the base station. Thus, the 4G LTE protocol may be said to have a substantially symmetrical architecture.
The above example shows the substantially symmetrical architecture of a communication system based on 4G LTE communication protocol. Other cellular communication protocols, such as GSM or CDMA, also may have substantially symmetrical architectures in which for the downlink and uplink, the complexity of signal processing at the base station is similar to that at the terminal device.
Note that various modifications can be made to the transmitter and the receiver and the base station and the terminal device, and the transmitter and receiver may be implemented using different electronic components. When we say that the 4G LTE protocol or other cellular communication protocol has a substantially symmetrical architecture, we mean that the signal processing specified by the LTE or other cellular communication protocol for the base station is comparable to that for the terminal device, even though the base station may have circuitry (including, e.g., data processor, memory, active/passive electronic components), software, and/or firmware different from those of the terminal device.
System Using Symmetric and Asymmetric Time-Reversal Communication
The following describes a time-reversal communication system that uses two types of time-reversal communication techniques, one symmetrical and the other asymmetrical.
FIG. 31A shows an exemplary downlink architecture 600 of an asymmetrical time-reversal communication system, and FIG. 31B shows an exemplary uplink architecture 610 of the asymmetrical time-reversal communication system. The downlink architecture 600 is a simplified version of that shown in FIG. 6, and the uplink architecture 610 is a simplified version of that shown in FIG. 7. In the asymmetrical time-reversal communication system, a base station may include a transmitter 602 of FIG. 31A and a receiver 612 of FIG. 31B, and a terminal device may include a receiver 604 of FIG. 31A and a transmitter 614 of FIG. 31B.
For the downlink, the transmitter 602 at the base station may transmit data to the receiver 604 at the terminal device. For the uplink, the transmitter 614 at the terminal device may transmit data to the receiver 612 at the base station. For the downlink, the transmitter 602 at the base station may perform waveform embedding, in which up-sampled sequences are embedded with a signature waveform (calculated based on the channel response obtained in a hand-shaking process). The receiver 604 at the terminal device may not need to perform signal processing that involves the signature waveform. For the uplink, the receiver 612 at the base station may perform information extraction, in which a convolution process may be performed between an input signal and the signature waveform. The transmitter 614 at the terminal device may not need to perform any signal processing that involves the signature waveform. Thus, the signal processing may be more complicated at the base station than at the terminal device, for both downlink and uplink.
FIG. 32 shows an exemplary symmetrical time-reversal communication system 620, in which a base station includes both a transmitter 602 and a receiver 604. The transmitter 602 and receiver 604 in FIG. 32 are the same as those in FIG. 31A. A terminal device may also include both a transmitter 602 and a receiver 604. For the downlink, the transmitter 602 of the base station may send data to the receiver 604 of the terminal device. For the uplink, the transmitter 602 of the terminal device may send data to the receiver 604 of the base station.
In the symmetric time-reversal downlink, the terminal device may first send a channel probing signal to the base station to allow the base station to obtain channel information and calculate a signature waveform for the terminal device. In the symmetric time-reversal uplink, the base station may first send a channel probing signal to the terminal device to allow the terminal device to obtain channel information and calculate a signature waveform for the base station. Thus, in the symmetrical time-reversal communication system, the signal processing performed by the base station and the terminal device may be similar.
In a symmetrical time-reversal communication system, when there are multiple terminal devices, the base station can use time division multiplexing to communicate with each of the terminal devices. Alternatively, the base station can use multiple antennas to receive focused signals from the multiple terminal devices.
Referring to FIG. 33, in some implementations, a time-reversal communication system 630 may include a base station 632 that is configured to communicate with terminal devices 634 using a symmetrical time-reversal architecture and an asymmetrical time-reversal architecture.
FIG. 34 shows an exemplary implementation of the architecture shown in FIG. 33. A base station 632 may include a transmitter 602 (FIG. 31A), a receiver 604, and a receiver 612 (FIG. 31B). A terminal device 634 may include a transmitter 602, a transmitter 614 (FIG. 31B), and a receiver 604. When operating in the asymmetrical time-reversal communication mode, for the downlink, the transmitter 602 of the base station 632 may send data to the receiver 604 of the terminal device 634. For the uplink, the transmitter 614 of the terminal device 634 may send data to the receiver 612 of the base station 632. When operating in the symmetric time-reversal communication mode, for the downlink, the transmitter 602 of the base station 632 may send data to the receiver 604 of the terminal device 634. For the uplink, the transmitter 602 of the terminal device 634 may send data to the receiver 604 of the base station 632.
The difference between symmetric time-reversal communication and asymmetric time-reversal communication may be in the uplink. In the asymmetric time-reversal communication, for the uplink, the terminal device may directly transmit the data, and the base station or access point may need to perform convolution to have virtual focusing and to extract uplink data. One advantage of the asymmetric time-reversal communication is that the complexity at the terminal device can be low.
By comparison, for the symmetric time-reversal communication, in the uplink, the terminal device may need to perform convolution of transmit data with a signature waveform before transmission, which may enable physical focusing of the transmit signal at the base station. Another device that is at a distance from the base station may not be able to receive the uplink data. Therefore, depending on different purposes, the terminal device may switch between different modes (i.e., asymmetric or symmetric time-reversal communication).
Each of the base station and terminal devices can include one or more processors and one or more computer-readable mediums (e.g., RAM, ROM, SDRAM, hard disk, optical disk, and flash memory). For example, the one or more processors can implement functions of the waveform embedding 176 (FIG. 6) and information extraction 200 (FIG. 7). The one or more processors can perform calculations based on one or more of Equations 1 to 32. The waveform embedding 176 and information extraction 200 can also be implemented using application-specific integrated circuits (ASICs). The term “computer-readable medium” refers to a medium that participates in providing instructions to a processor for execution, including without limitation, non-volatile media (e.g., optical or magnetic disks), and volatile media (e.g., memory) and transmission media. Transmission media includes, without limitation, coaxial cables, copper wire and fiber optics.
The features described above can be implemented advantageously in one or more computer programs that may be executable on a programmable system including at least one programmable processor coupled to receive data and instructions from, and to transmit data and instructions to, a data storage system, at least one input device, and at least one output device. A computer program may be a set of instructions that can be used, directly or indirectly, in a computer to perform a certain activity or bring about a certain result. A computer program can be written in any form of programming language (e.g., C, Java), including compiled or interpreted languages, and it can be deployed in any form, including as a stand-alone program or as a module, component, subroutine, a browser-based web application, or other unit suitable for use in a computing environment.
Suitable processors for the execution of a program of instructions include, e.g., both general and special purpose microprocessors, digital signal processors, and the sole processor or one of multiple processors or cores, of any kind of 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 executing instructions and one or more memories for storing instructions and data. Generally, a computer will also include, or be operatively coupled to communicate with, one or more mass storage devices for storing data files; such devices include magnetic disks, such as internal hard disks and removable disks; magneto-optical disks; and optical disks. Storage devices suitable for tangibly embodying computer program instructions and data include all forms of non-volatile memory, including by way of example semiconductor memory devices, such as EPROM, EEPROM, and flash memory devices; magnetic disks such as internal hard disks and removable disks; magneto-optical disks; and CD-ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, ASICs (application-specific integrated circuits). Flash memory and/or any type of known memory are within the scope of the invention.
While this specification contains many specific implementation details, these should not be construed as limitations on the scope of any inventions or of what may be claimed, but rather as descriptions of features specific to particular embodiments of particular inventions. Certain features that are described in this specification 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.
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. In certain circumstances, multitasking and parallel processing may be advantageous. Moreover, the separation of various system components in the embodiments described above should not be understood as requiring such separation in all embodiments, and it should be understood that the described program components and systems can generally be integrated together in a single software product or packaged into multiple software products.
Thus, particular embodiments of the subject matter have been described. Other embodiments are within the scope of the following claims. In some cases, the actions recited in the claims can be performed in a different order and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In certain implementations, multitasking and parallel processing may be advantageous.
A number of implementations of the invention have been described. Nevertheless, it will be understood that various modifications can be made without departing from the spirit and scope of the invention. For example, the system 150 can combine time division multiplexing with frequency division multiplexing to allow a base station to communicate with a first group of terminal devices at a first frequency and a second group of terminal devices at a second frequency. In this example, the base station can perform a hand-shaking process with two or more terminal devices at different frequencies. The base station 152 and terminal devices 154 can include more components that are not shown in the figures. For example, standard components such as analog-to-digital converters, digital-to-analog converters, and power amplifiers have been omitted from the figures.
Referring to FIG. 35, in some implementations, a mobile device 640 includes a first communication module 642 using a communication protocol having a substantially symmetric architecture, such as the Wi-Fi, Bluetooth, or cellular (e.g., GSM, CDMA, LTE) communication protocol. The mobile device 640 includes a second communication module 644 using a communication protocol having an asymmetric architecture, such as the asymmetric time-reversal communication protocol. The mobile device 640 uses the first communication module 642 to communicate with a base station 646 that uses the communication protocol having a substantially symmetric architecture. The mobile device 640 uses the second communication module 644 to communicate with devices 648 that use the communication protocol having an asymmetric architecture. For example, the mobile device can be, e.g., a mobile phone, a laptop computer, a tablet computer, a smart watch, or a head-mounted display. The base station 646 can be a Wi-Fi base station or a wireless telephone base station. The devices 648 can be, e.g., sensors mounted at various locations in a building or electronic devices worn on a user (e.g., body temperature sensor, heart rate monitor, insulin level monitor, hearing aid, or smart contact lens).
Referring to FIG. 36, in some implementations, a system 650 includes a group of robots 652 that communicate with devices 646 outside the group of robots using a communication protocol that has a substantially symmetrical architecture (e.g., Wi-Fi, Bluetooth, or cellular protocol, such as GSM, CDMA, or LTE). The device 646 can be, e.g., a Wi-Fi base station, a wireless telephone base station, a Bluetooth peripheral device, or another robot that communicates by Wi-Fi or Bluetooth. The communication within the group of robots 652 is performed using a communication protocol that has an asymmetrical architecture, such as asymmetrical time-reversal communication protocol. The robots 652 can be robotic arms operating at assembly lines in a factory. The robots can be autonomous vehicles, such as autonomous land-based vehicles, aquatic vehicles, or aerial vehicles. The robots 652 can be made small, having the sizes of insects, birds, or small animals, and used to monitor the environment. For example, one or more of the robots 652 can function as the base station, and the other robots 652 can function as terminal devices. Deploying a group of robots as shown in FIG. 36 may be useful for, e.g., search and rescue missions in collapsed buildings or tunnels where there may be severe multipath interference. For example, the group of robots 652 can be scouting robots used to explore mine shafts or underground caves.
In some implementations, components of a robot may communicate with one another using the time-reversal communication protocol. For example, a robot can have multiple arms each having multiple fingers, and each arm or finger may have multiple sensors and actuators. The sensors and actuators of the fingers and arms may communicate with one another or with a central controller using the asymmetrical time-reversal communication protocol, and the central controller may communicate with external devices using symmetrical communication protocols.
Referring to FIG. 37, in some implementations, a vehicle 660 may include a communication unit 662 that communicates with one or more devices 646 outside the vehicle 660 using a communication protocol that has a substantially symmetrical architecture, such as Wi-Fi, Bluetooth, or cellular communication protocol. The communication unit 662 communicates with one or more devices 664 within the vehicle 660 (or attached to the vehicle, or associated with the vehicle) using a communication protocol that has an asymmetrical architecture, such as the asymmetrical time-reversal communication protocol. For example, the communication unit 662 may communicate with a wireless telephone base station using GSM, CDMA, LTE, or LTE Advanced protocols. The communication unit 662 may communicate with a controller of an electric vehicle charging station using the Wi-Fi or Bluetooth protocol. The devices 664 can be sensors (e.g., temperature sensors, motion sensors, light sensors, acoustic sensors), information display devices, or audio/video devices (playback or recording devices). For example, car seats may be embedded with accelerometers and pressure sensors that measure the forces imparted to passengers and communicate with the communication unit 662 using a time-reversal communication protocol. A controller may use the information from the accelerometers and pressure sensors to control an active suspension system under the seats to provide a smoother ride to passengers. The vehicle 660 can be, e.g., a car, a bus, a ship, or an airplane. For example, a bus may have an electronic library that delivers digital content to e-readers provided to passengers using the time-reversal communication protocol.
In some implementations, the terminal device 154 (FIG. 4) can be a low-power wearable device, such as a body temperature sensor, a heart rate monitor, an electrocardiograph sensor, an insulin level monitor, a hearing aid, or a smart contact lens. The base station 152 can be a health monitor unit such as a electrocardiography recorder, a mobile phone, a tablet computer, a smart watch, a smart wristband, or other wearable fitness device having a data processor that executes a health and fitness program.
In some implementations, the base station 152 can be part of a control center of a building, and the terminal devices 154 can be sensors mounted at various locations in the building. The sensors can be, e.g., motion sensors, temperature sensors, humidity sensors, light sensors, or acoustic sensors. Because time-reversal communication requires little power consumption at the terminal devices, the sensors can be powered by photovoltaic devices and made at low costs. Many sensors can be placed at various locations in a building without the need to provide wiring to the sensors. For example, this allows a building manager to conveniently monitor various parameters in the building to help optimize utilization of various facilities in the building. For example, light sensors can be placed near conference room tables, and data from the light sensors can be used to control curtains, shutters, or the transmittance of window panels having variable transmittances, in order to optimize lighting conditions for users sitting around the tables. Various light reflectors and light guides can be adjusted based on data from light sensors to optimize the distribution of natural light in a building.
In some implementations, a mobile phone is configured to communicate with a wireless telephone base station using a cellular communication protocol, and communicate with a smart watch using either a second communication protocol that has an asymmetrical architecture or a third communication protocol that has a substantially symmetrical architecture. For example, the second communication protocol can be an asymmetrical time-reversal communication protocol, and the third communication protocol can be, e.g., the Wi-Fi or Bluetooth protocol.
In some implementations, a mobile device can be configured to communicate with a wireless telephone base station using a cellular communication protocol, and communicate with a network router using a communication protocol having an asymmetrical architecture, such as the asymmetrical time-reversal communication protocol. For example, this may be useful in a convention center where thousands of attendees may need to use mobile devices to access the Internet wirelessly. Conventional Wi-Fi based network routers may not be able to support thousands of users within a crowded space due to severe interference, whereas a network router using the time-reversal communication protocol may be able to support all the users by taking advantage of the spatial filtering characteristics of the conference rooms.
The time-reversal communication system may use a wider bandwidth than Wi-Fi, Bluetooth, and cellular communication systems. For example, the time-reversal communication system may use a 500 MHz bandwidth, an LTE system may use a 20 MHz bandwidth, and an LTE Advanced system may use a 100 MHz bandwidth.
In some implementations, a mobile device may communicate with a server using a first communication protocol having a substantially symmetric architecture, download a map and information about locations of access points that use a second communication protocol having an asymmetric architecture. When the mobile device moves to a location near the access points that use the second communication protocol, the mobile device uses the second communication protocol to establish links to the access points. For example, the access points can be network routers that allow the mobile device to access the Internet or a corporate intranet. The first communication protocol can be, e.g., the Wi-Fi, Bluetooth, or cellular communication protocol. The cellular communication protocol can be, e.g., GSM, CDMA, LTE, or LTE Advanced. The second communication protocol can be, e.g., the time-reversal communication protocol.
In this description, the terms used for various communication protocols are meant to encompass variants of the corresponding protocols. For example, the CDMA protocol is meant to include one or more of Wideband CDMA (WCDMA) and Time Division Synchronous CDMA (TD-SCDMA) protocols. The Wi-Fi protocol is meant to include one or more of IEEE 802.11, 802.11a, 802.11b, 802.11g, 802.11n, 802.11ac, and 802.11ad protocols.
It is envisioned that any and all of the implementations described herein may be used in what is commonly referred to as the “internet of things”. That is, communications may be between any device that is powered in any way, including by battery power, solar power, wind power, a fuel cell, a capacitor, and energy storage unit, a fly wheel, a power line, the mains and the like. Devices that may implement these protocols include phones, computers, robots, meters, lights, controllers, appliances, vehicles, displays, televisions, entertainment equipment, cooking equipment, heating/cooling equipment, safety systems, lighting systems, monitoring systems, tracking systems, and the like.
The term “time-reversal communication protocol” is meant to cover a communication protocol in which a base station uses information derived from a probe signal sent from a terminal device to generate a signature waveform for the terminal device, and embed the signature waveform into data signals transmitted to the terminal device so that the data signals, after propagating through a multipath channel, focus at the terminal device. The time-reversal communication protocol may not require the system to compute a time-reversed version of a received signal. A time-reversal multiuser downlink system may not necessarily compute a time-reversed version of a signal received from any user. In some implementations, the signature waveform is a time-reversed version of the channel response signal. In some implementations, the signature waveform can be a modified version of the time-reversed channel response signal. For example, the signature waveform can be designed to reduce interference. To generate the signature waveform, it may not be necessary to first generate a time-reversed signal of the channel response signal, for example, the time-reverse operation may be moved toward a later stage of the process, i.e., design the signature based on the channel response and then time-reverse the resulting signature. Another method of generating the signature waveform is to separate the time-reverse operation into many element-swap operations and distribute these operations in the signature waveform design algorithm.
Accordingly, other embodiments are within the scope of the following claims.

## Claims

1. An apparatus, comprising:
an input circuit configured to derive a channel response signal from a probe signal sent from a first device to the apparatus through multiple wireless propagation paths, and to receive a combined signals that includes a first signal from the first device and a second signal from a second device; and
a data processor configured to calculate a signature waveform for the first device based on the derived channel response signal, and to determine the first signal based on the combined signal and the calculated signature waveform;
wherein the data processor is further configured to generate a combined downlink signal from a combination of a first downlink waveform and a second downlink waveform, the first downlink waveform including data intended for the first device, the second downlink waveform including data intended for the second device.

2. The apparatus of claim 1 in which the apparatus comprises a base station that comprises the input circuit and the data processor.

3. The apparatus of claim 1 in which the data processor determines the first signal by performing a convolution operation, or an operation similar to the convolution operation, between the combined signal and the calculated signature waveform.

4. The apparatus of claim 1 in which the data processor determines the first signal by filtering the combined signals using the calculated signature waveform as a filtering parameter.

5. The apparatus of claim 1 in which the data processor calculates the signature waveform from the first for the first device by time-reversing the derived channel response signal.

6. The apparatus of claim 1 in which the data processor is further configured to determine the first downlink waveform based on the data intended for the first device and the calculated signature waveform.

7. The apparatus of claim 6 in which the data processor determines the first downlink waveform by performing a convolution, or an operation similar to the convolution, of the data intended for the first device and the calculated signature waveform.

8. The apparatus of claim 1 in which the combined downlink signal is configured to enable each of the first device and the second device to receive multipath signals that can be used to determine, by the first device, a first data signal intended for the first device, and to determine, by the second device, a second data signal intended for the second device, wherein the first and the second devices receiving the combined downlink signal through different propagation paths and respectively retrieving the first data signal and the second data signal.

9. The apparatus of claim 1 in which the input circuit is further configured to receive another combined signals that includes signals from three or more devices that communicate with the apparatus, and the data processor is further configured to calculate a corresponding signature waveform for each of the devices.

10. The apparatus of claim 1 in which the first device does not calculate the signature waveform.

11. The apparatus of claim 1 in which the first device does not perform any signal processing that uses the signature waveform.

12. The apparatus of claim 1 comprising a processor configured to reduce inter-symbol interference (ISI) or inter-user interference (IUI) or both.

13. The apparatus of claim 1 in which the probe signal comprises at least one of (i) substantially a pulse signal, (ii) a signal having a predesigned waveform, or (iii) a signal that includes data signals.

14. A system, comprising:
a base station; and
a first circuit configured to send a probe signal through multiple wireless propagation paths to the base station;
wherein the base station comprises an input circuit configured to receive a channel response signal derived from the probe signal,
wherein the base station is configured to receive combined signals that include a signal from the first circuit and a signal from a second circuit,
wherein the base station comprises a data processor configured to calculate a signature waveform for the first circuit based on the channel response signal, and determine the signal sent from the first circuit based on the combined signals and the signature waveform for the first circuit, and
wherein the data processor of the base station is configured to determine a combined downlink signal from a combination of a first downlink waveform and a second downlink waveform, the first downlink waveform including data intended for the apparatus, the second downlink waveform including another data intended for another apparatus.

15. The apparatus of claim 14 wherein the probe signal is substantially a pulse signal.

16. The apparatus of claim 14 wherein the probe signal is a predesigned waveform.

17. The apparatus of claim 14 wherein the probe signal includes data signals.

18. The apparatus of claim 14 further configured to receive data signals encoded on waveforms derived from a channel impulse waveform.

19. The apparatus of claim 14 in which the circuit is configured to receive the combined downlink signal from the base station and apply rate matching and channel decoding to the combined downlink signal.

20. A wireless communications system comprising:
a base station comprising a data processor,
at least two devices that communicate with the base station,
wherein the at least two devices each sends a probe signals to the base station,
wherein the base station detects and stores a channel response signals respectively derived from each of the received probe signal,
wherein the base station uses the stored channel response signals associated with each of the at least two devices to preferentially communicate with the each of the last least two devices; and
wherein the data processor of the base station is configured to generate a combined downlink signal from a combination of a first downlink waveform and a second downlink waveform, the first downlink waveform including data intended for one of the at least two devices, the second downlink waveform including data intended for another one of the at least two devices.

21. The wireless communications system of claim 20, comprising three or more devices.