Embodiments of the present invention relate to bearing determination. In particular; they relate to an apparatus, a method and a computer program for determining a bearing. Using, orthogonal frequency division multiplexed (OFDM) signals.
There are known techniques for determining a tearing using radio frequency (RF) signals. For example, an RF signal that is transmitted by a transmission apparatus may fee received at an antenna arrangement that comprises multiple antennas. A bearing from the antenna arrangement to the transmitting apparatus may be determined, for example, by measuring the received signal strength intensity (RSSI) at each of the multiple antennas.
According to various, but not necessarily all, embodiments of the invention there is provided art apparatus comprising: at least a first transformer configured to transform a first signal formed from a set of multiple orthogonal subcarriers and received via a first path, from a time domain to a frequency domain and produce for each of a plurality of the multiple orthogonal subcarriers a respective first coefficient; and a second transformer configured to transform a second signal formed from the set of multiple orthogonal subcarriers and received via a second path, different from the first path, from a time domain to a frequency domain and produce for each of a plurality of the multiple orthogonal subcarriers a respective second coefficient; and first coefficients and the plurality of second coefficients to determine a bearing for the apparatus.
According to various, but not necessarily all, embodiments of the invention there is provided a method, comprising: transforming a first signal formed from a set of multiple orthogonal subcarriers and received via a first path, from a time domain to a frequency domain, to produce for each of a plurality of the multiple orthogonal subcarriers a respective first coefficient; transforming a second signal formed from the set of multiple orthogonal subcarriers and received via a second path, different from the first path, from a time domain to a frequency domain, to produce for each of a plurality of the multiple orthogonal subcarriers a respective second coefficient; and processing the plurality of first coefficients and the plurality of second coefficients to determine a bearing for the apparatus.
According to various, but not necessarily all, embodiments of the invention there is provided an apparatus, comprising: means for transforming a first signal formed from a set of multiple orthogonal subcarriers and received via a first path, from a time domain to a frequency domain to produce for each of a plurality of the multiple orthogonal subcarriers a respective first coefficient; means for transforming a second signal formed from the set of multiple orthogonal subcarriers and received via a second path, different from the first path, from a time domain to a frequency domain to produce for each of a plurality of the multiple orthogonal subcarriers a respective second coefficient; and means for processing the plurality of first coefficients and the plurality of second coefficients to determine a bearing for the apparatus.
According to various, but not necessarily all, embodiments of the invention there is provided an apparatus, comprising: a diversity antenna arrangement for receiving a signal comprising multiple orthogonal subcarriers from a transmitter and comprising a first antenna at a first position for receiving the signal via a first path and a second antenna at a second position for receiving the signal via a second path, different to the first path; a first transformer configured to transform the signal received by the first antenna from a time domain to a frequency domain and produce for each of a plurality of the multiple orthogonal subcarriers a respective first coefficient; a second transformer configured to transform the signal received by the second antenna from a time domain to a frequency domain and produce for each of a plurality of the multiple orthogonal subcarriers a respective second coefficient; and processing circuitry configured to process the plurality of first coefficients and the plurality of second coefficients to determine a bearing for the transmitter.
According to various, but not necessarily all, embodiments of the invention there is provided an apparatus, comprising: processing circuitry configured: to determine a plurality of phase values from one or more received signals, each of the determined phase values falling within a predetermined range defined by a maximum value and a minimum value, to apply a scaling factor to the determined phase values to produce scaled phase values failing within the range, to determine an average scaled phase value by averaging the scaled-phase values, and to determine an average phase value for the plurality of phase values by applying the scaling factor to the average scaled phase value.
According to various, but not necessarily all, embodiments of the invention, there is provided a method, comprising: determining a plurality of phase values from one or more received signals, each of the determined phase values failing within a predetermined range defined by a maximum value and a minimum value; applying a scaling factor to the determined phase values to produce scaled phase values failing within the range; determining an average scaled phase value by averaging the scaled phase values; and determining an average phase value for the plurality of phase values by applying the scaling factor to the average scaled phase value.
For a better understanding of various examples of embodiments of the present invention reference will now be made by way of example only to the accompanying drawings in which;
The respective transformers Ti may, for example, be provided by different distinct hardware circuits arranged in parallel to operate simultaneously or may be provided by a single hardware circuit that is used in a time multiplexed manner so that it sequentially operates as each of the respective transformers Ti in turn.
Each of the N input signals si has traveled along a different signal path pi to arrive at the apparatus 2. A different path may result from having multiple receiver locations at the apparatus (receiver diversity) and/or may result from having multiple transmission locations for the received signals si (transmitter diversity).
Referring back to
For example, a received signal si may be represented by:
where the frequencies ωj are mutually orthogonal sub-carrier frequencies, and aij is a complex weighting coefficient that modulates the frequency ωj e.g. a received symbol.
The transmitted signal 12i originally transmitted along signal path pi (which is received as signal si) may be represented as
The complex weighting coefficient aij may be modelled as a multiplication of the original transmitted coefficient bij and a path-sensitive complex value xij.
The path-sensitive complex value xij can be represented as a magnitude Rij and a phase φij.
aij=bij*Rijexp(jφij)
where the frequencies ωj are mutually orthogonal sub-carrier frequencies, bij is a complex weighting coefficient e.g. a transmitted symbol and −π≦φij≦π.
The set of path-sensitive complex values xij record for each pairing of signal path pi and frequency ωj the effect that signal path pi has had in changing the originally transmitted data, the coefficient bij.
The transformer Ti in this example is configured to transform a received signal si from a time domain to a frequency domain and produce for each of the M multiple orthogonal sub-carriers ωj a respective coefficient aij to create a set Ai of coefficients {ai1, ai2, . . . aiM}. The coefficients are provided to processing circuitry 4.
The transform that is performed fey the transformer Ti may be, for example, an inverse discrete Fourier transform. A fast Fourier transform algorithm may in some embodiments be used to perform the discrete Fourier transform.
The set of coefficients {ai1, ai2, . . . aiM} are associated with a particular signal path pi. Each of the complex valued coefficients comprises a path-sensitive complex value xij that includes phase information φij and amplitude information Rij. This information is dependent upon the length of the signal path pi. For example, the difference between φij for path p1 and p2 depends upon the difference in position and orientation of the first transmission point 101 and the second transmission point 101 when there is transmitter diversity. For example, the difference between φij for path p1 and p2 depends upon the difference in position and orientation of the first antenna 201 and the second antenna 202 when there is receiver diversity. The information may also be subject to dispersion effects and have some dependence upon the dispersive properties of the signal path pi at the time the received signal si traveled along that path. The set of coefficients {ai1, ai2, . . . aiM} share in common that they have traveled along the same signal path pi but are each associated with a different frequency ωj and therefore potentially different dispersion effects.
The processing circuitry 4 may obtain phase information φij and amplitude information Rij using a known stored value of bij (e.g. the transmitted symbol) and the determined coefficient aij. For example, each of the determined coefficients aij may be compared with a known stored value of bij.
The processing circuitry 4, as illustrated by block 6, combines phase information φij and amplitude information Rij for the set of coefficients {ai1, ai2, . . . aiM} to obtain the representative phase information φi and representative amplitude information Ri for the signal path pi. This may compensate for dispersive effects and random noise effects.
The representative phase information φi and representative amplitude information Ri for the signal path pi may be the average phase information φi and average amplitude information Ri for the signal path pi.
For example, the average amplitude information Ri may be determined using the following equation:
The phase information φij that is determined using known stored values of bij and the determined coefficients aij relates to phase values that fall within a predetermined range defined by maximum and minimum values. For example, it may be that −π≦φij≦π.
In the event that the phase values φj for a particular signal path pi are situated close to the minimum and maximum boundaries of the range, it may not be possible to determine the average phase information φi correctly by simply calculating the arithmetic mean of the phase values φj.
A first method for determining the average phase information φi involves determining phase values φj for a particular signal path pi, and treating each phase value φj as a vector of arbitrary length. The vectors are summed and the angle of the resulting vector provides the average phase information φi.
A second method 100 for determining the average phase information φi, is illustrated in
At step 101 of
At step 102, a scaling factor φs is applied to each of the phase values φj. In this particular example, the scaling factor φs is subtracted from the phase values φj to produce a plurality of scaled phase values φj′:
φ′j=φj−φs
The scaling factor φs may be determined from the phase values φi. For example, the scaling factor φs may be set as:
φs=φ1−φm
where φm is the midpoint of the range for φij and φ1 is the first phase value that is determined for a particular signal path pi. For example, if −π≦φij≦π, then φm is zero, so φs=φ1.
If the variance of the phase values φj is small and the first phase value φ1 is located close to either the minimum or maximum boundary of the range, setting the scaling factor φs in accordance with the above equation will produce a plurality of scaled phase values φj′ that are located away from the minimum and maximum boundaries of the range.
In some embodiments of the Invention, the processing circuitry 4 determines whether the first phase value φ1 appears to be erroneous before it is used to set the scaling value φs. For example, it may do this by calculating the difference between the first phase value φ1 and some or all of the other phase values φi for that signal path pi. In the event that the average difference is above a threshold, the processing circuitry 4 may determine that the first phase value φ1 is erroneous, it may perform the same procedure for other phase values φj until an appropriate value is found for use in setting the scaling factor φs.
In order to ensure that all of the scaled phase values φj′ fall within the range −π≦φij′≦π, the scaling factor φs may be applied to the phase values φi using modular arithmetic (also known as modulo arithmetic). In this particular case, type of the modular arithmetic used is “modulo 2π” because the phase values φij are measured in radians and the range is 2π. If the phase values φij were measured in degrees, the type of modular arithmetic used would be “modulo 360”.
At step 103, an average scaled phase value φi′ is determined by calculating the arithmetic mean of the scaled phase values φi:
At step 104, the sealing factor φs is applied to the average scaled phase value φi′ to determine an average phase value φi. In this particular example, the scaling factor φs is added to the average scaled phase value φi′, because the scaling factor φs was subtracted from the phase values φj to produce a plurality of the scaled phase values φj′. Therefore:
φi=φ′i+φs
In order to ensure that the average phase value φi falls within the range −π≦φi≦π, the scaling factor φs may be applied to the average scaled phase value φi′ using modular arithmetic. As indicated above, in this particular case, the type of modular arithmetic used is “modulo 2π” because the phase values φij are measured in radians and the range is 2π.
The combination, for a signal path pi, of the phase information φij and amplitude information Rij for each of the sub-carrier frequencies ωj results in a fast and accurate assessment of the phase information φi and amplitude information Ri for the signal path pi without the need for frequency hopping.
If the transmitted coefficients bij are the same value b for all values of j then the set of coefficients ai1, ai2, . . . aiM} can be averaged and the average compared at block 8 to b to obtain the phase information φi and amplitude information Ri for the signal path pi.
If the transmitted coefficient bij are different for values of j then the set of coefficients ai1, ai2, . . . aiM} are compared at block 8 to the respective transmitted coefficient bij to obtain a set of path-sensitive complex values xij. The set of path-sensitive complex values xij are then averaged at block 6 to obtain the average phase information φi and amplitude information Ri for the signal path pi.
The transmitted coefficients bij may have a value b that is constant for all paths and frequencies, alternatively the transmitted coefficients bij may have a value bi that is constant for all frequencies of a particular signal path pi but changes to a different constant for different signal paths, alternatively the transmitted coefficients bij may have a value bij that may be different for different frequencies of a path and/or may be different for different paths.
At block 202 a received signal si is transformed from a time domain to a frequency domain. This produces for each of the M multiple orthogonal sub-carriers ωj a respective coefficient aij to create a set Ai of coefficients {ai1, ai2, . . . aiM}.
The set of coefficients {ai1, ai2, . . . aiM} are associated with a particular signal path pi. Each of the complex valued coefficients in a set comprise a path-sensitive complex value xij that includes phase information φij and amplitude information Rij. At block 204, phase information φij and amplitude information Rij for the set of coefficients {ai1, ai2, . . . aiM} is combined to obtain the representative phase information φj and representative amplitude information Ri for the signal path pi.
The combination, for a signal path pi, of the phase information φij and amplitude information Rij for each of the subcarrier frequencies ωj results in a fast and accurate assessment of the phase information φi and amplitude information Ri for the signal path pi without the need for frequency hopping.
Blocks 206 and 208, result in the repetition of blocks 202 and 204 for each of the N signal paths pi.
At block 208, the phase information φi and amplitude information Ri for all or some of the N signal paths is processed to determine the bearing of the apparatus 2.
At block 208, the processing circuitry may renormaiize the phase information φi and amplitude information Ri for the signal paths pi. For example, a particular signal path pr may be designated a reference path and phase information φr and amplitude information Rr for this path are used as a reference for the phase information φi and amplitude information Ri of the other paths.
The processing circuitry 4 may, for example, calculate φre=φi−φr for each i=1, 2 . . . N and calculate Rre=Ri/Rr for each i=1, 2 . . . N. φre represents the phase introduced by signal path pi relative to signal path pi. Rre represents the gain introduced by signal path pi relative to signal path pr. The set of pairs φre, Rre for each of the signal paths pi defines a bearing for the apparatus 2.
The processing circuitry 4 may access a lookup table that has measured reference pairs φre, Rre for each possible bearing (each possible combination of azimuth angle and elevation angle). A particular set of pairs φre, Rre may be used to look-up a bearing.
The lookup table may have measured reference pairs φre, Rre for a limited number of the possible bearings. A correlation may be performed between the measured pair φre, Rre and the set of reference pairs φre, Rre to identify a closest matching reference pair φre, Rre. The closest matching reference pair φre, Rre is then used to look-up a bearing in the lookup table.
The look-up table may be created by calibrating the apparatus 2. The values in the look-up table depend upon the relative displacement of the transmission points 101, 102 . . . if transmission diversity is used and depend upon the relative displacement of the antennas 201, 202 . . . if receiver diversity is used.
There are other mechanisms for calculating a bearing. For example, if for example, receiver diversity is used and the antennas 201, 202 . . . are arranged; along three orthogonal axis then the average phase information φi for each path in combination with straightforward trigonometry may be used by the processing circuitry 4 to determine the bearing.
The blocks illustrated in
In this example, the functionality of the transformers Ti are carried out in hardware by a single transform hardware circuit 40 and the functionality of the processing circuitry 4 is carried out in a processing hardware circuit 42. In some embodiments, the transform hardware circuit and the processing hardware circuit may be provided in a single component whereas in other embodiments they may fee provided as separate components. The transform hardware circuit sequentially operates as each of the transformers Ti. Switching circuitry 44 upstream of the transform hardware circuit 40 provides the received signals si in a time division multiplexed manner.
The transform hardware circuit 40 and the switching circuitry 44 may be part of a typical OFDM receiver. Although such a typical OFDM receiver may also have processing hardware circuitry it will not perform the operations of processing circuitry 4 described, for example, in relation to
Furthermore, it should be appreciated that the functionality of the transformers Ti and/or the processing circuitry 4 could, alternatively be implemented using a processor 50 and a memory 52 storing a computer program 54 as illustrated in
The processor 50 is configured to read from and write to the memory 52. The processor 50 may also comprise an output interface via which data and/or commands are output by the processor 50 and an input interface via which data and/or commands are input to the processor 50.
The memory 52 stores a computer program 54 comprising computer program instructions that control the operation of the apparatus 2 when loaded into the processor 50. The computer program instructions 54 provide the logic and routines that enables the apparatus 2 to perform one or more of the blocks illustrated in
The computer program 54 may arrive at the apparatus 2 via any suitable delivery mechanism 56. The delivery mechanism 58 may be, for example, a computer-readable storage medium, a computer program product, a memory device, a record medium such as a CD-ROM or DVD, an article of manufacture that tangibly embodies the computer program. The delivery mechanism may be a signal configured to reliably transfer the computer program. The apparatus 2 may in some implementations propagate or transmit the computer program 54 as a computer data signal.
Although the memory 52 is illustrated as a single component it may be implemented as one or more separate components some or all of which may be integrated/removable and/or may provide permanent/semi-permanent/dynamic/cached storage.
References to ‘computer-readable storage medium’, ‘computer program product’, ‘tangibly embodied computer program’ etc. or a ‘controller’, ‘computer’, ‘processor’ etc. should be understood to encompass not only computers having different architectures such as single/multi-processor architectures and sequential (e.g., Von Neumann)/parallel architectures but also specialized circuits such as field-programmable gate arrays (FPGA), application specific circuits (ASIC), signal processing devices and other devices. References to computer program, instructions, code etc. should be understood to encompass software for a programmable processor or firmware such as, for example, the programmable content of a hardware device whether instructions for a processor, or configuration settings for a fixed-function device, gate array or programmable logic device etc.
Although embodiments of the present invention have been described in the preceding paragraphs with reference to various examples, if should be appreciated that modifications to the examples given can be made without departing from the scope of the invention as claimed. For example, the method of determining an average phase value illustrated in
The description of the method of
Features described in the preceding description may be used in combinations other than the combinations explicitly described.
Although functions have been described with reference to certain features, those functions may be performable by other features whether described or not.
Although features have been described with reference to certain embodiments, those features may also be present in other embodiments whether described or not.
Whilst endeavoring in the foregoing specification to draw attention to those features of the invention believed to be of particular importance it should be understood that the Applicant claims protection in respect of any patentable feature or combination of features hereinbefore referred to and/or shown in the drawings whether or not particular emphasis has been placed thereon.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/EP2008/059666 | 7/23/2008 | WO | 00 | 3/25/2011 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2010/009763 | 1/28/2010 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
6791388 | Buchwald et al. | Sep 2004 | B2 |
7058002 | Kumagai et al. | Jun 2006 | B1 |
8184523 | Belotserkovsky et al. | May 2012 | B2 |
8374273 | Zhang et al. | Feb 2013 | B1 |
20090323865 | Bradley et al. | Dec 2009 | A1 |
20100002789 | Karabinis | Jan 2010 | A1 |
Number | Date | Country |
---|---|---|
1471703 | Oct 2004 | EP |
Entry |
---|
Yanxing Zeng et al., “Direct decoding of uplink space-time block coded multicarrier code division multiple access systems” Communications, 2004 IEEE International Conference on Paris, France Jun. 20-24, vol. 4, Jun. 20, 2004, pp. 2372-2376, XP010709687. |
Jian Zuo et al., “Development of TVA SuperPDC: Phasor applications, tools, and event replay”, Power and Energy Society General Meeting Conversion and Delivery of Electrical Energy in the 21st Century, 2008 IEEE, IEEE, Piscataway, NJ, USA, Jul. 20, 2008, pp. 1-8, XP031304064. |
Lin et al., “Applying Direction Finding Algorithms for Obtaining Processing Gains in OFDM System with Multiple Receiving Antennas” 2003 IEEE, pp. 248-251. |
Bossert et al., “Improved channel estimation with decision feedback for OFDM systems”, Electronics Letters, May 28, 1998, vol. 34, No. 11, pp. 1064-1065. |
Reede, “Ranging and Location for 802.22 WRANs”, Oct. 10, 2006, doc.: IEEE 802.22-06/0206r0. |
International Search Report for PCT/EP2008/059666 dated Aug. 4, 2009. |
Non-final Office Action in U.S. Appl. No. 13/534,185, dated Apr. 8, 2013. |
Non-Final Office Action dated Apr. 8, 2013 of related U.S. Appl. No. 13/534,185. |
Number | Date | Country | |
---|---|---|---|
20110194634 A1 | Aug 2011 | US |