# Patent

## Title
Determining channel coefficients in a multipath channel

## Abstract

A device, method, and computer-readable medium are provided. The device includes, but is not limited to, an antenna configured to receive a multipath signal and a processor operably coupled to the antenna to receive the multipath signal. The processor is configured to determine a training signal transmitted from a second device to create the received multipath signal; sample the received multipath signal; and determine channel coefficients based on the sampled multipath signal and the determined training signal using a mixed-norm convex optimization process, wherein the channel coefficients characterize a channel associated with the multipath signal.

## Background
REFERENCE TO GOVERNMENT RIGHTS
This invention was made with United States government support awarded by the following agencies: Navy/ONR N00014-06-1-1123. The United States government has certain rights in this invention.


BACKGROUND
Knowledge of channel state information (CSI) at a transmitter and/or a receiver can have a dramatic impact on communication over multipath channels in terms of power and spectral efficiency. In particular, coherent communication is generally far more efficient than non-coherent communication, but requires that the CSI be known at the receiver to maintain coherence. Training-based schemes, which involve probing the channel with known signaling waveforms and linear processing of the corresponding channel output, are commonly used to learn the CSI at the receiver. Recent measurement studies have shown that physical multipath channels tend to exhibit an approximately sparse multipath structure at high signal space dimension—time-bandwidth-antenna product. The CSI for such channels is characterized with significantly fewer dominant parameters compared to the maximum number dictated by the delay-Doppler-angle spread of the channel. Conventional training-based methods, often based on exhaustive probing coupled with a linear least squares approach or a non-linear parametric estimator, are ill-suited for exploiting the inherently low-dimensional nature of sparse or approximately sparse multipath channels.
SUMMARY
In an illustrative embodiment, a device is provided that includes, but is not limited to, an antenna configured to receive a multipath signal and a processor operably coupled to the antenna to receive the multipath signal. The processor is configured to determine a training signal transmitted from a second device to create the received multipath signal, sample the received multipath signal, and determine channel coefficients based on the sampled multipath signal and the determined training signal using a mixed-norm convex optimization process, wherein the channel coefficients characterize a channel associated with the multipath signal.
In another illustrative embodiment, a method of determining channel coefficients characterizing a channel associated with a multipath signal is provided. A multipath signal is received at a first device. Using a processor at the first device, a training signal transmitted from a second device to create the received multipath signal is determined, the received multipath signal is sampled, and the channel coefficients are determined based on the sampled multipath signal and the determined training signal using a mixed-norm convex optimization process.
In yet another illustrative embodiment, a computer-readable medium is provided. The computer-readable medium has computer-readable instructions stored thereon which, when executed by a processor, cause a device to perform the operations of the method of determining channel coefficients characterizing a channel associated with a multipath signal.
Other principal features and advantages of the invention will become apparent to those skilled in the art upon review of the following drawings, the detailed description, and the appended claims.



BRIEF DESCRIPTION OF THE DRAWINGS
Illustrative embodiments of the invention will hereafter be described with reference to the accompanying drawings, wherein like numerals denote like elements.
FIG. 1 shows a system in a scattering environment in accordance with an illustrative embodiment.
FIG. 2 shows a block diagram of a transceiver device of the system of FIG. 1 in accordance with an illustrative embodiment.
FIG. 3 depicts a flow diagram illustrating exemplary operations performed by the transceiver device of FIG. 2 in accordance with an illustrative embodiment.
FIG. 4 shows a schematic illustrating a virtual representation of a single-antenna, frequency-selective channel in accordance with an illustrative embodiment.
FIG. 5a shows example results derived using a least squares estimation method in accordance with an illustrative embodiment.
FIG. 5b shows example results derived using a lasso estimation method in accordance with an illustrative embodiment.
FIG. 6 shows a schematic illustrating a virtual representation of a single-antenna, doubly-selective channel in accordance with an illustrative embodiment.
FIG. 7 shows a schematic illustrating a virtual representation of an angle-delay-Doppler channel in accordance with an illustrative embodiment.



DETAILED DESCRIPTION
With reference to FIG. 1, a system 100 is shown in accordance with an illustrative embodiment. System 100 is located within a scattering environment 101 that may include, for example, a first scatterer 106a and a second scatterer 106b. There may be any number of scatterers in scattering environment 101 of the same or different types. Scatterers include any objects in the physical environment such as buildings, vehicles, vegetation, the terrain, weather, etc. capable of reflecting electromagnetic waves.
With reference to the illustrative embodiment of FIG. 1, system 100 is a communication system that may include a first transceiver device 102 and a second transceiver device 104. In an illustrative embodiment, first transceiver device 102 may include a base station of a wireless communication system and second transceiver device 104 may include a computer of any form factor, a personal data assistant, an instant messaging device, a phone, a radio, etc. In another illustrative embodiment, first transceiver device 102 may include a personal data assistant, a computer of any form factor, an instant messaging device, a phone, a radio, etc. In an alternative embodiment, system 100 may include first transceiver device 102 that is a sensor system, such as a radar or a sonar, to detect targets of various types within scattering environment 101.
With reference to the illustrative embodiment of FIG. 1, a wireless signal 108 transmitted from first transceiver device 102 to second transceiver device 104 is received at second transceiver device 104 as a delayed, attenuated, and, if either or both of first transceiver device 102 and second transceiver device 104 are moving, a Doppler shifted version of wireless signal 108. Wireless signal 108 also may be reflected and scattered by first scatterer 106a resulting in a second wireless signal 110 being received at second transceiver device 104. Second wireless signal 110 also is a delayed, attenuated, and possibly Doppler shifted version of wireless signal 108. Wireless signal 108 also may be reflected and scattered by second scatterer 106b resulting in a third wireless signal 112 being received at second transceiver device 104. Third wireless signal 112 also is a delayed, attenuated, and possibly Doppler shifted version of wireless signal 108. Wireless signal 108 also may be reflected and scattered by second scatterer 106b resulting in a fourth wireless signal 114 being received at first transceiver device 102. Fourth wireless signal 114 also is a delayed, attenuated, and possibly Doppler shifted version of wireless signal 108.
Second wireless signal 110, fourth wireless signal 114, and third wireless signal 112 are multipath signals. Multipath signal propagation results in multiple spatially distributed receive paths. Multipath propagation leads to signal fading—fluctuations in received signal strength—that can severely impact reliable communication and sensing functions. On the other hand, multipath is also a source of diversity based on the multiple statistically independent communication paths that can increase the rate and reliability of the communication and sensing functions. The impact of multipath fading versus diversity on performance depends on the amount of channel state information (CSI) available to the system. For example, knowledge of instantaneous CSI at the receiving device for coherent reception enables exploitation of diversity to combat fading. Additionally, gains in capacity and reliability are possible if even partial CSI is available at the transmitting device as well.
The statistical characteristics of a wireless channel or sensing environment depend on the interaction between scattering environment 101 and the signal space of first transceiver device 102 and second transceiver device 104. Signal space parameters include bandwidth (i.e. narrowband, wideband), the number of transceiver antennas, the antenna spacing, pulse duration, frequency, etc. Accurate channel modeling and characterization in time, frequency, and space, as a function of multipath and signal space characteristics enables exploitation of technological advances including frequency and bandwidth agility and reconfigurable antenna arrays and support coherent processing.
With reference to FIG. 2, a block diagram of first transceiver device 102 is shown in accordance with an illustrative embodiment. First transceiver device 102 may include an antenna 200, a transceiver 202, a processor 204, a computer-readable medium 206, and a channel estimation application 208. Second transceiver device 104 may include similar components. In another illustrative embodiment, first transceiver device 102 and/or second transceiver device 104 may include a separate transmitter and receiver. First transceiver device 102, of course, may include different and additional components such as a signal generator, a signal processor, etc. not explicitly shown in FIG. 2. First transceiver device 102 may include a plurality of antennas arranged in an array. The number and type of antennas of first transceiver device 102 may be different from the number and type of antennas of second transceiver device 104. The plurality of antennas may be arranged to form a uniform or a non-uniform linear array, a rectangular array, a circular array, a conformal array, etc. An antenna of the plurality of antennas may be a dipole antenna, a monopole antenna, a helical antenna, a microstrip antenna, a patch antenna, a fractal antenna, etc. The plurality of antennas may be reconfigurable spatially such that an antenna spacing between the plurality of antennas can be adjusted. Antenna 200 may be steerable.
Transceiver 202 processes electromagnetic signals received by antenna 200 and to be transmitted by antenna 200 under control of processor 204. First transceiver device 102 and/or second transceiver device 104 may include a plurality of transceivers that use the same or a different transmission/reception technology.
Processor 204 executes computer-readable instructions as known to those skilled in the art. The computer-readable instructions may be carried out by a special purpose computer, logic circuits, or hardware circuits. Thus, processor 204 may be implemented in hardware, firmware, software, or any combination of these methods. The term “execution” is the process of running an application or the carrying out of the operation called for by a computer-readable instruction. The computer-readable instructions may be written using one or more programming language, scripting language, assembly language, etc. Processor 204 executes a computer-readable instruction, meaning that it performs the operations called for by that computer-readable instruction. Processor 204 operably couples with transceiver 202 and computer-readable medium 206 to receive, to send, and to process information and/or to control the operations of first transceiver device 102. Processor 204 may retrieve a set of computer-readable instructions from a permanent memory device such as a read only memory device (ROM) and copy the instructions in an executable form to a temporary memory device that is generally some form of random access memory (RAM). First transceiver device 102 and/or second transceiver device 104 may include a plurality of processors that use the same or a different processing technology.
Computer-readable medium 206 is an electronic holding place or storage for information so that the information can be accessed by processor 204 as known to those skilled in the art. Computer-readable medium 206 can include, but is not limited to, any type of RAM, any type of ROM, any type of flash memory, etc. such as magnetic storage devices (e.g., hard disk, floppy disk, magnetic strips, . . . ), optical disks (e.g., compact disk (CD), digital versatile disk (DVD), . . . ), smart cards, flash memory devices, etc. First transceiver device 102 and/or second transceiver device 104 may have one or more computer-readable media that use the same or a different memory media technology. First transceiver device 102 and/or second transceiver device 104 also may have one or more drives that support the loading of a memory media such as a CD, a DVD, a flash memory card, etc.
Channel estimation application 208 includes operations which estimate and/or characterize the channel state within scattering environment 101 based on the signal space characteristics. The operations may be implemented using hardware, firmware, software, or any combination of these methods. With reference to the illustrative embodiment of FIG. 2, channel estimation application 208 is implemented in software, i.e. computer-readable instructions, stored in computer-readable medium 206 and accessible by processor 204 for execution of the computer-readable instructions that embody the operations of channel estimation application 208. The computer-readable instructions of channel estimation application 208 may be written using one or more programming languages, assembly languages, scripting languages, etc. The functionality provided by channel estimation application 208 may be distributed among one or more modules and across one or more devices. In an illustrative embodiment, first transceiver device 102 and/or second transceiver device 104 may characterize the channel state and send the channel state information back from the receiver to the transmitter to enable transceiver adaptation for enhanced spectral efficiency and interference management. For example, the channel state may be characterized by channel coefficients.
With reference to FIG. 3, illustrative operations associated with channel estimation application 208 of FIG. 2 are described. Additional, fewer, or different operations may be performed, depending on the embodiment. The order of presentation of the operations of FIG. 3 is not intended to be limiting. In an operation 300, a training signal is determined. In an operation 302, the determined training signal is transmitted using transceiver 202 and antenna 200. In an operation 304, a multipath signal is received. Depending on the embodiment, the training signal determination and transmission and the multipath signal reception may be performed at the same transceiver device. For example, if the transceiver device is associated with a radar or a sonar system. In an alternative embodiment, the training signal determination and transmission may be performed at first transceiver device 102 while the multipath signal reception may be performed at second transceiver device 104 or vice versa. In an illustrative embodiment, the training signal is determined at the transceiver device receiving the multipath signal as part of signaling from the transmitting transceiver device, an a priori determination as part of the system configuration, etc. so that the training signal is known at the transceiver device receiving the multipath signal and can thus be determined. In an illustrative embodiment, the training signal may not change once determined. In another illustrative embodiment, the training signal may be dynamically determined based, for example, on scattering environment 101 and/or on the current signal space characteristics for a given communication or sensing function.
In an operation 306, the received multipath signal is sampled. For example, the received multipath signal may be sampled in time, i.e. chip-rate sampling, projected onto a receive side signaling basis functions, i.e. multitone waveforms, to yield a discrete version of the received multipath signal, etc. In an operation 308, channel coefficients are determined based on the sampled multipath signal and the determined training signal using a mixed-norm convex optimization process. As used herein, a mixed-norm convex optimization process is a mathematical process for optimization of a convex objective function, possibly subject to convex equality and convex inequality constraints, such that the entire optimization setup (objective function and constraints, if they exist) contain at least two different vector norms, for example, l2 and l1 norms, l1 and l∞ norms, etc. In an illustrative embodiment, the mixed-norm convex optimization process is a Dantzig selector or a Lasso estimator. In an operation 310, a determination is made concerning whether or not to repeat the channel estimation process. For example, the channel estimation process may be repeated periodically, for each reception, at an initiation of a communication link, with a change in signaling parameters such as a frequency, bandwidth, or antenna spacing change, etc.
Illustrative channel estimation processes based on different signal space characteristics are described with additional detail below. A channel model is formed, a training signal is selected, and an estimator is applied for each example process.
In a first illustrative embodiment, first transceiver device 102 and second transceiver device 104 include single antennas and the channel is identified as frequency-selective because the channel impulse response remains constant over the time duration of interest and the channel maximum delay spread Tm is large relative to the inverse of the communication bandwidth W, such that TmW≧1. A “channel model” is a discrete approximation of the continuous time channel at the delay resolution dictated by the channel signaling parameters. The channel model captures the relationship between the clustering of physical paths and the sparsity of dominant channel coefficients. The complex baseband transmitted and received signals can be related as




  
    
      
        
          y
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              ∫
              0
              
                T
                m
              
            
            ⁢
            
              
                h
                ⁡
                
                  (
                  τ
                  )
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    τ
                  
                  )
                
              
              ⁢
              
                d
                τ
              
            
          
          +
          
            z
            ⁡
            
              (
              t
              )
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




