I. Field of the Disclosure
The technology of the disclosure relates generally to channel estimation that may be used, for example, for range finding.
II. Background
Computing devices abound in modern society, and more particularly, mobile communication devices have become increasingly common. The prevalence of these mobile communication devices is driven in part by the many functions that are now enabled on such devices. Increased processing capabilities in such devices means that mobile communication devices have evolved from pure communication tools into sophisticated mobile entertainment centers, thus enabling enhanced user experiences. Many services provided to mobile communication devices may depend on a location of the mobile communication device (i.e., so-called “location-based services”). Many techniques have been developed to find a distance between a transmitter and an object. The advent of ultra-wideband (UWB) communication has seen the development of real-time localization systems (RTLS), many of which are based on IEEE 802.15.4 (published by IEEE in 2015) and an amendment published in 2020, namely IEEE 802.15.4z. These standards leave open many implementation details and leave room for innovation.
Aspects disclosed in the detailed description include systems and methods for channel estimation. In a particular aspect, supercomplementary sets (SCSs) of sequences are defined and more optionally, but more particularly, supercomplementary zero-sum cross-correlation (SZC) sequence blocks are created from the SCSs and transmitted to a receiver. The receiver receives a signal that includes the SCSs or SZC sequence block(s) deformed by a channel frequency offset (CFO). The receiver may correlate the signal with the known SCS or SZC to create a signal with minimal or no side lobes from which the CFO may be removed and a channel estimate formed. Based on the channel estimate, a distance between a transmitter and the receiver may be determined and used for other purposes, such as the provision of location-based services.
In this regard, in one aspect, a source is disclosed. The source comprises an antenna. The source also comprises a transceiver coupled to the antenna. The source also comprises a control circuit coupled to the transceiver. The control circuit is configured to cause the transceiver to send a complementary low-sum cross-correlation (CLS) sequence block to a remote device for ranging.
In another aspect, a mobile device is disclosed. The mobile device comprises an antenna. The mobile device also comprises a receiver coupled to the antenna, the receiver comprising a correlator. The mobile device also comprises a control circuit coupled to the receiver. The control circuit is configured to cause the correlator to correlate a received signal with a CLS sequence block to derive a channel estimate.
In another aspect, a method of determining distance is disclosed. The method comprises transmitting a CLS sequence block from a source to a remote device. The method also comprises, at the remote device, correlating a received version of the CLS sequence block with a known CLS sequence block to derive a channel estimate. The method also comprises removing a carrier frequency offset from the channel estimate. The method also comprises sending information derived from the channel estimate from the remote device to the source.
The embodiments set forth below represent the necessary information to enable those skilled in the art to practice the embodiments and illustrate the best mode of practicing the embodiments. Upon reading the following description in light of the accompanying drawing figures, those skilled in the art will understand the concepts of the disclosure and will recognize applications of these concepts not particularly addressed herein. It should be understood that these concepts and applications fall within the scope of the disclosure and the accompanying claims.
It will be understood that although the terms first, second, etc. may be used herein to describe various elements, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first element could be termed a second element, and similarly, a second element could be termed a first element without departing from the scope of the present disclosure. As used herein, the term “and/or” includes any and all combinations of one or more of the associated listed items.
It will be understood that when an element such as a layer, region, or substrate is referred to as being “on” or extending “onto” another element, it can be directly on or extend directly onto the other element, or intervening elements may also be present. In contrast, when an element is referred to as being “directly on” or extending “directly onto” another element, no intervening elements are present. Likewise, it will be understood that when an element such as a layer, region, or substrate is referred to as being “over” or extending “over” another element, it can be directly over or extend directly over the other element, or intervening elements may also be present. In contrast, when an element is referred to as being “directly over” or extending “directly over” another element, no intervening elements are present. It will also be understood that when an element is referred to as being “connected” or “coupled” to another element, it can be directly connected or coupled to the other element, or intervening elements may be present. In contrast, when an element is referred to as being “directly connected” or “directly coupled” to another element, no intervening elements are present.
Relative terms such as “below” or “above” or “upper” or “lower” or “horizontal” or “vertical” may be used herein to describe a relationship of one element, layer, or region to another element, layer, or region as illustrated in the Figures. It will be understood that these terms and those discussed above are intended to encompass different orientations of the device in addition to the orientation depicted in the Figures.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure. As used herein, the singular forms “a,” “an,” and “the” are intended to include the plural forms as well unless the context clearly indicates otherwise. It will be further understood that the terms “comprises,” “comprising,” “includes,” and/or “including,” when used herein, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
Unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs. It will be further understood that terms used herein should be interpreted as having a meaning that is consistent with their meaning in the context of this specification and the relevant art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
Aspects in the detailed description include systems and methods for channel estimation. In a particular aspect, supercomplementary sets (SCSs) of sequences are defined and more optionally, but more particularly, supercomplementary zero-sum cross-correlation (SZC) sequence blocks are created from the SCSs and transmitted to a receiver. The receiver receives a signal that includes the SCSs or SZC sequence block(s) deformed by a channel frequency offset (CFO). The receiver may correlate the signal with the known SCS or SZC to create a signal with minimal or no side lobes from which the CFO may be removed and a channel estimate formed. Based on the channel estimate, a distance between a transmitter and the receiver may be determined and used for other purposes, such as the provision of location-based services.
Exemplary aspects of the present disclosure are well suited for real-time localization systems (RTLS) in impulse radio ultra-wideband (UWB) devices that use the IEEE 802.15.4 and 802.15.4z standards. These standards, and 802.15.4z in particular, require that the ranging signal from the source be transmitted without a guard interval. This requirement makes the use of certain channel impulse response (CIR) estimation processes non-compliant. Specifically, conventional CIRs may rely on scrambled timestamp sequences (STSs) consisting of aperiodic pseudorandom sequences located after a synchronization header packet. STS was introduced to create a CIR estimation field resilient to distance-reduction attacks. Conventional use of STS causes CFO deformation, which would normally be addressed by adding a guard interval longer than CIR. However, as noted, 802.15.4z prohibits such guard intervals.
Thus, exemplary aspects of the present disclosure contemplate the use of SZC formed from SCSs. This provides CIR estimation without any correlation artifacts with arbitrary sequence deformation consistent across a sequence block. This allows for transmission without a guard interval and thus allows for compliance with 802.15.4z.
In this regard,
The source 104 may further include a control circuit 112 coupled to a memory 114 and a transceiver 116. The source 104 may send and receive signals 118 through an antenna 120 using the transceiver 116. Operating software and/or signal symbols may be stored in the memory 114 as needed or desired.
The mobile device 106 may be any number of portable computing devices such as a smartphone, laptop, tablet, or the like and may include an antenna 122 through which signals may be received and sent.
As better illustrated in
The known symbol or sequence or block of sequences is shared with possible mobile devices (e.g., by being defined in a specification, defined by an application downloaded to the mobile device, being present in a widget that an application on the mobile device may access, or the like). The source 104 then transmits a symbol or sequence (block 304). During transmission, the signal path deforms the symbol or sequence in the signals 118 with a CFO (block 306).
The mobile device 106, and particularly the receiver 200, receives the deformed symbol or sequence (block 308). Note that the receiver 200 is considered to be accurate in frequency, so from the receiver's perspective, the transmission is deformed with the CFO. The received signal is processed (block 310), including optional steps of RF and BB processing (block 312) by the analog RF and BB processing circuit 206 and I/Q sampling (block 314) by the I/Q sampling circuit 208.
The correlator 210 then correlates the sampled signal with a known cross-correlation sequence (block 316) that was previously populated and stored in the mobile device, such as in the memory 204. This correlation will, because of the complementariness of the sequence, provide a low or no sideband signal from which removal of the CFO is readily accomplished (block 318). The signal is then accumulated (block 320) by the accumulator 214.
From the correlated signal, a channel estimate may be formed. From the channel estimate, a TOA and/or a distance may be calculated (block 322). In a first exemplary aspect, the channel estimate is calculated by the control circuit 202 and reported to the source 104, where the TOA and distance are determined. In a second exemplary aspect, the control circuit 202 also calculates the TOA and/or the distance and reports to the source 104. In still another exemplary aspect, the channel estimate is provided to a remote computing device (e.g., external device 108) which performs the calculations. From this calculation (wherever performed and to whatever reported), a decision about the provision of location-based services may be made.
While the hardware and basic process appear relatively straightforward, there is substantial work that must be done before the complementary sets used by the source 104 and the correlator 210 may be realized. The following section steps through the theory and the math that lies under the hardware and processes outlined above.
Accordingly, to assist in explaining the theory and math, a mathematical model of the receiver 200 with a single sample per preamble chip is used. However, it should be appreciated, in a real receiver, the sampling rate is higher than one sample per preamble chip; hence, a realistic model is up-sampled with regard to the model introduced here. Furthermore, a CFO also produces the effect of dilatation/expansion of the baseband signal in time since the same reference clock in the transmitter is used both for carrier and baseband signal generation. Therefore, baseband receiver processing also requires re-sampling. However, these considerations are omitted here for the sake of simplicity and without loss of generality.
As an initial foray into this discussion, a simplified model of the aperiodic preamble reception under 802.15.4 where the sampling rate is equal to the chip rate fp is provided. The preamble symbol duration is denoted Tsym=R/fp, where R is symbol duration in chips. The first part of the i-th preamble symbol for i≥0 consists of a chip sequence Si(n), non-zero for n=0, 1, . . . , N−1 transmitted at fp. The sequence si(n) goes through the channel with the impulse response (CIR) denoted h(n), non-zero for n=0, 1, . . . , M−1 with N+M−1≤R. Frequency offset between the transmitter and receiver is denoted Δf. The received sample sequence for i-th symbol is non-zero for 0≤n≤N+M−2:)
r
i(v,n)=ejvi
Here, “⊗” denotes aperiodic convolution, v=Δf/fp is the discrete-time frequency offset, and φ=2πvR is the phase offset between consecutive symbols, and {tilde over (s)}i(v,n)−si(n)e2πjvnis a version of si(n) deformed by CFO.
The correlator will correlate ri(v,n) with si(n):
c
i(v,n)=ri(v,n)⊗s*i(−n). (2)
Equation (2) can be rewritten as
c
i(v,n)=Xi,i(v,n)⊗h(n)ejvi. (3)
Here,
Xl,m(v,n)={tilde over (s)}l(v,n)⊗s*m(−n), (4)
represents cross-AF (CAF) of sl(n) and sm(n) at frequency v, non-zero for n=−N+1, . . . , 0, . . . , N−1. For m=1, Equation (4) represents an auto-ambiguity function (AF).
Frequency offset is removed after correlation and before accumulation by the carrier loop on a per-symbol basis. Therefore, the accumulator after k preamble symbols has a form
Note that since causality is omitted in correlation Equation (2), dk(v,n) has non-zero values for n=−N+1, . . . , 0, . . . , N+M−2. Furthermore, dk(v,n) represents an estimate of h(n) for n=0, 1, . . . ,M−1. Hence, the accumulator needs to contain only these indices.
From Equation (5), the quality of the channel estimate after k symbols depends on the shape of the mean AF:
that is, the level of its main lobe and its sidelobes.
Auto-AFs of an SCSS are complementary, and thus, the application of SCSSs in this scenario is straightforward. For k=I, where I is the size of the SCSS used
X
I(v,n)=δ(n)A(v), (7)
holds, where A(v) is the complex amplitude of the main lobe. Equations (5) and (7) yield
d
I(v,n)=IA(v)h(n), (8)
which represents a perfect channel estimation.
For unity-level (|si(n)|=1) sequences of length N from Equation (4) follows
Therefore, increasing CFO decreases the level of the channel estimation and rotates it but does not produce any correlation artifacts via sidelobes of
Consider a case in which R=N, that is, in which the preamble symbols are transmitted without any guard interval. The condition for the length of the channel is now M≤N. Due to the two adjacent symbols now overlapping with the current one in the channel, for I transmitted symbols, Equation (1) becomes
r
0(v,n)=[
r
i(v,n)=[
r
I−1(v,n)=[
for 0≤n≤N−1.
As noted previously, auto-AFs within an SCSS are complementary. However, this does not hold for the adjacent symbol CAFs in Equation (11), which generally do not complement. To overcome this problem, consider the following construct. Consider a case in which a block of I symbols is transmitted without a guard interval, as previously. However, now samples of 0-th symbol for n=−N+2, . . . , −1 are correlated with a sequence s−1(n) and accumulated into d(1), . . . , d(N−1). Similarly, samples of I−1-th symbol for n=N, . . . , 2N−2 are correlated with a sequence SI(n) and accumulated into d(0), . . . , d(N−2). Then, Equation (11) becomes
for 0≤n≤N−1.
Now, an SCSB gi(n) is constructed from si(n) of size I:
Sequences g−1(n) and g2I(n) are not transmitted but only are used in the receiver correlation, such as when deriving Equation (12); hence, the size of the transmitted sequence block is 2I.
From Equation (13), the block gi(n)
Xi,i+1(v, n)=−Xi+I,i+I+1(v,n), (14a)
Xi,i −1(v, n)=−Xi+I,i+I+1(v,n), (14b)
for i=0, . . . , I−1 holds. Equations (7), (12), and (14) yield
d
2I(v,n)=2IA(v)h(n), (15)
which represents a perfect channel estimation with the complex level reduction due to CFO, as in Equation (8).
If I is even, from Equation (14) follows that g−1(n)=g2I−1(n) and g2I(n)=g0(n). Therefore, SCSB gi(n) repeatedly transmitted L times without a guard interval produces d2I(v,n) (Equation 15) scaled L times, providing g−1(n) and g2I(n) are, as previously, used in the correlation of the precursors of 0-th transmitted symbol and the post-cursors of the last transmitted symbol, respectively.
Furthermore, the foregoing derivation holds if in {tilde over (s)}i(v,n), besides deformation by the frequency offset, other deformations, consistent from symbol to symbol, for example, video clipping, are included. Therefore, SCSB gi(n) is supercomplementary but now without a guard interval.
In a typical implementation, the device, according to the present disclosure, would have SCSS si(n) in the memory and generate SCSB gi(n) from si(n) and inversion logic based on the symbol counter−i, according to Equation (13).
When using an SCSB as previously described, the basic sequence block si(n) may be any SCSS constructed by any of a number of methods. For sequences of length N, minimal SCSS size is I=N. Therefore, si(n) and gi(n) minimally have N2 and 2N2 chips, respectively. In some scenarios, to minimize on-air time, it is useful to transmit only the si(n) block, that is, the first half of the gi(n) block without any guard interval. Therefore, it is worthwhile to construct an SCSS-based si(n) block with good correlation properties. This is done at v=0, as, generally, full complementing over deformations such as CFO cannot be expected. Correlation between Sl(n) and Sm(n) is denoted cl,m(n)=Xl,m(n,0). Correlation is done as in Equation (12) by choosing s−1(n)=sI−1(n) and sI(n)=s0(n). For v=0, Equation (12) becomes
for 0≤n≤N−1. With this choice, considering that from cm,l(n)=c*l,m(−n) holds from Equation (4), two sums of cross-correlations in Equation (16) are the following: Σi=0I−1ci,i−1(n−N) and Σi=0I−1ci,i+1(n+N) produce post-cursor and precursor cross-correlation sidelobes, respectively, that are conjugated mirror images of each other. Note that the same holds for the choice of S−1(n)=0 and sI(n)=0, which are not considered here.
Furthermore, since for SCSS the sidelobes of ci,i(n) fully complement, it is enough to observe only a single set of cross-correlation sidelobes:
Consider now that binary SCSS construction based on Walsh-Hadamard matrices can be easily expanded to polyphase or ternary SCSS designs. The search procedure of sidelobe minimization according to the present disclosure is as follows.
Construct normalized N×N normalized Hadamard matrix denoted HN with i-th row representing si(n) and then iteratively do the following.
The order of the steps 1 to 4 is unimportant and can be changed. These steps are consistent with the transformations that ensure HN remains an SCSS. Furthermore, these transforms represent elementary transforms of Hadamard matrices; that is, the resulting matrix always remains Hadamard.
If N is small enough and there is enough computing power, the search can be done deterministically instead of randomly.
For N=8, by using the foregoing procedure, Hadamard matrices with max{|sl(n)|}=0 for −(N−1)≤n≤−1, which are called ZCSSB, have been found. An example is Equation (18).
Furthermore, by modifying the foregoing search criterion to minimize max{|sl(n)|} for −(N−1)≤n≤(N−1), an N=8 ZCSSB, given in Equation (19), was found that has sl(n)=0 for −(N−1)≤n≤(N−1).
Then, for N=16, a ZCSSB was constructed from H81 as in Equation (20).
In Equation (20), H82 represents a version of H81 with each second row inverted, as given in Equation (21).
Direct search for N≥16 did not result in any ZCSSB, but only in LCSSB.
It is also noted that the operational steps described in any of the exemplary aspects herein are described to provide examples and discussion. The operations described may be performed in numerous different sequences other than the illustrated sequences. Furthermore, operations described in a single operational step may actually be performed in a number of different steps. Additionally, one or more operational steps discussed in the exemplary aspects may be combined. It is to be understood that the operational steps illustrated in the flowchart diagrams may be subject to numerous different modifications, as will be readily apparent to one of skill in the art. Those of skill in the art will also understand that information and signals may be represented using any of a variety of different technologies and techniques. For example, data, instructions, commands, information, signals, bits, symbols, and chips that may be referenced throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, optical fields or particles, or any combination thereof.
The previous description of the disclosure is provided to enable any person skilled in the art to make or use the disclosure. Various modifications to the disclosure will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other variations. Thus, the disclosure is not intended to be limited to the examples and designs described herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
The present application claims priority to U.S. Provisional Patent Application Ser. No. 63/319,977 filed on Mar. 15, 2022, and entitled “LOW-COMPLEXITY CHANNEL ESTIMATION USING SUPERCOMPLEMENTARY BLOCKS OF SEQUENCES,” the contents of which is incorporated herein by reference in its entirety. The present application also claims priority to U.S. Provisional Patent Application Ser. No. 63/373,902 filed on Aug. 30, 2022, and entitled “SYSTEMS AND METHODS FOR CHANNEL ESTIMATION,” the contents of which is incorporated herein by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
63373902 | Aug 2022 | US | |
63319977 | Mar 2022 | US |