US20020029110A1  System for determining the heading and/or attitude of a body  Google Patents
System for determining the heading and/or attitude of a body Download PDFInfo
 Publication number
 US20020029110A1 US20020029110A1 US09/870,628 US87062801A US2002029110A1 US 20020029110 A1 US20020029110 A1 US 20020029110A1 US 87062801 A US87062801 A US 87062801A US 2002029110 A1 US2002029110 A1 US 2002029110A1
 Authority
 US
 United States
 Prior art keywords
 attitude
 body
 determining
 antennas
 means
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Granted
Links
 239000000969 carrier Substances 0 abstract claims description 54
 239000011159 matrix materials Substances 0 description 6
 230000001131 transforming Effects 0 description 5
 239000011575 calcium Substances 0 description 4
 239000000203 mixtures Substances 0 description 3
 238000005096 rolling process Methods 0 description 3
 238000010276 construction Methods 0 description 2
 239000011295 pitch Substances 0 description 2
 239000000047 products Substances 0 description 2
 238000009825 accumulation Methods 0 description 1
 238000004422 calculation algorithm Methods 0 description 1
 238000004364 calculation methods Methods 0 description 1
 230000001427 coherent Effects 0 description 1
Images
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
 G01C21/00—Navigation; Navigational instruments not provided for in preceding groups G01C1/00G01C19/00
 G01C21/10—Navigation; Navigational instruments not provided for in preceding groups G01C1/00G01C19/00 by using measurements of speed or acceleration
 G01C21/12—Navigation; Navigational instruments not provided for in preceding groups G01C1/00G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
 G01C21/16—Navigation; Navigational instruments not provided for in preceding groups G01C1/00G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
 G01C21/165—Navigation; Navigational instruments not provided for in preceding groups G01C1/00G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with noninertial navigation instruments

 G—PHYSICS
 G01—MEASURING; TESTING
 G01S—RADIO DIRECTIONFINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCEDETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
 G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
 G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
 G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting timestamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
 G01S19/53—Determining attitude

 G—PHYSICS
 G01—MEASURING; TESTING
 G01S—RADIO DIRECTIONFINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCEDETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
 G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
 G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
 G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting timestamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
 G01S19/53—Determining attitude
 G01S19/54—Determining attitude using carrier phase measurements; using long or short baseline interferometry
 G01S19/55—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Leastsquares AMBiguity Decorrelation Adjustment] method