where h(τ) is the channel impulse response, x(t), y(t), and z(t) represent the transmitted, the received, and the additive white Gaussian noise (AWGN) waveforms, respectively.

Frequency-selective channels generate multiple delayed and attenuated copies of the transmitted waveform. For such “multipath” channels, h(τ) can be modeled as




  
    
      
        
          h
          ⁡
          
            (
            τ
            )
          
        
        =
        
          
            ∑
            
              i
              =
              1
            
            
              N
              p
            
          
          ⁢
          
            
              α
              i
            
            ⁢
            
              δ
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    i
                  
                
                )
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




and the transmitted and received waveforms related by





  
    
      
        
          y
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              ∑
              
                i
                =
                1
              
              
                N
                p
              
            
            ⁢
            
              
                α
                i
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    
                      τ
                      i
                    
                  
                  )
                
              
            
          
          +
          
            z
            ⁡
            
              (
              t
              )
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




which corresponds to signal propagation along Np physical paths, where αiεC and τiε[0, Tm] are the complex path gain and the delay associated with the i-th physical path, respectively.

The discrete path model of equation (2), while realistic, is difficult to analyze and identify due to nonlinear dependence on the real valued delay parameters {τi}. However, because the communication bandwidth W is limited, the continuous-time channel can be accurately approximated by a channel model with the aid of sampling theorems and/or power series expansions. The channel model provides a discrete approximation of frequency-selective channels by uniformly sampling the physical multipath environment in delay at a resolution commensurate with W, i.e.,




  
    
      
        
          y
          ⁡
          
            (
            t
            )
          
        
        ≈
        
          
            
              ∑
              
                i
                =
                0
              
              
                p
                -
                1
              
            
            ⁢
            
              
                β
                j
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    
                      j
                      /
                      W
                    
                  
                  )
                
              
            
          
          +
          
            z
            ⁡
            
              (
              t
              )
            
          
        
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          β
          j
        
        ≈
        
          
            ∑
            
              i
              ∈
              
                S
                
                  τ
                  ,
                  j
                
              
            
          
          ⁢
          
            
              α
              i
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              c
              ⁡
              
                (
                
                  j
                  -
                  
                    W
                    ⁢
                    
                      
                    
                    ⁢
                    
                      τ
                      i
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




where p=┌TmW┐+1, sinc(a)=sin(πa)/πa, and Sτ,j={i:τiε[j/W−1⁄2W,j/W+1⁄2W]}denotes the set of all physical paths whose delays lie within the delay resolution bin of width Δτ=1/W centered around the j-th resolvable virtual delay, {circumflex over (τ)}j=j/W·{βj}are termed as the channel coefficients in the delay space.

With reference to FIG. 4, a plurality of resolution bins 400 having a resolution of 1/W are shown for a time window [0, Tm]. The plurality of vertical arrows 402 indicate gains received from different paths at a delay τ. Equation (5) states that the coefficient βj for each of the plurality of resolution bins 400 is approximately the sum of the gains of all paths within the j-th delay resolution bin obtained at the center of each of the plurality of resolution bins 400. The accuracy of the approximation increases with the bandwidth W due to the resulting increase in delay resolution.
Based on the interarrival times between different multipath clusters within the delay spread, channels can be categorized as either “rich” or “sparse”. In a rich multipath channel, the interarrival times are smaller than the inverse of the communication bandwidth W. Sparse multipath channels, on the other hand, exhibit interarrival times that are larger than 1/W, for example as shown in FIG. 4, where not every resolution bin contains a multipath component. In particular, since a channel coefficient consists of the sum of the gains of all paths falling within its respective delay bin, sparse frequency-selective channels tend to have far fewer nonzero channel coefficients than the number of the plurality of resolution bins 400 at any fixed, but large enough bandwidth. Thus, the multipath channel may be considered sparse if the number of nonzero channel coefficients is less than the number of resolvable paths p that is p=┌TmW┐+1.
The continuous-time received signal y(t) can be chip-rate sampled at a rate of 1/W at the receiver to obtain an equivalent discrete time representation of equation (4) as

y=x*β+zy=Xβ+z  (6)

where * denotes the discrete convolution operator, βεCp is the vector of channel coefficients, and x, y, and z are the vectors of samples of x(t), y(t), and z(t), respectively. The dimensions of x, y, and z are dictated by the input signaling duration Ts, where x(t)=0 for all t∉[0,Ts]. n(=┌TsW┐) denotes the number of nonzero samples of x(t), resulting in xεCn and y,z εCn+p−1. The matrix X is an (n+p−1)×p Toeplitz-structured, convolution matrix formed from vector x.

In an illustrative embodiment, β is modeled as a deterministic, but unknown vector. A power constraint also may be applied on the transmitted signal x(t) that can be readily translated into an average power constraint on the entries of x. Without loss of generality, it can be assumed that maxk E[|xk|2]≦1 and that the entries of z correspond to an independent and identically distributed (iid) complex Gaussian white noise sequence.
The virtual representation of a frequency-selective channel captures its essential characteristics in terms of the channel coefficients {βj}, the determined (selected) training signal x, and the received multipath signal y. An estimator {circumflex over (β)}CS of the sparse channel can be formed by solving either an  penalized least squares problem using the “lasso” estimator, as described in R. Tibshirani, “Regression shrinkage and selection via the lasso,” J. Roy. Statist Soc. Ser. B, pp. 267-288, 1996 (hereinafter referred to as the lasso estimator), or a linear program called the “Dantzig selector” as described in E. Candés and T. Tao, “The Dantzig selector: Statistical estimation when p is much larger than n,” Ann. Statist, pp. 2313-2351, December 2007 (hereinafter referred to as the Dantzig selector).
As a first illustrative training signal used in a frequency-selective channel, a multi-carrier signal comprising an OFDM basis for the N0-dimensional signal space is generated via frequency shifts of a fixed prototype pulse can be used of the form

g(t):γk(t)=g(t)ej2πkt/T,kεS={0, 1, . . . , N0−1}  (7)

where N0=TW. In an illustrative embodiment, the pulse is assumed to have unit energy, ∫|g(t)|2dt=1. The training waveform may randomly dedicate Nr of the N0 basis elements as “pilot tones” such that





  
    
      x
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        
          ɛ
          
            N
            r
          
        
      
      ⁢
      
        
          ∑
          
            m
            ∈
            
              S
              r
            
          
        
        ⁢
        
          
            γ
            m
          
          ⁡
          
            (
            t
            )
          
        
      
    
  
  ,
  
    0
    ≤
    t
    <
    T
  




where the set of indices of pilot tones, Sr consists of Nr elements randomly selected from S and ε is the transmit energy budget available for training purposes. At the receiver, projecting the (noisy) received signal y(t) onto the OFDM basis waveforms yields




  
    
      
        y
        m
      
      =
      
        
          
            〈
            
              y
              ,
              
                γ
                m
              
            
            〉
          
          ≈
          
            
              
                
                  ɛ
                  
                    N
                    r
                  
                
              
              ⁢
              
                H
                m
              
            
            +
            
              
                z
                m
              
              ⁢
              m
            
          
        
        ∈
        
          S
          r
        
      
    
    ⁢
    
      
    
  




  
    
      
        where
        ⁢
        
          
        
        ⁢
        
          〈
          
            y
            ,
            
              γ
              m
            
          
          〉
        
      
      =
      
        ∫
        
          
            y
            ⁡
            
              (
              t
              )
            
          
          ⁢
          
            
              
                γ
                m
              
              ⁡
              
                (
                t
                )
              
            
            _
          
          ⁢
          
            d
            t
          
        
      
    
    ,
    
      {
      
        z
        m
      
      }
    
  



corresponds to an AWGN sequence and the channel coefficients are given by Hm≈H(f) f|f=m/T.

The virtual representation implies that




  
    H
    ⁡
    
      (
      f
      )
    
  
  ≈
  
    
      ∑
      
        l
        =
        0
      
      
        p
        -
        1
      
    
    ⁢
    
      
        h
        l
      
      ⁢
      
        
          e
          
            
              -
              j2π
            
            ⁢
            
              l
              W
            
            ⁢
            f
          
        
        .
      
    
  




As a result, the OFDM channel coefficients Hm can be written as





  
    H
    m
  
  =
  
    
      
        ∑
        
          l
          =
          0
        
        
          p
          -
          1
        
      
      ⁢
      
        
          h
          l
        
        ⁢
        
          e
          
            
              -
              j
            
            ⁢
            
              
            
            ⁢
            2
            ⁢
            π
            ⁢
            
              l
              
                N
                0
              
            
            ⁢
            m
          
        
      
    
    =
    
      
        u
        
          f
          ,
          m
        
        T
      
      ⁢
      h
    
  




where h is the p-dimensional vector of channel coefficients,




  
    
      
        u
        
          f
          ,
          m
        
      
      =
      
        
          
            [
            
              1
              ⁢
              
                
              
              ⁢
              
                ω
                
                  N
                  0
                
                
                  1
                  ⁢
                  m
                
              
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              
                ω
                
                  N
                  0
                
                
                  
                    (
                    
                      L
                      -
                      1
                    
                    )
                  
                  ⁢
                  m
                
              
            
            ]
          
          T
        
        ∈
        
          C
          L
        
      
    
    ⁢
    
      
    
  




  where




  
    
      
        ω
        
          N
          0
        
      
      =
      
        
          
            e
            
              
                -
                j
              
              ⁢
              
                
                  2
                  ⁢
                  π
                
                
                  N
                  o
                
              
            
          
          ⁢
          
            
          
          ⁢
          and
          ⁢
          
            
          
          ⁢
          n
        
        =
        0
      
    
    ,
    1
    ,
    ...
    ⁢
    
      
    
    ,
    
      
        N
        0
      
      -
      1.
    
  



y=Xh+z where the Nr sensing matrix X is comprised of {√{square root over (ε|Nr)}uf,mT:mεSr} as its rows and the AWGN vector z is distributed as CN(0Nr,INr).

The virtual model channel coefficients can be determined as a solution to the Dantzig selector as




  
    
      
        
          β
          ^
        
        =
        
          
            
              argmin
              
                
                  β
                  ~
                
                ∈
                
                  C
                  p
                
              
            
            ⁢
            
              
                
                
                  β
                  ~
                
                
              
              
                l
                1
              
            
            ⁢
            
              
            
            ⁢
            subject
            ⁢
            
              
            
            ⁢
            to
            ⁢
            
              
            
            ⁢
            
              
                
                
                  
                    X
                    H
                  
                  ⁢
                  r
                
                
              
              
                l
                ∞
              
            
          
          ≤
          λ
        
      
    
    
      
        (
        8
        )
      
    
  




where {circumflex over (β)} are the estimated channel coefficients, ∥∥l1 is an  norm, ∥∥l∞ is an  norm, p a maximum number of channel coefficients, Cp is a p-dimensional complex space, X is the sensing matrix {√{square root over (ε/Nr)}uf,mT: mεSr}, r is the Nr dimensional vector of residuals: r=y−X{tilde over (h)}, y is the sampled multipath signal, and λis a tuning parameter selected to control a number of channel coefficients to be determined. In an illustrative embodiment, λ(p,ε)=√{square root over (2ε(1+a)log p)} where ε is a transmit energy budget for the training signal and any a≧0. In a non-sparse environment, the same method can be used to determine the first const·Nr/log4 N0 largest channel coefficients. In another illustrative embodiment, a lasso estimator also may be used to determine the channel model coefficients.

As a second illustrative training signal used in a frequency-selective channel, a pseudo-random sequence {xk}k=0n−1, where n is the input probe duration and the xk's are iid realizations from a zero mean, unit variance distribution f(x). In an illustrative embodiment, {xk}k=0n−1 is a sequence of IID random variables drawn from a Rademacher distribution where n≧c1·log p·S2 where is any constant c1≧0, S is the number of desired coefficients to recover. The virtual model channel coefficients can be determined as a solution to the Dantzig selector as in equation (9), where the entries of {circumflex over (β)} are the estimated channel coefficients, ∥ ∥l1 is an  norm, ∥∥l∞ is an norm, p is a maximum number of channel coefficients, Cp is a p-dimensional complex space, X is a sensing matrix, r=y−X{tilde over (β)} is a vector of residuals, y is the sampled multipath signal, and λ is a tuning parameter selected to control a number of channel coefficients to be determined. Based on the first illustrative training signal, p=┌TmW┐+1 is the maximum number of nonzero channel coefficients, X is the (n+p−1)×p Toeplitz-structured, convolution matrix formed from the first illustrative training signal, n is the input probe duration, r is an (n+p−1)-dimensional vector of residuals, and λ can be chosen as √{square root over (2(1+a)log p)}·σ for any a≧0, and σ is a square root of the noise power. In an alternative embodiment, the lasso estimator can be used to determine the virtual model channel coefficients. In a non-sparse environment, the same method can be used to determine the first (n+p−1) largest channel coefficients.
With reference to FIGS. 5a and 5b, a comparison between use of a least squares estimator and use of a lasso estimator to determine the channel coefficients is shown. A training signal comprising an n=128 pseudo-random input probe to sense a p=128 length channel that has only ten channel coefficients 500 which were nonzero was used and are indicated as impulses in FIGS. 5a and 5b. The output of the channel was observed at a signal-to-noise ratio of −10 dB. Channel coefficients can be obtained by solving the least squares problem or a variant of it as

{circumflex over (β)}LS=(XHX)−1XHy.  (10)

The least squares estimates were obtained by pseudo-inverting the matrix X as in equation (10). The lasso estimates were obtained by executing a gradient projection for sparse reconstruction as described in Figueiredo, et al., “Gradient projection for sparse reconstruction: Application to compressed sensing and other inverse problems,” IEEE J. Select. Topics Signal Processing, pp. 586-597, December 2007. The results using the least squares estimator are shown in FIG. 5a. The results using the lasso estimator are shown in FIG. 5b. Using the lasso estimator, every nonzero channel coefficient 500 was identified and only one noise-only coefficient 502 was identified. In contrast, with reference to FIG. 5a, even a clairvoyantly thresholded least squares estimator would be unable to identify all of the nonzero channel coefficients 500 without also identifying a significant number of noise-only coefficients.
In a second illustrative embodiment, first transceiver device 102 and second transceiver device 104 include single antennas and the channel is identified as doubly-selective because the delay spread and Doppler spread of the channel are large relative to the inverse of the signaling bandwidth and duration, respectively, i.e., Wτmax≧1 and Tνmax≧1. W is the two-sided bandwidth. τmax is the delay spread of the channel and is defined as the maximum possible nonzero delay introduced by the channel. T is the packet duration. νmax is the maximum possible (one-sided) Doppler shift caused by the channel. In the second illustrative embodiment, the channel includes underspread channels, characterized by τmaxνmax<<1, and there is no interpacket interference in time and/or frequency, such that T>>τmax and W>>νmax. The dimension of the receiver signal space is N0≈TW, the time-bandwidth product.
The complex baseband transmitted and received signals in the absence of noise can be related as




  
    
      
        
          
            
              
                y
                ⁡
                
                  (
                  t
                  )
                
              
              =
              
              ⁢
              
                
                  ∫
                  0
                  
                    τ
                    max
                  
                
                ⁢
                
                  
                    h
                    ⁡
                    
                      (
                      
                        t
                        ,
                        τ
                      
                      )
                    
                  
                  ⁢
                  
                    x
                    ⁡
                    
                      (
                      
                        t
                        -
                        τ
                      
                      )
                    
                  
                  ⁢
                  
                    d
                    τ
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                ∫
                
                  
                    H
                    ⁡
                    
                      (
                      
                        t
                        ,
                        f
                      
                      )
                    
                  
                  ⁢
                  
                    X
                    ⁡
                    
                      (
                      f
                      )
                    
                  
                  ⁢
                  
                    e
                    
                      j2π
                      ⁢
                      
                        
                      
                      ⁢
                      ft
                    
                  
                  ⁢
                  
                    d
                    f
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  ∫
                  0
                  
                    τ
                    max
                  
                
                ⁢
                
                  
                    ∫
                    
                      -
                      
                        V
                        
                          max
                          /
                          2
                        
                      
                    
                    
                      V
                      
                        max
                        /
                        2
                      
                    
                  
                  ⁢
                  
                    
                      C
                      ⁡
                      
                        (
                        
                          v
                          ,
                          τ
                        
                        )
                      
                    
                    ⁢
                    
                      x
                      ⁡
                      
                        (
                        
                          t
                          -
                          τ
                        
                        )
                      
                    
                    ⁢
                    
                      e
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        2
                        ⁢
                        π
                        ⁢
                        
                          
                        
                        ⁢
                        vt
                      
                    
                    ⁢
                    
                      d
                      v
                    
                    ⁢
                    
                      d
                      τ
                    
                  
                
              
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




where x(t) and y(t) represent the transmitted and received waveforms, respectively, and X(f) is the Fourier transform of x(t). The channel is characterized by the time-varying impulse response, h(t,τ), or the time-varying frequency response, H(t,f), or the delay-Doppler spreading function, C(ν,τ) which are all three equivalent channel characterizations related to each other via Fourier transforms.

Doubly-selective channels generate multiple delayed, Doppler-shifted, and attenuated copies of the transmitted waveform. A discrete path model may be used to capture the characteristics of doubly-selective channels in terms of the physical propagation paths. In the discrete path model, the delay-Doppler spreading function of the channel can be expressed as




  
    
      
        
          C
          ⁡
          
            (
            
              v
              ,
              τ
            
            )
          
        
        =
        
          
            ∑
            
              i
              =
              1
            
            
              N
              p
            
          
          ⁢
          
            
              α
              i
            
            ⁢
            
              δ
              ⁡
              
                (
                
                  v
                  -
                  
                    v
                    i
                  
                
                )
              
            
            ⁢
            
              δ
              ⁡
              
                (
                
                  τ
                  -
                  
                    τ
                    i
                  
                
                )
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




and the transmitted and received waveforms related by





  
    
      
        
          y
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              i
              =
              1
            
            
              N
              p
            
          
          ⁢
          
            
              α
              i
            
            ⁢
            
              e
              
                j2π
                ⁢
                
                  
                
                ⁢
                
                  v
                  i
                
                ⁢
                t
              
            
            ⁢
            
              x
              ⁡
              
                (
                
                  t
                  -
                  
                    τ
                    i
                  
                
                )
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




which corresponds to signal propagation along Np physical paths, where αiεC, νiε[−νmax/2,νmax/2], and τiε[0,τmax] are the complex path gain, the Doppler shift, and the delay associated with the i-th physical path, respectively.

Though realistic, equation (13) is difficult to analyze and learn due to its nonlinear dependence on a potentially large number of physical parameters αi,νi,τi. However, because of the finite signaling duration and bandwidth, the discrete path model can be accurately approximated by the channel model with the aid of sampling theorems and/or power series expansions. The channel model provides a low dimensional approximation of the discrete path model by uniformly sampling the physical multipath environment in the delay-Doppler domain at a resolution commensurate with W and T (Δτ=1/W, Δν=1/T ). That is,




  
    
      
        
          y
          ⁡
          
            (
            t
            )
          
        
        ≈
        
          
            ∑
            
              l
              =
              0
            
            
              L
              -
              1
            
          
          ⁢
          
            
              ∑
              
                k
                =
                
                  -
                  K
                
              
              K
            
            ⁢
            
              
                h
                
                  l
                  ,
                  k
                
              
              ⁢
              
                e
                
                  j2π
                  ⁢
                  
                    k
                    T
                  
                  ⁢
                  t
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    t
                    -
                    
                      l
                      /
                      W
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        14
        )
      
    
  
  
    
      
        
          h
          
            l
            ,
            k
          
        
        ≈
        
          
            ∑
            
              i
              ∈
              
                
                  S
                  
                    τ
                    ,
                    l
                  
                
                ⋂
                
                  S
                  
                    v
                    ,
                    k
                  
                
              
            
          
          ⁢
          
            
              α
              i
            
            ⁢
            
              e
              
                -
                
                  jπ
                  ⁡
                  
                    (
                    
                      k
                      -
                      
                        Tv
                        i
                      
                    
                    )
                  
                
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              c
              ⁡
              
                (
                
                  k
                  -
                  
                    Tv
                    i
                  
                
                )
              
            
            ⁢
            sin
            ⁢
            
              
            
            ⁢
            
              c
              ⁡
              
                (
                
                  l
                  -
                  
                    W
                    ⁢
                    
                      
                    
                    ⁢
                    
                      τ
                      i
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




where sinc(a)=sin(πa)/πa and L=┌Wτmax┐+1 and K=┌Tνmax/2 denote the maximum number of resolvable delays and one-sided, Doppler shifts within the delay-Doppler spreading function, respectively. The set Sτ,l={i:τiε[/W−1⁄2W,l/W+1⁄2W)]} is the set of indices of all paths whose delays lie within the delay resolution bin of width Δτ=1/W centered around the l-th virtual delay, {circumflex over (τ)}l=/W, while SV,K={i:νiε[k/T−1⁄2T,k/T+1⁄2T)]} denotes the set of indices of all paths whose Doppler shifts lie within the Doppler resolution bin of width Δν=1/T centered around the k-th virtual Doppler shift, {circumflex over (ν)}k=k/T The parameters hl,k are identified as the channel coefficients in the delay-Doppler domain. Equation (15) states that the channel coefficient hl,k approximately consist of the sum of the gains of all paths whose delays and Doppler shifts lie within the (l, k)-th delay-Doppler resolution bin of size Δτ×Δν centered around the sampling point ({circumflex over (τ)}l,{circumflex over (ν)}k)=(/W,k/T) in the delay-Doppler domain, as illustrated in FIG. 6. In essence, the virtual representation of equation (14) effectively approximates a discrete path doubly-selective channel in terms of an N-dimensional parameter comprised of the channel coefficients hl,k, where N=L×(2K+1)=(┌Wτmax┐+1)·(2┌Tνmax/2┐+1)≈τmaxνmaxN0.

