The invention relates generally to precision GNSS location or attitude systems. More specifically, the invention relates to systems that utilize Global Navigation Satellite System (GNSS) signals to infer differential path length of carrier signals arriving at two or more antennas.
The determination of the orientation of an object's axes relative to a reference system is often of interest. Depending on the application, the orientation and reference system may be in two dimensions (2D) or in three dimensions (3D). In the case of two-dimensional systems, terms such as azimuth, heading, elevation, pitch, and inclination may be used in place of attitude.
There are many techniques in use to measure 2D and 3D attitude. Common techniques include using a magnetic compass to reference the object of interest to the local gravitational field, optical techniques to reference the object of interest to an earth-based or star-based reference frame, accelerometers to measure relative attitude or changes in attitude, and optical and mechanical gyroscopes for also measuring relative attitude. The merits of each technique are best judged according to the specific application or use. Likewise, each technique also exhibits disadvantages that may include accuracy, cost, and ease of use.
Recently, attitude determination using highly accurate space-based radio navigation systems has become possible. Such a radio navigation system is commonly referred to as a Global Navigation Satellite System (GNSS). A GNSS includes a network of satellites that broadcast radio signals, enabling a user to determine the location of a receiving antenna with a high degree of accuracy. To determine the attitude of an object, it is simply necessary to determine the position of two or more receiving antennas that have known placements relative to an object. Examples of GNSS systems include Navstar Global Positioning System (GPS), established by the United States; Globalnaya Navigatsionnay Sputnikovaya Sistema, or Global Orbiting Navigation Satellite System (GLONASS), established by the Russian Federation and similar in concept to GPS; and Galileo, also similar to GPS but created by the European Community and slated for full operational capacity in 2008.
Should it be necessary to improve the accuracy, reliability, or confidence level of an attitude or position determined through use of a GNSS, a Satellite-Based Augmentation System (SBAS) may be incorporated if one that is suitable is available. There are several public SBAS that work with GPS. These include Wide Area Augmentation System (WAAS), developed by the United States' Federal Aviation Authority, European Geostationary Navigation Overlay Service (EGNOS), developed by the European Community, as well as other public and private pay-for-service systems.
Currently the best-known of the available GNSS, GPS was developed by the United States government and has a constellation of 24 satellites in 6 orbital planes at an altitude of approximately 26,500 km. The first satellite was launched in February 1978. Initial Operational Capability (IOC) for the GPS was declared in December 1993. Each satellite continuously transmits microwave L-band radio signals in two frequency bands, L1 (1575.42 MHz) and L2 (1227.6 MHz). The L1 and L2 signals are phase shifted, or modulated, by one or more binary codes. These binary codes provide timing patterns relative to the satellite's onboard precision clock (synchronized to other satellites and to a ground reference through a ground-based control segment), in addition to a navigation message giving the precise orbital position of each satellite, clock correction information, and other system parameters.
The binary codes providing the timing information are called the C/A Code, or coarse acquisition code, and the P-code, or precise code. The C/A Code is a 1 MHz Pseudo Random Noise (PRN) code modulating the phase of the L1 signal and repeating every 1023 bits (one millisecond). The P-Code is also a PRN code, but modulates the phase of both the L1 and L2 signals and is a 10 MHz code repeating every seven days. These PRN codes are known patterns that can be compared to internal versions in the receiver. The GNSS receiver is able to compute an unambiguous range to each satellite by determining the time-shift necessary to align the internal code to the broadcast code. Since both the C/A Code and the P-Code have a relatively long “wavelength”—approximately 300 meters (or 1 microsecond) for the C/A Code and 30 meters (or 1/10 microsecond) for the P-Code, positions computed using them have a relatively coarse level of resolution.
To improve the positional accuracy provided by use of the C/A Code and the P-Code, a receiver may take advantage of the carrier component of the L1 or L2 signal. The term “carrier”, as used herein, refers to the dominant spectral component remaining in the radio signal after the spectral content resulting from the modulating PRN digital codes has been removed (e.g., from the C/A Code and the P-Code). The L1 and L2 carrier signals have wavelengths of about 19 centimeters and 24 centimeters, respectively. The GPS receiver is able to track these carrier signals and measure the carrier phase to a small fraction of a complete wavelength, permitting range measurement to an accuracy of less than a centimeter.
A technique to improve accuracy is realized by differencing GPS range measurements—known as Differential GPS (DGPS). The combination of DGPS with precise measurements of carrier phase leads to differential position accuracies of less than one centimeter root-mean-squared (i.e., centimeter-level positioning). Such accuracies are sufficient to determine the attitude of an object with 2 or more GPS GNSS antennas, typically spaced from 0.2 meters to 2 meters apart.
Accurate differential carrier phase is a primary concern for attitude determination or other precise GNSS positioning. Carrier phase data is available by tracking the carrier phase on either the L1 or L2 GPS signal. Navigation data is BPSK modulated onto both the L1 and the L2 carrier at a 50 Hz rate and, as such, the input carrier phase is subject to a 180 degree phase reversal every 20 milliseconds and the absolute phase can be inverted. The data modulation is removed from the carrier by means of a tracking loop known as a Costas loop.
A Costas loop results in a 180 degree phase ambiguity. That is, the Costas loop is just as likely to phase lock so that the binary 1's come out as binary 0's, and vice versa. The 180 degree phase ambiguity is of concern since it introduces an ambiguity of ½ of a carrier wavelength in the measured carrier phase. The wavelength of the L1 carrier is about 19 cm and the ½ cycle ambiguity is thus equivalent to 9.5 cm of measured phase.
Typically, in GPS receivers, the ½ cycle ambiguity is resolved by looking at certain data bits within the navigation message that are of a known value. If the bit is inverted over its known value, then the Costas loop is locked to the opposite phase and ½ cycle's worth of phase must be added to (or subtracted from) the measured carrier phase in order to maintain whole cycle phase alignment. It is of little consequence whether the ½ cycle is added or subtracted from the measured phase since, regardless, a whole cycle ambiguity is still present that must be removed by methods such as those described in U.S. Pat. No. 6,469,663 and/or U.S. patent application Ser. No. 11/243,112, entitled Attitude Determination Exploiting Geometry Constraints, to Whitehead et al., filed Oct. 4, 2005.
One problem that arises is that known data bits, such as those in the navigation message's preamble, arrive only so often. For example, the preamble itself is sent once every six seconds as the start of a 300 bit long sub-frame (see ICD-GPS-200). Other bits within the navigation message are known, or can be inferred from past data, but there are still many bits which are not known or which may not be predicted with 100% confidence. If the Costas loop is stressed (due to multipath fading, foliage attenuation, signal blockage, and the like) while the unknown bits are arriving, it may undergo a 180 degree phase shift that is not immediately detected. Hence, a ½ cycle error will arise in the measured phase that will persist until the next known data bit arrives. This may then reduce the accuracy in heading, pitch, or roll if the measured phase is used in an attitude determining device.
A further complication of a Costas loop is that the probability for a cycle slip is significantly higher when using a Costas loop as opposed to a conventional Phase Lock Loop (PLL). This is because the Costas loop is mathematically equivalent to a squaring loop that tracks the carrier phase at twice the carrier frequency. Phase tracking errors greater than 90 degrees may cause cycle slips in a Costas loop whereas phase errors of up to 180 degrees may be tolerated when using a PLL.
In an attitude system, such as that disclosed in Whitehead et al., the carrier phases arriving at two or more separate antennas are differenced with one another to create a differential carrier phase. The process of taking the difference cancels common mode errors such as satellite clock error and errors caused by propagation delays as the GNSS signal travels through the ionosphere and troposphere.
What is needed then is a Costas loop method that results in common ½ cycle ambiguity for carrier phases measured at both a master antenna and one or more slave antennas for a particular satellite. Being common, the ½ cycle ambiguity will cancel in the differential carrier phase.
Secondly, what is also needed is a method to make up for the loss of performance of a Costas loop over a conventional PLL. Again, since differencing is deployed, it is desirable that any cycle slips that arise on a carrier signal tracked from one antenna be present in the carrier signal tracked by a different antenna. Such cycle slips will cancel in the difference and will not affect the attitude or heading.
Thirdly, a method is needed that yields common, noise-induced effects in each individually tracked carrier phase so that the common-mode effects cancel in the differential carrier phase. A method with the aforementioned properties is applied to carrier tracking loops receiving data at two or more antennas that experience roughly similar motion or motion for which relative dynamic effects are low. The method further has the ability to track rapid clock-induced carrier phase when a common clock is employed.
Disclosed herein in an exemplary embodiment is a method of reducing Global Navigation Satellite System (GNSS) carrier tracking loop ambiguities comprising: receiving a plurality of GNSS satellite signals with a first antenna in operable communication with a first tracking device and a second antenna in communication with a second tracking device in at least one GNSS receiver; and sharing of data between the first tracking device and the second tracking device. The sharing is configured to facilitate a commonality in a carrier phase derived in the first and second tracking devices. The sharing also results in a cancellation of the commonality when a difference phase is formed between a carrier phase from the first tracking device and a carrier phase from the second tracking device.
Also disclosed herein in yet another exemplary embodiment is a system for reducing Global Navigation Satellite System (GNSS) carrier tracking loop ambiguities comprising: a first antenna in operable communication with a first tracking device configured to receive a plurality of GNSS satellite signals; and a second antenna in communication with a second tracking device configured to receive a plurality of GNSS satellite signals. The first tracking device and the second tracking device operable in at least one GNSS receiver. The first tracking device and the second tracking device are configured to share data therebetween to facilitate a commonality in a carrier phase derived in the first tracking device and the second tracking device. A difference phase is formed between a carrier phase from the first tracking device and another carrier phase from the second tracking device resulting in a cancellation of the commonality.
Further disclosed herein in another exemplary embodiment is a system for reducing Global Navigation Satellite System (GNSS) carrier tracking loop ambiguities comprising: means for receiving a plurality of GNSS satellite signals with a first antenna in operable communication with a first tracking device and a second antenna in communication with a second tracking device in at least one GNSS receiver. The system also includes means for sharing of data between the first tracking device and the second tracking device, the sharing configured to facilitate a commonality in a carrier phase derived in the first tracking device and the second tracking device. The sharing resulting in a cancellation in the commonality when a difference phase is formed between a carrier phase from the first tracking device and another carrier phase from the second tracking device.
Disclosed herein in yet another exemplary embodiment is a storage medium encoded with a machine-readable computer program code, the storage medium including instructions for causing a computing system to implement the abovementioned method for reducing Global Navigation Satellite System (GNSS) carrier tracking loop ambiguities.
Further, disclosed herein in yet another exemplary embodiment is a computer data signal, the computer data signal comprising code configured to cause a processor to implement the abovementioned method for reducing Global Navigation Satellite System (GNSS) carrier tracking loop ambiguities.
The invention also features a method for attitude determination or other differential GNSS positioning applications that is effective at reducing the occurrence of ½ ambiguities and cycle slips in the differential carrier phase. In an exemplary embodiment, this is accomplished by sharing information in the Costas tracking loops, designating one tracking device as a master while other slave tracking devices are driven by a carrier discriminator employing the master's in-phase in place of the slave's own in-phase data as is typically done. A related invention yields the benefit of partially canceling common noise-induced effects in differential carrier phase measured between two antennas.
Additional features, functions, and advantages associated with the disclosed methodology will be apparent from the detailed description which follows, particularly when reviewed in conjunction with the figures appended hereto.
To assist those of ordinary skill in the art in making and using the disclosed embodiments, reference is made to the appended figures, wherein like references are generally numbered alike in the several figures.
An exemplary embodiment of invention features a method and system for attitude determination or other GNSS positioning that significantly reduces the adverse effects of cycle slips and, simultaneously, problems resulting from ½ cycle phase ambiguities, such as may be encountered when employing a tracking loop, e.g., a phase tracking loop.
A preferred embodiment of the invention, by way of illustration, is described herein as it may be applied to attitude determination. While a preferred embodiment is shown and described by illustration and reference to attitude determination, it will be appreciated by those skilled in the art that the invention is not limited to attitude determination alone and may be applied to control systems, navigation, and the like, as well as combinations thereof. An embodiment of the invention as described herein may readily be applied to attitude determination as described in commonly assigned U.S. patent application Ser. No. 11/243,112, entitled Attitude Determination Exploiting Geometry Constraints, filed Oct. 4, 2005, the contents of which are incorporated by reference herein in their entirety.
It will further be appreciated that, while particular sensors, antennas, receivers and the like nomenclature associated therewith are enumerated to describe an exemplary embodiment, such terminology is utilized and described for illustration only and are not limiting. Numerous variations, substitutes, and equivalents will be apparent to those contemplating the disclosure herein.
In order to perform the prescribed functions and desired processing, as well as the computations therefor (e.g., the attitude determination processes, and the like), the receiver 108 may include, but not be limited to, a processor(s), computer(s), memory, storage, register(s), timing, interrupt(s), communication interface(s), and input/output signal interfaces, and the like, as well as combinations comprising at least one of the foregoing. For example, receiver 108 may include signal interfaces to enable accurate down-conversion and digitally sampling and tracking of satellite signals as needed to facilitate extracting the various timing and navigation information, including, but not limited to, the phase of the PRN code timing pattern. Additional features of the system 100, receiver 108 and the like, are thoroughly discussed herein.
One exemplary embodiment of an attitude determining system 100 shown in
A synchronization signal 214, denoted as sync, sent by the master tracking device 202 to the slave tracking devices 203 and 204 allows the master tracking device 202 and slave tracking devices 203, 204 to measure the code and carrier phase of each satellite signal 104 simultaneously. Furthermore, the RF down conversion within each tracking device 202, 203, and 204 and the sampling of data by each device is done using a common clock signal 220. When a single-difference phase observation is formed by subtracting the carrier (or code) phase measured by one tracking device with that measured by another tracking device for the same satellite, the portion of the phase due to the receiver's clock error is essentially eliminated in the difference.
The use of a single CPU connected to multiple shared synchronized tracking devices is advantageous. Data is easily shared among the tracking loops, and the carrier phase tracking loop may thus use a modified Costas loop. For any tracked satellite, e.g., 101, 102, 103, the in-phase reading from the Costas loop of a “master” tracking device, e.g., 202, dictates the sign of the quadrature-phase error discriminator not only for itself, but also for the discriminator of Costas loops in the remaining “slave” tracking devices, e.g., 203, 204. This forces the half-cycle phase ambiguity that is inherent in Costas tracking loops to be the same across all tracking devices for any given satellite so that it cancels when forming single-difference phase observations. Advantageously in this approach, the ambiguity is eliminated more quickly than methods that must wait to decode known navigation data bits that occur in headers or other parts of the navigation data.
In addition to driving the signals to baseband, the spreading Pseudo Random Noise (PN or PRN) sequence must be removed from the received signal (e.g., CA code spreads the L1 GPS signal). As depicted in the figure, this is accomplished by multipliers 414 and 415 which multiply the incoming signal by a replica model of the PN code that is generated by the Code NCO and PN generator 416. In the case of carrier tracking (as opposed to code tracking), the replica model is called the prompt PN code since it is steered by the code track-loop command 419 to align as closely as possible to the actual PN spreading code on the received signal (code tracking requires early and late combinations of a PN code, however this is not shown).
The data is then summed or low-pass filtered by summers 424 and 425, respectively. In the particular implementation shown, the summers perform sums over one millisecond which coincides with the length of the CA code's 1023 chip repetitive PN sequence (1023 chips with a 1.023 MHz chip rate yields a 1 millisecond span). Each summer accumulates and dumps its output. The in-phase output 426 is designated IP1ms and the quadrature phase output 427 is designated as Qp1ms where the superscript ‘1 ms’ indicates that summation is over one millisecond and the subscript ‘p’ indicates that this is the prompt channel (the channel relevant to carrier phase tracking).
It is often advantageous, at least for the purpose of recovering the navigation data bits, to increase the duration of the summations to span 20 ms which is the length of one navigation data bit. There should be no 180 degree phase shifts over the data bit's duration. Additional, but optional, summers 434 and 435 serve the purpose of continuing the sum for one data bit's duration and the Bit Synchronizer 436 controls the start and stop of the summers 434, 435 to coincide with the bit edge. The result is an in-phase bit aligned signal 436 denoted as IpBit and is quadrature counter part QpBit 437, respectively. The in-phase component, IpBit, is actually the demodulated navigation data bit. It is noteworthy to appreciate that in practice, this demodulated navigation data bit IpBit Bit can take on a sign opposite that of the true navigation data bit. However, the sign reversal is readily corrected with the parity algorithm specified in the ICD-GPS-200. The sign reversal, should it arise, is precisely a consequence of using a Costas loop with its potential for ½ cycle phase offset. The sign reversal remains for all arriving demodulated data bits until such a time that the track loop undergoes stress and slips back to the true ½ cycle phase alignment.
Continuing with
CD=Qpsign(Ip)
where the sign( ) function takes the arithmetic sign (+1 or −1) of the enclosed argument which, in this case, is the sign of Ip. Another common discriminator function is
CD=Tan−1(Qp/Ip)
Many discriminator functions are known and may readily be found in literature. Regardless of which is chosen, all serve the purpose of acting as a measure of phase tracking error. When the discriminator hovers around zero, the in-phase cosine signal 407 is properly aligned with the incoming carrier.
In
To facilitate further discussion and simplify the notation, henceforth, the superscripts denoting the duration of summation (such as bit and 1 ms) will be dropped with the understanding that a number of different durations are possible. The in-phase and quadrature phase data will be referred to simply as Ip and Qp hereafter.
Returning now to
The NCO's phase 460 is tapped off and sampled at a regular interval, say for example, every 1/20th of a second. Since the NCO phase is aligned to the arriving carrier phase by the track loop, the NCO phase provides a measure of the arriving carrier phase. This measure, termed the carrier phase observable, is used for carrier phase positioning, such as for example, an attitude system that calculates heading, pitch, and roll.
One thing that is evident, upon examination of the mathematical equations governing the carrier discriminators, is that if both Ip and Qp were to change sign, the discriminator output CD would remain unaffected since the change in sign cancels in the product (or quotient) of Ip and Qp. This is the desired effect of a Costas loop and results in an immunity of the Costas loop to the sign of the navigation data bit, which jointly effects the signs of both Ip and Qp. Otherwise, as for example if a simple PLL was used, every time the navigation data bit changed sign, the loop would react in an undesirable fashion by attempting to track the instantaneous 180 degree phase reversals, and likely the loop would behave erratically.
However, the benefit of a Costas loop of the prior art is also its disadvantage when one desires to recover the phase of the carrier without introducing a ½ cycle ambiguity. Returning to
phase error=(carrier phase of received signal)−(phase of NCO).
CD=Tan−1(Qp/Ip).
A conventional (non Costas) PLL is shown in
One way to achieve the benefits of a conventional PLL when employing a Costas loop is to know the sign of the modulating data in advance. It is then possible to use the discriminator
CD=Qpsign(Idata)
where Idata is defined as having the known sign of the modulating data. Since Idata is not a function of phase track error, the CD only depends on the term Qp and is thus similar to a conventional PLL discriminator in this respect. For example, the PLL discriminator shown in
It is worth noting that even if the sign of Idata is inverted so that it is opposite that of the true data, the Costas loop will simply undergo a onetime shift in phase of ½ cycle to negate the sign mismatch, but will then exhibit properties similar to those experienced when using the actual sign of the data. Consequently, the following discriminator also enables a Costas loop to behave like a conventional PLL, albeit with inverted data output.
CD=Qpsign(−Idata)
Referring now to
The slave track loop 650 has a subtle difference, however. Its own Ip data 652 is not fed to the salve carrier discriminator 654, but rather, Ip data 602 from the master track loop 600 is used in its place. The slave carrier discriminator 654 still makes use of its own Qp data 653. As with the master, the slave track filter 655 drives the carrier NCO 656 of the slave tracking loop 650.
Mathematically, the master and slave discriminators are
CDmaster=Qmastersign(Imaster)
CDslave=Qslavesign(Imaster)
Where we have adopted the following notation:
Now consider the slave's carrier discriminator
CDslave=Qslavesign(Imaster)
Clearly, CDslave depends only on the slave track loop through the quadrature-phase component, Qslave. For sign compensation, it no longer depends on its own in-phase component, Islave, but instead has an external dependence on the master's in-phase component, Imaster. As a consequence, the slave track loop 650 behaves more like a conventional PLL than a Costas loop. This is readily seen if you consider a carrier discriminator that uses the sign of the actual data
CDslave=Qslavesign(Idata).
If Imaster is taken as the true value of the data (Imaster=Idata), then the conventional PLL behavior is evident. The master's in-phase data does indeed maintain the same sign (or maintains consistent opposite sign) as the true data bits due to the nature of its Costas tracking loop. Furthermore, the data bits within the satellite broadcast signal arrive at master and slave antennas at substantially the same time as compared to the width of a data bit. When antennas are spaced apart by less than several meters (as they will be in a typical attitude determination device), the arrival time difference between master and slave data bits is less than ten nanoseconds, which is insignificant compared to the 20 millisecond data bit duration. Thus, as long as Imaster maintains a consistent sign or consistent opposite sign to the true data, the pull-in-zone for the slave is widened from λ/2 to λ.
In practice, the master track loop 600 will sometimes slip. If it slips ½ cycle (or any integer number of cycles plus ½ cycle), the slave control loop 650 will react by steering its NCO phase to an offset of ½ cycle to compensate for the sign reversal of Imaster. This is exactly the behavior that is desired for an attitude determining device, since ½ cycle phase offsets of the master e.g., 202, 600 and slave e.g., 203, 204, 650 are kept identical by the disclosed scheme. As such, any ½ cycle phase offsets will cancel in the difference between master and slave carrier phase. The carrier phase difference will only exhibit whole integer cycle phase ambiguities. Whole cycles are much easier to deal with than ½ cycles both in terms of estimating the ambiguities and when detecting cycle slips that may sometimes occur. Consequently, coupling the master's in-phase data 602 to the slave's carrier discriminator 654 yields a significant advantage.
In some attitude systems it is advantageous to not only deliver the in-phase data from the master tracking loop, e.g., 700 to the slave track loops, e.g., 701, 702, but share other track loop data as well. For example, when the master and slave antennas undergo motion that exhibits concurrent or similar dynamics, the tracking loops preferably need to react similarly and sharing of data between track loops is beneficial. In particular, in an attitude or heading system that experiences high translational accelerations (such as when mounted on an aircraft) but low rotational accelerations, it is likely that the dominant translational accelerations are seen by all antennas nearly equally. Track loops receiving data from each antenna should react nearly identically to track the common carrier phase accelerations.
In another situation, there may be minute physical accelerations experienced at the antennas, as for example, when a heading device is mounted on a large ship. Even if physical accelerations are small, the receiver's own oscillator undergoes frequency perturbations that result in apparent carrier phase accelerations. However, the oscillator is shared between master and slave tracking modules and the oscillator induced accelerations will be identical across each track loop.
As a solution of the aforementioned problem, in another exemplary embodiment, carrier phase tracking errors are summed (or averaged) and supplied to a responsive, high-bandwidth track loop filter that is shared among track loops. A lower bandwidth filter that is independent for each track loop, is utilized to track the slower but non-common dynamics using the individual phase errors of each system. A shared high-bandwidth track loop is employed to dominate the high frequency and noise response. As such, the noise induced into the individual carrier phase measurements will be similar and will cancel in the carrier phase differences that are formed as part of the attitude solution. This, advantageously, results in smoother differential carrier phase observations and ultimately attitude angles that are smoother as well.
Uniquely in this embodiment, phase tracking errors 805 and 855 of the master and slave, respectively, are each coupled to both track loops. Phase error 805 is fed to a summer 807 where it is added to phase error 855 that has been multiplied by a scale factor α using amplifier 856. Similarly, phase error 855 is fed to a summer 857 where it is added to phase error 805 that has also been multiplied by α using amplifier 806. The value α of amplifiers 806 and 856 is identical in this embodiment, however, other values may be employed. The summed signal produced by summer 807 is fed to the track filter 808 for the master tracking loop 801, which is either a proportional plus integral filter or some other track filter that results in a stable feedback control. The track filter 808 then delivers the NCO step 809 to the carrier NCO for the master tracking loop 801. Similarly, the summed signal produced by summer 857 is fed to the track filter 858 for the slave tracking loop 851, again, which is either a proportional plus integral filter or some other track filter that results in a stable feedback control. The track filter 858 then delivers the NCO step 859 to the carrier NCO for the slave tracking loop 851.
Continuing now to
Track errors 905 and 955 are summed at 920; the sum 921 is then delivered to the common track filter 922. The common track loop filter has its own set of parameters, such as proportional and integral gains. The output 923 of the common track loop is fed to the master NCO through summer 924 and to the slave NCO through summer 925.
The master's independent track loop filter 907 utilizes only the master's track error 905 which is scaled through the amplifier 906 having a gain of β. Similarly, the slave's independent track loop filter 957 utilizes only the slave's track error 955 which is scaled through amplifier 956 of gain β. In this embodiment, both track loops 907 and 957 employ identical filters (although it will be appreciated that they need not), which are typically adjusted to yield lower bandwidth response than the common track filter 922. The output 908 of the master's independent track filter is added to the common track filter output 923 and fed to the master's carrier NCO 910. The output 958 of the slave independent track filter is added to the common track filter output 923 and fed to the slave 's carrier NCO 960. Ultimate track loop response and degree of track loop independence is readily configurable by the designer.
To be consistent with the disclosure herein, the two systems need only ensure that the tracking loops of two or more receivers can share data at a rate consistent with the track loops update rate and that antennas that receive the carrier signals be sufficiently close together so that data bit transitions occur at roughly the same time (e.g., within 10%) relative to the in-phase and quadrature phase accumulation window start/stop times. It will be further appreciated that these constraints may be relaxed if the master and slave tracking loops run in non-real time (perhaps in post processing software) where methods can be applied to share and align data without concern for real-time processing constraints.
Furthermore, pipeline delays can be introduced into a real-time system's data flow so that even when antennas are spaced widely (or signals arrive at different times for other reasons such as hardware or filtering delays), data bit transitions can be aligned between track loops. In many circumstances, different bit arrival times are not an issue, for even if data accumulations occur over windows as small as one millisecond, antenna spacing would have to exceed 30 kilometers before signal travel delays induced a 10% misalignment of the accumulation windows.
It will be evident that there exist numerous numerical methodologies in the art for implementation of mathematical functions, in particular as referenced here, including, but not limited to, linearizations, least squares approximations, filters, Kalman filters, taking maximums, and summations. While many possible implementations exist, a particular method of implementation as employed to illustrate the exemplary embodiments should not be considered limiting.
The system and methodology described in the numerous embodiments hereinbefore provide a system and method of positioning or attitude determination that is effective and advantageously exploits Global Navigation Satellite System (GNSS) signals to infer differential path length of carrier signals arriving at two or more antennas. In particular, the described embodiments provide an improved tracking loop for carrier phase tracking. In addition, the disclosed invention may be embodied in the form of computer-implemented processes and apparatuses for practicing those processes. The present invention can also be embodied in the form of computer program code containing instructions embodied in tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other computer-readable storage medium, wherein when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus for practicing the invention. The present invention can also be embodied in the form of computer program code, for example, whether stored in a storage medium, loaded into and/or executed by a computer, or as data signal transmitted, whether a modulated carrier wave or not, over some transmission medium, such as over electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the computer program code is loaded into and executed by a computer, the computer becomes an apparatus for practicing the invention. When implemented on a general-purpose microprocessor, the computer program code segments configure the microprocessor to create specific logic circuits.
It will be appreciated that the use of “first” and “second” or other similar nomenclature for denoting similar items is not intended to specify or imply any particular order unless otherwise specifically stated. Likewise the use of “a” or “an” or other similar nomenclature is intended to mean “one or more”, unless otherwise specifically stated.
While the invention has been described with reference to an exemplary embodiment thereof, it will be understood by those skilled in the art that the present disclosure is not limited to such exemplary embodiments and that various changes may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In addition, a variety of modifications, enhancements, and/or variations may be made to adapt a particular situation or material to the teachings of the invention without departing from the essential spirit or scope thereof Therefore, it is intended that the invention not be limited to the particular embodiment disclosed as the best mode contemplated for carrying out this invention, but that the invention will include all embodiments falling within the scope of the appended claims.
Number | Name | Date | Kind |
---|---|---|---|
3727710 | Sanders et al. | Apr 1973 | A |
3899028 | Morris et al. | Aug 1975 | A |
4132272 | Holloway et al. | Jan 1979 | A |
4180133 | Collogan et al. | Dec 1979 | A |
4453614 | Allen et al. | Jun 1984 | A |
4637474 | Leonard | Jan 1987 | A |
4802545 | Nystuen et al. | Feb 1989 | A |
4858132 | Holmquist | Aug 1989 | A |
4918607 | Wible | Apr 1990 | A |
5031704 | Fleischer et al. | Jul 1991 | A |
5134407 | Lorenz et al. | Jul 1992 | A |
5152347 | Miller | Oct 1992 | A |
5156219 | Schmidt et al. | Oct 1992 | A |
5185610 | Ward et al. | Feb 1993 | A |
5202829 | Geier | Apr 1993 | A |
5207239 | Schwitaia | May 1993 | A |
5255756 | Follmer et al. | Oct 1993 | A |
5268695 | Dentinger et al. | Dec 1993 | A |
5296861 | Knight | Mar 1994 | A |
5343209 | Sennott et al. | Aug 1994 | A |
5369589 | Steiner | Nov 1994 | A |
5375059 | Krytsos et al. | Dec 1994 | A |
5390125 | Sennott et al. | Feb 1995 | A |
5416712 | Geier et al. | May 1995 | A |
5476147 | Fixemer | Dec 1995 | A |
5491636 | Robertson | Feb 1996 | A |
5511623 | Frasier | Apr 1996 | A |
5534875 | Diefes et al. | Jul 1996 | A |
5543804 | Buchler et al. | Aug 1996 | A |
5546093 | Gudat et al. | Aug 1996 | A |
5548293 | Cohen | Aug 1996 | A |
5561432 | Knight | Oct 1996 | A |
5592382 | Colley | Jan 1997 | A |
5610845 | Slabinski | Mar 1997 | A |
5612883 | Shaffer et al. | Mar 1997 | A |
5615116 | Gudat et al. | Mar 1997 | A |
5617317 | Ignagni | Apr 1997 | A |
5644139 | Allen et al. | Jul 1997 | A |
5564632 | Walker | Sep 1997 | A |
5664632 | Frasier | Sep 1997 | A |
5684696 | Rao et al. | Nov 1997 | A |
5691974 | Zehavi et al. | Nov 1997 | A |
5717593 | Gvili | Feb 1998 | A |
5725230 | Walkup | Mar 1998 | A |
5757316 | Buchler | May 1998 | A |
5765123 | Nimura et al. | Jun 1998 | A |
5777578 | Chang et al. | Jul 1998 | A |
5810095 | Orbach et al. | Sep 1998 | A |
5838562 | Gudat et al. | Nov 1998 | A |
5862501 | Talbot et al. | Jan 1999 | A |
5875408 | Bendett et al. | Feb 1999 | A |
5918558 | Susag | Jul 1999 | A |
5923270 | Sampo et al. | Jul 1999 | A |
5928309 | Korver et al. | Jul 1999 | A |
5933110 | Tang | Aug 1999 | A |
5935183 | Sahm et al. | Aug 1999 | A |
5940026 | Popeck | Aug 1999 | A |
5943008 | Dusseldorp | Aug 1999 | A |
5945917 | Harry | Aug 1999 | A |
5956250 | Gudat et al. | Sep 1999 | A |
5969670 | Kalafus et al. | Oct 1999 | A |
5987383 | Keller | Nov 1999 | A |
6014608 | Seo | Jan 2000 | A |
6052647 | Parkinson et al. | Apr 2000 | A |
6057800 | Yang et al. | May 2000 | A |
6061390 | Meehan et al. | May 2000 | A |
6062317 | Gharsalli | May 2000 | A |
6069583 | Silvestrin et al. | May 2000 | A |
6076612 | Carr et al. | Jun 2000 | A |
6122595 | Varley et al. | Sep 2000 | A |
6150978 | McBurney | Nov 2000 | A |
6191730 | Nelson, Jr. | Feb 2001 | B1 |
6191733 | Dizchavez | Feb 2001 | B1 |
6198992 | Winslow | Mar 2001 | B1 |
6199000 | Keller et al. | Mar 2001 | B1 |
6205401 | Pickhard et al. | Mar 2001 | B1 |
6229479 | Kozlov et al. | May 2001 | B1 |
6233511 | Berger et al. | May 2001 | B1 |
6236916 | Staub et al. | May 2001 | B1 |
6236924 | Motz | May 2001 | B1 |
6253160 | Hanseder | Jun 2001 | B1 |
6256583 | Sutton | Jul 2001 | B1 |
6259398 | Riley | Jul 2001 | B1 |
6285320 | Olster et al. | Sep 2001 | B1 |
6292132 | Wilson | Sep 2001 | B1 |
6313788 | Wilson | Nov 2001 | B1 |
6313789 | Zhodzishsky et al. | Nov 2001 | B1 |
6314348 | Winslow | Nov 2001 | B1 |
6325684 | Knight | Dec 2001 | B1 |
6345231 | Quincke | Feb 2002 | B2 |
6377889 | Soest | Apr 2002 | B1 |
6389345 | Phelps | May 2002 | B2 |
6397147 | Whitehead | May 2002 | B1 |
6421003 | Riley et al. | Jul 2002 | B1 |
6431576 | Viaud et al. | Aug 2002 | B1 |
6434462 | Bevly et al. | Aug 2002 | B1 |
6445983 | Dickson et al. | Sep 2002 | B1 |
6445990 | Manring | Sep 2002 | B1 |
6463374 | Keller et al. | Oct 2002 | B1 |
6466871 | Reisman et al. | Oct 2002 | B1 |
6469663 | Whitehead et al. | Oct 2002 | B1 |
6501422 | Nichols | Dec 2002 | B1 |
6515619 | McKay, Jr. | Feb 2003 | B1 |
6539303 | McClure et al. | Mar 2003 | B2 |
6542077 | Joao | Apr 2003 | B2 |
6553299 | Keller et al. | Apr 2003 | B1 |
6553311 | Aheam et al. | Apr 2003 | B2 |
6556942 | Smith | Apr 2003 | B1 |
6577952 | Geier et al. | Jun 2003 | B2 |
6587761 | Kumar | Jul 2003 | B2 |
6611228 | Toda et al. | Aug 2003 | B2 |
6622091 | Perlmutter et al. | Sep 2003 | B2 |
6631916 | Miller | Oct 2003 | B1 |
6643576 | O'Connor et al. | Nov 2003 | B1 |
6646603 | Dooley et al. | Nov 2003 | B2 |
6657585 | Kucik | Dec 2003 | B1 |
6671587 | Hrovat et al. | Dec 2003 | B2 |
6688403 | Bernhardt et al. | Feb 2004 | B2 |
6703973 | Nichols | Mar 2004 | B1 |
6711501 | McClure et al. | Mar 2004 | B2 |
6775802 | Gaal | Aug 2004 | B2 |
6792380 | Toda | Sep 2004 | B2 |
6865465 | McClure | Mar 2005 | B2 |
6922635 | Rorabaugh | Jul 2005 | B2 |
7027918 | Zimmerman et al. | Apr 2006 | B2 |
7031725 | Rorabaugh | Apr 2006 | B2 |
7162348 | McClure et al. | Jan 2007 | B2 |
20020029110 | Fukuda et al. | Mar 2002 | A1 |
20020038171 | Deguchi et al. | Mar 2002 | A1 |
20020072850 | McClure et al. | Jun 2002 | A1 |
20020165648 | Zeitler | Nov 2002 | A1 |
20020165669 | Pinto et al. | Nov 2002 | A1 |
20020169553 | Perimutter et al. | Nov 2002 | A1 |
20020175858 | Takahashi | Nov 2002 | A1 |
20030009282 | Upadhyaya et al. | Jan 2003 | A1 |
20030187560 | Keller et al. | Oct 2003 | A1 |
20040186644 | McClure et al. | Sep 2004 | A1 |
20040210357 | McKay et al. | Oct 2004 | A1 |
20040212533 | Whitehead | Oct 2004 | A1 |
20050116859 | Miller et al. | Jun 2005 | A1 |
20060251173 | Wang et al. | Nov 2006 | A1 |
Number | Date | Country |
---|---|---|
WO9836288 | Aug 1998 | WO |
WO0024239 | May 2000 | WO |
WO03019430 | Mar 2003 | WO |
WO2005119386 | Dec 2005 | WO |
Number | Date | Country | |
---|---|---|---|
20070085738 A1 | Apr 2007 | US |