Abstract
A system for determining the heading and/or attitude of a body receives radio waves from a plurality of positionfixing satellites using at least three antennas fixedly mounted at different positions of the body. To reliably obtain integer ambiguity solutions of carrier phases of the radio waves in a shorter time, the system directly determines integer ambiguities from attitude angle data obtained by an IMU attitude processing section when the integer ambiguities are to be redetermined in the event of an interruption of the received radio waves or a change in the combination of satellites to be used. This system provides a user with highly accurate uninterrupted heading and/or attitude angle information.
Description
 The present invention relates to a system for determining the heading and/or attitude of a body moving or stationary, such as an aircraft or a vessel.
 There is being developed an attitude determining system (ADS) which receives radio waves transmitted from a plurality of positionfixing satellites by a plurality of antennas disposed at specific positions of a body with one of the antennas being used as a reference antenna, measures the phases of carrier signals of the radio waves, determines the relative geometric relationship of the onboard antennas by determining the relative positions of the antennas with respect to the reference antenna, and thereby measures the threedimensional attitude of the body.
 Another example of a conventional system for determining the attitude of a body utilizes an inertial measurement unit (IMU) sensor like a rate gyro which detects angular velocity.
 The conventional ADS determines integer ambiguities involved in a double difference of carrier phase or in a single difference of carrier phase by observing the double difference of carrier phase or the single difference of carrier phase between the carrier phases of the radio waves transmitted from the positionfixing satellites to determine the positions of the antennas relative to the reference antenna. In the conventional integer ambiguity solution method, however, there arise errors in solution statistically and stochastically with the probability of finding a correct integer ambiguity solution being 85% to 95%, for example, wherein the time required for obtaining the integer ambiguity solution greatly depends on observation noise level and baseline lengths (i.e., the lengths between the reference antenna and the other antennas whose relative positions should be determined). Under an ordinary noise level, 30 seconds to a few minutes are required on average to obtain an integer ambiguity solution when the baseline lengths are approximately 1 m. When the noise level is high, it has occasionally be impossible to obtain an integer ambiguity solution even after 10 minutes of observation, for example.
 In practical situations, it becomes necessary to obtain integer ambiguity solutions when the combination of satellites used for position fix is altered or when a cycle slip (a situation in which information on integer components of a value in a counter for counting the carrier phases is lost due to inability of a coherent delay lock loop of a carrier reconstruction circuit to track) occurs due to noise, for example. Therefore, if the time for determining the integer ambiguity solution is long, an attitude output for a user would be frequently interrupted for impermissibly long periods of time. Also, if the rate of correct integer ambiguity solutions is low, the system would output unreliable attitude information. This has been a major bottleneck for practical application of the conventional ADS.
 In an attitude determining system using the aforementioned IMU sensor, on the other hand, a major problem has been the accumulation of errors which occurs due to the fact that the system determines the threedimensional attitude of a body by integrating its angular velocity components in the direction of each of three axes x, y, z.
 Previous efforts to overcome the aforementioned problems by integrating the conventional ADS and IMU system are found in U.S. Pat. No. 5,349,531, in which attitude information is interpolated with attitude determination results obtained from an output of an IMU sensor when an attitude data output from the ADS is interrupted. The attitude determining system of U.S. Pat. No. 4,754,280 basically determines the attitude of a body based on the output of the IMU sensor and outputs it as attitude information and utilizes the attitude determination results obtained by the IMU system just for estimating and removing IMU sensor errors.
 When a widely available lowcost IMU sensor is used in a system for determining the attitude of a body by means of the IMU sensor, however, significant errors would occur, making it impossible to perform accurate error modeling. Therefore, even if the IMU sensor error is estimated and compensated for, it can not be completely removed. Thus, when a maximum permissible error is about 1°, backup time ensured by the IMU is approximately 30 to 60 seconds at the maximum. This means that the ADS should restore an reliable attitude data output within this backup time, so that the conventional ADS/IMU integrated system has not been suited for practical use.
 Also, neither of the aforementioned U.S. Patents are intended to solve the aforementioned problems involved in obtaining an integer ambiguity solution of a single difference of carrier phase or a double difference of carrier phase in the ADS.
 In view of the foregoing, it is an object of the invention to provide an attitude determining system which can reliably determine an integer ambiguity solution in a shorter time using an output of an IMU sensor.
 Another object of the invention is to provide a measuring system which can obtain an integer ambiguity based on an output of an IMU sensor.
 The present invention makes it possible to determine integer ambiguities of single differences of carrier phase or double differences of carrier phase observed when determining relative positions of multiple antennas based on attitude angles of a body obtained by means of an IMU sensor to thereby shorten the time required for obtaining an integer ambiguity solution, so that attitude determining results can be continuously output to a user even if cycle slips frequently occur.
 According to one aspect of the present invention, a system for determining the heading and/or attitude of a body comprises
 at least three antennas which receive radio waves from a plurality of positionfixing satellites with one of the antennas being used as a reference antenna, the antennas being fixed at different positions of the body, means for determining the relative positions of the other antennas with respect to the reference antenna by observing single differences of carrier phase or double differences of carrier phase between carrier phases of the radio waves and by calculating integer ambiguities of the phase differences to thereby determine the heading and/or attitude of the body, an IMU sensor installed on the body, means for obtaining attitude angles of the body based on an output of the IMU sensor, and means for determining the integer ambiguities based on the attitude angles obtained when the integer ambiguities are to be redetermined.
 According to another aspect of the invention, a system for determining the heading and/or attitude of a body as described above has the said means for determining the integer ambiguities which comprises means for obtaining a baseline vector in the local coordinate system from the attitude angles obtained based on the output of the IMU sensor, means for obtaining the path difference with respect to the baseline vector, and means for obtaining the integer ambiguities based on the path difference with respect to the baseline vector and measured phase difference.
 According to another aspect of the invention, a system for determining the heading and/or attitude of a body comprises two antennas which receive radio waves from a plurality of positionfixing satellites, the antennas being fixed at different positions of the body with one of the antennas being used as a reference antenna, means for determining the relative positions of the other antennas with respect to the reference antenna by observing single differences of carrier phase or double differences of carrier phase between carrier phases of the radio waves and by calculating integer ambiguities to thereby determine the heading and/or attitude of the body, an IMU sensor installed on the body, means for obtaining attitude angles of the body based on an output of the IMU sensor, and means for determining the integer ambiguities based on the attitude angles obtained from the output of the IMU sensor when the integer ambiguities are to be redetermined.
 According to further aspect of the invention, a system for determining the heading and/or attitude of a body comprises a plurality of antennas which receive radio waves from a plurality of positionfixing satellites, the antennas being fixed at different positions of the body with one of the antennas being used as a reference antenna, means for determining the relative positions of the other antennas with respect to the reference antenna based on integer ambiguities of the phase differences to thereby determine the heading and/or attitude of the body, an IMU sensor installed on the body, means for obtaining attitude angles of the body from an output of the IMU sensor, and means for determining the integer ambiguities based on the attitude angles obtained from the output of the IMU sensor when the integer ambiguities are to be redetermined.
 Since the integer ambiguities of the single differences of carrier phase or the double differences of carrier phase observed when calculating the relative positions of the multiple antennas are directly determined from the attitude angles of the body obtained by means of the IMU sensor in this invention, it is possible to output attitude information to the user substantially without interruption even when the combination of satellites to be used is altered or a cycle slip has occurred. Furthermore, it becomes possible to output highly reliable attitude information because the rate of correct integer ambiguity solutions increases.
 FIGS. 1A and 1B are diagrams showing an example of an antenna layout depicted in a local coordinate system and an antenna coordinate system, respectively;
 FIG. 2 is a diagram showing the relationship between baseline vectors and single differences of carrier phase, etc.;
 FIG. 3 is a block diagram showing the construction of an attitude determining system for a body according to the present invention; and
 FIG. 4 is a flowchart showing operations performed by a reinitialization judgment section and an IMU/GPS integration processing section of an attitude data processing block shown in FIG. 3.
 First, an antenna coordinate system, a local coordinate system and coordinate transformation between these two coordinate systems are explained referring to FIGS. 1A and 1B.
 FIG. 1A shows an antenna layout depicted in the local coordinate system, in which the x, y and z axes represent north, east and vertical directions, respectively, with the position of a reference antenna (first antenna) taken as an origin, for instance. In FIG. 1A, a point with coordinates (x_{1} ^{n}, y_{1} ^{n}, z_{1} ^{n}) is the position of a second antenna, a point with coordinates (x_{2} ^{n}, y_{2} ^{n}, z_{2} ^{n}) is the position of a third antenna, L^{1} ^{n }is a baseline vector from the first antenna to the second antenna (i.e., the vector indicating the length and direction of a straight line connecting the reference antenna chosen from multiple antennas to one of the other antennas whose relative position should be determined), and L_{2} ^{n }is a baseline vector from the first antenna to the third antenna. The superscript “n” written to the right of the values of the coordinates and the baseline vectors indicates that these coordinates and baseline vectors are for the local coordinate system.
 FIG. 1B shows an antenna layout depicted in the antenna coordinate system, in which the x axis represents the heading (i.e., the direction in which a ship's bow or the nose of an aircraft is oriented), the y axis corresponds to a direction perpendicular to the x axis, and the z axis indicates a direction perpendicular to both the x and y axes, with the position of the reference antenna (first antenna) taken as an origin, for instance. In FIG. 1B, a point with coordinates (x_{1} ^{a}, y_{1} ^{a}, z_{1} ^{a}) is the position of the second antenna, a point with coordinates (x_{2} ^{a}, y_{2} _{a}, z_{2} ^{a}) is the position of the third antenna, L_{1} ^{a }is a baseline vector from the first antenna to the second antenna, and L_{2} ^{a }is a baseline vector from the first antenna to the third antenna. The superscript “a” written to the right of the values of the coordinates and the baseline vectors indicates that these coordinates and baseline vectors are for the antenna coordinate system.

 where L_{3} ^{n }is a cross product of L_{1} ^{n }and L_{2} ^{n}, and L_{3} ^{a }is a cross product of L_{1} ^{a }and L_{2} ^{a}.

 where θr is the angle of roll (i.e., the angle of rotation about the x axis of the antenna coordinate system) a body, θp is the angle of pitch (i.e., the angle of rotation about the y axis of the antenna coordinate system), and θy is the angle of yaw (i.e., the angle of rotation about the z axis of the antenna coordinate system). These angles are obtained respectively from the following equations:
 Since L_{1} ^{a}, L_{2} ^{a }and L_{3} ^{a }of equation (1) are known quantities, C_{a} ^{n }can be obtained if L_{1} ^{n}, L_{2} ^{n }and L_{3} ^{n }are determined from carrier phase difference observables. When the value of C_{a} ^{n }is obtained, attitude angles can be calculated by equations (3) to (5) above. Here, the angle of yaw is a direction referenced to north. Since measuring accuracy of the attitude angles is determined by the accuracy of calculation of the coordinates (x_{1} ^{n}, y_{1} ^{n}, z_{1} ^{n}) and (x_{2} ^{n}, y_{2} ^{n}, z_{2} ^{n}), an attitude determining system usually observes carrier phase differences and determines the relative positions of the antennas with high accuracy instead of performing position fixing by using a pseudorandom noise code. While either a single difference of carrier phase method or a double difference of carrier phase method may be used for position fixing, the double difference of carrier phase method is described in the following discussion:

 where
 m: satellite No.
 k: number of satellites
 λ: carrier wavelength of a received signal
 δρ_{j} ^{n}: double path difference of a baseline vector L_{j} ^{n }to a reference satellite S^{1 }and a satellite S^{m+1 }(m=1, 2, 3)
 δφ_{j} ^{m}: double difference of carrier phase observable between the reference satellite S^{1 }and the satellite S^{m+1 }
 δN_{j} ^{m}: double integer ambiguity difference between the reference satellite S^{1 }and the satellite S^{m+1 }
 δe_{j} ^{m}: double observable noise difference between the reference satellite S^{1 }and the satellite S^{m+1 }
 FIG. 2 is a graphical representation of the double difference of carrier phase, in which individual symbols have the following meanings:
 Δρ_{j} ^{m+1}: single path difference of the baseline vector L_{j} ^{n }to the satellite S^{m+1 }
 Δφ_{j} ^{m+1}: single difference of carrier phase observable of the baseline vector L_{j} ^{n }to the satellite S^{m+1 }
 ΔN_{j} ^{m+1}: single integer ambiguity difference of the baseline vector L_{j} ^{n }to the satellite S^{m+1 }
 Δe_{j} ^{m+1}: single observable noise difference of the baseline vector L_{j} ^{n }to the satellite S_{m+1 }
 Also, the individual symbols and the aforementioned variables have the following relationship:
 δρ_{j} ^{m}=Δρ_{j} ^{m+1}+Δρ_{j} ^{1 }
 δρ_{j} ^{m}=Δφ_{j} ^{m+1}+Δφ_{j} ^{1 }
 δN _{j} ^{m} =ΔN _{j} ^{m+1} +ΔN _{j} ^{1 }
 δe _{j} ^{m} =Δe _{j} ^{m+1} +Δe _{j} ^{1 }
 The double difference of carrier phase δρ_{j} ^{m }of equation (6) can be expressed by δρ_{j} ^{m}=δA_{j} ^{m}L_{j} ^{n }(j=1, 2), where both δA_{j} ^{m }and L_{j} ^{n }are functions of the coordinates (x_{1} ^{n}, y^{1} ^{n}, z_{1} ^{n}). Therefore, δA_{j} ^{m }and L_{j} ^{n}, or the coordinates (x_{1} ^{n}, y^{1} ^{n}, z_{1} ^{n}), are unknown quantities. Since there remain three unknown quantities if integer ambiguities are predetermined by one method or other, it is possible to determine the coordinates (x_{1} ^{n}, y_{1} ^{n}, z_{1} ^{n}) if at least three observables δφ_{j} ^{m }(i.e., the number of satellites≧4) are obtained.
 Expressing now coordinate vectors of a given baseline vector L in the local coordinate system and the antenna coordinate system as defined in FIGS. 1A and 1B by X^{n }and X^{a }respectively, the relationship between X^{n}, X^{a }and a coordinate transformation matrix C_{a} ^{n }(3×3) is expressed by equation (7) below:
 X ^{n} =C _{a} ^{n} X ^{a} (7)


 Where σ^{p} ^{2 }is the variancecovariance matrix of the baseline vector. σ_{θ} ^{2 }is the variancecovariance matrix of the attitude angles obtained from the IMU. H is a transformation matrix between the baseline vector and the attitude angles.
 The relationship between the standard deviation σ_{p }of errors of the baseline vector position and the value δρ obtained by projecting in the direction of viewing filed of a satellite is expressed by equation (10) below:
 δρ=δAσ_{p} (10)

 It is now assumed that antenna bases located at vertices of a triangle whose sides are 50 cm each are swung with a rolling angle of ±20° at rolling angle varying intervals of 60 seconds, a pitching angle of ±15° at pitching angle varying intervals of 70 seconds, and with a yawing angle of ±30° at yawing angle varying intervals of 180 seconds, and the direction of a satellite and its elevation angle are varied through 360° and between 0° and 80°, respectively. If errors of IMU sensors for sensing rolling, pitching and yawing angles are within a few degrees in terms of standard deviation under the conditions stated above, a range error calculated in the direction of viewing filed of the satellite is equal to or smaller than ¼ of the wavelength of the carrier of radio waves (i.e., approximately 5 cm) transmitted from the satellite.
 This means that it is possible to obtain errorfree integer ambiguities of carrier phases even when the ambiguities of the carrier phases are rounded to integers using the attitude angles obtained from outputs of the IMU sensors under the aforementioned conditions because the attitude angle error is usually a few degrees or less.
 Real number ambiguities are calculated as follows from the attitude angles.
 Since the position (x_{j} ^{a}, y_{j} ^{a}, z_{j} ^{a}) of a baseline vector L_{j} ^{a }(j=1, 2) in the antenna coordinate system is known, the position (x_{j} ^{n}, y_{j} ^{n}, z_{j} ^{n}) of the baseline vector L_{j} ^{n }(j=1, 2) in the local coordinate system is obtained from the attitude angles detected by the IMU sensors using equations (2) and (7).

 Since A_{j} ^{m }is a vector in the direction of viewing filed of a satellite m of the baseline vector L_{j} ^{n}, it can be determined from the position (x_{m} ^{n}, y_{m} ^{n}, z_{m} ^{n}) of each satellite and the position (x_{j} ^{n}, y_{j} ^{n}, z_{j} ^{n}) in the local coordinate system, where m is the number of the used satellite and j is the number of the baseline vector.

 where a value obtained by rounding ΔN_{j} ^{m }to an integer is the integer ambiguity (single integer ambiguity difference).
 In the case of the double difference of carrier phase method, it is possible to calculate a double difference of carrier phase real number ambiguity δN_{j} ^{m }from the aforementioned equation (6). Therefore, the invention is applicable to either the single difference of carrier phase method or the double difference of carrier phase method. It is to be noted that in the same way the invention can be applicable to one of any other numbers used for “j”.
 Now, the construction of the present attitude determining system for a body is described with reference to FIGS. 3 and 4.
 FIG. 3 is a block diagram of the attitude determining system. Designated by s1, s2, . . . sk in FIG. 3 are satellites of the Global Positioning System (GPS). GPS antennas receive radio waves from the multiple satellites, convert received signals into intermediate frequency (IF) signals, and amplify the IF signals. A GPS receiver calculates antenna positions and carrier phase differences from the signals received by the individual antennas and outputs these data together with information on the positions of the individual satellites to an attitude data processing block at specific time intervals (i.e., a few tens of milliseconds to a few seconds). Since the satellite positions obtained by position fixing, ephemeris information, etc. comply with the coordinate system of the GPS, these data are converted into the local coordinate system in the GPS receiver block or in the attitude data processing block. In the following discussion, it is assumed that the data are converted in the GPS receiver block.
 A satelliteuse planning section of the attitude data processing block selects satellites to be used based on the satellite ephemeris information and the antenna positions and gives the antenna positions and information on the selected satellites to a phase difference processing section. The phase difference processing section calculates observed quantities of single (or double) phase differences from carrier phase difference signals received from the GPS receiver block and delivers the resultant data to an integer ambiguity determination section. The integer ambiguity determining section determines integer ambiguities of the single (or double) phase differences of the relevant satellites, evaluates the goodness of the integer ambiguities, and gives integer ambiguity information to a GPS attitude processing section. The GPS attitude processing section calculates baseline vectors in the local coordinate system from the integer ambiguities and carrier phase difference observables output from the phase difference processing section using a known method. Attitude angles are calculated from the baseline vectors for the local coordinate system thus obtained and the known baseline vectors for the antenna coordinate system from equations (1) to (5). The attitude angles obtained by the GPS are hereinafter referred to as “GPS attitude angles.”
 Referring to FIG. 3, an IMU sensor block includes a rate gyro for detecting angular velocities around three perpendicular axes as well as an xaxis amplifier, yaxis amplifier and zaxis amplifier for amplifying angular velocity outputs. This IMU sensor block outputs data on the angular velocities around the individual axes to an IMU attitude processing section of the attitude data processing block. More specifically, the IMU sensor block outputs angular velocities osis about axes of an IMU sensor coordinate system referenced to a righthanded inertial coordinate system whose x axis coincides with the line of intersection of the earth's equatorial plane and a meridian plane passing through the vernal equinox and y axis coincides with an eastwarddirected line perpendicular to the x axis. What are needed here are angular velocities ω^{s} _{ns }about coordinate axes of the IMU sensors in the local coordinate system. The relationship between ω^{s} _{is }and ω^{s} _{ns }is given by ω^{s} _{ns}=ω^{s} _{is}−ωs_{in}, where ω^{s} _{in }is an angular velocity obtained when an angular velocity in the local coordinate system referenced to the inertial coordinate system is measured from the IMU sensor coordinate system. Since ω^{s} _{is}>>ω^{s} _{in }in the case of a lowspeed body like a ship, ω^{s} _{is }may be disregarded. In the case of a highspeed body like an aircraft, ω^{s} _{ns }is calculated from the position of the body using a known method.
 The IMU attitude processing section calculates the attitude angles using Euler's equation or Gilmore's algorithm which are known. The attitude angles thus obtained are hereinafter referred to as “IMU attitude angles.”
 An IMU/GPS integration processing section combines the IMU attitude angles and the GPS attitude angles and outputs individual attitude angles including the angles of roll, pitch and yaw.
 FIG. 4 is a flowchart showing operations performed by a reinitialization judgment section and the IMU/GPS integration processing section of the attitude data processing block shown in FIG. 3. First, the attitude data processing block selects satellites to be used for position fixing and takes in observation data obtained by receiving radio waves from the selected satellites (step n1). If integer ambiguities are not determined yet under initial conditions, baseline vectors are estimated (real number solutions are determined) for the individual antennas. A Kalman filter or the least squares method, for instance, is used for this estimation. Subsequently, the integer ambiguities are determined by rounding the real number solutions to integers using a known method (steps n2→n3→n5→n6→n5 . . . ).
 Next, the goodness of the estimated integer ambiguities is verified by conducting a χ^{2}test or by judging whether baseline lengths coincide with actual values (step n7). If the result of verification is unsatisfactory, the aforementioned operation for determining the integer ambiguities is reexecuted (steps n8→n5 . . . ). If the result of verification is satisfactory, baseline vectors are calculated using the obtained integer ambiguities as known quantities, and GPS attitude angles are calculated from relative positions of the individual antennas (step n9). Then, IMU attitude angles are calculated from angular velocities around the three axes detected by the IMU sensors at the same time that the GPS attitude angles were calculated (step n10).
 Next, the GPS attitude angles and the IMU attitude angles are combined. Specifically, the IMU attitude angles are converted to attitude angles defined in the same local coordinate system as for the GPS attitude angles by using the GPS attitude angles as initial values of a strapdown of the IMU sensor block. Also, IMU sensor errors are estimated from differences between the GPS attitude angles and the IMU attitude angles (step n11).
 The GPS attitude angles are attitude angles of the axes of the antenna coordinate system referenced to the local coordinate system, and the IMU attitude angles are attitude angles of the axes of the IMU sensor coordinate system referenced to the local coordinate system. Although it is desirable to install the antennas and IMU sensors on the body such that the axes of the antenna coordinate system and the axes of the IMU sensor coordinate system exactly match each other, this is difficult in practice. Accordingly, a deviation of the coordinate axes is compensated for by a mathematical operation involving rotation of the coordinate axes.
 The aforementioned estimation of the IMU sensor errors is made by estimating ambiguity and drift errors which constitute the IMU sensor errors using information on the differences between the GPS attitude angles and the IMU attitude angles. The IMU sensor errors thus estimated are removed from the IMU attitude angles.
 Subsequently, the reinitialization judgment section of the attitude data processing block judges whether any cycle slip has occurred. If no cycle slip has occurred, the attitude data processing block outputs the obtained attitude angles (steps n12→n13→n15).
 When the occurrence of a cycle slip has been detected, the integer ambiguities are directly determined from IMU attitude angle data by the aforementioned method and the baseline vectors are determined from the integer ambiguities (steps n13→n14→n4). This estimation of the baseline vectors in step n4 and the estimation of the baseline vectors in step n3 differ in that the number of states in the baseline vector estimation in step n3 is three (x, y, z) plus m number of ambiguities whereas the number of states in the baseline vector estimation in step n4 is three only (x, y, z) in a case where a Kalman filter is used, for instance. Therefore, computational complexity of the baseline vector estimation in step n4 is rather small so that the baseline vectors can be estimated in an extremely short period of time.
 From this time onward, the attitude data processing block carries out the aforementioned verification of the goodness of the estimated integer ambiguities in step n7 and the later operations (steps n7→n8 . . . ) Since no error usually occurs in the integer ambiguities determined in step N14, the operation flow may proceed from step n4 to step n9, skipping steps n7 and n8.
 Since the integer ambiguities are directly determined from the IMU attitude angle data as stated above when system reinitialization has become necessary due to cycle slips or a change in the combination of satellites to be used, for instance, it is not necessary to spend a long period of time for determining the integer ambiguities. As a consequence, it is possible to output the attitude angles of the body to be used for determining the integer ambiguities substantially without interruption and thereby provide a user with highly accurate attitude angle information continuously.
Claims (4)
1. A system for determining the heading and/or attitude of a body, said system comprising:
at least three antennas which receive radio waves from a plurality of positionfixing satellites with one of the antennas being used as a reference antenna, the antennas being fixed at different positions of the body;
means for determining the relative positions of the other antennas with respect to the reference antenna by observing single differences of carrier phase or double differences of carrier phase between carrier phases of the radio waves and by calculating integer ambiguities of the phase differences to thereby determine the heading and/or attitude of the body;
an IMU sensor installed on the body;
means for obtaining attitude angles of the body based on an output of the IMU sensor; and
means for determining the integer ambiguities based on the attitude angles obtained when the integer ambiguities are to be redetermined.
2. A system for determining the heading and/or attitude of a body as defined in claim 1 wherein said means for determining the integer ambiguities comprises means for obtaining a baseline vector in the local coordinate system from the attitude angles obtained based on the output of the IMU sensor, means for obtaining the path difference with respect to the baseline vector, and means for obtaining the integer ambiguities based on the path difference with respect to the baseline vector and measured phase difference.
3. A system for determining the heading and/or attitude of a body, said system comprising:
two antennas which receive radio waves from a plurality of positionfixing satellites, the antennas being fixed at different positions of the body with one of the antennas being used as a reference antenna;
means for determining the relative positions of the other antennas with respect to the reference antenna by observing single differences of carrier phase or double differences of carrier phase between carrier phases of the radio waves and by calculating integer ambiguities to thereby determine the heading and/or attitude of the body;
an IMU sensor installed on the body;
means for obtaining attitude angles of the body based on an output of the IMU sensor; and
means for determining the integer ambiguities based on the attitude angles obtained from the output of the IMU sensor when the integer ambiguities are to be redetermined.
4. A system for determining the heading and/or attitude of a body, said system comprising:
a plurality of antennas which receive radio waves from a plurality of positionfixing satellites, the antennas being fixed at different positions of the body with one of the antennas being used as a reference antenna;
means for determining the relative positions of the other antennas with respect to the reference antenna based on integer ambiguities of the phase differences to thereby determine the heading and/or attitude of the body;
an IMU sensor installed on the body;
means for obtaining attitude angles of the body from an output of the IMU sensor; and
means for determining the integer ambiguities based on the attitude angles obtained from the output of the IMU sensor when the integer ambiguities are to be redetermined.
Priority Applications (2)
Application Number  Priority Date  Filing Date  Title 

JP2000164684  20000601  
JP2000164684  20000601 
Publications (2)
Publication Number  Publication Date 

US20020029110A1 true US20020029110A1 (en)  20020307 
US6424915B1 US6424915B1 (en)  20020723 
Family
ID=18668270
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

US09/870,628 Active US6424915B1 (en)  20000601  20010601  System for determining the heading and/or attitude of a body 
Country Status (3)
Country  Link 

US (1)  US6424915B1 (en) 
AU (1)  AU775676B2 (en) 
GB (1)  GB2368218B (en) 
Cited By (17)
Publication number  Priority date  Publication date  Assignee  Title 

US20040212533A1 (en) *  20030423  20041028  Whitehead Michael L.  Method and system for satellite based phase measurements for relative positioning of fixed or slow moving points in close proximity 
WO2005071431A1 (en) *  20040123  20050804  Novatel Inc.  Inertial gps navigation system with modified kalman filter 
US20070016371A1 (en) *  20050713  20070118  Honeywell International Inc.  Methods and systems of relative navigation for shipboard landings 
US20080004796A1 (en) *  20060630  20080103  Wolfgang Hans Schott  Apparatus and method for measuring the accurate position of moving objects in an indoor environment 
US20080221794A1 (en) *  20041207  20080911  Sagem Defense Securite  Hybrid Inertial Navigation System Based on A Kinematic Model 
US20110231133A1 (en) *  20100317  20110922  Qualcomm Incorporated  Methods and Systems for Wireless Platform Attitude Determination 
WO2012075464A1 (en) *  20101203  20120607  Qualcomm Incorporated  Inertial sensor aided heading and positioning for gnss vehicle navigation 
CN104569914A (en) *  20141229  20150429  南京航空航天大学  Motion attitude sensing method based on electromagnetic wave polarization threedimensional electric/magnetic signal stationary characteristics 
EP3012585A1 (en) *  20141020  20160427  Honeywell International Inc.  Systems and methods for attitude fault detection in one or more inertial measurement units 
US20170192103A1 (en) *  20160106  20170706  Honeywell International Inc.  Systems and methods for vehicle attitude determination 
US9719788B2 (en)  20100524  20170801  Robert J. Wellington  Determining spatial orientation information of a body from multiple electromagnetic signals 
WO2017203108A1 (en)  20160524  20171130  Robert Schegerin  Device and method for piloting an aircraft 
US9880562B2 (en)  20030320  20180130  Agjunction Llc  GNSS and optical guidance and machine control 
US10038239B2 (en) *  20160125  20180731  Sky Perfect Jsat Corporation  Antenna adjusting apparatus and antenna adjusting method 
USRE47101E1 (en)  20030320  20181030  Agjunction Llc  Control for dispensing material from vehicle 
US10197681B2 (en)  20141226  20190205  Furuno Electric Co., Ltd.  State calculating device, method of calculating state, and state calculating program 
US10514469B2 (en)  20141226  20191224  Furuno Electric Co., Ltd.  Attitude angle calculating device, method of calculating attitude angle, and attitude angle calculating program 
Families Citing this family (47)
Publication number  Priority date  Publication date  Assignee  Title 

JP4446569B2 (en) *  20000724  20100407  古野電気株式会社  Carrier phase relative positioning device 
JP3726884B2 (en) *  20010425  20051214  学校法人日本大学  Attitude estimation apparatus and method using inertial measurement apparatus, and program 
JP4116792B2 (en)  20011219  20080709  古野電気株式会社  Carrier phase relative positioning device 
US6697736B2 (en) *  20020206  20040224  American Gnc Corporation  Positioning and navigation method and system thereof 
US7885745B2 (en)  20021211  20110208  Hemisphere Gps Llc  GNSS control system and method 
JP3892401B2 (en) *  20030120  20070314  Ｈｏｙａ株式会社  Manufacturing method of disk substrate for perpendicular magnetic recording medium, and manufacturing method of perpendicular magnetic recording disk 
US8140223B2 (en)  20030320  20120320  Hemisphere Gps Llc  Multipleantenna GNSS control system and method 
US8686900B2 (en) *  20030320  20140401  Hemisphere GNSS, Inc.  Multiantenna GNSS positioning method and system 
US7689354B2 (en)  20030320  20100330  Hemisphere Gps Llc  Adaptive guidance system and method 
US8271194B2 (en)  20040319  20120918  Hemisphere Gps Llc  Method and system using GNSS phase measurements for relative positioning 
US8583315B2 (en)  20040319  20131112  Agjunction Llc  Multiantenna GNSS control system and method 
US8214111B2 (en)  20050719  20120703  Hemisphere Gps Llc  Adaptive machine control system and method 
US8138970B2 (en)  20030320  20120320  Hemisphere Gps Llc  GNSSbased tracking of fixed or slowmoving structures 
US8594879B2 (en)  20030320  20131126  Agjunction Llc  GNSS guidance and machine control 
US8190337B2 (en)  20030320  20120529  Hemisphere GPS, LLC  Satellite based vehicle guidance control in straight and contour modes 
US8265826B2 (en) *  20030320  20120911  Hemisphere GPS, LLC  Combined GNSS gyroscope control system and method 
US8954251B2 (en)  20041005  20150210  Vision Works Ip Corporation  Absolute acceleration sensor for use within moving vehicles 
US9878693B2 (en)  20041005  20180130  Vision Works Ip Corporation  Absolute acceleration sensor for use within moving vehicles 
US9327726B2 (en)  20041005  20160503  Vision Works Ip Corporation  Absolute acceleration sensor for use within moving vehicles 
US8437935B2 (en)  20041005  20130507  Vision Works Ip Corporation  Absolute acceleration sensor for use within moving vehicles 
US7239953B2 (en) *  20041005  20070703  Vision Works, Llc  Absolute acceleration sensor for use within moving vehicles 
US7310062B1 (en) *  20050728  20071218  Rockwell Collins, Inc.  Dual antenna diversity method to detect GPS signal tampering 
GB2429128B (en) *  20050808  20090311  Furuno Electric Co  Apparatus and method for carrier phasebased relative positioning 
US7835832B2 (en)  20070105  20101116  Hemisphere Gps Llc  Vehicle control system 
US8000381B2 (en)  20070227  20110816  Hemisphere Gps Llc  Unbiased code phase discriminator 
JP4976948B2 (en) *  20070725  20120718  古野電気株式会社  Attitude measurement device 
US7948769B2 (en)  20070927  20110524  Hemisphere Gps Llc  Tightlycoupled PCB GNSS circuit and manufacturing method 
US7808428B2 (en)  20071008  20101005  Hemisphere Gps Llc  GNSS receiver and external storage device system and GNSS data processing method 
WO2009100463A1 (en)  20080210  20090813  Hemisphere Gps Llc  Visual, gnss and gyro autosteering control 
US8018376B2 (en)  20080408  20110913  Hemisphere Gps Llc  GNSSbased mobile communication system and method 
US8217833B2 (en)  20081211  20120710  Hemisphere Gps Llc  GNSS superband ASIC with simultaneous multifrequency down conversion 
US8386129B2 (en)  20090117  20130226  Hemipshere GPS, LLC  Rasterbased contour swathing for guidance and variablerate chemical application 
US8085196B2 (en)  20090311  20111227  Hemisphere Gps Llc  Removing biases in dual frequency GNSS receivers using SBAS 
US8185261B2 (en) *  20090707  20120522  The Aerospace Corporation  Systems and methods for attitude propagation for a slewing angular rate vector 
US8311696B2 (en)  20090717  20121113  Hemisphere Gps Llc  Optical tracking vehicle control system and method 
US8401704B2 (en)  20090722  20130319  Hemisphere GPS, LLC  GNSS control system and method for irrigation and related applications 
US8174437B2 (en)  20090729  20120508  Hemisphere Gps Llc  System and method for augmenting DGNSS with internallygenerated differential correction 
US8334804B2 (en)  20090904  20121218  Hemisphere Gps Llc  Multifrequency GNSS receiver baseband DSP 
US8649930B2 (en)  20090917  20140211  Agjunction Llc  GNSS integrated multisensor control system and method 
US8130142B2 (en)  20090921  20120306  Appareo Systems, Llc  GNSS ultrashort baseline heading determination system and method 
US8548649B2 (en)  20091019  20131001  Agjunction Llc  GNSS optimized aircraft control system and method 
US8583326B2 (en)  20100209  20131112  Agjunction Llc  GNSS contour guidance path selection 
US9829582B2 (en) *  20110919  20171128  Raytheon Company  Method and apparatus for differential global positioning system (DGPS)based real time attitude determination (RTAD) 
US9371002B2 (en)  20130828  20160621  Vision Works Ip Corporation  Absolute acceleration sensor for use within moving vehicles 
US9855986B2 (en)  20130828  20180102  Vision Works Ip Corporation  Absolute acceleration sensor for use within moving vehicles 
US10114126B2 (en)  20150430  20181030  Raytheon Company  Sensor installation monitoring 
WO2019231345A1 (en) *  20180529  20191205  Limited Liability Company "Topcon Positioning Systems"  Gnssbased attitude determination algorithm and tripleantenna gnss receiver for its implementation 
Family Cites Families (13)
Publication number  Priority date  Publication date  Assignee  Title 

US5657025A (en) *  19950807  19970812  Litton Systems, Inc.  Integrated GPS/inertial navigation apparatus providing improved heading estimates 
AU2111797A (en) *  19960118  19970811  Monfort, Inc.  Method and apparatus for determining meat lean/fat ratio 
FR2743892B1 (en) *  19960119  19980213  Sextant Avionique  aid system to the aircraft was flying with a headup display 
US6088653A (en) *  19961231  20000711  Sheikh; Suneel I.  Attitude determination method and system 
US5757316A (en) *  19970201  19980526  Litton Systems, Inc.  Attitude determination utilizing an inertial measurement unit and a plurality of satellite transmitters 
US6154691A (en) *  19970902  20001128  Honeywell International Inc.  Orienting a satellite with controlled momentum gyros 
US6234799B1 (en) *  19980406  20010522  American Gnc Corporation  Realtime IMU simulator 
US6023988A (en) *  19980810  20000215  Case Corporation  Calibration system for a power shift transmission 
US6167347A (en) *  19981104  20001226  Lin; ChingFang  Vehicle positioning method and system thereof 
US6246960B1 (en) *  19981106  20010612  ChingFang Lin  Enhanced integrated positioning method and system thereof for vehicle 
US6178363B1 (en) *  19981222  20010123  The Boeing Company  Inertially augmented GPS landing system 
US6331835B1 (en) *  19990202  20011218  The Charles Stark Draper Laboratory, Inc.  Deeplyintegrated adaptive GPSbased navigator with extendedrange code tracking 
US6317688B1 (en) *  20000131  20011113  Rockwell Collins  Method and apparatus for achieving sole means navigation from global navigation satelite systems 

2001
 20010531 GB GB0113251A patent/GB2368218B/en active Active
 20010531 AU AU50082/01A patent/AU775676B2/en active Active
 20010601 US US09/870,628 patent/US6424915B1/en active Active
Cited By (32)
Publication number  Priority date  Publication date  Assignee  Title 

US10168714B2 (en)  20030320  20190101  Agjunction Llc  GNSS and optical guidance and machine control 
USRE47101E1 (en)  20030320  20181030  Agjunction Llc  Control for dispensing material from vehicle 
US9880562B2 (en)  20030320  20180130  Agjunction Llc  GNSS and optical guidance and machine control 
US9886038B2 (en)  20030320  20180206  Agjunction Llc  GNSS and optical guidance and machine control 
US20040212533A1 (en) *  20030423  20041028  Whitehead Michael L.  Method and system for satellite based phase measurements for relative positioning of fixed or slow moving points in close proximity 
WO2005071431A1 (en) *  20040123  20050804  Novatel Inc.  Inertial gps navigation system with modified kalman filter 
US8165795B2 (en) *  20041207  20120424  Sagem Defense Securite  Hybrid inertial navigation system based on a kinematic model 
US20080221794A1 (en) *  20041207  20080911  Sagem Defense Securite  Hybrid Inertial Navigation System Based on A Kinematic Model 
US7474962B2 (en) *  20050713  20090106  Honeywell International Inc.  Methods and systems of relative navigation for shipboard landings 
US20070016371A1 (en) *  20050713  20070118  Honeywell International Inc.  Methods and systems of relative navigation for shipboard landings 
US7761233B2 (en) *  20060630  20100720  International Business Machines Corporation  Apparatus and method for measuring the accurate position of moving objects in an indoor environment 
US20080004796A1 (en) *  20060630  20080103  Wolfgang Hans Schott  Apparatus and method for measuring the accurate position of moving objects in an indoor environment 
US8321166B2 (en)  20100317  20121127  Qualcomm Incorporated  Methods and systems for wireless platform attitude determination 
US20110231133A1 (en) *  20100317  20110922  Qualcomm Incorporated  Methods and Systems for Wireless Platform Attitude Determination 
WO2011116119A1 (en) *  20100317  20110922  Qualcomm Incorporated  Methods and systems for wireless platform attitude determination 
US9719788B2 (en)  20100524  20170801  Robert J. Wellington  Determining spatial orientation information of a body from multiple electromagnetic signals 
CN103339473A (en) *  20101203  20131002  高通股份有限公司  Inertial sensor aided heading and positioning for GNSS vehicle navigation 
WO2012075464A1 (en) *  20101203  20120607  Qualcomm Incorporated  Inertial sensor aided heading and positioning for gnss vehicle navigation 
US9803983B2 (en)  20101203  20171031  Qualcomm Incorporated  Inertial sensor aided heading and positioning for GNSS vehicle navigation 
US9816818B2 (en)  20101203  20171114  Qualcomm Incorporated  Inertial sensor aided heading and positioning for GNSS vehicle navigation 
US9891054B2 (en)  20101203  20180213  Qualcomm Incorporated  Inertial sensor aided heading and positioning for GNSS vehicle navigation 
EP3012585A1 (en) *  20141020  20160427  Honeywell International Inc.  Systems and methods for attitude fault detection in one or more inertial measurement units 
US9880021B2 (en) *  20141020  20180130  Honeywell International Inc.  Systems and methods for attitude fault detection in one or more inertial measurement units 
US20160290825A1 (en) *  20141020  20161006  Honeywell International Inc.  Systems and methods for attitude fault detection in one or more inertial measurement units 
US10197681B2 (en)  20141226  20190205  Furuno Electric Co., Ltd.  State calculating device, method of calculating state, and state calculating program 
US10514469B2 (en)  20141226  20191224  Furuno Electric Co., Ltd.  Attitude angle calculating device, method of calculating attitude angle, and attitude angle calculating program 
CN104569914A (en) *  20141229  20150429  南京航空航天大学  Motion attitude sensing method based on electromagnetic wave polarization threedimensional electric/magnetic signal stationary characteristics 
US10261194B2 (en) *  20160106  20190416  Honeywell International Inc.  Systems and methods for vehicle attitude determination 
US20170192103A1 (en) *  20160106  20170706  Honeywell International Inc.  Systems and methods for vehicle attitude determination 
US10038239B2 (en) *  20160125  20180731  Sky Perfect Jsat Corporation  Antenna adjusting apparatus and antenna adjusting method 
FR3051900A1 (en) *  20160524  20171201  Robert Schegerin  Device and method for controlling an aircraft 
WO2017203108A1 (en)  20160524  20171130  Robert Schegerin  Device and method for piloting an aircraft 
Also Published As
Publication number  Publication date 

AU5008201A (en)  20011206 
GB2368218B (en)  20040428 
GB0113251D0 (en)  20010725 
US6424915B1 (en)  20020723 
AU775676B2 (en)  20040812 
GB2368218A (en)  20020424 
Similar Documents
Publication  Publication Date  Title 

US6768959B2 (en)  Apparatus and method for accurate pipeline surveying  
US6427122B1 (en)  Positioning and data integrating method and system thereof  
US6240367B1 (en)  Full fusion positioning method for vehicle  
US8326532B2 (en)  Arrangement for and method of two dimensional and three dimensional precision location and orientation determination  
US7146740B2 (en)  Methods and apparatus for automatic magnetic compensation  
US4405986A (en)  GSP/Doppler sensor velocity derived attitude reference system  
US7346452B2 (en)  Inertial GPS navigation system using injected alignment data for the inertial system  
EP0870175B1 (en)  A zero motion detection system for improved vehicle navigation system  
KR19980070731A (en)  Method and apparatus for determining attitude using inertial measurement unit and multiple satellite transmitters  
RU2363013C2 (en)  Method of combined use of kinematic mode in real time systems and regional, wide area or carrier phase global positioning system  
US6915205B2 (en)  Apparatus for detecting location of movable body in navigation system and method thereof  
EP1839070B1 (en)  Improved integrity and continuity satellite positioning receiver  
EP0870173B1 (en)  Improved vehicle navigation system and method  
AU2004203383B2 (en)  Inertial GPS navigation system with modified kalman filter  
DE60031868T2 (en)  Navigation system and method for following the position of an object  
US5583774A (en)  Assuredintegrity monitoredextrapolation navigation apparatus  
AU2002315595B2 (en)  An inertial/GPS navigation system  
JP2005031082A (en)  Test method of integrity of gps measuring, error detection method in special vehicle, mapping method and system of gps multipath level, and system installed in vehicle for testing integrity of gps measuring  
US7961143B2 (en)  Partial search carrierphase integer ambiguity resolution  
CA2664994C (en)  Hybrid ins/gnss system with integrity monitoring and method for integrity monitoring  
US6029111A (en)  Vehicle navigation system and method using GPS velocities  
US6408245B1 (en)  Filtering mechanization method of integrating global positioning system receiver with inertial measurement unit  
US6760663B2 (en)  Solution separation method and apparatus for groundaugmented global positioning system  
Hare  Depth and position error budgets for mulitbeam echosounding  
US7937190B2 (en)  Apparatus for an automated aerial refueling boom using multiple types of sensors 
Legal Events
Date  Code  Title  Description 

AS  Assignment 
Owner name: FURUNO ELECTRIC COMPANY, LIMITED, JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:FUKUDA, MASARU;SATO, SADAO;TODA, HIROYUKI;REEL/FRAME:012206/0939 Effective date: 20010823 

STCF  Information on status: patent grant 
Free format text: PATENTED CASE 

CC  Certificate of correction  
FPAY  Fee payment 
Year of fee payment: 4 

FPAY  Fee payment 
Year of fee payment: 8 

FPAY  Fee payment 
Year of fee payment: 12 