With reference to FIG. 6, a plurality of resolution bins 600 having a resolution of 1/W in the delay domain space, and 1/T in the Doppler spread domain space are shown for a time window [0, τmax]. The plurality of dots 602 indicate gains received from different paths at a delay τ and Doppler shift ν.
Based on the interspacings between the multipath clusters within the delay-Doppler domain, doubly-selective channels can be characterized as either “rich” or “sparse”. In a rich multipath channel, the interspacings are smaller than Δτ=1/W in delay and Δν=1/T in Doppler shift. Sparse multipath channels, on the other hand, exhibit interspacings that are larger than Δτ and/or Δν, for example, as shown in FIG. 6, where not every resolution bin contains a multipath component or “dot” indicating a multipath component. In particular, since a channel coefficient consists of the sum of the gains of all paths falling within its respective delay-Doppler resolution bin, sparse doubly-selective channels tend to have far fewer than N nonzero channel coefficients at any fixed, but large enough signaling duration and/or bandwidth. Thus, the doubly selective channel may be considered sparse if D<<N where D=|{(k):hl,k>0}| and N=L×(2K+1)≈τmaxνmaxN0 is the total number of resolvable delays and Doppler shifts (channel coefficients) within the delay-Doppler spread.
As a first illustrative training signal used in a doubly-selective channel, a single-carrier spread spectrum waveform corresponding to a particular spreading code with binary phase shift keying as the modulation scheme can be used. The resulting training waveform x(t) can be represented as




  
    
      
        
          
            x
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              ∑
              
                n
                =
                0
              
              
                
                  N
                  0
                
                -
                1
              
            
            ⁢
            
              
                x
                n
              
              ⁢
              
                
                  I
                  
                    [
                    
                      0
                      ,
                      
                        T
                        c
                      
                    
                    )
                  
                
                ⁡
                
                  (
                  
                    t
                    -
                    
                      nT
                      C
                    
                  
                  )
                
              
            
          
        
        ,
        
          0
          ≤
          t
          <
          T
        
      
    
    
      
        (
        16
        )
      
    
  




where I[0,Tc)(t) is the chip waveform, Tc≈1/W is the chip duration, and {xnε is the spreading code corresponding to the training signal waveform. The output of the channel corresponding to x(t) is given by





  
    
      
        
          y
          ⁡
          
            (
            t
            )
          
        
        ≈
        
          
            
              ∑
              
                l
                =
                0
              
              
                L
                -
                1
              
            
            ⁢
            
              
                ∑
                
                  k
                  =
                  
                    -
                    K
                  
                
                K
              
              ⁢
              
                
                  h
                  
                    l
                    ,
                    k
                  
                
                ⁢
                
                  e
                  
                    j2π
                    ⁢
                    
                      k
                      T
                    
                    ⁢
                    t
                  
                
                ⁢
                
                  x
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        l
                        /
                        W
                      
                    
                    )
                  
                
              
            
          
          +
          
            
              z
              ⁡
              
                (
                t
                )
              
            
            ⁢
            
              
            
            ⁢
            0
          
        
        ≤
        t
        <
        
          T
          +
          
            τ
            max
          
        
      
    
    
      
        (
        17
        )
      
    
  




where z(t) is a zero-mean, circularly symmetric, complex AWGN waveform. For spread spectrum waveforms, chip-rate sampling of y(t) at the receiver yields an equivalent discrete-time representation





  
    
      
        
          
            y
            n
          
          =
          
            
              
                ∑
                
                  l
                  =
                  0
                
                
                  L
                  -
                  1
                
              
              ⁢
              
                
                  ∑
                  
                    k
                    =
                    
                      -
                      K
                    
                  
                  K
                
                ⁢
                
                  
                    h
                    
                      l
                      ,
                      k
                    
                  
                  ⁢
                  
                    e
                    
                      j2π
                      ⁢
                      
                        k
                        
                          N
                          0
                        
                      
                      ⁢
                      n
                    
                  
                  ⁢
                  
                    x
                    
                      n
                      -
                      l
                    
                  
                
              
            
            +
            
              z
              n
            
          
        
        ,
        
          n
          =
          0
        
        ,
        1
        ,
        ...
        ⁢
        
          
        
        ,
        
          
            N
            0
          
          +
          L
          -
          2
        
      
    
    
      
        (
        18
        )
      
    
  




where {zn} corresponds to a zero-mean, circularly symmetric, complex AWGN sequence and N0≈TW is the dimension of the receiver signal space.

If Ñ0=N0+L−1 and an Ñ0-length sequence of vectors {xnεCL} is defined comprising the spreading code {xn} as xn=[xn xn−1 . . . xn−(L−1)] for n=0, 1, . . . , Ñ0−1 where the notational understanding is that xi=0 for i∉{0, 1, . . . N0−1}, and if




  
    
      
        H
        =
        
          [
          
            
              
                
                  h
                  
                    0
                    ,
                    
                      -
                      K
                    
                  
                
              
              
                
                  h
                  
                    0
                    ,
                    
                      
                        -
                        K
                      
                      +
                      1
                    
                  
                
              
              
                ...
              
              
                
                  h
                  
                    0
                    ,
                    K
                  
                
              
            
            
              
                
                  h
                  
                    1
                    ,
                    
                      -
                      K
                    
                  
                
              
              
                
                  h
                  
                    1
                    ,
                    
                      
                        -
                        K
                      
                      +
                      1
                    
                  
                
              
              
                ...
              
              
                
                  h
                  
                    1
                    ,
                    K
                  
                
              
            
            
              
                ⋮
              
              
                ⋮
              
              
                ...
              
              
                ⋮
              
            
            
              
                
                  h
                  
                    
                      L
                      -
                      1
                    
                    ,
                    
                      -
                      K
                    
                  
                
              
              
                
                  h
                  
                    
                      L
                      -
                      1
                    
                    ,
                    
                      
                        -
                        K
                      
                      +
                      1
                    
                  
                
              
              
                ...
              
              
                
                  h
                  
                    
                      L
                      -
                      1
                    
                    ,
                    K
                  
                
              
            
          
          ]
        
      
    
    
      
        (
        19
        )
      
    
  




is the L×(2K+1) matrix of channel coefficients, and {unεC2K+1} is an Ñ0-length sequence of phase vectors given by

un=[ωN0KnωN0(K−1)n . . . ωN0−Kn]T  (20)

where





  
    ω
    
      N
      0
    
  
  =
  
    e
    
      
        -
        j
      
      ⁢
      
        
          2
          ⁢
          π
        
        
          N
          o
        
      
    
  




and n=0, 1, . . . , Ñ0−1, the sequence {yn} in equation (18) can be written as





  
    
      
        
          
            
              
                y
                n
              
              =
              
              ⁢
              
                
                  
                    x
                    n
                    T
                  
                  ⁢
                  
                    Hu
                    n
                  
                
                +
                
                  z
                  n
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    (
                    
                      
                        u
                        n
                        T
                      
                      ⊗
                      
                        x
                        n
                        T
                      
                    
                    )
                  
                  ⁢
                  
                    vec
                    ⁡
                    
                      (
                      H
                      )
                    
                  
                
                +
                
                  z
                  n
                
              
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  
                    
                      (
                      
                        
                          u
                          n
                          T
                        
                        ⊗
                        
                          x
                          n
                          T
                        
                      
                      )
                    
                    ⁢
                    h
                  
                  +
                  
                    z
                    n
                  
                
              
              ,
              
                n
                =
                0
              
              ,
              1
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  
                    N
                    ~
                  
                  o
                
                -
                1
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




If h=vec(H)εCN is the vector of channel coefficients, and the yn's are stacked into an Ñ0 dimensional vector y, the following system of equations results

y=Xh+z  (22)

where the Ñ0×N sensing matrix X is comprised of {unTxnT} as its rows: X=└u0x0 . . . uÑ0−1xÑ0−1┘T. In an illustrative embodiment, the spreading code {xn} is generated from a Rademacher distribution, i.e., the xn's independently take values of √{square root over (ε/N0)} or −√{square root over (ε/N0)} with probability of 0.5 each where ε is a transmit energy budget for the training signal.

The channel model coefficients can be determined as a solution to the selector




  
    
      
        
          h
          ^
        
        =
        
          
            
              argmin
              
                
                  h
                  ~
                
                ∈
                
                  C
                  N
                
              
            
            ⁢
            
              
                
                
                  h
                  ~
                
                
              
              
                l
                ⁢
                
                  
                
                ⁢
                1
              
            
            ⁢
            
              
            
            ⁢
            subject
            ⁢
            
              
            
            ⁢
            to
            ⁢
            
              
            
            ⁢
            
              
                
                
                  
                    X
                    H
                  
                  ⁢
                  r
                
                
              
              
                l
                ∞
              
            
          
          ≤
          λ
        
      
    
    
      
        (
        23
        )
      
    
  




where ĥ are the estimated channel coefficients, ∥∥l1 is an  norm, ∥∥l∞ is an l∞ norm, N=L×(2K+1)=(┌Wτmax┐+1)·(2┌Tνmax/2┐+1)≈τmaxνmaxN0 a maximum number of non-zero channel coefficients, CN is a set of non-zero channel coefficients, X is the sensing matrix X=└u0x0 . . . uÑ0−1xÑ0−1┘T, r is the Ñ0 dimensional vector of residuals: r=y−X{tilde over (h)}, y is the sampled multipath signal, and λ is a tuning parameter selected to control a number of channel coefficients to be determined. In an illustrative embodiment, λ(N,ε)=√{square root over (2ε(1+a)log N)} where εis a transmit energy budget for the training signal and any a≧0. In another illustrative embodiment, a lasso estimator may be used to determine the channel model coefficients. In a non-sparse environment, the same method can be used to determine the first S=const·√{square root over (N0/log p)} largest channel coefficients.

As a second illustrative training signal used in a doubly-selective channel, a multi-carrier signal comprising a complete orthogonal short-time Fourier (STF) basis for the N0-dimensional signal space is generated via time and frequency shifts of a fixed prototype pulse g(t):γl,k(t)=g(t−lTo)ej2πkWot,(l,k)εS={0, 1, . . . , Nt−1}×{0,1, . . . , Nf−1} where Nt=T/T0 and Nf=W/W0. In an illustrative embodiment, the prototype pulse is assumed to have unit energy, ∫|g(t)|2dt=1, and completeness of {γl,k} stems from the underlying assumption that T0W0=1, which results in a total of NtNf=TW/T0W0=N0 basis elements. The training waveform may randomly dedicate Nr of the N0 STF basis elements as “pilot tones” such that




  
    
      
        
          
            x
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                ɛ
                
                  N
                  r
                
              
            
            ⁢
            
              
                ∑
                
                  
                    (
                    
                      n
                      ,
                      m
                    
                    )
                  
                  ∈
                  
                    S
                    r
                  
                
              
              ⁢
              
                
                  γ
                  
                    n
                    ,
                    m
                  
                
                ⁡
                
                  (
                  t
                  )
                
              
            
          
        
        ,
        
          0
          ≤
          t
          <
          T
        
      
    
    
      
        (
        24
        )
      
    
  




where the set of indices of pilot tones, Sr consists of Nr elements randomly selected from S and ε is the transmit energy budget available for training purposes. At the receiver, assuming that the basis parameters T0 and W0 are matched to the channel parameters τmax and νmax so that γl,k's serve as approximate eigenfunctions for sufficiently underspread channels, projecting the (noisy) received signal y(t) onto the STF basis waveforms yields





  
    
      
        
          y
          
            n
            ,
            m
          
        
        =
        
          
            
              〈
              
                y
                ,
                
                  γ
                  
                    n
                    ,
                    m
                  
                
              
              〉
            
            ≈
            
              
                
                  
                    ɛ
                    
                      N
                      r
                    
                  
                
                ⁢
                
                  H
                  
                    n
                    ,
                    m
                  
                
              
              +
              
                
                  z
                  
                    n
                    ,
                    m
                  
                
                ⁢
                
                  
                
                ⁢
                
                  (
                  
                    n
                    ,
                    m
                  
                  )
                
              
            
          
          ∈
          
            S
            r
          
        
      
    
    
      
        (
        25
        )
      
    
  




where y, γn,m=∫y(t)γn,m(t)dt, {zn,m} corresponds to an AWGN sequence and the STF channel coefficients are given by Hn,m≈H(t,f)|(t,f)=(nT0,mW0).

The virtual representation of a doubly-selective channel implies that




  
    H
    ⁡
    
      (
      
        t
        ,
        f
      
      )
    
  
  ≈
  
    
      ∑
      
        l
        =
        0
      
      
        L
        -
        1
      
    
    ⁢
    
      
        ∑
        
          k
          =
          
            -
            K
          
        
        K
      
      ⁢
      
        
          h
          
            l
            ,
            k
          
        
        ⁢
        
          e
          
            j
            ⁢
            
              
            
            ⁢
            2
            ⁢
            
              
            
            ⁢
            π
            ⁢
            
              k
              T
            
            ⁢
            t
          
        
        ⁢
        
          
            e
            
              
                -
                j
              
              ⁢
              
                
              
              ⁢
              2
              ⁢
              
                
              
              ⁢
              π
              ⁢
              
                l
                W
              
              ⁢
              f
            
          
          .
        
      
    
  




As a result, the STF channel coefficients Hn,m can be written as





  
    
      
        
          
            
              
                H
                
                  n
                  ,
                  m
                
              
              =
              
              ⁢
              
                
                  ∑
                  
                    l
                    =
                    0
                  
                  
                    L
                    -
                    1
                  
                
                ⁢
                
                  
                    ∑
                    
                      k
                      =
                      
                        -
                        K
                      
                    
                    K
                  
                  ⁢
                  
                    
                      h
                      
                        l
                        ,
                        k
                      
                    
                    ⁢
                    
                      e
                      
                        j
                        ⁢
                        
                          
                        
                        ⁢
                        2
                        ⁢
                        
                          
                        
                        ⁢
                        π
                        ⁢
                        
                          k
                          
                            N
                            t
                          
                        
                        ⁢
                        n
                      
                    
                    ⁢
                    
                      e
                      
                        
                          -
                          j
                        
                        ⁢
                        
                          
                        
                        ⁢
                        2
                        ⁢
                        
                          
                        
                        ⁢
                        π
                        ⁢
                        
                          l
                          
                            N
                            f
                          
                        
                        ⁢
                        m
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  u
                  
                    f
                    ,
                    m
                  
                  ′
                
                ⁢
                
                  Hu
                  
                    t
                    ,
                    n
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  (
                  
                    
                      u
                      
                        t
                        ,
                        n
                      
                      ′
                    
                    ⊗
                    
                      u
                      
                        f
                        ,
                        m
                      
                      ′
                    
                  
                  )
                
                ⁢
                
                  vec
                  ⁡
                  
                    (
                    H
                    )
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  (
                  
                    
                      u
                      
                        t
                        ,
                        n
                      
                      ′
                    
                    ⊗
                    
                      u
                      
                        f
                        ,
                        m
                      
                      ′
                    
                  
                  )
                
                ⁢
                h
              
            
          
        
      
    
    
      
        (
        26
        )
      
    
  




where H is the L×(2K+1) matrix of channel coefficients as in equation (19), h=vec(H)εCN, uf,m=[1 ωNf1m . . . ωNf(L−1)m]TεCL, and ut,n=[ωNtKn ωNt(K−1)n . . . ωNt−Kn]TεC2K+1. Stacking the received training symbols {yn,m} into an Nr-dimensional vector y yields the following system of equations

y=Xh+z  (27)

where the Nr×N sensing matrix X is comprised of {√{square root over (ε/Nr)}(ut,nTuf,mT):(n,m)εSr} as its rows and the AWGN vector z is distributed as CN(0Nr, INr).

The channel model coefficients can be determined as a solution to the Dantzig selector as




  
    
      
        
          h
          ^
        
        =
        
          
            
              argmin
              
                
                  h
                  ~
                
                ∈
                
                  C
                  N
                
              
            
            ⁢
            
              
                
                
                  h
                  ~
                
                
              
              
                l
                1
              
            
            ⁢
            
              
            
            ⁢
            subject
            ⁢
            
              
            
            ⁢
            to
            ⁢
            
              
            
            ⁢
            
              
                
                
                  
                    X
                    H
                  
                  ⁢
                  r
                
                
              
              
                l
                ∞
              
            
          
          ≤
          λ
        
      
    
    
      
        (
        28
        )
      
    
  




where ĥ are the estimated channel coefficients, ∥∥l1 is an l1 norm, ∥∥l∞ is an l∞ norm, N=L×(2K+1)=(┌Wτmax┐+1)·(2┌Tνmax/2┐+1)≈τmaxνmaxN0 a maximum number of channel coefficients, CN is a N-dimensional complex space, X is the sensing matrix {√{square root over (ε/Nr)}(xt,nTxf,mT):(n,m)εSr}, r is the Nr dimensional vector of residuals: r=y−X{tilde over (h)}, y is the sampled multipath signal, and λ is a tuning parameter selected to control a number of channel coefficients to be determined. In an illustrative embodiment, λ(N,ε)=√{square root over (2ε(1+a)log N)} where ε is a transmit energy budget for the training signal and any a≧0. In a non-sparse environment, the same method can be used to determine the first S=const·Nr/log4N0 largest channel coefficients. In another illustrative embodiment, a lasso estimator also may be used to determine the channel model coefficients.

In a third illustrative embodiment, first transceiver device 102 and second transceiver device 104 include uniform linear arrays of NT transmit antennas and NR receive antennas forming a multi-antenna (MIMO) channel. In the absence of noise, the corresponding baseband transmitted and received signal can be related as




  
    
      
        
          
            y
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              ∫
              
                
                  -
                  W
                
                /
                2
              
              
                W
                /
                2
              
            
            ⁢
            
              
                H
                ⁡
                
                  (
                  
                    t
                    ,
                    f
                  
                  )
                
              
              ⁢
              
                X
                ⁡
                
                  (
                  f
                  )
                
              
              ⁢
              
                e
                
                  j
                  ⁢
                  
                    
                  
                  ⁢
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  π
                  ⁢
                  
                    
                  
                  ⁢
                  f
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
              ⁢
              
                
              
              ⁢
              
                d
                f
              
            
          
        
        ,
        
          0
          ≤
          t
          ≤
          T
        
      
    
    
      
        (
        29
        )
      
    
  




where W is the two-sided bandwidth, T is the packet duration, y(t) is the NR dimensional received signal, X(f) is the (element-wise) Fourier transform of the NT-dimensional transmitted signal x(t), and H(t,f) is the NR×NT time-varying frequency response matrix of the channel. A MIMO channel can be accurately modeled in terms of the multiple spatially distributed physical paths as





  
    
      
        
          H
          ⁡
          
            (
            
              t
              ,
              f
            
            )
          
        
        =
        
          
            ∑
            
              n
              =
              1
            
            
              N
              p
            
          
          ⁢
          
            
              β
              n
            
            ⁢
            
              
                a
                R
              
              ⁡
              
                (
                
                  θ
                  
                    R
                    ,
                    n
                  
                
                )
              
            
            ⁢
            
              
                a
                T
                H
              
              ⁡
              
                (
                
                  θ
                  
                    T
                    ,
                    n
                  
                
                )
              
            
            ⁢
            
              e
              
                j
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                
                  
                
                ⁢
                π
                ⁢
                
                  
                
                ⁢
                
                  v
                  n
                
                ⁢
                t
              
            
            ⁢
            
              e
              
                
                  -
                  j
                
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                
                  
                
                ⁢
                π
                ⁢
                
                  
                
                ⁢
                
                  t
                  n
                
                ⁢
                f
              
            
          
        
      
    
    
      
        (
        30
        )
      
    
  




which represents signal propagation over Np paths. βn denotes the complex path gain, θR,n the angle of arrival (AoA) at the receiver, θT,n the angle of departure (AoD) at the transmitter, τn the relative delay, and νn the Doppler shift associated with the n-th path. The NT×1 vector aT(θT) and the NR×1 vector aR(θR) denote the array steering and response vectors, respectively, for transmitting/receiving a signal in the direction θT, θR and are periodic in θ with unit period where θ=d sin(φ)/λ, d is the antenna spacing, φ is the physical angle measured with respect to the array broadside, and λ is the wavelength of propagation.

While the physical model of equation (30) is highly accurate, it is difficult to analyze and learn due to its nonlinear dependence on a potentially large number of physical parameters. However, because of the finite transmit and receive array apertures, signaling duration, and bandwidth, the physical model can be well-approximated by a channel model with the aid of a Fourier series expansion. The channel model provides a low-dimensional approximation of equation (30) by uniformly sampling the multipath environment in the angle-delay-Doppler domain at a resolution commensurate with the signal space parameters: (ΔθR, ΔθT, Δτ, Δν)=(1/NR, 1/NT, 1/W, 1/T) That is,




  
    
      
        
          H
          (
          
            
          
          ⁢
          
            t
            ,
            
              
            
            ⁢
            f
          
          )
        
        ≈
        
          
        
        ⁢
        
          
            ∑
            
              i
              =
              1
            
            
              N
              R
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              
                N
                T
              
            
            ⁢
            
              
                ∑
                
                  l
                  =
                  0
                
                
                  L
                  -
                  1
                
              
              ⁢
              
                
                  ∑
                  
                    m
                    =
                    
                      -
                      M
                    
                  
                  M
                
                ⁢
                
                  
                    
                      H
                      v
                    
                    (
                    
                      
                    
                    ⁢
                    
                      i
                      ,
                      
                        
                      
                      ⁢
                      k
                      ,
                      
                        
                      
                      ⁢
                      l
                      ,
                      
                        
                      
                      ⁢
                      m
                    
                    )
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      a
                      R
                    
                    (
                    
                      
                    
                    ⁢
                    
                      i
                      
                        N
                        R
                      
                    
                    )
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      a
                      T
                      H
                    
                    (
                    
                      
                    
                    ⁢
                    
                      k
                      
                        N
                        T
                      
                    
                    )
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    e
                    
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      π
                      ⁢
                      
                        m
                        T
                      
                      ⁢
                      t
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    e
                    
                      j
                      ⁢
                      
                        
                      
                      ⁢
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      π
                      ⁢
                      
                        l
                        W
                      
                      ⁢
                      f
                    
                  
                
              
            
          
        
      
    
    
      
        (
        31
        )
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            
              H
              v
            
            ⁡
            
              (
              
                i
                ,
                k
                ,
                l
                ,
                m
              
              )
            
          
          ≈
          
            
              ∑
              
                n
                ∈
                
                  
                    S
                    
                      R
                      ,
                      i
                    
                  
                  ⋂
                  
                    S
                    
                      T
                      ,
                      k
                    
                  
                  ⋂
                  
                    S
                    
                      τ
                      ,
                      l
                    
                  
                  ⋂
                  
                    S
                    
                      v
                      ,
                      m
                    
                  
                
              
            
            ⁢
            
              β
              n
            
          
        
      
    
    
      
        (
        32
        )
      
    
  




where a phase and an attenuation factor have been absorbed in the βn's in equation (32). In equation (31), NR, NT, L=┌Wτmax┐+1 and M=┌Tνmax/2┐ denote the maximum number of resolvable AoAs, AoDs, delays, and one-sided Doppler shifts within the channel angle-delay-Doppler spread, respectively. Due to the fixed angle-delay-Doppler sampling of equation (30), which defines the fixed basis functions in equation (31), the virtual representation is a linear channel representation completely characterized by the channel coefficients Hν(i,k,m). Additionally, equation (32) states that each Hν(i,k,m) is approximately equal to the sum of the complex gains of all physical paths whose angles, delays, and Doppler shifts lie within an angle-delay-Doppler resolution bin of size ΔθR×ΔθT×Δτ×Δν centered around the virtual sample point ({circumflex over (θ)}R,i, {circumflex over (θ)}T,k, {circumflex over (τ)}l, {circumflex over (ν)}m)=(i/NR, k/NT, l/W, m/T) in the angle-delay-Doppler domain as shown with reference to FIG. 7. The virtual representation of equation (30) effectively approximates a physical MIMO channel in terms of a Dmax-dimensional parameter comprising the channel coefficients Hν(i,k,l,m), where Dmax=NRNTL(2M+1).

Based on the interspacings between the multipath clusters within the angle-delay-Doppler domain, MIMO channels can be characterized as either “rich” or “sparse”. In a rich MIMO channel, the interspacings are smaller than ΔθR,ΔθT,Δτ,Δν). Sparse MIMO channels, on the other hand, exhibit interspacings that are larger than ΔθR, ΔθT, Δτ, and/or Δν such that not every angle-delay-Doppler bin of size ΔθR×ΔθT×Δτ×Δν contains a physical path as shown with reference to FIG. 7. In particular, since a channel coefficient consists of the sum of the gains of all paths falling within its respective angle-delay-Doppler resolution bin, sparse MIMO channels tend to have far fewer than Dmax nonzero channel coefficients at any fixed. but large enough number of antennas, signaling duration, and bandwidth. A MIMO channel is D-sparse if D=|SD|<<Dmax, where Dmax=NRNTL(2M+1) is the total number of resolvable paths (channel coefficients) within the angle-delay-Doppler spread of the channel and SD={(i,k,m):|Hν(i,k,l,m)|>0}.
If the signaling parameters are chosen so that the channel is non-selective




  
    
      
        H
        =
        
          
            
              ∑
              n
            
            ⁢
            
              
                β
                n
              
              ⁢
              
                
                  a
                  R
                
                ⁡
                
                  (
                  
                    θ
                    
                      R
                      ,
                      n
                    
                  
                  )
                
              
              ⁢
              
                
                  a
                  T
                  H
                
                ⁡
                
                  (
                  
                    θ
                    
                      T
                      ,
                      n
                    
                  
                  )
                
              
            
          
          ≈
          
            
              A
              R
            
            ⁢
            
              H
              v
            
            ⁢
            
              A
              T
              H
            
          
        
      
    
    
      
        (
        33
        )
      
    
  




where AR and AT are NR×NR and NT×NT unitary discrete Fourier transform (DFT) matrices, respectively. The NR×NT beamspace matrix Hν couples the virtual AoAs and AoDs, and its entries are given by the Dmax=NRNT channel coefficients {Hν(i,k)}. The channel is assumed to be D-sparse in the angular domain (D=|SD|<<Dmax) and it remains constant over the packet signaling duration T with a blockfading assumption corresponding to Tνmax<<1.

To learn the NR X NT (antenna domain) matrix H, part of the packet duration T is used to transmit a training signal to the receiver. Stacking the Mtr=TtrW received vector-valued training signals {y(n),n=1, . . . , Mtr} into an Mtr×NR matrix Y yields the following system of equations




  
    
      
        Y
        =
        
          
            
              
                ɛ
                
                  M
                  tr
                
              
            
            ⁢
            
              XH
              T
            
          
          +
          W
        
      
    
    
      
        (
        34
        )
      
    
  




where ε is the total transmit energy budget available for training, X is the collection of Mtr training signal vectors {x(n), n=1, . . . , Mtr} stacked row-wise into an Mtr×NR matrix with the constraint that ∥X∥F2=Mtr, and W is an Mtr×NR matrix of unit-variance AWGN.

Given receiver processing in the antenna domain, let χ={1, . . . , NR}×{1, . . . , NT} be the set of indices of elements within the antenna domain matrix H, and let the number of receive signal space dimensions dedicated to training be Ntr≧c2·log5Dmax·D for some constant C2>0, and choose χtr to be a set of Ntr ordered pairs sampled uniformly at random from χ. The number of transmitted training signals may be defined as Mtr=|{k:(i,k)εχtr}| resulting in ytr=√{square root over (ε/Mtr)}{H(i,k)}χtr+w where ytr is the vector of received training data and w ̃CN(0Ntr,INtr).
If vec−1 denotes the inverse of the vec operator defined as stacking the columns of a matrix into a vector, and λ=(2{tilde over (ε)}(1+a)log Dmax)1/2 for any a≧0 where {tilde over (ε)}=εNtr/DmaxMtr, aT,iT and aR,iT denote the i-th row of AT* and AR, respectively, and considering an Ntr×Dmax sensing matrix Utr that is comprised of {(aT,kTaR,iT):(i,k)εχtr} as its rows, the estimate of the channel coefficients is obtained from ytr as follows:

Ĥ=AR(vec−1(DS(√{square root over (ε/Mtr)}Utr,ytr,λ)))ATH  (35)

The virtual channel model coefficients in equation (35) can be determined as a solution to the Dantzig selector, DS, defined as




  
    
      
        
          h
          ^
        
        =
        
          
            
              argmin
              
                
                  h
                  ~
                
                ∈
                
                  C
                  
                    D
                    max
                  
                
              
            
            ⁢
            
              
                
                
                  h
                  ~
                
                
              
              
                l
                1
              
            
            ⁢
            
              
            
            ⁢
            subject
            ⁢
            
              
            
            ⁢
            to
            ⁢
            
              
            
            ⁢
            
              
                
                
                  
                    X
                    H
                  
                  ⁢
                  r
                
                
              
              
                l
                ∞
              
            
          
          ≤
          λ
        
      
    
    
      
        (
        36
        )
      
    
  




where ĥ are the estimated channel coefficients, ∥∥l1 is an l1 norm, ∥∥l∞ is an ∞ norm, Dmax=NRNT a maximum number of non-zero channel coefficients, CDmax is a Dmax-dimensional complex space, X is the sensing matrix √{square root over (ε/Mtr)}{(aT,kTaR,iT):(i,k)εχtr}, r is the D dimensional vector of residuals: r=(ytr−Xtrĥ), y is the sampled multipath signal, and λ is a tuning parameter selected to control a number of channel coefficients to be determined such as λ=(2{tilde over (ε)}(1+a)log Dmax)1/2. In a non-sparse environment, the same method can be used to determine the first D largest channel coefficients. In another illustrative embodiment, a lasso estimator may be used to determine the channel model coefficients.

Given receiver processing in the beamspace domain, let hν,iT denote the i-th row of Hνand define Di to be the number of nonzero virtual coefficients in each row of Hν; that is, Di=∥hν,i∥0. The number of training vectors may be chosen as Mtr≧c6·log(NT/maxiDi)·maxiDi for some constant c6>0 and each training vector x(n),n=1 . . . , Mtr is an iid vector of binary random variables taking values +1/√{square root over (NT)} or −1/√{square root over (NT)} with probability 0.5 each. At the receiver, the received training signals {y(n), n=, . . . , Mtr} are stacked into an Mtr×NR matrix Y to yield equation (34).
To determine the channel coefficients in an illustrative embodiment, define Yν=YAR* and let yν,i denote the i-th column of Yν, select a λ=(2ε(1+a)(log Dmax)/NT)1/2 for an a≧0, and define ĥν,i=DS(√{square root over (ε/Mtr)}XAT*,yν,i,λ) for i=1, . . . , NR where DS denotes the Dantzig selector. The estimate of H is obtained as:

Ĥ=AR[ĥν,1 . . . ĥν,NR]TATH  (36)

In a non-sparse environment, the same method can be used to determine the first const·Mtr/log NT largest nonzero channel coefficients per row of channel matrix Hν. In another illustrative embodiment, a lasso estimator may be used to determine the channel model coefficients.

If the signaling parameters are chosen so that the MIMO channel is frequency-selective such that Wτmax≧1, limited to blockfading corresponding to Tνmax<<1, and assuming that the communication packet is comprised of N0≈T/(Tf+τmax)≧NT OFDM vector-valued symbols with each OFDM symbol consisting of Q=TfW≧┌Wτmax┐+1 tones where Tf<T denotes the OFDM symbol duration, the channel model of equation (31) reduces to




  
    
      
        
          H
          ⁡
          
            (
            f
            )
          
        
        ≈
        
          
            ∑
            
              l
              =
              0
            
            
              L
              -
              1
            
          
          ⁢
          
            
              
                H
                ~
              
              ⁡
              
                (
                l
                )
              
            
            ⁢
            
              e
              
                
                  -
                  j
                
                ⁢
                
                  
                
                ⁢
                2
                ⁢
                
                  
                
                ⁢
                π
                ⁢
                
                  l
                  W
                
                ⁢
                f
              
            
          
        
      
    
    
      
        (
        37
        )
      
    
  




where the antenna domain matrices {tilde over (H)}'s are defined in terms of the beamspace matrices: {tilde over (H)}≈ARHVATH. The channel frequency response H(f) is characterized by Dmax=NRNTL channel coefficients {Hν(i,k,}, out of which only D<<Dmax coefficients are assumed to be nonzero. To learn this Dmax-dimensional channel, Mtr of N0Q OFDM tones may be selected as “pilot tones” to transmit known vector-valued training signals to the receiver over the selected tones. The transmitted and received training signals in this case are related to each other as





  
    
      
        
          
            y
            ⁡
            
              (
              
                n
                ,
                q
              
              )
            
          
          =
          
            
              
                
                  ɛ
                  
                    M
                    tr
                  
                
              
              ⁢
              
                H
                ⁡
                
                  (
                  q
                  )
                
              
              ⁢
              
                x
                ⁡
                
                  (
                  
                    n
                    ,
                    q
                  
                  )
                
              
            
            +
            
              w
              ⁡
              
                (
                
                  n
                  ,
                  q
                
                )
              
            
          
        
        ,
        
          
            (
            
              n
              ,
              q
            
            )
          
          ∈
          
            P
            tr
          
        
      
    
    
      
        (
        38
        )
      
    
  




where ε is the transmit energy budget available for training, the matrix H(q)=H(f)|f=q/T∫,w(n,q) ̃CN(0NR,INR),Ptr⊂{1, . . . , N0}×{0, . . . , Q−1}:|Ptr|=Mtr, and the set of transmit training vectors x(n,q) is designed such that ΣPtr∥x(n,q)∥22=Mtr.

For an illustrative training signal used in a MIMO channel, let hν,iT denote the i-th row of Hν and define Hν,i=[hν,i(0) . . . hν,i(L−1)]. Further, define Di to be the number of nonzero virtual coefficients in




  
    H
    
      v
      ,
      i
    
  
  (
  
    
      D
      i
    
    =
    
      
        ∑
        
          l
          =
          0
        
        
          L
          -
          1
        
      
      ⁢
      
        
          
          
            
              h
              
                v
                ,
                i
              
            
            ⁡
            
              (
              l
              )
            
          
          
        
        0
      
    
  
  )




and let the number of pilot tones be Mtr≧c9·log5(NTQ)maxiDi for some constant c9>0. Define the pilot tones Ptr to be a set of Mtr ordered pairs sampled uniformly at random from {1, . . . , NT}×{0, . . . , Q−1}, and the corresponding training vectors as {x(n,q)=en:(n,q)εPtr}. At the receiver, the received signals {y(n,q)}Ptr are stacked row-wise to yield an Mtr×NR matrix Y consisting of y(n,q)'s as its rows.

To determine the channel coefficients, define an Mtr×NTL matrix Utr that is comprised of {(uf,qTaT,nT):(n,q)εPtr} as its rows where,




  
    u
    
      f
      ,
      q
    
    T
  
  =
  
    [
    
      
        e
        
          
            -
            j
          
          ⁢
          
            
          
          ⁢
          2
          ⁢
          
            
          
          ⁢
          π
          ⁢
          
            q
            Q
          
          ⁢
          0
        
      
      ⁢
      
        
      
      ⁢
      ...
      ⁢
      
        
      
      ⁢
      
        e
        
          
            -
            j
          
          ⁢
          
            
          
          ⁢
          2
          ⁢
          
            
          
          ⁢
          π
          ⁢
          
            q
            Q
          
          ⁢
          
            (
            
              L
              -
              1
            
            )
          
        
      
    
    ]
  




and aT,nT denotes the n-th row of AT*, choose λ=(2ε(1+a)(log Dmax)/NT)1/2 for any a≧0, define Yν=YAR*, and let yν,i denote the i-th column of Yν. The estimate of Hν,i is given as follows:

Ĥν,i=vec−1(DS(√{square root over (ε/Mtr)}Utr,yν,i,λ))  (39)

where DS denotes the Dantzig selector. In a non-sparse environment, the same method can be used to determine the first const·Mtr/log5(NTQ) largest nonzero channel coefficients for Hν,i corresponding to an output of each receive antenna. In another illustrative embodiment, a lasso estimator may be used to determine the channel model coefficients.

In another exemplary embodiment, a training signal comprising one OFDM symbol may be used. The described processes can further be applied in a MIMO, doubly selective signal space that employs STF signaling.
The word “illustrative” is used herein to mean serving as an example, instance, or illustration. Any aspect or design described herein as “illustrative” is not necessarily to be construed as preferred or advantageous over other aspects or designs. Further, for the purposes of this disclosure and unless otherwise specified, “a” or “an” means “one or more”. Further, for the purposes of this disclosure, “and” or “or” means “and/or”. The illustrative embodiments may be implemented as a method, machine, or article of manufacture using standard programming and/or engineering techniques to produce software, firmware, hardware, and/or any combination thereof to control a device to implement the disclosed embodiments.
The foregoing description of illustrative embodiments of the invention have been presented for purposes of illustration and of description. It is not intended to be exhaustive or to limit the invention to the precise form disclosed, and modifications and variations are possible in light of the above teachings or may be acquired from practice of the invention. The embodiments were chosen and described in order to explain the principles of the invention and as practical applications of the invention to enable one skilled in the art to utilize the invention in various embodiments and with various modifications as suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims appended hereto and their equivalents.

## Claims

1. A device comprising:
an antenna configured to receive a multipath signal; and
a processor operably coupled to the antenna to receive the multipath signal, the processor configured to
determine a training signal transmitted from a second device to create the received multipath signal;
sample the received multipath signal; and
determine channel coefficients based on the sampled multipath signal and the determined training signal using a mixed-norm convex optimization process, wherein the channel coefficients characterize a channel associated with the multipath signal.

2. The device of claim 1, wherein the device is the second device.

3. The device of claim 2, wherein the device is a sensor.

4. The device of claim 2, wherein the antenna is further configured to transmit the training signal.

5. The device of claim 1, wherein the antenna is further configured to transmit the determined channel coefficients to the second device.

6. A method of determining channel coefficients characterizing a channel associated with a multipath signal, the method comprising:
receiving a multipath signal at a first device;
using a processor at the first device, determining a training signal transmitted from a second device to create the received multipath signal;
sampling the received multipath signal; and
determining channel coefficients based on the sampled multipath signal and the determined training signal using a mixed-norm convex optimization process, wherein the channel coefficients characterize a channel associated with the multipath signal.

7. The method of claim 6, further comprising determining a resolution value based on a signaling characteristic of the determined training signal.

8. The method of claim 7, wherein the received multipath signal is sampled uniformly at the determined resolution value.

9. The method of claim 7, wherein the channel coefficients are determined using a discrete approximation of a continuous time channel based on the determined resolution value.

10. The method of claim 6, wherein the mixed-norm convex optimization process comprises a Dantzig selector.

11. The method of claim 10, wherein the mixed-norm convex optimization




  
    β
    ^
  
  =
  
    
      argmin
      
        
          β
          ~
        
        ∈
        
          C
          P
        
      
    
    ⁢
    
      
        
        
          β
          ~
        
        
      
      
        l
        1
      
    
  




subject to ∥XHr∥l∞≦λ where {circumflex over (β)} are the determined channel coefficients, ∥∥l1 is an l1 norm, ∥∥l∞ is an l∞ norm, p is a maximum number of channel coefficients, Cp is a p-dimensional complex space, X is a sensing matrix, r=y−X{tilde over (β)} is a vector of residuals, y is the sampled multipath signal, and λ is a tuning parameter selected to control a number of channel coefficients determined.

12. The method of claim 11, wherein X comprises an (n+p−1)×p Toeplitz-structured, convolution matrix formed from the training signal, n is ┌TSW┐ where TS is a signaling duration of the training signal and W is a bandwidth of the training signal.

13. The method of claim 12, wherein the training signal comprises a pseudo-random sequence randomly selected from a zero mean, unit variance distribution.

14. The method of claim 11, wherein X comprises



  
    
      X
      =
      
        
          ⌊
          
            
              
                u
                0
              
              ⊗
              
                x
                0
              
            
            ⁢
            
              
            
            ⁢
            ...
            ⁢
            
              
            
            ⁢
            
              
                u
                
                  
                    
                      N
                      ~
                    
                    0
                  
                  -
                  1
                
              
              ⊗
              
                x
                
                  
                    
                      N
                      ~
                    
                    0
                  
                  -
                  1
                
              
            
          
          ⌋
        
        T
      
    
    ⁢
    
      
    
  




  
    
      
        
          where
          ⁢
          
            
          
          ⁢
          
            u
            n
          
        
        =
        
          
            [
            
              
                ω
                
                  N
                  o
                
                Kn
              
              ⁢
              
                
              
              ⁢
              
                ω
                
                  N
                  o
                
                
                  
                    (
                    
                      K
                      -
                      1
                    
                    )
                  
                  ⁢
                  n
                
              
              ⁢
              
                
              
              ⁢
              ...
              ⁢
              
                
              
              ⁢
              
                ω
                
                  N
                  o
                
                
                  -
                  Kn
                
              
            
            ]
          
          T
        
      
      ,
      
        

      
      ⁢
      
        
          ω
          
            N
            0
          
        
        =
        
          e
          
            
              -
              j
            
            ⁢
            
              
                2
                ⁢
                π
              
              
                N
                o
              
            
          
        
      
    
    ⁢
    
      
    
  




  
    
      
        and
        ⁢
        
          
        
        ⁢
        n
      
      =
      0
    
    ,
    1
    ,
    ...
    ⁢
    
      
    
    ,
    
      
        
          N
          ~
        
        0
      
      -
      1
    
    ,
    
      
        
          N
          ~
        
        0
      
      =
      
        
          N
          0
        
        +
        L
        -
        1
      
    
    ,
    
      
        N
        0
      
      ≈
      TW
    
    ,
  



T is a packet duration of the training signal, W is a bandwidth of the training signal, L=┌Wτmax┐+1 τmax is the delay spread of the channel, K=┌Tvmax/2┐, and vmax is a maximum, one-sided Doppler shift caused by the channel.

15. The method of claim 14, wherein the training signal comprises a single-carrier spread spectrum waveform.

16. The method of claim 15, wherein a spreading code, {xn}, of the single-carrier spread spectrum waveform is generated from a Rademacher distribution, where the xn's independently take values of √{square root over (ε/N0)} or −√{square root over (ε/N0)} with probability of 0.5 each where ε is a transmit energy budget for the training signal.

17. The method of claim 6, wherein the sampled multipath signal is formed by projecting the received multipath signal onto a plurality of basis elements.

18. The method of claim 11, wherein the mixed-norm convex optimization process is performed for each receive antenna of a plurality of receive antennas.

19. The method of claim 6, wherein the mixed-norm convex optimization process comprises a Lasso estimator.

20. A non-transitory computer-readable medium having computer-readable instructions stored thereon which, when executed by a processor, cause a device to:
determine a training signal transmitted from a second device to create a multipath signal received at the device;
sample the received multipath signal; and
determine channel coefficients based on the sampled multipath signal and the determined training signal using a mixed-norm convex optimization process, wherein the channel coefficients characterize a channel associated with the multipath signal.