Avionics Navigation Systems Second Edition
Avionics Navigation Systems Second Edition
10 9
CONTENTS
Preface xvii
Acknowledgments xxi
List of Contributors xxiii
1 Introduction 1
Myron Kayton
1.1 Definitions
1.2 Guidance versus Navigation
1.3 Categories of Navigation 2
1.4 The Vehicle 3
1.4.1 Civil Aircraft 3
1.4.2 Military Aircraft 5
1.5 Phases of Flight 7
1.5.1 Takeoff 7
1.5.2 Terminal Area 7
1.5.3 En route 7
1.5.4 Approach 8
1.5.5 Landing 8
1.5.6 Missed Approach 9
1.5.7 Surface 9
1.5.8 Weather 9
1.6 Design Trade-offs 9
1.7 Evolution of Air Navigation II
1.8 Integrated Avionics 15
1.8.1 All Aircraft 15
1.8.2 Military Avionics 16
1.8.3 Architecture 17
1.9 Human Navigator 19
l.l Introduction 21
2.2 Geometry of the Earth 23
v
vi CONTENTS
3.1 Introduction 55
3.2 Inertial System Characteristics 57
3.3 An Integrated Stellar-Inertial System 61
3.4 Integrated Doppler-Inertial Systems 64
3.5 An Airspeed-Damped Inertial System 67
3.6 An Integrated Stellar-Inertial-Doppler System 68
3.7 Position Update of an Inertial System 69
3.8 Noninertial GPS Multisensor Navigation Systems 69
3.9 Filtering of Measurements 70
3.9.1 Single Sensor, Stationary Vehicle 70
3. 9.2 Multiple Sensors, Stationary Vehicle 71
3.9.3 Multiple Sensors, Moving Vehicle 72
3.10 Kalman Filter Basics 72
3.1 0.1 The Process and Measurement Models 73
3.10.2 The Error Covariance Matrix 75
3.10.3 The Recursive Filter 75
3.11 Open-Loop Kalman Filter Mechanization 77
3.12 Closed-Loop Kalman Filter Mechanization 79
3.13 GPS-INS Mechanization 81
CONTENTS vii
7 .I Introduction 313
7.2 The System 314
7.3 Instruments 317
7 .3.1 Accelerometers 317
7.3.2 Gyroscopes 324
7.3.3 Optical Gyroscopes 326
7.3.4 Mechanical Gyroscopes 342
7.3 .5 Future Inertial Instruments 34 7
7.4 Platforms 348
7 .4.1 Analytic Platform (Strapdown) 348
7.4.2 Gimballed Platform 361
7 .4.3 Inertial Specifications 364
7.5 Mechanization Equations 365
7.5 .1 Coordinate Frames 365
7.5.2 Horizontal Mechanization 368
7.5.3 Vertical Mechanization 373
7.6 Error Analysis 376
7.6.1 Purpose 376
x CONTENTS
References 705
Index 741
PREFACE
The purpose of this book is to present a unified treatment of the principles and
practices of modern navigation sensors and systems. This second edition is a
total rewrite of the first edition.
During the 28 years since the first edition was published, there have been
tremendous changes in the science and practice of navigation: the introduction
of navigation satellites that provide, for the first time in history, global, con-
tinuous precise navigation; an enormous increase in the speed and memory of
digital computers, accompanied by a sharp decrease in their size and cost; the
invention of clever algorithms, based primarily on Kalman filters, that mix the
outputs of several sensors to produce a best estimate of position, velocity and,
sometimes, of time; and the proliferation of avionics on aircraft, interconnected
by digital data buses, so that navigation is only one of several avionic subsys-
tems.
This book was written for the navigation system engineer, whether user or
designer, who is concerned with the practical application of newly developed
technology, and for the technical specialist who wishes to learn about adja-
cent specialties. It is an engineer-oriented text that will serve a wide spectrum
of readers, from the systems analyst who writes mathematical models to oper-
ations personnel who want to learn about the avionics equipment in their air-
craft. This book applies to civil and military aircraft, helicopters, and unmanned
aerial vehicles. It covers the speed range from hovering helicopters to hyper-
sonic transports. For all those vehicles, it discusses the state-of-the-art and the
development of new systems that are likely to be introduced in the future.
Each chapter first presents basic functions and fundamental principles. It then
discusses design characteristics, equipment configurations, sources of error, and
typical performance levels. It closes with a projection of future trends. Topics
such as comparative performance levels, weights, and costs of equipment are
covered wherever possible. Most chapters assume a knowledge of undergradu-
ate physics and mathematics; some assume a knowledge of electronic circuits.
References are collected at the end of the book, chapter by chapter, for the
interested reader to use as background reading and to pursue the subject in
more depth. The index is comprehensive enough to allow readers to find topics
outside their area of specialty. It includes a glossary of acronyms. Chapters 2
through 15 conclude with illustrative problems that clarify points in the text
and lead the reader into new areas. These problems will be useful to university
instructors who use the text as part of a course in avionics, guidance and control,
xvii
xviii PREFACE
Seymour Schoen, T. J. Thomas, Carl Wiley and Willis Wing. The Acknowledg-
ments explain which first-edition material was re-used.
Due to increasing specialization, the second edition was written by a much
larger, more diverse team, whose members are the foremost current experts in
their fields. We wish to thank them for their generous contributions of time
in preparing drafts, editing, and re-editing in order to give you, the reader, a
coherent, unified book.
While the art and science of navigation is hundreds of years old, the last
50 years have produced exciting new sensors and systems that permit an accu-
racy and level of safety never before seen on moving vehicles. We hope that
this second edition presents the fundamentals and enough details to stimulate
innovation and the development of ever-improving systems of navigation.
MYRON KA YTON
Santa Monica, California
WALTER R. FRIED
Santa Ana, Calijbrnia
January 1997
ACKNOWLEDGMENTS
xxi
xxii ACKNOWLEDGMENTS
WALTER R. FRIED, M.S. (Editor; lead author, Chapters 6 and 10), Consultant,
Hughes Aircraft Company, Santa Ana, CA
JAMES R. HUDDLE, Ph.D. (Lead author, Chapter 3), Chief Scientist and Head
of Advanced System Engineering, Litton Guidance and Control Division,
Woodland Hills, CA
JOHN G. MARK, Ph.D. (Co-author, Chapter 7), Chief Scientist, Litton Guid-
ance and Control Division, Woodland Hills, CA
RICHARD H. MCFARLAND, Ph.D., P.E. (Co-author, Chapter 13), Director,
Emeritus, Avionics Engineering Center. Ohio University, Athens, OH
CLYDE A. MILLER, Ph.D. (Lead author, Chapter 14), Program Director for
Research, Federal Aviation Administration, Washington, DC
PETER MORRIS (Co-author, Chapter 4 ), The Analytical Sciences Corporation,
Reading, MA
STEPHEN S. OSDER (Author, Chapter 8), Consultant, formerly McDonnell-
Douglas Fellow, Scottsdale, AZ
JACK 0. PEARSON, Ph.D. (Lead author, Chapter 11), Vice President, Radar
and Communication Systems, Hughes Aircraft Company, El Segundo, CA
JOHN A. SCARDINA, Ph.D. (Co-author, Chapter 14), Team Leader for Air Traf-
fic Management, Federal Aviation Administration, Washington, DC
CARY R. SPITZER (Author, Chapter 15), President, AvioniCon, formerly,
National Aeronautics and Space Administration, Williamsburg, VA
DANIEL A. T AZARTES (Lead author, Chapter 7), Senior Member of Technical
Staff, Litton Guidance and Control Division, Woodland Hills, CA
BAHAR UTTAM, Ph.D. (Lead author, Chapter 4), President, Synetics Corpo-
ration, Wakefield, MA
A. J. VAN DIERENDONCK, Ph.D. (Author, Chapter 5), AJ Systems, Los Altos,
CA
D. B. VICKERS, M.S. (Lead author, Chapter 13), Technical Director, Avionics
Engineering Center. Ohio University, Athens, OH
WILLIAM M. WATERS, Ph.D. (Co-author, Chapter 13), Senior Consultant,
Naval Research Laboratory, Washington, DC
EDGAR A. WESTBROOK, (Co-author, Chapter 6), Technical Staff, retired, The
MITRE Corporation, Bedford, MA
WILLIS G. WING, (Co-author, Chapter 9), Sperry Gyroscope Company,
retired, Glen Head, NY
1 Introduction
1.1 DEFINITIONS
The term "guidance" has two meanings, both of which are different from "nav-
igation":
C- Band SATCOM
itary usage. The avionics bay is below the cockpit in the space between the
radome and nose wheel well (in many civil aircraft, the avionics bay is aft of
the nosewheel). Avionics and air-data sensors are located in the bay. Access is
beneath the aircraft.
Trans-Pacific hypersonic aircraft may be developed in the twenty-first cen-
tury that will navigate as does the Space Shuttle: inertial boost, GPS or celestial
midcourse, and GPS or other radio approach. They will compete with electronic
mail and teleconferences.
shock and vibration loads. By the year 2000, most military aircraft will carry
GPS receivers.
1.5.1 Takeoff
The takeoff phase begins upon taxiing onto the runway and ends when climb-
out is established on the projected runway centerline. The aircraft is guided
along the centerline by hand-flying or a coupled autopilot based on steering
signals (from an ILS localizer since 1945). Two important speed measurements
are made on the runway. The highest ground speed at which an aborted takeoff
is possible is precomputed and compared, during the takeoff run, to the actual
ground speed as displayed by the navigation system. The airspeed at which the
nose is lifted ("rotation") is precalculated and compared to the actual airspeed
as displayed by the air-data system. Barometric altitude rate or GPS-derived
altitude rate (inertially smoothed) is measured and monitored.
1.5.3 En Route
The en-route phase leads from the origin to the destination and alternate des-
tinations (an alternate destination is required of civil aircraft operating under
instrument flight rules). From the 1930s to the 1990s, airways were defined
by navigation aids over land and by latitude-longitude fixes over water. The
width of airways and their lateral separation depended on the quality of the
defining navaids and the distance between them. The introduction of inertial
navigation systems and DME in the 1970s caused aviation authorities to cre-
ate "area-navigation" airways (RNAV) that do not always interconnect VOR
navaids [7: AC-90-45A] (see Section 2.7.4).
Beginning in the 1990s, GPS has allowed precise navigation anywhere, not
just on airways. Given the extensive use of on-board collision-avoidance equip-
8 INTRODUCTION
ment and the trend toward reducing government budgets, "free-flight" is being
introduced in controlled airspace. Each aircraft would agree on a route before
takeoff and then be free to change the route, after interaircraft communication
verified that the risk of collision is sufficiently low. En-route surveillance by
independent ground-based radars may disappear or be replaced by position fixes
and reports via com-nav satellites. Busy terminal areas and airport surfaces are
likely to remain under central, positive control.
In the United States in 1996, the en-route navigation error must be less
than 2.8 nmi over land and 12 nmi over oceans (2-sigma) [7: AC-20-130]. As
regional maps become available in digital form with aeronautical annotation
(e.g., minimum en-route altitude), aircraft in undeveloped areas will use GPS
for en-route navigation and nonprecision approaches.
1.5.4 Approach
The approach phase begins at acquisition of the landing aid and continues until
the airport is in sight or the aircraft is on the runway, depending on the capa-
bilities of the landing aid (Chapter 13 ).
During an approach, the decision height (DH) is the altitude above the run-
way at which the approach must be aborted if the runway is not in sight. The
better the landing aids, the lower the decision height. Decision heights are pub-
lished for each runway at each airport (Chapter 13). The decision height for a
Category III landing is I 00 ft or less. By law, an approach may not even be
attempted unless the horizontal visibility, measured by a runway visual range
(RVR) instrument, exceeds a threshold that ranges from zero (Category IIIC,
not approved anywhere in 1996) to 800 meters (Category I). A nonprecision
approach has electronic guidance only in the horizontal direction. An aircraft
executing a nonprecision approach must abort if the runway is not visible at
the minimum descent altitude, which is typically 700 ft above the runway.
In 1996, civil aircraft outside the ex-Soviet bloc used the Instrument Land-
ing System for low-ceiling, low-visibility approaches. A Microwave Landing
System had been approved by the International Civil Aviation Organization for
precision approaches and was being installed at major international airports,
especially in Europe. In the United States, Loran and GPS had been approved
for nonprecision approaches at many airports. (Landing aids are discussed in
Chapter 13.)
1.5.5 Landing
The landing phase begins at the decision height (when the runway is in sight)
and ends when the aircraft exits the runway. Navigation during flare and decrab
may be visual or the navigation set's electrical output may be coupled to an
autopilot. A radio altimeter measures the height of the main landing gear above
the runway for guiding the flare. The rollout is guided by the landing aid (e.g.,
the ILS localizer). Landing navigation is described in Chapter 13.
DESIGN TRADE-OFFS 9
1.5.7 Surface
Aircraft movement from the runway to gates, hangars, or revetments is a major
limit on airport capacity in instrument meteorological conditions. Surface navi-
gation is visual on the part of the crew, whereas the ground controllers observe
aircraft visually or with a surface surveillance radar. No matter how good
the surface navigation, collision avoidance among aircraft and ground vehi-
cles requires central guidance, typically provided by a human controller with
computer assistance. Position reports (e.g., via GPS) from aircraft that are con-
cealed in radar shadows reduce the risk of collision and help keep unwanted
aircraft off active runways.
1.5.8 Weather
Instrument meteorological conditions (IMC) are weather conditions in which
visibility is restricted, typically less than 3 miles as defined by law. Aircraft
operating in IMC are supposed to fly under instrument flight rules (IFR), defined
by law in each country (Chapters 13 and 14).
The navigation-system designer conducts trade-offs for each aircraft and mis-
sion to determine which navigation systems to use. Trade-offs consider the fol-
lowing attributes:
The earliest aircraft were navigated visually. Pilots had an anemometer for air-
speed, a barometer for altitude, and a magnetic compass for heading. Artificial
horizons and turn-and-bank indicators allowed pilots to hold attitude and head-
ing in clouds, hence motivating the installation of navigation aids. Lighted bea-
cons were installed across the United States in the 1920s to mark airmail routes.
Starting in I 929, four-course radio beacons were also added to the lighted air-
ways to guide aircraft. Four-course beacons were installed in France, South
America, and North Africa. In the 1930s, aircraft were equipped with medium-
frequency and high-frequency direction finders (MF/DF and HF jDF) that mea-
12 INTRODUCTION
sured the bearing of broadcast stations relative to the axis of the aircraft. A fix
was obtained by plotting the direction toward two or more stations. Beacons
near an airport allowed aircraft to fly a "nonprecision approach" to the runway
(Chapter 13). Vertical beacons at 75 MHz, called z-beacons or marker beacons,
were installed along the four-course airways and along approaches to runways
to give a positive indication of position (Chapters 13 and 14).
Air-traffic control was procedural, following the precedent of railroad
"block" clearances. Overland airways connected radio beacons; overwater air-
ways were defined on a map, hundreds of miles apart. The airways were divided
into longitudinal blocks of 20- to 30-minutes flying time. The air-traffic con-
troller relied on the pilots' report of position, allowed only one aircraft at a time
to enter a block, and kept the block free of other traffic until the pilot reported
leaving. The size of the block was commensurate with the uncertainties in nav-
igation at the time.
During World War II, meteorologists learned to route aircraft to take advan-
tage of the cyclonic winds that eire le around high- and low-pressure regions at
mid- and high-latitudes. Bellamy [ 12] states that the transatlantic flying time
was reduced an average of I Oo/c compared to a great-circle track, with occa-
sional savings exceeding 25%, by taking advantage of cyclonic tail winds.
These pressure-pattern routes were plotted graphically in the 1940s-1960s but
are now computed routinely in airline and military dispatch offices.
Crosswinds cause an aircraft to "drift" perpendicularly to its longitudinal
axis. From the 1930s to the 1960s, drift angle was measured in flight with a
downward-looking telescope that observed the direction of movement of the
ground, when it was visible. From the 1940s to the 1960s, drift was estimated
over oceans by observing trends in the difference, D, between the readings
of the radio altimeter and pressure altimeter. Bellamy showed that in cyclonic
winds, drift is proportional to the horizontal gradient of D [12]. The introduc-
tion of Doppler and inertial navigators in the 1960s and 1970s allowed drift
to be observed directly. The Doppler navigator measures the direction of the
ground-speed vector relative to the aircraft's centerline. The inertial navigator
subtracts in-flight-measured airspeed from the measured ground velocity to cal-
culate wind, hence lateral drift.
After World War II, VOR stations (Chapter 4) and Instrument Landing Sys-
tems (ILS, Chapter 13) were installed. VOR/DME and ILS have been the basis
of navigation in western countries ever since. During the 1960s, air-traffic con-
trollers came to rely on surveillance radar in densely populated airspace (Chap-
ter 14). The controller identified the aircraft on his screen, hence eliminating
the need for a position report from the crew. Radar surveillance of air traffic
is called "positive control," which, in 1996, existed in the United States, most
of Canada, western Europe, and Japan. In the late 1990s, the automatic report-
ing of on-board-derived position began to supplement (perhaps eventually to
replace) radar surveillance.
The former Soviet republics have ICAO navigation aids and ILS at about 50
international airports and on corridors connecting them to the borders. Over-
EVOLUTION OF AIR NAVIGATION 13
flying western aircraft navigate inertially and with Omega, GPS, and nondirec-
tional beacons. Since the late 1960s, domestic civil and military aircraft have
used an L-band range-angle system known by its Russian acronym, RSBN, and
not standardized by ICAO. It has 176 channels between 873 and 1000 MHz.
Domestic airports guide landing aircraft with ground-based precision approach
radar (PAR) using verbal commands to the crew. At international airports, PARs
monitor aircraft on ILS approaches. In the 1990s, the former Soviet republics
were purchasing western avionics equipment.
The People's Republic of China depended on imported Russian nondirec-
tional beacons and PARs until the late 1970s, when it began to install western
radars, ILS, VOR, and DME. In the 1990s, China installed VHF air-to-ground
radio relays throughout most of the nation [ 13]. In 1996, western air-traffic con-
trol and navigation equipment was being installed throughout Southeast Asia
and Indonesia.
Outside the developed world, major cities and some airways had
VOR/DME-based procedural traffic control, so aircraft filing flight plans could
be separated from each other by human controllers. Polar areas, the South
Atlantic Ocean, and much of the Pacific and Indian Oceans had no navaids and
no control whatsoever. Most of the rest of the world was divided into Flight
Information Regions that advised crews of weather conditions and the status of
airports and navigation aids but did not separate traffic. Position reports over
oceans and in remote areas are mostly by HF radio but, beginning in the 1990s,
were being made via satellite (e.g., North Pacific and Atlantic Oceans). In 1996,
a few airlines were transmitting GPS-inertial position over digital data links via
geostationary communication satellites over the Pacific Ocean, a system called
Automatic Dependent Surveillance, the first step in the Future Air Navigation
System (FANS, Chapter 14). Outside the United States and Canada, most air-
craft pay directly for traffic control services.
Until the 1970s, precise absolute time could not economically be measured
on a vehicle. Hence, radio navigation aids were built that measured the dif-
ference in time of arrival of radio signals from ground stations. The earliest
(hyperbolic Loran and Decca, some military systems) date from the 1940s. As
airborne clocks became more stable in the 1970s, "passive" or "one-way" rang-
ing systems could solve for position and the absolute clock offset by processing
precisely timed signals from several stations. Direct-ranging Loran and Omega
(as distinguished from hyperbolic Loran and Omega, all discussed in Chapter
4), GPS and GLONASS (Chapter 5), and JTIDS (Chapter 6) are examples of
such one-way ranging systems. As airborne clocks become more accurate in
the twenty-first century, absolute time of arrival will be directly measurable
and clock offsets will become negligible.
GPS and GLONASS are based on one-way passive range measurements
to several stations, most of which are spacecraft (Figure 1.2). A few stations
are ground-based pseudolites whose transmissions mimic those of spacecraft.
Chapter 5 describes the GPS and GLONASS systems. The receiver in the air-
plane computes position, velocity, the offset in the airborne clock, and, in some
14 INTRODUCTION
Figure 1.2 Global Positioning System Spacecraft, Block IIF (courtesy of Rockwell).
L-band antenna array, S-band control antenna, and solar array are v isible.
receivers, the ionospheric delay (Chapter 5). In 1996, the military modes of
GPS achieved 20-meter (2dnns) accuracy anywhere in the world , while the civil
mode could achieve 40-meter accuracy but was intentionally degraded to I00-
meters, a handicap to civil navigation that may be discontinued before the year
2000. The United States and Russia have announced that GPS and GLONASS
wi II be avai Iable worldwide , free of charge. for a least 15 years and there-
after with 6 years ' warning of the end of service. Nevertheless, worldwide civil
authorities are reluctant to rely on military -controlled navigation aids that might
be switched oil or degraded during hostilities. The advent of continuous GPS
allows the use of AHRS-quality inertialjattitude-ret'erence systems (Chapters
7 and 9) in all but the most demanding military applications. The undetected
loss of a navigation signal or the failure of a receiver could be catastrophic.
especially during a landing at low decision height.
A widespread method of improving GPS accuracy and monitoring the signals
is to install a ground station that receives GPS signals and transmits position
errors or ranging errors and satellite failure status on a radio link to nearby air-
craft. This diffi' rentiill CPS (DGPS) can achieve centimeter accuracy for fixed
observers and 1- to 5-meter accuracy on aircraft that can solve at tens of itera-
tions per second or whose velocity calculations are smoothed by an inertial nav-
igator. The United States was experimenting with a nationwide DGPS system
INTEGRATED AVIONICS 15
3. Engine control. This is the electronic control of engine thrust, often called
throttle management. Afterburner and thrust reversers may be controlled
manually, perhaps via a thrust-by-wire control system.
4. Flight management. This subsystem stores the coordinates of en-route
waypoints and calculates the steering signals to fly toward them. It cal-
culates climb and descent profiles that may be followed with or without
constraints on the time at which designated fixes and altitudes are crossed.
Crossing fixes at predetermined times and altitudes is sometimes called
four-dimensional navigation; it requires that the flight management sub-
system control engine thrust. In 1996 all flight management subsystems
stored waypoints in digital form. By the year 2000 many will store dig-
ital maps of the en-route airspace, standard approaches (called STARs in
the United States), standard departures (called S!Ds in the United States),
approach plates, and checkli~;ts (see Chapter 14 ).
5. Subsystem monitoring and control. Faults in all subsystems are displayed,
as are recommended actions to be taken. This subsystem includes wired
logic and software for the automatic reconfiguration of faults in time-criti-
cal subsystems (e.g., flight control, where a fault can destroy the aircraft
in less than three seconds). Quick-responses to safety-critical faults were
automated in flight-control systems by the 1980s. In the 1990s, the trend
was to automate the responses to slower-acting faults, thus reducing the
workload in one-pilot and two-pilot aircraft. The failure-monitoring sub-
system may include an on-board maintenance recorder, the radio trans-
mission of faults to reduce repair time, and an accident recorder whose
data survive a crash (required by law on many aircraft).
6. Collision avoidance. This subsystem predicts impending collisions with
other aircraft or the ground and recommends an avoidance maneuver
(Chapter 14 ).
7. Weather detection. This sub~.ystem observes weather ahead of the aircraft
so that the route of flight can be altered to avoid thunderstorms and areas
of high wind-shear. The sensors are usually radars (Chapter II) and lasers
(Chapter 8).
8. Emergency locator transmitter ( ELT). This subsystem is triggered auto-
matically on high-g impact or manually. In 1996, ELTs emit distinc-
tive tones on 121.5, 243, and 406 MHz [8]. These frequencies (and per-
haps soon 1.6 GHz) are monitored by search-and-rescue aircraft and by
SARSAT-COSPAS satellites.
I. Radar, infrared, and other target sensors. These may have their own dis-
plays and controls or may share multipurpose devices.
2. Weapon management
• Fire control. Calculates lead angle for aiming guns and unguided rockets
at other aircraft and at ground targets.
• Stores management, that initializes and launches guided weapons: mis-
siles and bombs.
3. Electronic countermeasures. This subsystem detects, locates, and iden-
tifies enemy emitters of electromagnetic radiation. It may also generate
jamming signals. In 1996, electronic countermeasures were often so com-
plex that they were installed in an externally carried pod on specially
equipped aircraft.
4. Mission planning. Pre-flight mission planning is usually done at the air-
base by a computer that prepares coordinated flight plans for an entire
squadron. On-board software replans routes through enemy defenses
based on en-route observations. En-route replanning requires on-board
digital maps of the terrain and the real-time detection of enemy radars.
5. Formation flight. This subsystem maintains formation flight in instrument
meteorological conditions. It once consisted of beacons, transponders, and
communication links but is being replaced by relative GPS.
1.8.3 Architecture
Before the 1960s, electrical and electronic systems on aircraft consisted of inde-
pendent subsystems, each with its own sensors, analog computers, displays, and
controls. The appearance of airborne digital computers in the 1960s created the
first integrated avionic systems. The interconnectivity of airborne electronics is
called architecture. It involves six aspects:
l. Displays. They present information from the avionics to the pilots (Chap-
ters 9 and 15). The information consists of vertical and horizontal navi-
gation data, flight-control data (e.g., speed and angle of attack), and com-
munication data (radio frequencies). The displays show the status of all
subsystems including their faults. Displays consist of dedicated gauges,
dedicated glass displays, multipurpose glass displays, and the support-
ing symbol generators. In 1996, flat-panel vertical- and horizontal-situa-
tion displays were displacing cathode-ray tubes as "glass displays." Mul-
tipurpose displays of text and block diagrams are flat-panel matrices sur-
rounded by buttons whose labels change as the displays change. On-board
digital terrain data, used for mission planning, can be displayed on the
horizontal situation display or on a head-up display.
2. Controls. The means of inputting information from the pilots to the
avionics. The flight controls traditionally consist of rudder pedals and a
18 INTRODUCTION
wires and by using fiberoptic data buses. Aircraft constructed with a con-
tinuous metal skin have an added layer of Faraday shielding. Neverthe-
less, direct lightning strikes on antennas destroy input circuits and may
damage feed cables. A nearby strike may induce enough current to do the
same. Composite airframe structures can be transparent to radiation, thus
exposing the avionics and power systems to external fields.
7. Standards. The signals in space created by navaids are standardized by
the International Civil Aviation Organization (ICAO), Montreal, a United
Nations agency [3]. These standards are written by committees that con-
sist of representatives of the member governments. Interfaces among
airborne subsystems, within the aircraft, are standardized by ARINC
(Aeronautical Radio, Inc.), Annapolis, Maryland, a nonprofit organiza-
tion owned by member airlines [1]. Other requirements are imposed
on airborne equipment by two nonprofit organizations supported by
member entities (mostly airframe and avionics manufacturers and gov-
ernment agencies). In the United States, RTCA, Inc. (Formerly Radio
Technical Commission for Aeronautics), Washington D.C., defines the
environmental specifications and test procedures for airborne hardware
and software, and writes performance specifications for airborne equip-
ment [5]. In Europe, EUROCAE (European Organisation for Civil Avi-
ation Equipment), Paris, produces specifications for airborne equipment,
some of which are in conjunction with RTCA [II]. Government agen-
cies in all major nations define rules governing the usage of naviga-
tion equipment in flight, weather minimums, traffic separation, ground
equipment required, pilot training requirements, and so on [6-1 0] (see
Chapters I 3, I 4). Some of these rules are standardized internationally
by ICAO. U.S. military organizations once issued their own standards for
airborne circuit boards but have accepted civil standards since the early
1990s.
Large aircraft often had (and a few still had in 1996) a third crew member, the
flight engineer, whose duties were to operate engines and aircraft subsystems
such as air conditioning and hydraulics. Aircraft operating over oceans once
carried a human navigator who used celestial fixes, whatever radio aids were
available, and dead reckoning to plot the aircraft's course on a paper chart (some
military aircraft still do). Those navigators were trained in celestial observato-
ries to recognize stars, take fixes, compute position, and plot the fixes.
The navigator's crew station disappeared in civil aircraft in the 1970s,
because inertial, Doppler, and radio equipment came into use that automati-
cally selected stations, calculated position, calculated waypoint steering, and
accommodated failures. Hence, instead of requiring a skilled navigator on each
20 INTRODUCTION
2.1 INTRODUCTION
The navigation equations describe how the sensor outputs are processed in the
on-board computer in order to calculate the position, velocity, and attitude of the
aircraft. The navigation equations contain instructions and data and are part of
the airborne software that also includes moding, display drivers, failure detec-
tion, and an operating system, for example. The instructions and invariant data
are usually stored in a read-only memory (ROM) at the time of manufactur-
ing. Mission-dependent data (e.g., waypoints) are either loaded from a cock-
pit keyboard or from a cartridge, sometimes called a data-entry device, into
random-access memory (RAM). Waypoints are often precomputed in a ground-
based dispatch or mission-planning computer and transferred to the flight com-
puters.
Figure 2.1 is the block diagram of an aircraft navigation system. The system
utilizes three types of sensor information (as explained in Chapter 1):
1. Absolute position data from radio aids, radar checkpoints, and satellites
(based on range or differential range measurements).
2. Dead-reckoning data, obtained from inertial, Doppler, or air-data sen-
sors, as a means of extrapolating present position. A heading reference is
required in order to resolve the measured velocities into the computational
coordinates.
3. Line-of-sight directions to stars, which measure a combination of position
and attitude errors (as explained in Chapters I and 12).
Heading
attitude
~--------------------------------------~
I
To cockpit displays
pointing sensors
'?-
I i I
I
I
Dead-
Inertial }
-.:::.. reckoning
----------· I
I
I
air data I
computations I
I
Doppler I
I
Position I Waypoints
I
velocity I
attitude · ·
Most
l
I
I
t ~ Time to go
Star line Celestial __... Position
probable Course
of sight
...::.
equations
position
-= computations ~ Range, bearing to displays, FMS
1---1
computations ~ Steering signals to autopilot
Position in ..., To map display
sensors ~
..:::.. To weapon computers
}
Radio (VOR, Velocity
Loran, Orne Positioning
Satellite (G
~~ computations J Pc
Radar data
where !1 is the inertial angular velocity of the Earth (15.04107 degjhr) and
R is the radius vector from the mass center of the Earth to a point where the
field is to be computed. The direction of g is the "plumb bob," or "astronomic"
vertical [ 10].
In cooling from a molten mass, the Earth has assumed a shape whose surface
is a gravity equipotential and is nearly perpendicular to g everywhere (i.e., no
horizontal stresses exist at the surface). For navigational purposes, the Earth's
surface can be represented by an ellipsoid of rotation around the Earth's spin
axis. The size and shape of the best-fitting ellipsoid are chosen to match the
sea-level equipotential surface. Mathematically, the center of the ellipsoid is
at the mass center of the Earth, and the ellipsoid is chosen so as to minimize
the mean-square deviation between the direction of gravity and the normal to
the ellipsoid, when integrated over the entire surface. National ellipsoids have
been chosen to represent the Earth in localized areas, but they are not always
good worldwide approximations. The centers of these national ellipsoids are
not exactly coincident and do not exactly coincide with the mass center of the
Earth [9]. In 1996, the World Geodetic System (WGS-84, [20]) was the best
approximation to the geoid, based on gravimetric and satellite observations.
Reference [20] contains transformation equations that convert between WGS-84
and various national ellipsoids. The differences are typically hundreds of feet,
though some isolated island grids are displaced as much as a mile from WGS-
84. The navigator does not ask that the Earth be mapped onto the optimum
ellipsoid. Any ellipsoid is satisfactory for worldwide navigation if all points on
Earth are mapped onto it.
The geometry of the ellipsoid is defined by a meridian section whose semi-
major axis is the equatorial radius a and whose semiminor axis is the polar
radius b, as shown in Figure 2.2. The eccentricity of the elliptic section is
defined as e= Va 2 - b 2 /a and the ellipticity, or flattening, asf =(a- b)/a.
The radius vector R makes an angle F c with the equatorial plane, where F c
is the geocentric latitude; R and F c are not directly measurable, but they are
sometimes used in mechanizing dead-reckoning equations.
The geodetic latitude F T of a point is the angle between the normal to the ref-
erence ellipsoid and the equatorial plane. Geodetic latitude is our usual under-
standing of map latitude. The term "geographic latitude" is sometimes used syn-
onymously with "geodetic" but should refer to geodetic latitude on a worldwide
ellipsoid.
The radii of curvature of the ellipsoid are of fundamental importance to dead-
~ Polar axis of Earth
and axis of ellipsoid
Surface of Earth
tA = astronomic latitude of P
tT = geodetic latitude of P
Nav.1on;ar. 4>c = geocentric latitude of P
=
PC= h height above
reference ellipsoid
OE ==a= semimajor axis
=
OD = b semiminor axis
Center of ellipsoid
and mass
center of Earth tc
0
AvI
Figure 2.2 Meridian section of the Earth, showing the reference ellipsoid and gravity field.
GEOMETRY OF THE EARTH 25
The prime radius of curvature, Rp, is the radius of the best-fitting circle to a
vertical east-west section of the ellipsoid:
Rp = a "" a [ I + -d- sm
. 2 F T] (2.3)
(l-e2sin 2 FT) 112 2
The Gaussian radius r~f curvature is the radius of the best-fitting sphere to the
ellipsoid at any point:
(2.4)
The radii of curvature are important, because they relate the horizontal com-
ponents of velocity to angular coordinates, such as latitude and longitude; for
example,
· Veast
lcos FT = - - - (2.5)
Rp +h
where h is the aircraft's altitude above the reference ellipsoid, measured along
the normal to the ellipsoid (nearly along the direction of gravity), and 1 is its
longitude, measured positively east.
For numerical work a= 6378.137 km = 3443.918 nmi,f = 1/298.2572, and
d- = f(2 -f) [18]. (One nautical mile = 1852 meters exactly, or 6076.11549
ft.) The angle between the gravity vector and the normal to the ellipsoid, the
deflection qf the vertical, is commonly less than 10 seconds of arc and is rarely
greater than 30 seconds of arc [18]. The magnitude of gravity at sea level is
within 0.02 emf sec 2 . It decreases 10- 6 g for each 10-ft increase in altitude above
sea level [ 161.
26 THE NAVIGATION EQUATIONS
The position, velocity, and attitude of the aircraft must be expressed in a coordi-
nate frame. Paragraph l below describes the rectangular Earth-centered, Earth-
fixed (ECEF) coordinate frame, y; Paragraph 2 describes the Earth-centered
inertial (ECI) coordinate frame, Xi, which simplifies the computations for iner-
tial and stellar sensors. Other Earth-referenced orthogonal coordinates, called
Z;, can simplify navigation computations for some navaids and displays. Para-
graphs 3-5 describe coordinates commonly used in inertial navigation systems.
Paragraph 6 describes coordinates used in land navigators or in military air-
craft that support ground troops. Paragraphs 7-9 describe coordinates that were
important before powerful airborne digital computers existed.
Y3
North Pole
Geodetic wander-azimuth
coordinates
the small relative motions of stars ("proper motion") and have defined
an "average" ECI reference frame [11]. To an accuracy of w-s degjhr,
an ECI frame can be chosen with its x 3 -axis along the mean polar axis
of the Earth and with its x 1- and x 2 -axes pointing to convenient stars (as
explained in Chapter 12). ECI coordinates have three navigational func-
tions. First, Newton's laws are valid in any ECI coordinate frame. Sec-
ond, the angular coordinates of stars are conventionally tabulated in ECI.
Third, they are used in mechanizing inertial navigators, Section 7.5.1.
3. Geodetic spherical coordinates. These are the spherical coordinates of
the normal to the reference ellipsoid (Figure 2.2). The symbol z 1 repre-
sents longitude 1; z2 is geodetic latitude F r, and z3 is altitude h above
the reference ellipsoid. Geodetic coordinates are used on maps and in the
mechanization of dead-reckoning and radio-navigation systems. Transfor-
mations from ECEF to geodetic spherical coordinates are given in [9] and
[23].
4. Geodetic wander azimuth. These coordinates are locally level to the ref-
erence ellipsoid. 23 is vertically up and z2 points at an angle, a, west of
true north (Figure 2.3). The wander-azimuth unit vectors, z 1 and z2, are
in the level plane but do not point east and north. Wander azimuth is the
most commonly used coordinate frame for worldwide inertial navigation
and is discussed below and in Section 7.5.1.
5. Direction cosines. The orientation of any z-coordinate frame (e.g., navi-
gation coordinates or body axes) can be described by its direction cosines
relative to ECEF y-axes. Any vector V can be resolved into either the y-
or z-coordinate frame. The y and z components of V are related by the
equation
(2.7)
where
The navigation computer calculates in terms of the Cu, which are usable
28 THE NAVIGATION EQUATIONS
sin F = C33,
or
tan a = - -
c13 (2.9)
C23
wherever they converge. In polar regions, where a and l are not mean-
ingful, the navigation system operates correctly on the basis of the C;i.
Section 7.5.1 describes an inertial mechanization in direction cosines. If
the z-coordinate frame has a north-pointing axis, a = 0.
6. Map-grid coordinates. The navigation computer can calculate position in
map-grid coordinates such as Lambert conformal or transverse Mercator
xy-coordinates [13]. Grid coordinates are used in local areas (e.g., on mil-
itary battlefields or in cities) but are not convenient for long-range naviga-
tion. A particular grid, Universal Transverse Mercator (UTM), is widely
used by army vehicles of the western nations. The U.S. Military Grid
Reference System (MGRS) consists of UTM charts worldwide except,
in polar regions, polar stereographic charts [ 13]. The latter are projected
onto a plane tangent to the Earth at the pole, from a point at the opposite
pole.
7. Geocentric spherical coordinates. These are the spherical coordinates of
the radius vector R (Figure 2.2). The symbol z1 represents longitude
l; Z2 is geocentric latitude F c, and z3 is the radius. Geocentric coordi-
nates are sometimes mechanized in short-range dead-reckoning systems
using a spherical-Earth approximation. Initialization requires knowledge
of the direction toward the mass center of the Earth, a direction that is
not directly observable.
8. Transverse-pole spherical coordinates. These coordinates are analogous
to geocentric spherical coordinates except that their poles are deliberately
placed on the Earth's equator. The symbol z1 represents the transverse
longitude; Z2, the transverse latitude; and Z3, the radius. They permit non-
singular operation near the north or south poles, by placing the transverse
pole on the true equator. Transverse-polar coordinates involve only three
z; variables instead of nine direction cosines. However, they cannot be
used for precise navigation, since the transverse equator is elliptical, com-
DEAD-RECKONING COMPUTATIONS 29
Dead reckoning (often called DR) is the technique of calculating position from
measurements of velocity. It is the means of navigation in the absence of posi-
tion fixes and consists in calculating the position (the Zi-coordinates) of a vehi-
cle by extrapolating (integrating) estimated or measured ground speed. Prior to
GPS, dead-reckoning computations were the heart of every automatic naviga-
tor. They gave continuous navigation information between discrete fixes. In its
simplest form, neglecting wind, dead reckoning can calculate the position of a
vehicle on the surface of a flat Earth from measurements of ground speed V R
and true heading wy:
where x- x 0 and y- y0 are the east and north distances traveled during the mea-
surement interval, respectively. Notice that a simple integration of unresolved
ground speed would give curvilinear distance traveled but would be of little
use for determining position.
Aircraft heading (best-available true heading) is measured using the quanti-
ties defined in Figure 2.4. With a magnetic compass, for example, the best avail-
able true heading is the algebraic sum of magnetic heading and east variation.
With a gimballed inertial system, the best available true heading is platform
heading (relative to the Zi computational coordinates) plus the wander angle
a (Section 7.5.2). When navigating manually in polar regions, dead-reckoned
velocity is resolved through best available grid heading.
30 THE NAVIGATION EQUATIONS
True
north Grid
1' north Magnet1c
j1. north
(D1rection of 1
:2b~:~:r-t --- - ~~
1, Tr = true track
'Pr =true headmg
f \ ,J Gri?tion .,./
Vanation I j ~ J'b
~_j_/ ~ ~'I$'? V w = wind speed
I o~
I / Tr ~,-4i
I. . / j I
..:.~
Dntt
angle, o
True f -- i 7-Y':r \ Sideslip
angle,~
beanng, BT 1 1
1
True east
------ ---
( a)
V E = Earth speed
vg = ground speed = horizontal
(b) component of Vr;
YA =air-mass flight
path angle
YT = Earth-referenced
flight path angle
In the presence of a crosswind the ground-speed vector does not lie along
the aircraft's center line but makes an angle with it (Figure 2.4). The true-o
track angle T 1 , the angle from true north to the ground-speed vector, is pre-
ferred for dead-reckoning calculations whe n it is available. The drift angle 0 can
be measured with a Doppler radar or a drift sight (a downward-pointing tele-
DEAD-RECKONING COMPUTATIONS 31
scope whose reticle can be rotated by the navigator to align with the moving
ground).
In a moving air mass
where 8 is the pitch angle, VTAS is the true airspeed and (3 is the sideslip angle.
On a flat Earth, the north and east (or grid north and grid east) distances traveled
are found by integrating the two components of velocity with respect to time.
On a curved Earth, the position coordinates are not linear distances but angular
coordinates. Equations 2.5 show a method for transforming linear velocities
to angular coordinates. The accuracy of airspeed data is limited by errors in
predicted windspeed and by errors in measuring airspeed and drift angle.
The dead-reckoning computer can process Doppler velocity. If the Doppler
radar measures ground speed V., and drift angle o,
Vnorth = Vg cos (1/;T + o)
Veast = Vg sin (lj;T + o) (2.12)
For example, in a spherical coordinate frame whose z3 -axis lies along the posi-
tion vector R:
dRI =dR
- - R+(wyzXR)
A
dt y dt
where the first term is the rate of change of radius, along the radius vector,
32 THE NAVIGATION EQUATIONS
and Wyz is the angular velocity of the Z;-coordinate frame relative to y;; R is
the unit vector in the direction of R. In direction-cosine mechanizations, the
C;J are related to the C;J by Equation 7.40, where w; - ll; of that equation is
identical to Wyz of this one.
2.5 POSITIONING
Y3
Radio station i
Surface of earth
Yl
Figure 2.5 Light-of-sight distance.
34 THE NAVIGATION EQUATIONS
v · (R,.;- Ro)
sin E = -----,------,-- (2.14)
IR.,;- Rol
If n is the unit north-pointing horizontal vector at the aircraft, the azimuth angle
of the line of sight is
These vector equations can be resolved into any coordinate frame. For exam-
ple, in ECEF,
(2.16)
In one-way ranging systems, where the clock offset and range are to be cal-
culated, the measured pseudorange vector from the ith station, R;m, is R 0 -
R,;, corrected for the unknown offset of the airborne clock and for propagation
delays in the atmosphere, expressed in distance units. The magnitude of the
pseudorange in any coordinate frame (e.g., ECEF) is
R;m = ryc(TOA)
Rim(k) = [(xk- x,.;) 2 + (yk- y,;) 2 + (Zk- z.,;) 2 ] 1/ 2 -ryctk (2.17)
where
is the measured pseudorange from the aircraft to the ith
station
Rim(k) is the calculated pseudorange in the kth iteration.
TOA; is the time of arrival of the signal from the ith station
relative to the expected time of arrival as measured by
the aircraft's clock
c is the speed of light in vacuum = 2.99792458x 108 m/sec
POSITIONING 35
The stations may be moving (e.g., satellites) or stationary (e.g., GPS pseudolites).
Four pseudorange measurements are needed to solve for the four unknowns in
Equations 2.17: aircraft position and clock offset. When more than four measure-
ments are made, the equations are overdetermined so that a solution requires a
model of the ranging errors, for example, using a Kalman filter (Chapter 3). The
airborne computer usually solves for its position by assuming a position and clock
offset, calculating the pseudoranges to four stations from Equation 2.17, compar-
ing to the measured pseudoranges (with respect to its own clock), and iterating
until the calculated and measured pseudoranges are close enough. The next iter-
ation is chosen as follows: In the kth iteration, the assumed position is xk Yk Zk
whose range to the ith station differs from the measured range R;m by LlXk =
R;m(k)- R,;. The components of LlXk in the navigation coordinates are LlXk, LlYb
and LlZk. The sensitivity of R;111 to position is
(2.18)
where dR;m/dX1 are the direction cosines between the line of sight to the ith
station and the jth coordinate axis and Llt k is the error in estimating clock offset.
If the assumed position and clock offset were correct, LlXk. LlYk, LlZk, and Lltk
would be zero and LlR;k would also be zero. But if the assumed position were
misestimated by LlXk, the error along the line of sight would be the dot product
of the unit vector along the line of sight with LlXk. Thus after computing Rimk
from Equation 2.17, the next iteration is Rim(k + l) = Rim(k)+LlXk. Iterations cease
when the difference between the calculated and measured pseudorange is within
the desired accuracy. A recursive filter allows a new calculation of position
and clock offset after each measurement (Chapters 3 and 5). In Equation 2.17,
Earth-based line-of-sight navaids use an average index of refraction YJ (Chapter
4), whereas satellite-based navaids, whose signals propagate mostly in vacuum,
assume that YJ = I and correct for the atmosphere with a model resident in the
receiver's software (Section 5.4).
of arrival of the navigation signal from two or more stations and compute its
own position as follows:
Three or more stations are needed if the aircraft's clock is not synchronized
to the Loran or Omega transmitters. If a receiver incorporated a sufficiently sta-
ble clock, only two stations would be needed for a direct-ranging fix. Two-sta-
tion fixes with a synchronized clock or three-station fixes with an asynchronous
clock result in two position solutions, at the intersections of two circular, each
having a vortex at a transmitter. The correct position can be found by receiving
an additional station or from a priori knowledge.
the intervening Earth's surface and the presence of the sunlight terminator
between the aircraft and the station.
• Subtract the times to two stations to calculate the predicted difference in
propagation time.
• Measure the difference in time of arrival of the signals from the two sta-
tions.
• Subtract the measured and predicted time differences to the two stations.
• Calculate the time-difference gradients from which is calculated the most
probable position of the aircraft after the measurements (see Section 2.5.2
of the First Edition of this book and Chapter 4 of this book).
• Iterate until the residual is smaller than the allowed error.
These navigation systems obtain occasional updates when the aircraft overflies
a patch of a few square miles, chosen for its unique profile [5]. A digital map
of altitude above sea level, h,, is stored for several parallel tracks; see Figure
2.6. For example, if 0.1-nmi accuracy is desired, h,(t) must be stored in 200-ft
squares sampled every 0.2 sec at 600 knots.
The aircraft measures the height of the terrain above sea level as the differ-
ence between barometric altitude (Chapter 8) and radar altitude (Chapter 10);
see Figure 2.7. Each pair of height measurements and the dead-reckoning posi-
tion are recorded and time-tagged.
After passing over the patch, the aircraft uses its measured velocity to calcu-
late the profile as a function of distance along track, hm(x), and calculates the
cross-correlation function between the measured and stored profiles:
Stored track + 1
-------~-
Nominal track - -~
------------~------~~~--
Stored track -1 ~
-- --
-------~- ~------~
hradar
Sea level
Figure 2.7 Measurement of terrain altitude.
ilA
r:f>ms(7) =
f
0
hm(x)h,(x- 7) dx (2.19)
where the map patch has a length A. The integration is long enough (n > l) to
ensure that the patch is sampled, even with the expected along-track error. The
computer selects the track whose cross-correlation is largest as the most prob-
able track. The computer selects the x-shift of maximum correlation 7 as the
along-track correction to the dead-reckoned position. Heading drift is usually
so small that correlations are not required in azimuth. The algorithm accom-
modates offsets in barometric altitude caused by an unknown sea-level setting.
The width of the patch depends on the growth rate of azimuth errors in the
dead-reckoning system. Simpler algorithms have been used ("mean absolute
differences") and more complex Kalman filters have been used [5].
Terrain correlators are built under the names TERCOM, SITAN, and TER-
PROM. They are usually used on unmanned aircraft (cruise missiles) and can
achieve errors less than 100 ft [ 1, 8]. The feasibility of this navigation aid
depends on the existence of unique terrain patches along the flight path and
on the availability of digital maps of terrain heights above sea level. The U.S.
Defense Mapping Agency produces TERCOM maps for landfalls and mid-
course updates in three-arcsec grids.
radio beacons. The computation begins with the best-estimate of the present
position of the aircraft and ends by delivering computed range and bearing to
other vehicle subsystems (Figure 2.1 ). Waypoints may be loaded before depar-
ture or inserted en route. The navigation computer, mission computer, or flight-
management computer performs the steering calculations.
Range and bearing to a destination can be calculated by using either the
spherical or the plane triangle of Figure 2.8. If flat-Earth approximations are
satisfactory, the xy coordinates of the aircraft are computed using the dead-
reckoning Equation 2.1 0; x 1 and y 1 of the targets are loaded from a cassette or
from a keyboard. Then, range D and bearing BT to the target, measured from
true north, are
1 X- Xt
BT =tan (2.20)
Y- Yt
The crew will want a display of relative bearing (BR = BT - 1/JT) or relative
track (TR = BT- T T). BT is the true bearing of the target. Relative bearing BR
is the horizontal angle from the longitudinal axis of the aircraft to the target,
and relative track T R is the horizontal angle from the ground track of the aircraft
to the target (Figure 2.4 ).
1
If .::1}.. and Llcll are less than radian, the plane triangle solution exceeds the
Ya
North pole
Waypoint
(2.21)
D
cos ~ = sin cp sin cp 1 + cos cp cos cp 1 cos (/\ - 1\ 1)
Rc
(2.22)
Rc is the Gaussian radius of curvature, Equation 2.4. At long range, where the
absolute and percent errors are largest, they are usually least significant. Within
I 00 nmi of the aircraft, the Earth can be assumed flat within an error of 0.3
nm1.
Steering and range-bearing computations can be performed directly from the
Z; navigation coordinates or from the direction cosines Ci to prevent singular-
ities near the north and south poles.
Knowing the measured or computed ground speed, an aircraft can be steered
in such a manner that the ground speed vector-not the longitudinal axis of
the aircraft-tracks toward the de~:ired waypoint (relative track is the steering
command that is nulled). The difference between heading toward the target
and tracking toward the target is significant only for helicopters; the vehicle
eventually arrives there in either case, by slightly different paths.
Two general kinds of steering to a destination are commonly used: (I) steer-
ing directly from the present position to the destination and (2) steering along
a preplanned airway or route. The former results in area navigation (Section
2.7.4) using the shortest (though not necessarily the fastest) route to the des-
tination, whereas the latter is representative of flying along assigned airways
(Chapter 14). Either steering method may be solved by the plane-triangle or
spherical-triangle calculation.
The rhumb line is used by ships and simple aircraft. It is defined by flying
at a constant true heading to the local meridians. The resulting flight path is a
straight line on a Mercator chart. A1,rcraft sometimes divide a complex route into
rhumb-line segments so that each segment can be flown at constant heading.
More often, the continually changing heading toward the next waypoint is
recomputed and fed to the autopilot. Since the great circle maps into a near-
straight line on a Lambert conformal chart, the crew can monitor the flight
path by manual plotting, if desired. In the twenty-first century, electronic map
displays will show the moving aircraft on charts.
COURSE COMPUTATION 41
(2.23)
(2.24)
(Figure 2.8). The lateral speed V 2 is the magnitude of the dot product of the air-
craft's velocity with this unit vector. The range to go from R 3 to R 2 is given by
Equation (2.20) or (2.22). Time to go is calculated from the proposed velocity
schedule.
The fastest route is neither the great circle nor the airway because of winds,
especially because of the stratospheric jet stream whose speed often exceeds I 00
knots. Thus where high-altitude aircraft are not confined to airways, they follow
preplanned "pressure-pattern" routes that take advantage of cyclonic tail winds.
(2.25)
IR3 x Rzl
(2.26)
IR311Rzl
(2.27)
p2
_... ---- --,
--- Pt
''
'' '
''
I
I ''
I
I
I
I
I
I
I
I
.'
VI and v3. The position PI is found from the triangle PI VI V3. The aircraft's
position must be known well enough to exclude the false solution at P 2 . An
artificial airway is defined by the points A 1 and A 2 • D and L are usually found
iterati vel y:
5. The next estimate of p 1 is along the vector I::J.pk in Figure 2.9, whose
components along p 1 and P3 are I::J.p 1 and I::J.p3.
6. Repeat step 2 and iterate until I::J.p; are acceptably small.
44 THE NAVIGATION EQUATIONS
After determining the aircraft's position P 1 , the distance-to-go and lateral dis-
placement are calculated as in Section 2.7. Lis sent to the autopilot, as explained
in Section 2.7.3, and D is used to calculate time-to-go.
In the 1990s, civil aircraft were being allowed the freedom to leave RNAV
airways (Section 1.5.3) using GPS, inertial, Omega, and Loran navigation, none
of which constrain aircraft to airways.
Deterministic errors are added algebraically and statistical errors are root
sum squared. The total is sometimes called total system error.
The two- or three-dimensional vector position error r can be defined as indi-
cated minus actual position. A series of measurements taken on one navigation
system or on any sample of navigation systems will yield a series of position
measurements that are all different but that cluster around the actual position.
If the properties of the navigation systems do not change appreciably with age,
if the factory is neither improving nor degrading its quality control, and if all
systems are used under the same conditions, then the statistics of the series
of measurements taken on any one system are the same as those taken for a
sample ("ensemble") of systems. Mathematically it is said that the statistics are
ergodic and stationary.
If the position errors are plotted in two dimensions, as shown in Figure 2.1 0,
NAVIGATION ERRORS 45
Arbitrary y axis
Principal x axis
• •
Arbitrary x axis
•
Actual position _or
Figure 2.10 Two-dimensional navigation errors. In principal axes, the x andy statistics
are independent.
it will generally be found that the average position error r = Et:.r;/N) is not zero
(t:.ri are the individual position errors; N is the total number of measurements).
Two measures of performance are the mean error and the circular error probabil-
ity (CEP) (also known as the circular probable error, CPE). The CEP is usually
considered to be the radius of a circle, centered at the actual position (but more
properly centered at the mean position of a group of measurements) that encloses
50% of the measurements. The mean error and the CEP may be suitable as crude
acceptance tests or specifications, but they yield little engineering information.
More rigorously, the horizontal position error should be considered as a
bivariate (two-dimensional) distribution. The mean error r and the directions
of the principal axes x andy, for which errors in x and y are uncorrelated, must
be found. To find the principal axes, a convenient orthogonal coordinate system
(x',y') is established with its origin at r. Then the standard deviation (or rms)
a in each axis and the correlation coefficient p between axes are calculated:
L,(x;)2
N
L,(y;)2
N
~ ''
L...XiYi
p=-- (2.28)
Gx'G:/
46 THE NAVIGATION EQUATIONS
From these quantities are determined a new set of coordinates, xy, which are
rotated () from x' y':
2ax'av'
tan 2() = 2
·
2
p (2.29)
ax'~ ay'
The origin of the new xy-coordinates coincides with the origin of x' y'. The com-
ponents of the position errors along x and y are uncorrelated and can be consid-
ered separately. In inertial systems, the principal axes are usually the instrument
axes. In Doppler systems, the principal axes are along the velocity vector and
normal to it or along the aircraft axis and normal to it if the antenna is body-
fixed (as was usually the case in 1996). In Loran systems, the principal axes
are found by diagonalizing the covariance matrix. The rms errors in the new
coordinates can be calculated anew or can be found from the errors in x'y' from
[7, p. 598].
The one-dimensional statistics along either of the principal axes will now be
discussed. First, the mean and standard deviation are found in each independent
axis. The errors along each axis are plotted separately as cumulative distribution
curves, which show the fraction of errors less than x versus x. This curve allows
all properties of the statistics to be determined. In many systems, experimental
cumulative distributions will fit a Gaussian curve.
If the one-dimensional errors are indeed Gaussian, their statistics have the
following properties:
Mean x
Standard deviation (rms) a
50% of the errors lie within x ± 0.675a (probable error)
68.3% of the errors lie within x±a
95.4% of the errors lie within x±2a
99.7% of the errors lie within x±3a
Navigation systems, civil and military, are often specified by the fraction of
navigation errors that fall within a circle of radius E p, centered on the mean.
Figure 2.11 shows this probability if ax and a y are Gaussian distributed and
uncorrelated. For example, if ax = 0.4 nmi and a y = 1.0 nmi, b = ax/ ay = 0.4.
If the radius of the circle is E p = 1.5 nmi, then a = E p /a y = 1.5, and from the
graph, P = 0.85. Thus, 85% of the errors fall within a circle of radius 1.5 nmi.
Weapons often inflict damage in a circular pattern. Hence, military tactical
navigation systems are sometimes specified by the CEP, the radius of the circle
that encompasses 50% of the navigation errors (which are inherently elliptically
distributed):
. av
CEP = 0.59(ax + ay) ± 3% 1f - · < ax < 3a y (2.30)
3
When ax= ay =a, the CEP = 1.18a.1 and, from Figure 2.11, 95% of the
navigation errors lie within a circle of radius 2.45a.
o..0.8
"'
"''2
'0
~0.7~--~--~~~~--~----4----4----~
0
Q)
u
~0.6~--~--~~~4----4----4----4----~
c:
.r.
·~ 0.5 1-----+-+->f----t+-1--4------t----+-----+----+-l
Ep
~
a=-
Uy
0
~ 04 f-----++--+-11----4-------'u. - - - x<Uy--H
b= u;
~03r--~~-+~----+------+----+-----+----+-l
..6
"'
..0
0
0: 0. 2 f--1--ft-J'f-----1-----+------1----+------t-----t-1
2 3
a
Figure 2.11 Probability of an error lying within a circle of radius tp. (ax and av are
the uncorrelated standard deviations in x andy.) From unpublished paper by Bacon and
Sondberg.
48 THE NAVIGATION EQUATIONS
(2.31)
(2.32)
(PDOP) 2 = (2.33)
(HDOP) 2 = (2.34)
where TDOP is the time dilution of precision, the contribution of clock error
to the error in pseudorange. Equations for GDOP, PDOP, and HDOP, when the
standard deviations in range to each station are different, are provided in [ 12].
Dilution of precision plays an important role in radio-ranging computations,
especially for Loran (Chapter 4) and GPS (Chapter 5). Detailed GDOP equa-
tions for Loran are in Section 4.5.1, for Omega in Section 4.5.2, and for GPS
in Section 5.5.2. Receivers usually flag a PDOP or HDOP greater than approx-
imately 6 as an indication of poor geometry of the radio stations, hence a poor
fix.
Traditional aeronautical charts are printed on paper. They are of three kinds:
Military targeting charts show the expected location of defenses, the initial
50 THE NAVIGATION EQUATIONS
approach fix, the direction of approach, a visual picture of the target in season
(e.g., snow covered) and the preplanned escape route.
Since World War II, experiments have been made with analog charts driven
by automatic navigation equipment. Paper charts were unrolled or scanned onto
CRTs, while an aircraft "bug" was driven by the navigation computations. The
systems were limited by cost, reliability, and the need for wide swaths of chart
to allow for diversion. Their use was confined to some helicopters and experi-
mental military aircraft.
In 1996, digital maps were well··established in surveying data bases, the cen-
sus, automotive navigation, and other specialized uses. Manufacturers of nav-
igation sets created their own data base of navaids and airports or purchased
one. Small digital data bases were included in the navaid's ROM whereas large
data bases, especially those that included terrain, were usually delivered to cus-
tomers on CD-ROM. National cartographic services in the developed world
were all converting from paper maps to digital data bases. The U.S. Defense
Mapping Agency (DMA) issued a standard for topographic maps on CD-ROM
[19], and several other nations' cartographic agencies were doing the same. The
U.S. Defense Mapping Agency produces separate data bases for terrain eleva-
tion and cultural features. They can be stored separately and superimposed on
an airborne display. A U.S. National Imagery Transmission Format was created
to send and store digital data. The GRASS language was widely used in the
United States to manipulate DMA data [6]. Private companies were produc-
ing remarkably diverse Geograph:tc Information System (GIS) data bases. In
1996, at least one company (Jeppesen) was producing digital approach plates
on CD-ROM (4]. RTCA published a guide to aeronautical data bases [24a] as
did ARINC [26a].
The technical challenges have been (I) to standardize the medium (e.g., CD-
ROM), (2) to standardize the format of data stored on the medium so that any
disc could be loaded onto any aircraft, just as any chart can be carried on any
aircraft, and (3) to develop on-board software that displays sections of the chart
across which the aircraft seems to move (moving map or moving bug, or both)
and orient the chart properly (north up or velocity vector up). As the aircraft
nears the edge of the chart, the software must move to a new section while
avoiding hysteresis when flying near the edge of a chart. The expectation is
that CD-ROM en-route and approach charts will be readily available to military
users before the year 2000.
Digital chart displays have provisions for weather or terrain overlays (from
airborne radar or from uplinked data), and provisions for traffic overlays (from
on-board TCAS, ground uplink on Mode-S, or position broadcasts from other
aircraft, Chapter 14). Civil airlines, driven by cost considerations, may gradually
abandon their practice of purchasing charts and distributing them to the crews
in hard copy. Instead, they may at first print charts on demand in the dispatch
room from a central data base and, later, distribute portable digital charts on
CD-ROM or via radio uplink to be loaded into the aircraft avionics when the
crew boards.
SOFTWARE DEVELOPMENT 51
loaded into the flight computer via the on-board data bus. Revisions of
flight software may be issued from time to time.
7. A copy of the flight software is usually delivered to a training facility,
where it is used to check out crews in a ground simulator. The simulator
may be a part-task computer-based trainer (CBT); a terminal that emu-
lates the navigation keyboards, on-board computer, and displays; or it may
be a high-fidelity emulation of the cockpit and avionics. A high-fidelity
simulator may incorporate a flight computer that contains the navigation
software or may rely on a scientific computer, programmed to emulate the
flight computer. In a CBT, sensor inputs are simulated; in a high-fidelity
simulator, they may come from real or simulated hardware. Simulator
training is cheaper and often more effective than flight training.
8. The final task in the preparation of the navigation software is the evalua-
tion of its performance during flight versus the specification. This is done
by the aircraft manufacturer or operator.
The increasing capability of airborne digital computers will permit more com-
plex algorithms to be solved. Companies that specialize in aircraft navigation
will continue to build libraries of proprietary algorithms that they incorporate
into their products. Crew interfaces will become more graphical to reduce work-
load and reduce errors in loading data. Direct loads from the ground via Mode-S
links and other data links (some via satellite) will be commonplace.
By the year 2000, on-board CD-ROM readers will display charts and flight-
manual data on military aircraft. The civil aviation industry may prefer to print
up-to-date paper charts for each flight in the dispatch rooms, downloaded from
a central data base, as an alternative to procuring and distributing them.
The software verification costs assigned to each aircraft will be substan-
tial, because they are amortized over a few hundred units, even with standard
libraries of routines.
PROBLEMS
2.1. The direction cosine matrix [C] transforms the Earth-centered inertial
coordinates y; into the locally-level navigation coordinates Zj· Let a = 0
when the aircraft is on the equator, and let the initial matrix be
[C] = [ ~ 0
0 !l
Let a= A..
PROBLEMS 53
(a) If the aircraft flies 90° due east, show the direction cosine matrix.
0I 00 0I]
Ans.
[0 I 0
(b) If the aircraft flies 90° due north from its original position, show the
direction cosine matrix.
Ans:
(c) If the aircraft flies 30° due east on the equator at 600 knots from its
original position, what are the C;i and the C;/ Use Equation 7.40.
Ans:
. 5
[C)=-- --J3
I -J3
I
--f3]1 hr
-]
57.3 [ I _-J3 O
2.2. An aircraft flies 3 hrs east then 2 hrs north at 300 knots at an altitude of
3 nmi, starting at 40° north latitude.
(a) Find its position using the flat-Earth dead-reckoning equations.
Ans. x = 900 nmi, y = 600 nmi.
(b) Find the final latitude-longitude using spherical-Earth equations with
the Gaussian radius of curvature.
Ans. <P = 49.98°, .1A = 19.54°.
(c) Find the final latitude-longitude using the ellipsoidal-Earth equations,
(2.5).
Ans. <P = 49.99°, .11-. = 19.50°.
(d) Find the distance from the start to the destination of case a using the
flat-Earth range equation, 2.20.
Ans. 1081.7 nmi.
(e) Find the distance to the destination of case b using the spherical-Earth
approximation with the Gaussian radius of curvature.
Ans. 1022.7 nmi.
(f) Find the distance to the destination of case b using the flat-Earth equa-
tions plus the correction of Equation 2.21.
Ans. 1027.0 nmi.
54 THE NAVIGATION EQUATIONS
2.4. Derive Equation 2.21. Hint: solve the spherical triangle for the cosine of
the range angle and express the coordinates of the waypoint as the present
position plus a small increment. Expand the sines to third order and the
cosines to second order.
2.5. Verify the calculations on page 46 for the probability of falling inside a
rectangle of edge 1.2 nmi (x = 0.6 nmi) by 2.4 nmi (y = 1.2 nmi). Let ax
= 0.4 nmi and a y = 2.4 nmi.
3 Multisensor Navigation
Systems
3.1 INTRODUCTION
• Linear and angular acceleration and rate in body coordinates for vehicle
control purposes.
• Vehicle state relative to the air mass including orientation-angle of attack,
sideslip, and airspeed, again for vehicle control purposes.
Two fundamental error sources affect the error behavior of an inertial system.
These are the errors in the measurements of force made by the accelerometers
and the errors in the measurement of angular change in orientation with respect
to inertial space made by the gyroscopes. The basic mechanization of an iner-
tial system, which is described in detail in Chapter 7, is depicted schematically
in Figure 3.1. In this figure it is seen that the force measurements made by the
accelerometers are first transformed to a selected navigation coordinate frame
that is typically the local geodetic coordinates of east, north, and local verti-
cal. These measurements are then compensated for the force of gravity with
a mathematical model, such that vehicle acceleration with respect to inertial
space is obtained. The resulting variable, after correction for Coriolis accelera-
tion, is then integrated once into velocity and a second time into position change
with respect to the Earth. Additionally, the gyroscopic measurements of angular
change with respect to inertial space are modified using the system computed
velocity and the Earth's rotation rate vector to reflect the rotation of the local
vertical due to earth rotation and the vehicle change in position as it travels over
the surface of the Earth. In this manner the orientation of the accelerometer
axes relative to the Earth at the present position of the vehicle is continuously
computed. The result is that there are three sources of change in orientation
error of the accelerometers with respect to an Earth-fixed reference coordinate
frame: (I) integrated gyro drift rate; (2) integrated error in system computed
velocity which results from error in the measurement of acceleration-due to
accelerometer measurement errors, imperfect knowledge of the local force of
gravity and the current error in the know ledge of orientation of the accelerome-
ter sensing axes which causes a misresolution of any accelerometer force mea-
surement including that of the gravity vector; and (3) error in the orientation
of the navigation coordinate axes which changes as they rotate with respect to
inertial space.
The result of this interaction of error effects is that the error characteristic of
an inertial system for the computation of velocity, position, and instrument axes
orientation is described by the sinusoidal Schuler oscillation that has a period of
approximately 84.4 minutes. Due to this oscillatory characteristic, the position
error response to a step of constant accelerometer measurement error is not a
quadratic in time but a bounded Schuler oscillation as shown in Figure 3.2.
Note that the error in position and the "tilt," which is the error in orientation of
the accelerometer sensing axes with respect to the local level plane, are equal
for any acceleration measurement error.
The velocity error response to a step of constant gyro drift rate error, as
shown in Figure 3.3, is a bounded Schuler oscillation that has a constant error
dictated by the magnitude of the gyro drift rate. The error in position is also
characterized by a Schuler oscillation but diverges in time in proportion to the
integrated velocity error. The tilt due to the gyro drift rate is a bounded Schuler
oscillation with zero mean. This occurs because any tilt results in a counter-
U1
00 CORIOLIS ACCELERATION
COMPUTATION
ACCELEROMETER ~-----------,
FORCE TRANSFORMATION
I
(}
v
I
MEASUREMENTS TO NAVIGATION 1 (}
- '- COORDINATE R
FRAME I •-
MATHEMATICAL
MODEL FOR
MATHEMATICAL
MODEL FOR
I
EARTH'S EARTH'S
GRAVITY VECTOR ROTATION VECTOR
GYRO
MEASUREMENTS
l
I "·' > I
I
I Velocity
I
I Error
" I
I ov I
R. I f-......1'-------. </> Platform
Tilt
tiA
Acceleration M
Measurement Error Gyro Drift
Rate
I 00 Position
Error
t: j t----+-'----~
0.5
~j 0 f------lr-----r------1
V - Step of Acceleration
Measurement Error
1
-0.5 0.5
-1~--~----~~~~--~
0 2 4 8 2 4 6 8
6
§.] §.-~
Ul Figure 3.2 Inertial system error response to a step of acceleration measurement error.
\C
Constant acceleration measurement error induces a zero-mean Schuler oscillation in velocity
error and an identical nonzero mean oscillation in tilt and position error.
""
0 r-----------------Dg ----------------,
I
I
I
I
I g Gravity Magnitude I
I R Radius of Curvature of Earth I
I w, ~ jf, (}Schuler Frequency )I
I Velocity 1
I I I I E~or I * ~ I I ! o¢Pl~~~rm
M
Acceleration
Measurement Error
I <ld
Gyro Drift
Rate
I
b - Step of Gyro Drift Rate
I
I I • o8 Position
Error
[£:]
[R b~''t \ /
t/>(1) = sin (w)) ~
t~3 °
5113feet per
-1-: .... ~
1 ~
-0.5 / ..~ -6
oO(I) = [£:]
@ 84.4 minutes
.~
[sio (w,l) - (w))]
-1 0 2 4 6 8 2
' 4 6 8 -8
0 2 4 6 8
COMPUTED DIRECTION TO
STAR FROM UNTILTED TRUE DIRECTION TO STAR (AT oo)
PLATFORM AT ERRONEOUS
POSITION -------------
POINTING DIRECTION TO
STAR AT TRUE POSITION
FROM TILTED PLATFORM
8
=
IF <I> 08, THEN NO ERROR IN 1jJ IS
OBSERVED, ALTHOUGH THERE EXISTS
• ERROR IN COMPUTED BOTH ORIENTATION ERROR <j>, AND
POSITION OF 68 ONLY POSITION ERROR 08
n
(NO TIL RESULTS IN
A POINTING ERROR
TO STAR 0Fc56
Figure 3.4 Angular position and orientation errors observed with a star-tracker with
inertial stabilization. A stellar sensor observes the error 1/;, the difference between the
position error (M) and the orientation error (1>) of the inertial system.
2. Error in system computed angular position on the surface of the Earth (j().
This causes an error in pointing to the star even if the reference coordinate
frame provided by the inertial system is perfectly aligned with respect to
the Earth or inertial space at its actual position. This situation is indicated
at the computed position in Figure 3.4.
3. Error in system time that results in an orientation error of an Earth-refer-
enced frame with respect to the stellar background due to Earth rotation.
This error is usually minimal and is not considered further here.
The implications of this situation is obtained by referring to the error block dia-
gram of Figure 3.5 where the error 1/; is shown. Since with an inertial system
any error in the accelerometer measurement results in an equal error in both
position (j() and orientation ¢, then for such an error source 1/; will be zero and
will not be detected with a star observation. On the other hand, a constant gyro
drift rate produces an oscillatory diverging error in the angular position (j() and
a bounded oscillatory error in orientation ¢, so the resulting error 1/; is a ramp
AN INTEGRATED STELLAR-INERTIAL SYSTEM 63
r--------1 -g
GYRO DRIFT RATE
t---+--rp TILT
ACCELERATION
MEASUREMENT ERROR
1j1 = [<j>--b9] = f bd dt
Figure 3.5 Illustration of the 1/.; vector observed with a star-tracker. An inertially sta-
bilized stellar sensor is used to measure the effects of gyro drift rate but cannot observe
the effect of acceleration measurement errors.
proportional to the gyro drift rate. Since this diverging error is detectable with
the inertially stabilized star-tracker, then clearly the error in position can be
corrected or "reset" when a star shot is taken.
Most stellar-inertial systems have been mechanized with inertial systems
using stabilized gimbal control loops (as described in Chapter 7). This mech-
anization permits a telescope and inertial instrument assembly to be manufac-
tured that has only two low-bandwidth pointing loops that maintain the tele-
scope line of sight to the star being observed, usually for the elevation and
bearing angles to the selected star. High-bandwidth pitch, roll, and heading sta-
bilization control loops employing gyro measurements of small angular dis-
turbances isolate the inertial instruments and tracker assembly from vehicle
angular motion. When such a gimballed mechanization is employed, the sys-
tem gyros are rotated slowly with respect to inertial space to maintain a desired
orientation of the instrument assembly with respect to the Earth, usually that
of the local-level mechanization. Because the gyro sensing axes change slowly
with respect to inertial space, the detected components of the 1/; vector, which
is the integral of gyro drifts projected onto inertially fixed axes, can be used to
calibrate the gyros. This is possible because the detected components of the 1/;
vector can be correlated with specific gyro-sensing axes.
Since the component of the 1/; error vector along the line of sight to an indi-
vidual star cannot be observed, two star shots are generally required to com-
pletely correct the system. Ideally these lines of sight are orthogonal to enhance
observability. It is possible to obtain a complete calibration of the system gyro
drift rate vector, when the following two conditions occur:
2. The gyro sensing axes have not changed orientation with respect to iner-
tial space by a great amount since the last point in time when a pair of
stars was observed.
1
As a rule of thumb, the uncompensated effect of the deflection of the vertical results in an inertial
system position error divergence of about 0.1 nmi in one hour, thereafter being in proportion to
the square root of time.
INTEGRATED DOPPLER-INERTIAL SYSTEMS 65
ACCELERATION
ljJ
BIAS STEP
0V R e-----+j
Figure 3.6 Integration of a reference speed sensor with an inertial system. A reference speed
sensor is used to observe the Schuler oscillations induced by acceleration measurement error
and gyro drift rate and the velocity bias error induced by constant gyro drift rate.
AN AIRSPEED-DAMPED INERTIAL SYSTEM 67
Figure 3.7 Block diagram of a fixed-gain (KD, Kp) speed sensor damped-inertial sys-
tem mechanization.
When a Doppler radar (or air data sensor) is combined with a stellar-inertial
system, significant benefits are obtained. Stellar measurements are used to cali-
brate the inertial system gyro drift rates, thereby reducing any mean error in the
inertial system computed velocity components. Consequently, when the stellar-
inertial system velocity components are compared with commensurate velocity
components from a reference velocity sensor, the mean errors in the reference
velocity sensor components can be calibrated. Effects of the errors in the mea-
surement of acceleration (accelerometer measurement error and uncertainty in
the gravity vector) can now be more effectively removed using the calibrated
reference velocity sensor. Recall that the errors induced by an error in the vehi-
cle acceleration measurement cause divergent velocity, position and orientation
errors in a stellar-inertial system, not observable with the star tracker.
The most attractive attribute of a stellar-inertial system when combined with
a reference velocity sensor is that the resulting errors in all the navigation vari-
ables are bounded. The residual error levels in these variables are of course a
function of the accuracy of the individual sensors employed. Accuracies on the
order of a hundred feet in position, a fraction of a foot per second in velocity
and subarcsecond in orientation are achievable with such a system. An added
attraction of such a system is that it is fully autonomous.
NONINERTIAL GPS MULTISENSOR NAVIGATION SYSTEMS 69
Although most integrated navigation systems are inertially based, there are
exceptions. In particular, GPS receiver module cards have been developed that
can be embedded in the chassis of other sensors such as a Doppler radar [ 12] or
Loran [ll, 14, 15]. In the case of Loran, the integration was motivated by the
fact that the GPS satellites may not be sufficient in number at all times to pro-
vide positioning data with guaranteed integrity worldwide. A further attraction
of Loran is that it does not suffer from the line-of-sight limitation of GPS. It
employs a low-frequency signal that propagates along the surface of the Earth
(the ground wave), and it is time-synchronized using atomic clock controlled
70 MULTISENSOR NAVIGATION SYSTEMS
transmitters. In the case of Doppler radar, the motivation for GPS integration is
to significantly reduce the effect of Doppler associated error sources on navi-
gation performance. These include speed measurement errors over water, scale
factor error and bias, and the errors in a low-accuracy inertial attitude and head-
ing reference system (AHRS). Of course, the Doppler-AHRS system provides
the ability to dead-reckon position during GPS outages.
(3.1)
where
Xa is the actual position of the stationary vehicle
Ed; is a deterministic error in the ith measurement
~; is a random error in the ith measurement
(3.2)
The best estimate of position in this case is the average of the measurements
(3.3)
(3.4)
is the standard deviation of the random error a~ divided by the square root of
FILTERING OF MEASUREMENTS 71
(3.5)
(3.6)
where only the (N - I )th estimate and the Nth measurement are required.
(3.7)
(3.8)
The minimum error in this least squares estimate is the weighted sum of the
independent errors ox In and 8x2N:
(3.9)
72 MULTlSENSOR NAVIGATION SYSTEMS
(3.10)
(3 .11)
where it is evident that the variance of the error in the combined estimate oxN
is less than the errors obtained from either of the two independent position
sensors.
The Kalman filter requires that all error states are modelable as zero mean noise
processes with known variances, power spectral densities, and time correlation
parameters. Thus, the various error quantities to be estimated and the associ-
KALMAN FILTER BASICS 73
ated measurement noises are all random processes whose correlation structure is
assumed to be known. The Kalman filter then obtains estimates of the states of
these stochastic processes, which are described by a linear or linearized mathe-
matical model. It accomplishes this goal by capitalizing on the known correla-
tion structure of the various processes involved and the measurements of linear
combinations of the error states. To do this, both the error propagation in time
and the measurement processes are expressed in vector form. This provides
a convenient way with linear matrix algebra to keep track of relatively com-
plex relationships among all the quantities of interest. This is one of the main
features that distinguishes a Kalman filter from most digital signal processing
applications. The navigation system integration problem is multiple-input and
multiple-output in nature, and matrix algebra is essential in keeping track of
the relationships among the variables.
In the multisensor navigation system application, the behavior of the sen-
sor error states are described by linear differential or equivalent finite dif-
ference equations. Comparison of measurements between navigation sensors
are described by linear combinations of these error states. The filter accom-
plishes the task of error state estimation by a more complex recursive procedure,
rather than forming a simple weighted sum of the individual measurements as
described above.
Under the assumption of Gaussian noise distributions, the Kalman filter min-
imizes the mean square error in its estimates of the modeled state variables
[5]. The Gaussian assumption is usually a reasonable assumption in navigation
applications because the noise effects that arise in the measurements are often
due to a summation of many smaller random contributions. Thus, by the cen-
tral limit theorem of statistics, there is then a tendency toward the Gaussian
distribution regardless of the distribution of the individual contributions. Many
multisensor navigation systems have successfully implemented Kalman filters
while implicitly making this assumption.
(3.12)
where
xk is the n x I process state vector at time tb being errors in position,
velocity, and attitude and sensor errors
74 MULTISENSOR NAVIGATION SYSTEMS
(3 .13)
(3 .14)
where
Zk is the m x 1 vector measurement at time tk
Ilk is am X 1 measurement white noise vector at time tk
Hk is the m X n matrix giving the ideal relationship between Zk and xk
when no noise is present
(3.15)
Note that there are now four known matrices that describe the process and mea-
surement models: ¢b Qb Hb and Rk. Thus, the design of a Kalman filter
requires considerable a priori knowledge about the dynamics and statistics of
the various processes to be modeled. The state transition matrix r:f>k describes
how the state vector Xk would propagate from one step to the next in the absence
of a driving function. The Qk parameter tells us something about the noise in
the Xk process. If the elements of ~b are large, this means that a large amount of
randomness is inserted into the process with each step. The Hk matrix describes
the linear relationship between sensor measurements Zk and the error states to
be estimated xk. Note that relatively complicated relationships can be accom-
modated as long as they are known and linear. Finally, Rk describes the mean-
square measurement noise errors. Generally, large values in the Rk matrices
means poor measurements. Note that all four of the key matrices are permitted
to vary with time (i.e., with the index k).
Perhaps the most challenging aspect of applying Kalman filtering in multi-
sensor navigation applications is that of establishing the mathematical equations
that describe the physical situation at hand and casting the equations in the form
dictated by Equations 3.12 and 3 . 14. Fortunately, this can be done with a rea-
KALMAN FILTER BASICS 75
(3.16)
The diagonal elements of the error covariance matrix are the variances of the
error in the filter estimate of the navigation system error state vector. The off-
diagonal elements of the P matrix are co variances between different error states
in the vector, and they contain important information as to the degree of corre-
lation of one error state with another. The use of such correlation information
in the gain computation is what distinguishes the Kalman filter from simpler
mechanizations. In Kalman filtering, a measurement is used not only to update
estimates of navigation error variables directly involved in the observation (e.g.,
position) but also to update estimates of error variables not directly involved
(velocity, sensor errors, etc.).
Pk = [1-Kk Hk ]Pk
w General Form for an Arbitrary
Gain
T
Uncorrected
Inertial Position,
Velocity &
Platform Orientation
Corrected Inertial System Outputs
Inertial
Navigation
System
Knertial System
.Error Estimate Prediction of
Propagation Aiding Sensor
[ Measurements
Kalman
Filter Observable
[ Design Measurement
Error
Uncorrected
Aiding Sensor
Q iding Sensor
·ror Estimate
l'ropagation
Corrected Aiding
Sensor
Aiding Sensors Measurements Measurements
such as GPS,
Doppler or Loran
Kalman filter theory. The details as to exactly how this linearization is gener-
ally done are given in references [2, 6]. The basic idea is to choose a reference
trajectory in state space that is close to the true trajectory and then write equa-
tions for the perturbation variables that represent the difference between the true
and reference trajectory for use in the filter equations. The equations in terms
of the perturbation variables are linear when higher-order terms are neglected.
In the block diagram of Figure 3.9, the INS-corrected output is used for the
reference trajectory as a matter of convenience, since it provides all the basic
navigation quantities of vehicle position, velocity, acceleration, and orientation
in a continuous manner. None of the other navigation sensors, when consid-
ered individually, provide the same complement of information in a continuous
manner. Thus, the INS is the logical choice for the reference, even though its
accuracy may be less than some of the other sensors at discrete instants in times
when their measurements are available.
The final reason for choosing the INS-based integration methodology has
to do with maintaining high dynamic response in the position, velocity, and
attitude state variables available from the inertial system without filtering. The
usual price associated with filtering is time delay or sluggish response. This
CLOSED-LOOP KALMAN FILTER MECHANIZATION 79
In the mechanization of the last section, the INS is not corrected internally
throughout the time span of the mission. Clearly, if the internally computed
inertial system navigation variables diverge too far from their true values, the
linearization assumption becomes suspect, and the associated modeling inaccu-
80 MULTISENSOR NAVIGATION SYSTEMS
racy can lead to difficulties. In the early days of Kalman filtering, system engi-
neers discovered that such divergence could be avoided by feeding the filter
error estimates back to correct internally the inertial system at each time point
where an aiding sensor measurement was available. This will then reduce the
difference between the corrected real-time computed navigation variables and
the true values, provided that the Kalman filter is producing good error esti-
mates. When the filter error estimates are fed back in this manner, the mecha-
nization is called the closed-loop Kalman filter mechanization [2]. This method
has been used extensively and successfully in a variety of actual navigation
applications. It is especially important to use the closed-loop Kalman mecha-
nization in applications where the mission length is relatively long and the error
model on which the filter equations are based are a simplification of the actual
linear error model of the system. Such simplifications are made to reduce the
computational burden on the real-time computer.
Figure 3.10 illustrates schematically the mechanization of the integrated INS
system for the closed-loop mechanization. This diagram and that of Figure 3.9
are both conceptual in that the INS correction that takes place is usually just a
matter of correcting certain numerical values in the inertial system computer. If
the error modeling assumptions are valid, the performance of the open-loop and
closed-loop Kalman filter mechanizations become essentially equivalent. Both
mechanizations of the Kalman filter have been used extensively, so the deci-
sion as to which should be used wi II depend on the particular situation at hand.
The closed-loop mechanization is generally used except when it is desirable to
preserve the pure INS output undisturbed by internal corrections. In this case,
-..
Inertial Corrected Inertial System Outputs
Navigation
System
T
I nertial
system
j~
,, Prediction of
Aiding Sensor
Measurements
Cor rections
l-
Kalman
Filter
--
Aiding Design
Se nsor
Cor rections , 1 Observable
Measurement j~
Error
Aiding
Sensors
Carrected Aiding Sensor Measurements
-....
Figure 3.10 Closed-loop Kalman filter architecture.
GPS-INS MECHANIZATION 81
One of the most important multisensor navigation systems is that which inte-
grates an inertial navigation system (INS) with a Global Positioning System
(GPS) receiver. GPS (Section 5.5) is a very accurate, worldwide satellite nav-
igation system using ranging that will be operating into the twenty-first cen-
tury. The GPS-INS integration is nearly always done using a Kalman filter. A
simplified Kalman filter error model for accomplishing such an integration is
presented in Section 13.3.
y y
DMEStation
(Xt.Yt)
p/
0
Aircraft Z
(x,y)
X
0
/~_:
(x, y)•-LSX-1 Ay
(x,* y*)
(3.17)
From Figure 3.11 b, p is the measured range within some measurement error
and p * is the predicted range based on the aircraft's assumed reference position
which is slightly incorrect. Note that there is a linear relationship between Llp
and the perturbations Llx and Lly. Thus, using Llp as the observation (difference
of the commensurate measurements of range from the DME and estimated range
from the INS) processed by the Kalman filter and choosing the incremental
quantities Llx and Lly as error state variables rather than the total position states
x andy, the required linear measurement relationship for Kalman filter theory
is satisfied.
Referring again to Figures 3.9 and 3.1 0, when the corrected INS position
computation is used to estimate the range to the DME station, the differencing
operation produces Llp using the corrected range measurement from the aid-
ing sensor with the addition of some measurement noise. This example accom-
plishes exactly what small-perturbation theory says is necessary to linearize the
measurement process. Note the coefficients of Llx and Lly are just the direction
cosines between the line of sight t:o the DME station and the respective x- and
y-axis. The GPS measurement situation is similar in this regard and the direc-
tion cosines also appear in the modeling details discussed below.
ub White Noise
Ud ··--.t·lr--J----,1
White Noise ...__ ___,
d •6\-·--~·1 J
L------'
b
•
Figure 3.12 GPS clock error model.
is called clock bias. This bias in seconds can be scaled by the speed of light
to obtain a range bias in meters. Similarly, the rate of change of clock bias is
usually given in meters per second.
It is difficult to model crystal clock errors with a first-order state model; see
reference [8] or Chapter I 0 of reference [6]. The 2 state model shown in Figure
3.12 used in most Kalman filter implementations in GPS receivers allows for
the estimation of both clock bias and drift. Numerical values for the spectral
densities of the white noise forcing functions ud and ub depend on the quality
of the crystal clock. The references show how to convert from the usual clock
specifications to the Kalman filter parameters. This 2-state model is embedded
in the larger ll-state error model for the whole GPS/INS system. (For additional
discussion of clock characteristics, see Section 5.3.2.)
""'"
Aximuth
Gyro White
Noise
Ez
-~----- I I·
az +g --
ay,z are the Vehicle Accelerations Along the y and z Axes
Wx is the Platform Angular Rate About the x Axis
R is the Radius of the Earth
Figure 3.14 Block diagram of a simplified inertial system y-axis error model for the 11-statc
GPS-INS Kalman filter.
oe
Ul
86 MULTISENSOR NAVIGATION SYSTEMS
ACCELERATION
WHITE NOISE
bAz t--+--•Az
Vertical Position
Error
Vertical Velocity
Error
(3.19)
A ••
'-lZ = o''A ,- + ( -2g ) A~ (3.21)
R '-1.(
rf>x = Ez (3.22)
GPS-lNS MECHANIZATION 87
The accelerometer measurement errors oA,, oAy, and oA, and the gyro drift
rates Ex, E_p and Ez are mutually independent white noise processes with known
spectral densities. Finally, the differential equations for the clock errors are from
the block diagram of Figure 3.12.
For GPS clock errors,
b = d + U!J
d= u" (3.23)
Once the state variables are defined, it is a routine matter to put the differential
equations into state-space form:
x = Fx + t (3.24)
where F is the system state vector dynamics matrix describing the dynamic
coupling between the (error) states. In expanded form Equation 3.24 is
88 MULTISENSOR NAVIGATION SYSTEMS
0 0 0 0 0 0 0 0 0 0
0 0 -(g + Az) 0 0 0 0 0 Ay 0 0
X] 0 0 0 0 0 0 0 Wx 0 0 X!
R
X2 X2
X3 0 0 0 0 0 0 0 0 0 0 X3
X4 0 0 0 0 0 (g + Az) 0 0 -Ax 0 0 X4
xs xs
I
X6 0 0 0 0 --- 0 0 0 -Wy 0 0 X6
X7 R X7
Xg 0 0 0 0 0 0 0 0 0 0 Xg
Xg Xg
XJO
XJj
0 0 0 0 0 0 2( !) 0 0 0 0 XJ()
Xjj
0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0
F
0
oAx
Ey
0
oAy
+ Ex (3.25)
0
OAz
Ez
Uh
UJ
The final step in determining the process model is to specify the ¢k and Qk
matrices. lf the update interval /11 is relatively small. ¢k can be approximated
with just the first two terms of the Taylor series expansion of eFt:.r:
This approximation must be used with care. For example, if the aircraft expe-
riences high dynamics, then !:lt will have to be quite small or else additional
terms must be included [6].
Computation of the Qk matrix in a high-dimensional system is usually an
onerous task, especially if there is nontrivial coupling among the various state
GPS-INS MECHANIZATION 89
variables as in this example. However, once the F matrix, the spectral densities
of the forcing functions, and the Llt propagation interval are specified, then Qk
is numerical computable [6, 7 J.
Four pseudorange measurements (one per satellite) are sufficient for a stand-
alone GPS solution. Each pseudorange measurement is of the general form
z = Hx + v (3.27)
X4
0 0 0 0 0 0
n
hi7
{i]
hi4
[ h, xs
h2I 0 0 h24 0 0 h27 0 0
z= 0 0 0 0 0 0
X6
h3I h34 h37
X7
h4I 0 0 h44 0 0 h47 0 0
Xg
H X9
xw
XII
(3.28)
90 MULTISENSOR NAVIGATION SYSTEMS
where the hiJ are the respective direction cosines between lines of sight to the
various satellites and the navigation coordinate axes. The R matrix reflects the
additive white measurement noise components in which in this model would
be a 4 x 4 matrix:
(3.29)
R would usually be specified as diagonal with all terms along the diagonal being
equal. The numerical values of the terms are chosen to match the expected
variances of the white pseudorange measurement errors. As a practical matter,
these terms are usually specified to be larger than the expected error variances
to compensate for the inaccuracy in modeling the measurement errors as white
noise. The GPS measurement frequency (usually about I Hz) is very high rel-
ative to the characteristic Schuler time constant of the inertial system errors,
so implementing an R matrix with large values does not significantly degrade
performance due to the filtering provided by the INS.
Inertial Error States The inertial system error states that are generally the
next in importance are the "biases" associated with the accelerometers and
gyros. These are the random forcing functions that are calibrated at the start of a
mission but then slowly wander away from their initial values during the course
of the flight. It is especially important to keep calibrating these forcing func-
tions continuously if it is anticipated that the INS might have to dead-reckon in
the free-inertial mode for a significant time period without GPS measurements.
If the biases are estimated in flight, the system errors grow more slowly during
the free inertial period. Usually a Jirst-order Markov or a random-walk process
is used to represent an inertial instrument error source, so only one error state is
added to the model for each source. Such an elaboration of the inertial system
error model would add an additional three error states for three gyro drift rate
errors and an additional three error states for three accelerometer measurement
errors.
The primary trend is to employ low-cost, low-to-medium performance iner-
tial components to bridge the periods of GPS outage in a cost-effective manner.
These instruments will likely employ micro-machined silicon technology for
the accelerometers and fiber optics for the gyros. In 1996 these technologies
PRACTICAL CONSIDERATIONS 91
did not yet obtain the high-accuracy performance levels of the more expen-
sive technologies used in the past. These instruments require more elaborate
error modeling, including such additional states as scale factor error, mutual
mechanical misalignments of instrument sensing axes, and bias changes as a
function of measured environmental variations to obtain good calibration.
In addition to the instrument error states discussed above, the effects of
unknown variations in the gravity vector should be considered as error states in
some Kalman filter designs that include an inertial system. This consideration is
important because the gravity disturbance vector introduces errors in the force
measurements made by the accelerometers that can significantly affect some
applications.
GPS Error States Besides the receiver clock error model, filter designers have
also been concerned with the error in the pseudorange measurement and the
error in the Doppler or integrated Doppler (delta-range) measurement to each
satellite due to residual satellite clock and orbit errors and transmission path
effects. Inclusion of an error state for each of these measurements can obtain a
calibration of the measurement under certain conditions. For example, precise
knowledge of the vehicle location and velocity at the measurement time can
provide such a condition. Further, measurements from satellites currently being
tracked can be used to calibrate the measurements from a new satellite when
an initial track is established.
Other refinements of the basic !!-state filter are also possible. Refinements
in the error dyanmics of the model itself may be required in some applications.
The example used a north-oriented coordinate frame, but a different frame can
be used (e.g., the azimuth-wander frame discussed in Chapter 7). The !!-state
model includes all the basic quantities that need to be estimated since the filter
estimates position, velocity, instrument or platform coordinate frame orienta-
tion, and GPS receiver clock errors. All of these quantities are observable in
the !!-state filter once measurements to four or more satellites have been made
available. The observability implies that the usual problems of platform level-
ing, gyrocompassing, and damping the Schuler oscillation are all taken care of
automatically by the Kalman filter without any special ad hoc procedures. In
the past, the principal constraint on the elaboration of the error model has been
limited computer resources. However, in the 1990s, Kalman filter designs based
on models with several tens of error states were being implemented.
• Measurement synchronization
92 MULTISENSOR NAVIGATION SYSTEMS
• Measurement editing
• Tuning parameter adjustment
• Filter equation implementation
(3.30)
(3.31)
only guarantees this condition when the gain is optimal. Note that even
though the equation implemented in the computer is for the optimal gain
(3.32)
\ Raw Sensor
Measure nent Inputs
•
•
........
......
"Master''
Kalman Filter
Estimates of Vehicle
• ~
Position, Velocity &
I
Orientation
L--..-----
""
VI
96 MULTISENSOR NAVIGATION SYSTEMS
The future will see significant activity focused on integration of the naviga-
tion sensors with other avionics sensors on the aircraft. In military systems, the
information collected by many sensors is being fused into a central data base to
ensure their optimal use. Civil aircraft systems are likely to see a slower pace
of multisensor integration, largely because of the procurement practices of the
airlines. They purchase equipment from different suppliers based on compet-
itive prices and a desire to interchange "black boxes." Furthermore, reliabil-
ity and safety considerations may lead to a preference for functinoal redun-
dancy and a partially federated approach. In 1995, the civil aviation industry
began to investigate centrally fused data in a program called Integrated Modu-
lar Avionics (IMA), which is somewhat equivalent to the functional integration
and resource-sharing programs represented by the military ICNIA and ICNIS
programs. These efforts are likely to continue and accelerate in the twenty-first
century.
High-powered, low-cost digital processors will facilitate large-scale integra-
tion and enable redundant sensors w achieve a high level of fault tolerance. The
self-contained black-box subsystems that drove the federated avionics systems
architecture in the past will eventually disappear as the more information-rich
centralized system architecture evolves. Significant savings should be realized
in cost, size, weight, and power. The one major disadvantage of the federated
system is that the central-computer software engineers must be expert in each
sensor. Libraries of standard sensor modules may appear near the turn of the
twenty-first century.
Because of its significant cost/performance benefits, the combination of a
low-cost, low-to-moderate performance inertial measurement unit and a GPS
receiver will be a widely used multisensor system for many types of air vehicles
for years to come.
PROBLEMS
3.1. Refer to the simplified inertial system error model shown in Figure 3.14,
and make the following assumptions:
PROBLEMS 97
Compute the steady state values of (east) tilt ¢x and azimuth misalign-
ment ¢, that will result from observing (north) velocity error Ll.v during
the process of initial alignment of the inertial system
Hint: Also see Chapter 7. Ans.: ¢, = 200 arcsec ¢x = 5 an·sec
3.2. Assume an inertial system has been initially aligned as in Problem 2.1
(equilibrium conditions obtained), with a variance in (east) tilt of [5
arcsecf and a variance in azimuth misalignment of [200 arc-sec] 2 , and,
additionally, there exists an independent error in (north) position, with a
variance of [500 feet] 2 . What are the variances in tilt, azimuth misalign-
ment, and position error when a star tracker is integrated with the inertial
system and a star is observed directly overhead (along the local vertical),
if the tracker measurement error is zero and processed by a Kalman filter?
Ans.:a;, = [3.5 arc sec ]2 , ai,
= [353 .ft] 2 and a;,
= [200 arcsecf
What are these variances if a star is observed on the horizon along the
local north line?
Ans.: a;, = [5 arcsec] 2 , ai, = [500 .ftf and a;, = 0
3.3. Consider an inertial system initially aligned as in Problem 3.1 but inte-
grated with a reference speed sensor (e.g., a Doppler radar), where the
vehicle accelerates in the eastward direction at Ax = I 0 fps 2 for I 0 sec-
onds to a velocity of Vx = I 00 fps. What are the steady-state values of the
tilt and azimuth misalignment variances if Doppler speed measurements
with no error are processed by a Kalman filter?
Ans.: a;, = [5 arcsec] 2 and a;, = [200 arcsecf
What are the steady-state values of tilt and azimuth misalignment variances
if y (north) position change measurements (e.g., as obtained from a GPS
receiver) with no error are processed by Kalman filter?
Ans.: a;, = [5 arcsecf and a~, = 0
3.4. For an integrated reference speed sensor and inertial system, where the
speed sensor has a constant bias error with variance of II fpsf and the
inertial system has a constant gyro drift rate with variance of [0.0 I o /hr] 2 ,
compute the steady-state variance of the system velocity error, a~v·
Hint: Refer to Fip,ure 3.6. Ans.: azv = [.707 fpsf
What is the steady-state variance of this velocity error if, in addition, a
98 MULTISENSOR NAVIGATION SYSTEMS
4.1 INTRODUCTION
This chapter discusses the basic principles of terrestrial radio navigation sys-
tems, the radio propagation and noise characteristics and the major system per-
formance parameters. The systems described in detail include all the impor-
tant point source systems, such as direction finders, VOR, DME, and Tacan;
and the hyperbolic systems, such as Loran-C, Decca, and Omega. All of these
systems have been used worldwide and have provided accurate and reliable
positioning and navigation in one or two dimensions for many years. In 1996
hundreds of thousands of civil and military aircraft throughout the world were
equipped with these systems, and many of these will be used for years to come.
Satellite radio nevigation systems (e.g., GPS) are discussed in Chapter 5 and
military integrated radio communication-navigation systems are described in
Chapter 6.
Display or
Modulator Transmiuer - R=iver - Processor - Data Bus
Interface
Abbreviation
Name Frequency Wavelength
Very low frequency VLF 3 to 30kHz 100 to 10 km
Low frequency LF 30 to 300kHz 10 to I km
Medium frequency MF 300 to 3000 kHz I km to 100m
High frequency HF 3 to 30 MHz 100 to 10m
Very high frequency VHF 30 to 300 MHz 10 to I em
Ultrahigh frequency UHF 300 to 3000 MHz I m to IO em
Superhigh frequency SHF 3 to 30 GHz IO to 1 em
Extremely high
frequency EHF 30 to 300 GHz 10 to I mm
At the higher frequencies, the following letter designations for certain fre-
quency bands have been widely accepted, although they do not have official
status (they are frequently related to standard wave-guide sizes) [35]:
GENERAL PRINCIPLES 101
With progress constantly being made at still higher frequencies, other sys-
tems of nomenclature will be required. However, radio navigation, as defined in
this chapter, is confined primarily to the bands lying between VLF and UHF,
where the above nomenclature is likely to remain in use. Regardless of fre-
quency, the following general rules apply in free space:
where R is the range between antennas in the same units as those for the
antenna area.
3. Multiple antennas may be used at both ends of the path to increase the
effective antenna area. Such increases in area produce an increase in
directivity or gain and result in more of the transmitted power reach-
ing the receiver. The gain G of an antenna (in the direction of maxi-
mum response) is equal to its directivity D times its efficiency. The maxi-
mum effective aperture or effective area of an antenna is equal to D'A/47r.
It is defined as the ratio of the power in the terminating impedance to
the power density of the incident wave, when the antenna is oriented
for maximum response and under conditions of maximum power trans-
fer [ 19j. It is also defined as the physical area times the antenna aper-
ture efficiency (or absorption ratio) [19, 36]. The directivity or gain of
an antenna is usually expressed as a ratio with respect to either a hypo-
thetical isotropic radiator or a half-wave dipole. A dipole has an effective
area of about 0.13 times the square of the wavelength [36]. A transmitter
102 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
Thus, for fixed effective antenna areas, the power transferred from trans-
mitter to receiver increases as the square of the frequency. However, this
is accompanied by a corresponding increase in directivity. Such direc-
tivity is of no concern in fixed point-to-point service and is of advan-
tage in reducing external noise pickup. In many moving-platform appli-
cations, such as aircraft, a high level of directivity is a distinct disad-
vantage. However, when the use of tuned dipole antennas is assumed,
the power transferred decreases as the square of the frequency. This is
seen from Equation 4.1 where the receiver antenna area for a dipole, (i.e.,
0.13 A2 ) is substituted in the numerator.
4. The minimum power that a receiver can detect is referred to as its sensitivity.
Where unlimited amplification is possible, sensitivity is limited by the noise
existing at the input of the receiver. Such noise is of two main types:
a. External. Due to other unwanted transmitters, electrical-machinery
interference, atmospheric noise, cosmic noise, and the like.
b. Internal. Depending on the state of the art and approaching, as a lower
limit, the thermal noise across the input impedance of the receiver,
which is given by
(4.3)
where
NP is the noise power (in watts)
k is the Boltzmann's constant (1.38 X w- 23 JoulesjKelvin)
T is the temperature (in Kelvin)
A.f is the bandwidth (in Hertz)
No
NF= - - - - (4.4)
kTob.fGa
where No is the noise power out of a practical receiver and N 1 is the noise
power out of an ideal receiver at standard temperature T 0 , of available
gain Ga, and of bandwidth D.f.
5. The minimum bandwidth occupied by the system is proportional to the
information rate. For most navigational purposes, the necessary informa-
tion rate is quite low. For instance, to navigate in a given direction to an
accuracy of 500 ft with an aircraft that cannot change its position more
than 500 ft in that direction in any one second, new information is needed
only once per second. However, most practical systems have employed
many times this minimum bandwidth. The reasons include (a) the need for
other services, such as communications on the same channel, (b) the use
of pulse techniques to aid in resolving multiple targets and to reduce the
effects of multipath transmission, and (c) the use of spectrum-spreading
techniques to improve signal-to-noise ratio (S/N), accuracy of range mea-
surements, reduction of effects due to interference of site errors. (Spread-
ing the spectrum beyond that needed by the information rate itself has the
same effect as increasing the power, provided that optimum techniques
are used at each end of the link [2, 27].)
where
PT is the transmitter power
PN is the noise power in receiver
(SjN)REQ is the required signal-to-noise ratio in receiver
NF is the receiver noise figure
FN is the noise improvement factor due to modulation method and
bandwidth spreading (e.g., frequency modulation)
is the transmitter antenna gain
104 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
10 37 63 37 62
100 58 99 57 92
500 87 195 71 125
1000 110 245 79 145
2000 132 86 165
than the inverse square, as would occur in free space. In addition, further losses,
increasing with frequency, are encountered. Table 4.1 gives some typical exam-
ples. At low frequencies, ranges up to 5000 miles or more are obtainable, if suf-
ficient power can be generated to overcome atmospheric noise, path attenuation
and to compensate for low antenna efficiencies [42].
Second, at low frequencies, it is physically difficult to construct a vertical
transmitting antenna large enough to be half a wavelength (or its electrical
equivalent, i.e., a quarter-wave antenna above a perfectly conducting plane).
Therefore, the antenna is generally much shorter than the ideal and is res-
onated to the operating frequency by external series inductance of the lowest
possible losses. The result, despite the best engineering practice, is the radiation
of considerably less power over a very narrow bandwidth than that generated
by the transmitter. Nevertheless, ground-wave service is sufficiently attractive
for many applications so that low efficiencies are tolerated in some lower-fre-
quency applications.
Third, in most parts of the world and at most times of the year, atmospheric
noise at low frequencies is so much greater than receiver noise that additional
transmitter power must be used. This noise is generated mostly by lightning
flashes. As shown in Table 4.2, at the latitude of the United States, atmospheric
Sky Waves In a region lying between 50 and 500 km above the Earth's sur-
face, radiation from the sun produces a set of ionized layers called the iono-
sphere [38, 42]. The location and density of these layers depends on the time of
day and, to a lesser extent, the season and the !!-year sunspot cycle. The iono-
sphere acts as a refractive medium; when the refractive index is high enough
in relation to the frequency of a radio wave, it bends the radio wave and will,
under favorable conditions, return the wave back to Earth.
Figure 4.2 shows a simplified picture of the geometry involved. At A, the
radio wave strikes the refractive layer at too steep an angle and, although it is
bent, is not sufficiently affected to return to Earth; it continues out into space
(unless it encounters a more heavily ionized layer further out). At B, the radio
wave strikes the refractive layer a1: a more oblique angle, is bent sufficiently to
travel somewhat parallel to Earth, and is finally bent sufficiently to return to
Earth. At C, the wave arrives at the refractive layer with glancing incidence and
immediately returns to Earth. At D, the refractive index is too low in relation
to frequency to seriously deflect the radio wave, which then travels on out to
space; generally, this happens at frequencies above 30 MHz.
From this geometry it is evident that return to Earth occurs only at some
minimum distance for a given frequency and degree of ionization. This is called
the maximum usable frequency for that distance. Signals at higher frequencies,
GENERAL PRINCIPLES 107
·: :...
·.·.·.' :: . ... ... .
. . . . .. .
. . .' .
: ·.·.· .·
if returned at all, will be returned only at greater distances. This critical distance
is known as skip distance; inside it there is no return to Earth at the particular
operating frequency. If more than one ionizing layer are present, there may be
various skip distances for the same frequency.
At those frequencies and distances where ionospheric reflection occurs, the
attenuation of the radio signal is only that due to the spreading out of the power
over the surface of the Earth and is, consequently, proportional to distance.
Conversely, as indicated in the previous paragraphs, ground-wave attenuation
is very much greater, except at the lowest frequencies. At frequencies of around
I MHz, the signal level produced at the receiver by the two types of transmis-
sions is likely to look like that shown in Figure 4.3.
As the frequency increases, the ground-wave curve will move to the left and
the sky-wave curve to the right, leaving a gap (due to skip distance) where nei-
ther wave produces a usable received signal. In the region where the ground
wave and sky wave are about equal, severe fading will occur due to the ran-
domly varying phase of the sky-wave signal with respect to the ground-wave
signal. Even when sky-wave signal strength is adequate, serious distortion of
its modulation may occur due to the different paths simultaneously traveled by
the signal between transmitter and receiver. These are called multipath effects.
The differential time delay between these paths may reach several milliseconds,
thus preventing faithful reproduction of modulation frequencies above a few
hundred Hertz.
Therefore, sky-wave transmission is quite variable, and its efficacy is highly
dependent on the distance to the receiver, the frequency used, and the time of
day. For these reasons, the general practice in the 3- to 30-MHz communica-
tion band has been to use receivers and transmitters that would readily tune
over the whole band and to change frequencies from hour to hour, depending
on the distances required and on the condition of the ionosphere. Much work
108 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
"'c::
Oil
v;
-o
"'
>
iii
u
"'
a::
has gone into the creation of charts predicting maximum usable frequencies and
skip distances [42]. More recent developments include propagation-frequency
evaluators that quickly evaluate (by frequency-scanning techniques) the best
path to be used for communication between two points at a given time. By
use of such techniques, ionospheric reflection has been a major long-distance
communication aid. Until the advent of wideband submarine cables and com-
munication satellites, these frequencies were the mainstay of the transoceanic
telephone and radio networks. The highly variable characteristics of the iono-
sphere, which cause different frequencies to travel by different paths, led to the
development of many ingenious sreerable antenna systems for this service. The
use of ionospheric reflection for navigation systems has been confined almost
exclusively to ground-based direction finders.
Conversely, sky-wave transmission is considered a handicap, rather than an
aid, to those navigation systems that depend on groundwaves. In such systems,
the almost direct, reasonably predictable ground wave is contaminated by sky-
wave energy that has arrived by a devious path. The mixture of the two often
produces serious errors not only in distance but also in bearing measurements,
since the effective reflecting point is not necessarily on the vertical plane join-
ing the transmitting and receiving stations. Methods for reducing such sky-wave
contamination include (I) the use of tall antenna structures for improved ver-
tical directivity resulting in transmission fields being concentrated along the
ground and less toward the sky, and (2) the use of only the leading edge of
pulse transmission, since this edge arrives sooner by ground wave then by sky
wave and is, therefore, uncontaminated (this usually requires greater bandwidth
than that required by the information rate).
but from approximately 100 MHz to 3 GHz, the transmission path is highly pre-
dictable and is unaffected by time of day, season, precipitation, or atmospher-
ics. Above 3 GHz, absorption and scattering by precipitation and by the atmo-
sphere begin to be noticed, and they become limiting factors above I 0 GHz.
Furthermore, above that frequency, atmospheric absorption does not increase in
a smooth manner but rather is characterized by narrow peak-absorption bands
and by narrow "windows" of relatively reduced absorption. Tables 4.3 and
4.4 [ 14] show attenuation effects due to fog and rain (in addition to free-
space loss).
Because of absorption above I 0 GHz, transmission at such frequencies is
severely limited within the Earth's atmosphere. High-flying aircraft and space
vehicles of course are under no such restrictions.
In designing antenna systems for line-of-sight frequencies, it often happens
that due to the relatively short wavelengths, the antenna is spaced away from a
reflecting object such as the ground by a critical number of wavelengths, which
has a marked effect on the overall antenna pattern. For instance, at 1 GHz, the
wavelength is about I ft. If this practice were used at 1 GHz, a quarter-wave
structure might be built with its base on the ground. Since even nearby blades
of grass would seriously mar its performance (not to mention persons and vehi-
cles moving about nearby), this would obviously be impractical. Instead, such
an antenna would likely be mounted on a pole, say, 10 ft high, so as to clear the
10 3 15 37 100
10 4 7 12 20 43
10-5 3 6 9 20
10 6 3 4 8
10-7 4
I. Raising the antenna high enough above the ground, m relation to fre-
quency, to make the null structure so fine that the slightest irregularities
in the reflecting surface will break up the null pattern. At heights of I 00
wavelengths or more, the problem can usually be ignored.
2. Placing a horizontal counterpoise immediately below the antenna so as
to make its effective height quite small. The counterpoise shortens the
path of the reflected ray, thereby raising the angle at which cancellation
occurs. To be effective, such a counterpoise must be many wavelengths
in diameter.
3. Using high vertical antenna directivity (either by antenna arrays or by
parabolic reflectors) and then pointing the resulting narrow antenna beam
slightly above the horizon. The uptilt reduces the energy striking the
ground and, therefore, reduces the reflected wave.
4. Making the null on one frequency occur at the same time that a maximum
is occurring on another by using frequency diversity or wide-spectrum
modulation which allows several frequencies to be used simultaneously
[I].
5. Introducing vertical diversity up to an appreciable vertical angle via two
antenna systems, one at half the height of the other with the null of one
corresponding to the maximum of the other. However, the frequencies fed
to these two antenna systems cannot be coherent; otherwise a new null
pattern will appear at angles where the signals are otherwise equal. This
method is, in general, limited to receiving systems where two separate
receivers can be used.
6. Attempting other forms of diversity that make use of two or more paths
simultaneously. The term "diversity" in radio propagation refers to this
use of paths with different frequencies, polarizations, and so on in order
to make reflections occur at different points on each path.
R in nautical miles
h in feet
Figure 4.6 Line-of-sight range.
+------
Direction of arrival
of signal
+------
Receiver B
------
------
Figure 4.7 Direction-finder principle.
SYSTEM DESIGN CONSIDERATIONS 113
"'N
I
I
(a) (b)
(c)
BD- AD= k on any one line
Three pairs of stations provide a fix
(d)
Figure 4.8 Common geometric position fixing schemes: (a) Rho-theta (p-{}); (b) theta-
theta({}-{}); (c) rho-rho (p-p); (d) hyperbolic.
Rho-theta systems provide a unique fix from a single station, and the LOPs
always cross at right angles. Theta-theta systems provide a unique fix from
two stations. The geometric accuracy is highest when the lines cross at right
angles and is poor on a line connecting the stations. Rho-rho systems provide
an ambiguous fix from two stations and a unique fix from three stations. Geo-
metric accuracy is greatest within the triangle formed by the three stations and
gradually decreases as the vehicle moves outside and away from the triangle.
The hyperbolic system uses LOPs that each define a constant difference in
distance to two stations. Such systems operate under conditions where the deter-
mination of absolute distance to the station is impractical. Three pairs of stations
are needed for a unique fix; however, for many practical applications, two pairs
suffice. Geometric accuracy is very much a function of the relative station loca-
tions. Poor geometry leads to a property frequently called geometric dilution of
precision (GDOP).
Another method, called pseudoranging, has been developed and is used in
certain modern radio navigation systems, such as GPS (Chapter 5) and one
114 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
mode of JTIDS-RelNav (Chapter 6). In this method the user receiver and the
reference station(s) are assumed not to be synchronized in time. By measur-
ing several (in general, at least four) such pseudoranges (versus true ranges
when time synchronization does exist), the user's three-dimensional position
and its time offset (from the transmitter or system time) can be determined
(Section 2.5). Three such pseudoranges are sufficient to determine the user's
two-dimensional position (provided that the user's altitude is known). Finally,
the term direct ranging (rho-rho and rho-rho-rho) has been applied to a hyper-
bolic system used in a true ranging mode by achieving some form of time syn-
chronization (see Sections 4.5.1 and 4.5.2).
The measurement of distance and direction by radio gives accurate results
only if the radio path between the points being measured is direct and the propa-
gation speed is known. In practice, the path between a transmitter and a receiver
on the Earth's surface may be quite devious and may, in fact, be a combination
of many paths. The problems of reducing such multipath transmissions and of
recognizing the direct path are major reasons for the multiplicity of radio nav-
igation systems that have been proposed or are in use.
identifies two or more possible positions of the vehicle, with the same
set of measurements, with no indication of which is the most nearly cor-
rect position. (This is not a problem with Loran-C, since the ambiguous
fix is a great distance from the desired fix. Ambiguous lines of position
(LOP) occur in the Omega system, since there is no means to identify
particular points of contact phase (lanes) that recur throughout the cover-
age area. Because of this ambiguity, Omega receivers must be initialized
to a known position, and the lanes counted as they are crossed . )
6. Capacity. Capacity is the number of users that the radio-navigation system
can accommodate simultaneously. For example, there is no restriction on
the number of receivers that may use Loran-C, Omega, or a VOR station
simultaneously; on the other hand, DME and Tacan are currently limited
to about II 0 users for traffic handling.
4.4.1 Direction-Finding
Direction-finding represents the earliest use of radio for navigational purposes;
it continues to perform a useful function, particularly in those parts of the world
that have not yet adopted the more specialized navigation aids. Its chief attrac-
tion lies in the fact that, with the proper receiving equipment, the direction
of a transmitter can be found. Such transmitters do not necessarily have to be
specially designed for direction-finding; they can be broadcast stations, com-
munication stations, navigation stations, or any other kind of radiating system.
The chief drawback of direction-finding is that quite elaborate receiving
equipment must be used if the best accuracy is to be obtained. Most aircraft
are unable to accommodate such equipment. Direction-finders for aircraft nav-
igation may, therefore, be grouped into two broad classes:
1. Fixed loop. Intended for flying radial courses to and from the ground
station by orienting the aircraft for minimum signal
2. Rotatable loop. Hand-operated systems that were abandoned because of
the work load they imposed on the pilot
3. Rotating loop. Driven by a motor and forming part of a servo system
that automatically rotates the loop until a null is found and then stops,
sometimes referred to as a radio compass. Early loops were about nine
inches in diameter and were housed in teardrop-shaped plastic enclosures
about one foot away from the aircraft skin.
4. Fixed, crossed loops, with a motor-driven goniometer. Forming part of a
servo system that automatically displays bearing in the cockpit. The prime
advantage of this system over those using the physically rotating loop is
that all moving parts (except the indicator) are in the radio-receiver box.
Antenna projection from the aircraft with such a system in as low as one
inch, with horizontal dimensions of about one foot. Typical airline-type
equipment weigh less than 20 lb.
c/
Direction of
transmitter
~~~--~
Figure 4.9 Direction-finding loop.
Combined
Pattern of
/ .. ------ ......... ,/patterns
omnidirectional
antenna alone \
\
\
I
Plane I
of loop :n
I
I
I
I
Pattern of /
I
------ _,..,"
loop alone
Figure 4.10 Loop and sense antenna patterns: A = left-hand loop pattern; B = right-
hand loop pattern, 180° out of phase with A; C = omnidirectional pattern, 180° out of
phase with A; D = C + B - A.
POINT SOURCE SYSTEMS 119
away from it. A sense antenna can be added when the signal ambiguity must
be resolved. The sense antenna adds an additional 90° phase shift. As the loop
changes direction, its phase will vary with respect to the constant sense antenna
voltage resulting in the cardioid pattern shown in Figure 4.1 0. The combined
pattern has only one null position. Since the omnidirectional antenna, and its
phase and amplitude relation to the loop, are less precisely definable than the
loop itself, it is customary to use the loop alone for precise directional measure-
ment. The sense of the bearing can then be determined by coupling the vertical
or sense antenna to the loop and rotating the loop 90° in a specified direction,
noting whether the signal increases or decreases as the loop antenna is rotated.
nal crossover direction is found by turning the aircraft toward the transmitting
station. Equipment designed for both direction finding and homing uses a rotat-
ing antenna that generates a similar pair of cardioid patterns, whose equisignal
crossover direction is found. Accuracy is about sc along the axis of the aircraft
but reaches 30° broadside. The direction-finding attachment is carried by many
U.S. military aircraft and is useful for air-to-air direction finding and homing
during rendezvous and refueling. It is also of value in locating downed flyers
who carry small UHF rescue beacons.
Cone of silence
I
I I
I I
wgggyc~,--
Figure 4.11 Nondirectional beacon, vertical pattern.
POINT SOURCE SYSTEMS 121
popularity because (l) they are inexpensive, (2) they are omnidirectional, and
(3) they place responsibility for accuracy entirely on the airborne receiver.
Nondirectional beacons are probably the least expensive way by which a
government can claim that it has equipped its airways with "radio aids to nav-
igation." In 1996, many thousands were in service around the world, and the
United States maintains approximately 177,000 nondirectional beacons for civil
aviation use. This number is expected to increase by about 7000 a year for the
next ten years [48].
Airway To four-course
-----------~range, VOR,
or ILS
Fan-marker pattern,
viewed from above
Figure 4.12 Fan-marker pattern.
122 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
Antenna patterns
Along Across
airway airway
Direction
of airway
cned by capacitive loading, recessed into the aircraft, and covered by a dielec-
tric sheet. Streamlined antenna packages are 3 ~ by 6 in. weighing 18 oz. The
marker beacons are gradually being phased out as an en-route aid in view of the
implementation of area-coverage systems, such as VOR/DME RNAV, Loran-C,
and GPS. However, along instrument landing approaches, the 75 MHz marker
remains a standard piece of equipment (Chapter 13).
consistent and is limited by only two major factors: (1) propagation effects,
including vertical pattern effects, and site and terrain errors, and (2) instrument
errors in reading 30-Hz phase differences in the airborne equipment.
Antenna patterns
Feed point
8 " '" '
View from above
Figure 4.14
m
Alford loop.
Vort'"'
124 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
Bndge Bridge
Continuous wave,
modulated by tone
wheel, voice, 1020 Hz
Identity tone
9960Hz
FM ±480Hz
at 30 Hz Modulation
eliminator
Tone wheel
Continuous wave
northwest-southeast pair of Alford loops and then, 90" later, to the northeast-
southwest pair of Alford loops. When combined with the modulated energy
applied simultaneously to all loops, this variation generates a rotating cardioid.
Each pair of loops is fed via a balanced bridge network. Each bridge has three
arms that are each abou( one-quarter wavelength long, the fourth arm being half
a wavelength longer. Energy fed into one corner of the bridge does not appear at
the diagonally opposite corner. The bridge, therefore, allows the mixing of two
signals and application of the result of two loads without the loads affecting
each other and without the signal sources affecting each other. The phasing
between tone wheel and goniometer and the physical placement of the Alford
loops are such that the two 30-Hz signals are exactly in phase when viewed
from magnetic north.
This seemingly elaborate arrangement serves two main purposes:
POINT SOURCE SYSTEMS 125
The four Alford loops are arranged in a tight square and then placed half
a wavelength above a metal-mesh counterpoise about 39 ft in diameter. This
counterpoise also acts as the roof of the transmitter house. The loops are pro-
tected from the weather by a plastic randome, often hemispherical in shape. If
a Tacan antenna is collocated with the VOR, the randome is conical in shape,
somewhat resembling an Indian tepee.
Audio output
Receiver
108-118 MHz
Bearing
reference frequency. After filtering, this is compared with the variable phase.
Several grades of receivers are currently in use.
The airline type of equipment uses a remotely tuned crystal-controlled super-
heterodyne receiver and has at leasr two types of display. One display compares
one 30-Hz sine wave with the other 30-Hz sine waves, the two signals being
brought into phase by a motor-driven phase shifter forming part of a servo loop.
The shaft position of this motor, therefore, displays bearing directly and may
be remoted by selsyns to other parts of the aircraft and to auto-pilots. Another
display shows (on a vertical left-right needle) the phase difference between
one 30-Hz signal and a manually phase-shifted 30-Hz signal representing the
desired bearing. The sensitivity of the vertical needle is usually arranged for
a full-scale deflection of ±I oc around the manually selected bearing and thus
shows angular deviation from the desired track.
The simplest types of receivers use manual tuning and only the left-right
type of display around a manually selected bearing. Both types of receivers are
commonly arranged to also receive the 108- to 112-MHz instrument-landing-
system localizer signals. Typical receivers weigh 20 lb for the airline type and
5 lb for the simplest type, exclusive of antenna. Over 200,000 airborne sets
have been installed, about half of them for light aircraft.
It was previously mentioned that one of the problems of the VOR is the
difficulty of accurately measuring phase shifts at 30Hz. Much circuit refinement
has taken place for the better grades of receiver. This includes, for instance,
the use of identical circuits for both 30-Hz signal paths wherever possible so
that temperature effects will be common to both. The result is that instrument
accuracy of better than 1o is achieved in airline-type equipment.
480 wavelengths per second at VOR radio frequencies. The 9960-Hz frequency
difference is consequently varied by ±480 Hz at a 30-Hz rate, with a phase
dependent on the bearing of the receiver.
In the receiver, the output of the amplitude-modulation detector contains all
the signals present with the conventional VOR. Phase comparison between the
two 30-Hz sine waves is performed as before, the only difference being that the
30-Hz amplitude-modulated signal is the reference and the 30-Hz frequency-
modulated signal is the variable. Since the instrumentation is concerned only
with the difference between the two, normal operation results with a standard
VOR receiver.
However, since the aperture of the ground antenna is approximately five
wavelengths, as compared with less than half a wavelength with the four Alford
loops in a standard VOR ground station, a tenfold reduction in site error is the-
oretically possible. Actual measurements at formerly "impossible" sites verify
this. At a good site, maximum deviations measured during a 20-mi orbital flight
were reduced from 2.8° with a standard VOR to 0.4° with a Doppler VOR [3].
Residual errors can probably be reduced to 0.1 o.
The importance of the Doppler VOR lies in the improvement it provides
without any change being made to the airborne equipment. Every airborne set
can benefit from it.
Distance
reading
Airborne
interrogator
Aircraft skin
Ground
transponder
or beacon
100 being a more typical number. The pulse-repetition rate of the interrogators
is deliberately made somewhat unstable, and the interrogator is designed to
recognize only those replies whose pulse-repetition rate and phase are exactly
the same as its own.
In any line-of-sight geographical area, there is the possibility of providing
136 beacons, each handling I 00 or more aircraft. Since each beacon's duty cycle
is only 2% under these conditions., room exists to expand the system to handle
heavier traffic. Modern techniques permit the airborne interrogation rate to be
decreased substantially.
To and from
antenna
Beacon
identrty
code
Tacan
bearmg
01stance
drsplay
I
130 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
1--:r~
noise
I I
150 .;ec~~~
~ t'"'-less than beacon recovery time
I I
I I
~.,...,...~~~~.,..,__,..
~:
Own Own~:
interrogations (5) replies ( 4) 1
Many different forms of circuit have been devised for these functions. They
all depend on the sequence of wave forms shown in Figure 4.20. This figure
shows five consecutive snapshots of an imaginary oscilloscope whose sweep is
started by the airborne interrogation from a single aircraft and whose deflection
circuit is connected to the output of the receiver in that aircraft.
In this instance, if one assumes a maximum desired range of 200 nm, the
sweep is 2400 p.sec long. Since the ground beacon is transmitting an average of
3000 pulse-pairs per sec, each sweep will display, on the average, about seven
pulses. These will be quite randomly spaced, except those generated in response
to our own interrogation. At an interrogation rate of 30 per sec, even the fastest
aircraft does not move by as much as a pulse width from one interrogation to
the next. The desired replies therefore occupy an almost fixed position on the
oscilloscope display, whereas those intended for other aircraft move in a random
manner. The dotted line shows the fixed (or slowly changing) position of the
desired reply. On scan 3, the desired reply is missing; this is because the beacon
has just replied to another aircraft and has not yet "recovered." Recovery time is
typically on the order of I 00 p.sec. Desired replies may also be missing because
of other random effects. However, all airborne DME ranging circuits are based
on the principle that, within a given time slot, many more desired replies will
be received than undesired replies.
The basic objective of all DME ranging circuits is to locate the time slot
in which the desired replies are actually occurring. This is the search process,
and it is usually conducted at the highest permissible pulse-repetition rate ( 150
pulse-pairs per second) in order to save time, which, depending on the tech-
nique, may vary from I to 20 sec. Once this time slot has been found, the track
mode commences and can be conducted at a much lower pulse-repetition rate,
usually between 5 and 25 pulse-pairs per second.
Search is typically performed as follows: A gate is generated I 0 p.sec wide,
POINT SOURCE SYSTEMS 131
vides 3-kw peak output, together with an antenna gain of 9 dB. Otherwise, the
ground circuits follow the same principles as the airborne ones, with a 63-MHz
intermediate-frequency receiver amplifier being used. The number of aircraft
that a beacon can handle is usually based on the assumption that 95% of the
aircraft will be in the track mode at not over 25 interrogations per sec; 5% are
in the search mode, at not over !50 interrogations per sec. For I 00 aircraft, this
means about 3000 pulse-pairs per sec.
The duty cycle of the ground transmitter is therefore much greater than that
of the airborne equipment, and the average power consumption is also greater.
Most beacons are operated on the constant-duty-cycle principle, whereby
receiver gain is increased until 3000 pps appear at the output of the receiver. In
the absence of interrogation, these pulses will all be due to receiver noise; with
interrogation from less than I 00 aircraft, they are a mixture of noise and inter-
rogations; with interrogation from more than I 00 aircraft, they are the inter-
rogations from the 100 nearest aircraft. After the 3000-pulse limit is passed,
the gain is automatically reduced. This constant-duty cycle has the following
advantages:
For the simplest, most reliable circuitry, the beacon is arranged not to
receive while transmitting (self-oscillation could otherwise result); furthermore,
to reduce interrogation by multiputh echoes of strong interrogation pulses, it is
desirable to reduce receiver gain for a short while after each genuine interroga-
tion. Some interrogations are consequently lost; the amount of this countdown
is typically on the order of 20%. Thus, an airborne equipment interrogating at
25 pps receives only 20 replies. Airborne tracking circuits are, however,
designed to operate at this reduced rate.
The delay between transmission and reception is nominally 50 fJ.Sec. For
greatest accuracy, this must be maintained constant; considerable circuit refine-
ment is used to retain this value, independent of interrogation strength and envi-
ronmental effects. Typical en-route-type beacons exhibit a total variation for
±0.5 fJ.Sec, corresponding to a distance error of ±0.04 mi [30]. Beacons asso-
ciated with instrument-landing systems may be designed to be more accurate,
due to the smaller spread of interrogation-signal levels. The ICAO requires an
overall system accuracy of 0.5 mi or 3%, whichever is greater.
Under the control of an external keyer, usually common to the associated
VOR, the beacon transmits an identify signal. Typically, this occurs for about
3 sec every 37 sec. During this time the random pulses are replaced by regu-
POINT SOURCE SYSTEMS 133
larly spaced pulses at 1350 pulse-pairs per sec. These activate a 1350-Hz tuned
circuit in the aircraft and are keyed with a three-letter Morse code, ~ sec per
dot and ~ sec per dash. During this time, the airborne ranging circuit is in the
memory condition.
Since the DME system, unlike the VOR system, is not a passive system, it
has an inherent capacity limitation. The value generally quoted is 110 aircraft
per beacon.
I. Because of its higher frequency (960 to 1215 MHz versus 108 to 118
MHz), the Tacan beacon antenna can be smaller, it is therefore more suit-
able for shipboard and mobile use.
2. The multilobe principle, to enhance bearing accuracy, is built into all
equipment, ground based and airborne.
3. Both distance and bearing are obtained via the same radio-frequency
channel, providing certain equipment economies.
The system is in general use by the U.S. Navy and Air Force, and by NATO
military forces. In 1996, over 800 facilities were maintained for the U.S. DoD
with a DoD user population of 13,000 [48].
Plastic cylinders
rotating at 15 rps ~----- ------ Central radiator, stationary
~- One parasitic element, rotating
--........JIUL--~
A
Nine parasitic elements, rotating
Speed
control
Receiver
the 960- to 1215-MHz range. Equipment has been built with from 1 to 11 ver-
tical elements, depending on the kind of site for which the set is intended.
All transmission and reception is by this central radiator. At a radius of about
3 in. and usually mounted on a plastic cylinder is the 15-Hz parasitic rotat-
ing element. At a radius of about 18 in. is another plastic cylinder on which
are mounted nine parasitic elements, 40° apart. These superimpose a 135-Hz
amplitude modulation on the tram:mitted signal. Depth of modulation is about
20% for each of these signals. On the same shaft that rotates the parasitic ele-
ments are three reference-pulse disks. These generate I, 9, and 90 low-level
pulses per revolution, respectively, by varying the magnetic inductance of a
solenoid. These pulses are fed down to the transponder. The motor that rotates
this whole assembly is usually of ac type, its speed controlled to better than
1% by a servo system in which the reference-pulse frequency is compared to
a frequency standard, such as a tuning fork.
When installed aboard a ship, the Tacan antenna is stabilized in two planes.
In the horizontal plane, compensation is provided to ensure that the reference
pulses do not shift with the heading of the ship but remain oriented to north.
In the vertical plane, compensation is provided for the roll of the ship. (Early
systems also provided for pitch compensation, but this was subsequently found
to be unnecessary.)
The transponder is a constant-duty-cycle DME beacon to which the bearing-
reference pulses have been added. Once per revolution, coincident with the
maximum of the antenna pattern pointing east, a so-called north reference pulse
POINT SOURCE SYSTEMS 135
code is emitted. This comprises 24 pulses, the spacing between pulses being
alternately 12 and 18 p.sec. When these pulses are decoded in the airborne
equipment, they become 12 pulses. spaced 30 p.sec apart. This pulse train is
initiated by the one-per-revolution reference from the antenna.
Eight times per revolution, the 135Hz reference pulse group is emitted. (The
ninth group coincides with the north pulse and is intentionally omitted.) This
comprises 12 pulses spaced 12 p.sec apart. The circuitry of the transponder is
arranged in such a way that the reference pulse groups take priority over the
normal constant-duty-cycle pulses. The overall transmitted pulse envelope is
shown in Figure 4.22.
The 1350-Hz identity tone. transmitted every 30 sec, is derived from the 90
pulses-per-revolution disk on the antenna shaft. thus producing phase coherence
between identity and reference pulses and allowing each to be received without
interference from the other. The identity code comprises 1350 groups per sec,
each composed of four pulses spaced 12. 100. and 12 p.sec, respectively. The
reason for the I 00 p.sec spacing between the 12 p.sec pairs is that this combi-
nation produces the least bearing error during identity transmissions. reducing
the necessity for bearing memory circuits in the airborne equipment.
The DME interrogations are amplitude modulated by the rotating antenna.
reducing the effective sensitivity of the Tacan beacon about 3 dB below that of
an ordinary DME beacon. Although the use of a separate, nonmodulated receiv-
ing antenna would avoid this loss. such an arrangement has not been found
necessary in actual practice.
the amplitude modulation of the pulses over the required range of expected
signal strengths. This is usually taken to vary from minimum usable signal up
to about I mw of signal at the receiving antenna.
Figure 4.23 is a generic block diagram of the airborne Tacan bearing circuit.
Following decoding, the amplitude-modulated signal is filtered into two sine
waves, one at 15 and one at 135 Hz. The "north" pulse activates a 33.3 kHz
ringing circuit, whereas the 135 Hz reference pulse group activates an 83.3
kHz ringing circuit. These reference pulses are continually compared with the
two sine waves and actuate two motor-driven servo systems, geared together
9: I. Whenever the 135-Hz signal is present and the 15-Hz signal is within
±20" of its correct position, the 135-Hz signal controls the servo. In effect, the
bearing accuracy is determined by the nine-lobe antenna pattern of the ground
beacon, with the one-lobe pattern used to resolve ambiguity, which otherwise
would occur every 40". As with DME, both static and velocity memories have
been applied to airborne bearing circuits to carry them through short-term signal
dropouts. Solid-state airborne equipment typically weights 20 lb and occupies
about 1ft3 . Modern receivers incorporate digital implementations of some of
the receiver functions depicted in Figure 4.23.
i
135-Hz Phase
filter shifter f-- Comparator f---
L________ _ _
From DME
nterrogator
Peak
rider
~ I
9: 1 .gear
y
I
Servoamplifier
I
15-Hz Phase I
f-.-- Comparator 1--- I
filter shifter I
f-----'
0
I
Bearing
display
North
L-..;.. burst
decoder
l
-
VJ
-...!
Figure 4.23 Airborne Tacan bearing circuit.
138 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
4.4.8 VORTAC
Since Tacan beacons can be more readily installed on ships and at tactical sites
than VOR beacons, large numbers of military aircraft are equipped with Tacan.
To save these aircraft the cost of carrying additional equipment for navigating
the ICAO VOR/DME airways, several countries, including the United States,
use the VORT AC system. In this system each VOR station, instead of being
collocated with a DME, is collocated with a Tacan beacon (which also provides
DME service) to provide rho-theta navigation to both civil and military aircraft.
Civil aircraft read distance from the Tacan beacon and bearing from the VOR
beacon. Military aircraft read both distance and bearing from the Tacan beacon.
Thus each type of aircraft fits into the same air-traffic management system,
regardless of which type of airborne equipment it carries. In 1996, it is estimated
that there are more than 200,000 users in the United States alone.
At the ground station, the VOR central antenna is housed in a plastic cone
that supports the Tacan antenna. Leads to the Tacan antenna pass through the
middle of the VOR antenna, along its line of minimum radiation, and do not dis-
turb the VOR pattern. In the case of Doppler VOR (Section 4.4.5), the antennas
are arranged in a circle outside the cone.
4.5.1 Loran
Loran (long-range navigation) is a hyperbolic radio-navigation system that has
evolved over a period of years, beginning just before the outbreak of World
War II in Europe. The Loran-C ~·ystem [8, 11, 44] has benefited greatly from
analysis of the shortcomings of previous systems, It uses ground waves at low
frequencies, thereby securing an operating range of over I 000 mi, independent
of line of sight. Second, it uses pulse techniques to avoid sky-wave contami-
nation. Third, being a hyperbolic system, it is not subject to the site errors of
point-source systems. Fourth, it uses a form of cycle (phase) measurement to
improve precision. It inherently provides a fine-coarse readout of low inherent
ambiguity. All modern Loran systems are of the Loran-C variety. (Loran-A and
Loran-D configurations no longer exist.)
HYPERBOLIC SYSTEMS 139
Loran-C users fall into the two general categories: navigation users and pre-
cise time and time interval (PTTI) users. By far the larger population of direct
users is in the navigation category. An even larger group of indirect users ben-
efits from a PTTI application of Loran-C, in which digital switching, signaling,
and timing of the nation's telephone system is accomplished using Loran-C.
Every telephone subscriber in the United States is an indirect beneficiary of the
Loran-C system.
signal format calls for simultaneous transmission in both chains. This function
is called blanking.
Wave Form and Signals in Space Each station transmits signals that have
standard pulse leading-edge characteristics. Each pulse consists of a 100-kHz
carrier that rapidly increases in amplitude in a prescribed manner and then
decays at a rate that depends upon the particular transmitter and transmitting
antenna characteristics. The leading edge of the standard Loran-e pulse antenna
current wave form, against which the actual antenna current wave form is com-
pared, is defined as i(t):
-2(t-T)}
i(t) = A(t-T) 2 exp { sin (0.27rt+¢) for T < t < 65+7, (4.6)
65
where
Sampling point
50% amplitude
j,,. .
30 ~o~sec
because it affects the system phase accuracy. The standard sampling point for
a Loran-e receiver is the positive going zero crossing of the phase decoded
pulse on its third cycle (approximately 30 f.J-Sec) after the arrival of the ground
wave. This tracking is accomplished by a phase locked loop. In addition, it
affects the apparent signal-to-noise ratio as seen by the receiver and, therefore,
the available receiver accuracy at a given averaging time. Amplitude stability
is important, because it affects the EeD of a transmitted Loran-e pulse and
thereby affects the ability of a receiver to lock-on and track the correct cycle.
Ionosphere
Groundwave
Transmitter~~ - Ground
/
/
A /, Signal Propagation Mechanisms
Space Wave
/
,.,,....,.
---- ................
........
........
', /
/
,.,."'
/ .........
...__ ........
Earth
Figure 4.25 Ground-wave and sky-wave modes of propagation.
Over the chain coverage area, the propagation anomalies exhibit both spatial
and temporal variations.
The temporal variations fall into two primary categories: diurnal and sea-
sonal. The diurnal variations are short-term propagation effects caused primar-
ily by local weather changes and day /night transitions along the signal path.
Variations in the refractive index of the atmosphere versus height from the
ground (vertical lapse rate) contribute to the short-term propagation errors. The
diurnal time-difference (TD) variations tend to be relatively small, on the order
of tens of nanoseconds.
The larger category of temporal variations are the seasonal effects, which
are most pronounced over land paths. These long-term errors tend to be peri-
HYPERBOLIC SYSTEMS 143
where cP; is the half-angle subtended by station i and r is the correlation coeffi-
cient between two LOPs, which i~: taken to be 0.5. Note that ¢ 1 +ct>2 is equal to
the crossing angle of the LOPs (8). The relationship between 2drms and GDOP
IS
2drms = 2a-.-
Ko { ~--
1
+ + -- 2rcos 8
--- (4.8)
sm e sin 2 cP! sin 2 ¢2 sin cP! sin ¢2
where a is the timing error and K 0 is the constant 500 ft/ f.J-Sec.
500r-------------------------------------------------------~
400
300
u 200
5l0 100
c:: 0
~
~-100
~ -200
-300
-400
-500 0 100 150 200 250 300 350
50
Day
Figure 4.26 Seasonal variation of repeatable accuracy.
ity is very good. The 1994 Federal Radio Navigation Plan [48] cites a Loran-C
repeatability error range at 60-300 ft 2drms.
• TD out of tolerance
• ECD out of tolerance
• Improper phase code or GRI
• Master or secondary station output power or master station off-air
Direct Ranging There are some Loran-C users who do not employ Loran-C in
the hyperbolic mode but rather in the direct range rho-rho-rho mode or the rho-
rho mode. The rho-rho-rho process involves a minimum of three transmitters
and use of an iterative computation to obtain a fix. Direct range rho-rho mode
146 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
Antenna
Switch
To
Antenna
Network
,.....,...._
piing networks, and output networks. The TOPeO and display units contain
built-in monitoring and fault detection circuitry, and if a fault is detected the
TOPeO automatically switches coupling and output networks. The TOPeO and
display units serve as centralized alarm panels and status displays for the trans-
mitting station.
y_
RF
Front End
RF
Limiter
-- ir
I I Envelope II Limiter 1...
I
r
I
Sum J:
I I____.._
I r --Delay I,
Data
Latch
.
Data
• .
Microprocessor
Controller
I ,.. ..,I Input/Output
~
Latch
+ I Timing And
Strobe ,..
I --~
search for master and limited range search for secondaries. When the correla-
tion reaches a certain threshold, after ensuring no false locks have occurred on
secondary signals, acquisition is complete.
PGTR identification is the process of ensuring that the receiver is operating
on the ground wave of the signals. Ensuring operation of the ground wave,
sometimes called guard sampling or ground-wave location, operates on the
principle that the ground-wave signal from a station always arrives at a receiver
before the sky wave because of the longer sky-wave path. It is necessary to find
the ground wave because its timing, and hence its position-locating qualities, is
stable, while the sky wave is not. Typically, the acquisition process locates the
sky wave because of its much larger amplitude. Ground-wave location proceeds
by using signal detection algorithms at the signal 30 to 60 f-1-Sec ahead of the
receiver's reference time. If signals are found, the receiver timing is advanced
and the process repeated. This continues until no signals are found at two or
more successive locations. Often, multiple independent tests are made after no
signal is detected to account for the possibility of the ground-wave and sky-
wave signals summing out of phase and creating a null that might otherwise be
presumed to be the start of the ground wave.
Tracking is the process of maintaining a constant timing relationship between
the receiver's time reference and the PGTR for each signal being tracked. In an
automatic tracking receiver, circuits within the equipment automatically adjust
the time reference and update the display to provide continuous readings. These
receivers also provide alarms or warnings advising the operator of undesirable
signal conditions or transmitter blinking.
Loran-e signal reception can be impaired by interference from other signals
broadcast on slightly different frequencies (generally low-frequency commu-
nications). To avoid the degradation in S/N associated with these interfering
sources, Loran-e sets are equipped with notch filters that can be used to atten-
uate the interfering signal. Some receivers are equipped with preset notch filters,
others with adjustable notch filters, and yet others that automatically search for
interfering signals near the Loran-e band and dynamically notch out any inter-
ference.
152 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
M-X
AREA OF FIX
UNCERTAINTY
M-Y
M-Y
Area of Fix
Uncertainty
4.5.2 Omega
Principles and System Configuration In 1996, the Omega VLF radio-naviga-
tion system comprised eight transmitting stations located throughout the world.
At each station, continuous-wave (CW) signals are transmitted on four com-
mon frequencies and one station-unique frequency. The signal frequencies are
time-shared among the stations so that a given frequency is transmitted by only
one station at any given time.
To support medium accuracy navigation, the signal transmissions from all
stations are phase-synchronized to about l p,sec. For purposes of time transfer
and to facilitate the systemwide synchronization procedure, the signal timing is
maintained to within an accuracy of about 0.5 p,sec with respect to coordinated
universal time (UTC).
Omega signals are subionospheric; that is, they are propagated between the
Earth's surface and the D-region of the ionosphere. Because VLF signal attenu-
ation is low, the signals are propagated to great ranges, typically 5000 to 15,000
nmi. Signals with amplitudes as low as l 0 p, V/meter can often be detected and
used for navigation. Of primary interest to navigation users is the signal phase
which provides a measure of transmitter-receiver distance. The fractional part
of a cycle (or lane, which is the equivalent distance measure) is generally the
156 TERRESTRIAL RADIO- NAVIGATION SYSTEMS
only measurable component of· the signal phase. thus leading to lane ambiguity.
However, the lane ambiguity problem is reduced through the use of multipl e
frequencies and is resolved for navigation through a process of continuous lane
count.
When used as a stand-alone system for navigation, an Omega receiver pro-
vides an accuracy of 2 to 4 nmi 95 % of the time 143, 48J. In the differential
mode of operation, where a receiver utilizes Omega signal phase corrections
transmilled from a nearby monitor station, a position accuracy of about 500
meters can the attained. Because Omega is a continuous VLF phase-measur-
ing system, it has been appropriately integrated with noncontinuous , high-accu-
racy sensors. The resulting system has an accuracy that is comparable to the
high-accuracy navigation aid and degrades relatively slowly in time when the
high-accuracy aid is unavailable. As commonly used in overocean civil air-
line configurations, an Omega receiver is combined with an inertial navigation
system, so that the Omega system error effectively "bounds'' the error of the
inertial system.
The signals from the eight Omega transmitting stations shown in Figure 4.34
provide continuous signal coverage over most of the globe. The suite of elec-
tronics equipment (mainly signal generation, control , and amplification units) is
virtually the same for all stations in the system , but the station antennas differ
substantially. Because they radiate long-wavelength VLF signals. the antennas
are the largest physical structures at the stations. Three types of antennas are
employed in the Omega system: (I) grounded tower, (2) insulated tower, and
(3) valley-span. Each has an associated signal monitoring facility about 20 to
50 km from the effective phase center of the antenna. These unmanned facilities
perform several functions, including monitoring the performance of the asso-
II
;
i
, :1'' ··.- , __;
·; ',·
' tJ']~pan
i,Hawaii
•
i .• .... .
Liberia i
.
·.:11 /]' ~
i ; • La Reunion 1·
Australia
() ·~· ;:_,:
' '
{!Argentina
L __ _ __ _ ___;;____;_, ____ __ __ _ _
Segments
Station
Norway A
ubena B
Hawaii c
Nonh Dakota D
La Reun1on E
Argentina
Australia G
Japan H
Wave Form and Signal in Space The principles of Omega navigation usage
depend almost entirely on the assumed relationship between the signal phase
HYPERBOLIC SYSTEMS 159
In free space, the wave number is given by 1 k 0 =f/c, wheref is the fre-
quency of the signal and c is the speed of light. The nominal wave number is
given by knom = 0.9974ko, which is chosen as an intermediate value between
observed night and day wave numbers on seawater paths. The exact value of
the nominal wave number is not critical; it is only important that the value be
near the average over all time and space conditions.
The Omega PPC may be thought of as the variation of the "true" Omega
signal phase (¢) from the nominal phase (¢nom):
where r:f>nom = knomr, and where r is the distance between a transmitter and a
receiver over the surface of the Earth. The assumption of a fixed wave number
that relates nominal phase and distance over the surface of the Earth is the basis
for the so-called nominal model of Omega signal phase/distance relationships.
The calculation of the PPC is based on a semiempirical model of phase vari-
ation as a function of the electromagnetic characteristics of a signal path from
transmitter to receiver [22-23, 24-25].
Thus, by Equation (4.9), the PPC and the nominal model together determine
the predicted phase for a given station, signal frequency, position, and time.
In this relation the nominal phase is the "dominant" term in the sense that it
accounts for approximately 99% of the cumulative phase from the signal source,
that is, the distance between the transmitting station and the receiver in units
of wavelength. Measured in cycles of nominal signal wavelength (somewhat
larger than a free-space wavelength) at I 0.2 kHz, the nominal phase is 100 to
500 for typical paths, whereas the PPC is usually between -3.00 and +3.00
cycles, with a resolution of 0.01 cycle (a unit referred to as a centicycle). The
predicted phase has a typical diurnal variation of 0.5 to 2 cycles, amounting to
about 0.2 to 2% of the nominal phase.
Figure 4.36 shows a typical diurnal observed phase profile (measured with
respect to a precise time or frequency standard) in which the path illumina-
tion conditions, nominal phase, and two sample PPC values are identified. The
figure illustrates the higher (retarded) phase during path night and the lower
(advanced) phase in path day, with a total diurnal shift of about 0.65 cycle.
Since the phase is a function of effective ionospheric height which varies with
the relative sun angle (solar zenith angle), the observed phase exhibits a "bowl-
shaped" profile during the clay with less variation at night. The phase profile
changes from day to night behavior during path transition when the sunset ter-
minator cuts the path. The figure illustrates the time-independence of the nom-
inal phase and the consequent time-dependence of the PPC values.
1An alternative definition of the free-space wave number used in many texts is ko = 27r.f /c.
HYPERBOLIC SYSTEMS 161
120
iii'
Ill
u>- eo
~
Nominal Phase
'EIll 60
.!:!..
Ill
en 40
~
.c:
~
20
Transition -+--- Night ---+-- Transition _...,.,___ _ Day ---~
0
0 02 04 06 08 10 12 14 16 18 20 22 24
UT
Figure 4.36 Typical diurnal phase behavior.
I. The lower boundary (the Earth"s surface) has a relatively high conduc-
tivity (greater than 10- 3 mhojm over most areas of the Earth). so waves
do not readily penetrate the surface.
2. The Earth's atmosphere (at altitudes between 0 and 70 km) has an
extremely low concentration of charged particles and thus acts as a vac-
uum to VLF waves.
3. The D and E regions of the ionosphere (70 to II 0 km) have low average
conductivity (about I o-s mhojm) but have a steep conductivity gradient
between 70 and I 00 km which serves to reflect VLF waves.
The above conditions also lead to low attenuation of the propagated signal; for
example, over a range of about 1000 km the signal amplitude is reduced (on
average) by a factor of two.
Factors effecting Omega signal propagation include the action of the Earth's
magnetic field, the structure of the ionosphere, solar control and the effects of
the !!-year sunspot cycle, and the presence of two or more propagation modes.
The Earth's magnetic field introduces an anisotropy into the behavior of
VLF waves interacting with the ionosphere. That is to say, signal propagation
depends upon the direction of propagation. This anisotropy is strongest on paths
perpendicular to the geomagnetic field (east-west paths). The presence of two
or more signal propagation modes with comparable amplitude will cause the
phase to become a strongly oscillatory function of distance, thus rendering the
signals unusable for navigation/positioning.
162 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
The ionosphere is quite sensitive to the net incident solar illumination. Dur-
ing the day, solar photoionization maintains a small, but stable ionized compo-
nent which is not present in nighttime regions. Solar control of the ionosphere
introduces a strong diurnal dependence on Omega signal propagation.
Station A Station B
(before moving distance Llr) is zero. However, after moving distance Llr toward
the station B along the baseline, the Jane phase difference is given by
Since the wave number is given by k = ljf.., where f.. is the signal wave-
length, it follows from the above that Llr = l/(2k) = f../2. At points away from
the baseline, the Jane width increases with the diverging hyperbolic curves as
f../(2 sin (1/; As/2)), where 1/; AB is the angle subtended by the two stations at the
receiver location.
A hyperbolic Jane is actually the family of all hyperbolic curves with phase
differences between 0 and I cycle that lie within the Jane boundaries. A hyper-
bolic curve normally has two branches, corresponding to the positive and nega-
tive values of the range difference. In the case of Omega, the sign of the phase
difference is known from measurement, which limits the receiver location to
a single branch. In the vicinity of the receiver (near the baseline), the hyper-
bolic curve resembles a planar hyperbola with foci located at the two associated
transmitting stations. On scales of 1 to I 0 nmi, these hyperbolic curves are well-
approximated by straight lines. At points well away from the baseline. however,
the spherical shape of the Earth causes the hyperbolic curve to close on itself in
a quasi-elliptical shape, with one of the stations at one focus and the antipode
of the other station at the other focus.
If the appropriate Jane is shown for a second pair of stations (which may
include a station common to the first pair), then a phase-difference measurement
with respect to these two stations establishes a second hyperbolic curve, whose
intersection with the first curve determines the receiver position (see Figure
4.38). It is possible that these two curves could intersect in two locations, but the
correct intersection is easily resolved for one or more of the following reasons:
I. One of the two intersections is usually relatively far from the known
approximate location of the receiver.
2. An independent third pair of stations (if available) provides another hyper-
bolic curve that passes near one of the two intersections.
3. A moving receiver shows successive fixes consistent with vehicle speed
for one of the intersections and inconsistent for the other.
In cases where more than two independent hyperbolic curves are available
(i.e., more than four usable Omega station signals are accessible), the multiple
curves do not, in general, intersect at a single point due to the effects of noise
164 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
Station
''' Station
A ', B
---$------------
Baseline
------~---
''
'
Rho-Rho Method The rho-rho technique requires only two range measure-
ments for a fix. As in the hyperbolic case, it is assumed that the correct lane is
initially known and successive measurements are processed over small enough
distance/time intervals so that lane changes are readily tracked. To obtain an
accurate estimate of the distance traveled based on successive phase measure-
ments, the processor must have access to a frequency /time reference (clock)
of sufficient stability so that the reference is effectively synchronized to the
Omega station during the period between precision updates.
The change in station-receiver distance, obtained from two successive phase
measurements of the station signal, places the receiver's new position on a cir-
HYPERBOLIC SYSTEMS 165
cular curve (centered on the station location) within the appropriate lane. Since
the receiver's previous position is assumed known, some points on the new
circle are more likely candidates than others for the new position, based on
platform velocity and maneuvering limits. However, the new position is accu-
rately determined only when a distance change to a second station is obtained
from successive phase measurements. In this case, a second circle is established
that intersects the first at the new receiver location. Although two intersections
are possible, the correct intersection can be resolved using methods listed above
for the hyperbolic case. This method corresponds to a system of two equations
and two unknowns. In cases in which more than two usable station signals are
available, a least-squares technique can be used to resolve the multiple inter-
sections that arise as a result of measurement noise or phase prediction error.
nearly constant but often differs significantly from the average observed phase.
This difference is referred to as the PPC bias error which varies in magnitude
from 0 to 30 centicycles (cecs). Also occurring in this 15-30 day period at a
fixed hour are random (nonsysternatic) day-to-day variations in the observed
phase on the order of I to 5 cec;. Since these random variations (which are
due to ionospheric fluctuations) are not reflected in the PPCs, they make up the
random component of PPC error.
When converting phase measurements to position, the bias and random com-
ponents of phase error produce corresponding bias and random components
of position error. The transformation of the phase error components to posi-
tion error components depends upon the individual phase errors of all signals
received and the geometrical configuration of the receiver and stations corre-
sponding to the received signals. If the magnitude of the random phase errors
is assumed to be the same for all ~ ignals received and the bias error is assumed
to be zero, then the radial position error standard deviation (a r) can be obtained
by multiplying the phase error standard deviation (aq,) by a scalar factor known
as the geometric dilution of precision (GDOP). For a least-squares method of
position determination, used when multiple redundant signals are present, the
following form 2 of GDOP [24, App. B] is obtained:
1
GDOP=
2
q I 'I 1/2
L ~~ 2
sin (((3;- (3i )/2)
q-2 q-1 q
L L L sin 2
(((3k- ~~j )/2) sin 2 (((3;- (3k)j2) sin 2 (((3j- (3;)/2)
i= I j=i+ I k=j+ I
(4.11)
where q is the number of usable signals received and (3; is the bearing to the
ith station (corresponding to the ith usable signal). The GDOP becomes very
large whenever at least q - I stations have bearings which are nearly equal.
Another property of the GDOP is that the GDOP for q station signals is never
greater than the GDOP for any subset (> 3 stations) of q. This means that for
least-squares position processing, the use of additional (usable) signals does not
degrade, and typically improves, the resulting position accuracy.
For moving vehicles performing navigation, the bias error is effectively
removed at initialization, leaving only phase error due to noise (typically less
than 1 cec ). However, the paths from the station to the receiver eventually
change (both in space and time) enough so as to become decorrelated with
the original configuration of station signal paths to the receiver and the initial
correction no longer applies. From this decorrelation time until the next preci-
sion update, the Omega receiver is subject to PPC bias and random errors and
the effect of GDOP. Omega-only accuracies have been reported for aircraft of
2.7 to 3.3 nmi 95% of the time l32, 33].
that signal phase from different stations cannot be compared (in an absolute
sense) to determine position. Because the received VLF signal is generally sta-
ble in time, VLF navigation requires an initial calibration to permit subsequent
phase tracking of the signals from selected VLF stations. Accurate phase track-
ing requires an on-board precise frequency standard or a correction based on
an estimate of the frequency /time offset of the receiver's internal clock. This
estimate is usually obtained from Omega signal processing in the rho-rho-rho
mode.
In addition to internal differences in signal processing, signals from the two
systems are processed differently regarding external information. For example,
all known OmegajVLF receivers use externally supplied PPCs to correct the
measured Omega phase prior to navigation use, whereas few, if any, currently
operational receivers correct VLF signal phase measurements. This means that,
for most receivers, the received VLF signal phase is not accurately related to
distance over the Earth's surface, a problem that is not necessarily amelio-
rated by redundant measurements. External deselection data regarding modal
and long-path signals are available for Omega but not for VLF. Failure to de-
select VLF modal signals is potentially a more serious problem for navigation
than the lack of VLF PPCs, since modal phase excursions can be large and
sudden, often resulting in cycle slips or advances.
As a result of the signal-processing differences, due to the internal and exter-
nal information bases, receiver-processing algorithms treat Omega and VLF sig-
nals differently. Once acquired and initialized, Omega signal-processing alone is
robust and will fail only under unusual circumstances (e.g., cycle shifts or fewer
than three signals above the minimum SjN). VLF signal-processing schemes
generally rely on the presence of Omega signals and other aids in the receiver's
navigation filter. In most receivers VLF signals are closely monitored with
frequent cross-consistency checks. Normally, OmegajVLF receivers are pro-
grammed to exclude initialization with VLF signals alone, since this repre-
sents a "degraded" mode. Current FAA certification procedures require than an
OmegajVLF receiver system operate satisfactorily with Omega signals alone.
The principal functions of the timing and control subsystem are signal gener-
ation and phase control. The signal source is a precision cesium beam frequency
standard of 9.193 GHz with a stability of 5 parts in 10 12 . Three cesium stan-
dards are used for frequency drift comparison and control, and are maintained
HYPERBOLIC SYSTEMS 169
as reserves in the event of failure of the on-line standard. Phase control is main-
tained by comparing the RF signal phase to the phase of the antenna current
reference signal fed back from the antenna tuning subsystem. The signal phase
is advanced or retarded to insure that its phasing at the antenna coincides with
the appropriate UTC epoch.
The transmitter subsystem consists of those devices that amplify the signal
generated in the timing and control subsystem. The RF signal from the timing
and control subsystem is first raised to a level of 160-V RMS by the input
amplifier. The driver amplifier further raises the signal level to a nominal 520-
v RMS and the final amplification is performed by the power amplifier that
boosts the signal voltage and current to a peak power of 150 kw. Following
this final amplification stage, the signal is fed to the antenna tuning subsystem.
The antenna-tuning subsystem is designed to tune the antenna at the RF sig-
nal frequency by impedance matching the antenna to the input circuit. This
ensures the maximum effective radiated power at the antenna for a given input
signal power. Based on the long keying pulses from the timing and control
subsystem and the current samples received from the current transformer, the
antenna-tuning control first implements fine inductive tuning through the var-
iometers. The antenna-tuning control signal activates a mechanical drive that
moves the variometer coil to the appropriate position for matching impedance.
The long keying pulses activate antenna relays that connect the appropriate var-
iometer into the main antenna circuit. The RF signal is then transferred to the
"helix," a large helical coil that acts as a coarse tuning device for the antenna.
The helix is equipped with separate taps for each signal frequency transmitted.
Finally, the RF signal is conducted to the antenna structure itself from which the
signal is radiated. The structural feature which principally differentiates Omega
stations is the antenna structure. Two basic designs are utilized: tower and val-
ley span. The tower antennas are further classified as either the grounded or
insulated type.
the input signal in units of clock cycles. 3 A digital phase tracking technique
used in many Omega receivers is the phase lock loop in which the reference
phase is shifted by an amount that depends on the previous phase measure-
ment and the time-averaged phase computed at the previous measurement time.
The time-averaging refers to a moving average that differs from the average at
the previous loop cycle by a weighted value of the previous measurement. A
second-order phase lock loop is designed to track the time rate of change of
phase in a manner similar to that of the first-order loop.
Like most signal-tracking circuitry. the basic function of the phase lock loop
is to reduce the effective bandwidth (inversely proportional to the effective time
constant) so as to best reproduce the desired signal. For aircraft receivers. time
constants typically range from I 00 to 200 sec. Shorter time constants do not
provide sufficient averaging or noise rejection. and longer time constants may
exceed the time required for aircraft maneuvers. such as sharp turns. Since the
duty cycle for each of the common frequency Omega signals is I 0%. the e,ffec-
tive phase measurement time comtant is I 0 to 20 sec. Using standard assump-
tions [24]. these time constants correspond to noise equivalent bandwidths of
0.025 to 0.013 Hz. When compared to the input bandwidth of I 00 Hz. these
narrow output bandwidths correspond to gains of better than 35 dB. Thus. sig-
nals with SjN as low as -20 to -30 dB in the 100Hz receiver input bandwidth
can be effectively utilized in aircraft Omega receivers.
After the signal phase measurements are made. PPCs are computed using an
appropriate modeljalgorithm and added to the measured phase to produce an
"idealized" phase value that can be readily used in the subsequent positioning
calculations. Although the PPCs require receiver position as an input. the PPCs
are not sensitive to precise position since they vary less than 0.05 cycle over
ranges of 50 to I 00 km. Thus, the PPCs can be accurately computed from only
approximate knowledge of position.
Before determining position, the (idealized) signal phases are usually
weighted based on the expected relative accuracy of the phase measurement.
This accuracy is most commonly determined by the estimated SjN. which. for
phase lock loop receivers. is clo~:e!y related to the rms loop error. If the esti-
mated Sj N is below a preset threshold. such as -30 dB in a I 00-Hz bandwidth.
the signal phase is usually excluded (given zero weight) in the position solu-
tion. In addition to these weighting and exclusion procedures based on internally
derived SjN data. the signal phases are edited by invoking external information
concerning the signals. External information usually refers to signal deselection
data that are generally extracted from known coverage information. including
modal "maps'' and data on the occurrence of long-path signals. The resulting
signals that arc not deselected or excluded are further screened for acceptable
geometry. In some receiver mechanizations. all common frequency signals from
a station must be acceptable to be used in the position fix; in others. only a sin-
0The reference clock/oscillator commonly has a frequency of I to 5 MHz but may be converted
to a lower frequency.
HYPERBOLIC SYSTEMS 171
gle acceptable signal frequency from a station is necessary for inclusion in the
fix algorithm. Position change estimates are then formed from the weighted
and edited Omega signal phase data at the common Omega signal frequencies.
The estimates are computed by means of a least-squares or Kalman estimation
technique (see Chapter 3). The Omega-based calculation of position change
is frequently combined with the aircraft-supplied true airspeed and heading or
inertial system information to furnish the best position estimate.
An airborne Omega receiver block diagram is shown in Figure 4.39; a pho-
tograph of an Omega/VLF receiver for commercial aviation applications is
depicted in Figure 4.40.
4.5.3 Decca
Hyperbolic systems other than Loran and Omega exist and are used for nav-
igation. One such example is the Decca system [31] developed by the British
and used extensively during the later stages of World War II. In 1996, its major
area of implementation is in northwestern Europe where it is primarily used by
shipping companies.
Decca is based on the measurement of differential arrival times (at the vehic-
ular receiver) of transmissions from two or more synchronized stations (typi-
cally 70 mi apart). As an illustration. consider two stations (A and B) I 0 mi apart
and each radiating synchronized radio-frequency carriers of I 00 kHz. Assume
that there is some way by which each station can be identified. The wavelength
at this frequency is 3000 meters, or about 2 mi. On a line between the stations.
the movement of a vehicle D one mile toward one station and one mile away
from the other station will cause the vehicle to traverse one cycle of differen-
tial radio-frequency phase. There wilL therefore, be I 0 places along the line AB
where the signals from the two stations will be in phase. As the vehicle moves
laterally away from this line, isophase LOPs can be formed (each line being a
hyperbola) with the stations as foci and BD- AD as a constant for each LOP.
Site error virtually vanishes in such a system. and the accuracy depends
entirely on the constancy of propagation between the stations and the vehicle.
In an effort to avoid line-of-sight limitations, Decca uses a low frequency (70
to 130 kHz), which is subject to sky-wave contamination. and uses continuous
waves. which preclude the separation of ground waves from sky waves. Thus,
despite the low frequency (whose ground-wave range is on the order of I000
mi), practical Decca coverage is limited to areas where sky-wave strength does
not exceed about 50% of ground-wave strength. This is typically 200 mi.
A typical Decca chain consists of a master station and three slave stations.
A typical station has a 2-kw crystal-controlled transmitter feeding a 300-ft
antenna. The slave stations are referred to by the color of the phase meter asso-
ciated with each at the receiver. Each station transmits a stable continuous wave
frequency that bears a fixed relationship to the frequencies of the other three
stations. Phase comparison therefore produces a family of hyperbolic LOPs of
.......
-...l
N
Atmosphenc
Noise
Internal Noise
(POSitiOn. limel
CDU
Signal Form Instruments
Cond1t10n1ng Phase
POSitiOn Autopilot
(Ga1n L1m1ting, Track
Loop Estimate
Filtenng) Miss1on Com-
SwitChing
puter
Reference
Phase
L_____ w~~~-~-----J
Editing/ 1 I
constant phase. The spaces between these lines are called lanes. The intersect ion
of two LOPs provide a position fix.
4.5.4 Chayka
ehayka (meaning "sea gull") is a pulse-phase radio-navigation system similar
to the Loran-e system. It is used in Russia and surrounding territories and seas.
By using ground waves at low frequencies, the operating range is over 1000 mi;
by using pulse techniques, sky-wave contamination can be avoided. The system
is designed to provide both a means of determining an accurate use r positi on
and a source of hig h-accuracy time sig nals. The system can support a n unlim-
ited number of users since the computations are pe rformed at the user receiver
and position determination is possible at any time of day or year, regardl ess of
meteorological conditions.
Each of the stations in the Chayka networks transmit pulses with standard
characteristics. The pulse consists of a I 00-kHz carrier wave that increases from
zero to a maximum and then decays at a specific rate to form the envelope of
the signal. All slave stations transmit signals in packets of eight pulses; the
masters emit a ninth pulse for identification. The interval between pulse onsets
is 1.0 ± 0.05 J-tSec. In addition, in order to provide the possibility for automatic
detection and identification of signals and to reduce the influence of multiply
reflected signals, the signals are phase-coded with the slaves all having the same
phase (i.e., oo or 180°) and the master phase differing by exactly 180°.
The repetition periods of the radio signals are selected based on a trade-off:
maximizing the average power of the signals at the receivers while preventing
any signal overlap within the network operating region. Since all slave stations
transmit signals with identical phase codes, each slave station transmits with
its own specific code delay relative to the master signal. The magnitude of the
code delays are selected such that the order of reception of slave station signals
is identical everywhere within the network operating region.
Wave Form and Signals in Space Each station transmits signals with stan-
dard pulse modulation charactcris1:ics. Each pulse consists of a 100-kHz carrier
wave modulated by an envelope that depends on the specific transmitting sta-
tion equipment. Two types of radio transmitter (RT) stations are currently in
use: those with vacuum tubes and those using impact excitation of the linear
output circuit. The envelope of the vacuum tube RTs can be approximated by
where U 111 is the pulse amplitude and t 111 is the time interval from the onset to
the peak of the pulse.
The envelope of the impact excitation RTs can be approximated by
Coverage In 1996 the Chayka system consisted of only two networks, Euro-
pean and Eastern. Coverage from the European network is centered near
Moscow and includes most of the area between 5c and 50' East longitude and
FUTURE TRENDS 175
40° and 65° North latitude (e.g .. Eastern Europe. western portion of the former
Soviet Union. the Black Sea, and part of the Caspian Sea). Coverage from the
Eastern network includes most of the area between 135° and 160° East lon-
gitude, and 35° and 65° North latitude (e.g .. eastern shoreline of the former
USSR. portions of Japan. and the surrounding areas of the Pacific Ocean).
Terrestrial radio-navigation systems will continue to play a major role for air-
craft navigation throughout the world for many years. Since the U.S. satellite-
based GPS had achieved full operational capability (FOC) in 1995. followed
in 1996 by the Russian GLONASS. there had been expectations that these
satellite systems would quickly replace the terrestrial systems such as VOR,
VOR/DME. Loran, and Omega. However, this was not the case and is not
likely to occur in the near future. The reasons for this include (I) the widespread
implementation of equipment by aircraft owners and the cost of replacement by
satellite receivers. (2) the lack of available air-traffic management operational
procedures compatible with satellite-based systems. (3) the absence of full sole-
means navigation system status of GPS. and (4) the fact that issues involving
system accuracy, integrity. availability. and continuity of service of the satel-
lite systems had not been fully resolved (Chapters 5. 13, and 14 ). Therefore.
the terrestrial radio-navigation systems will continue to be used for many years
on a global basis. In the more distant future, some of these systems will be
decommissioned when their utility will have been fully replaced by that of the
satellite systems.
By 1993, the U.S. Coast Guard had implemented full coverage of the conti-
nental United States by Loran-C chains and the FAA had authorized Loran-C
for supplemental navigation for en-route and nonprecision approaches. At least
ten U.S. airports had received approval for Loran-C approaches. As a result.
there was extensive use of airborne Loran-C receivers on U.S. General Aviation
aircraft and that usage is likely to continue for some time until GPS receivers
are widely implemented on General Aviation aircraft.
Since 1990. a number of major studies have been conducted and published
that show the advantages and discuss techniques of combining data from Loran
and GPS for aircraft navigation [45. 46, 47]. Among the major advantages are
the mitigation of the effects of GPS coverage outages caused by satellite shut-
downs or poor geometry and. conversely. that of Loran coverage outages due to
ground station shutdowns. high atmospheric noise levels. or precipitation static.
In addition Loran data could provide on-board fault detection and isolation of
GPS satellites, in connection with GPS Receiver Autonomous Integrity Mon-
itoring (RAIM. Section 5.7.2). The combining of GPS and Loran data (e.g.,
with a Kalman filter) can be at the pseudorange level and mutual time synchro-
nization can also be included [46]. Therefore. research and development on the
integration of Loran and GPS is likely to continue in the future.
176 TERRESTRIAL RADIO-NAVIGATION SYSTEMS
PROBLEMS
4.7. What are the factors that impact the accuracy of a Loran-e fix?
Ans.: Signal-to-noise ratio at the receiver, crossing angle of the
Loran-C lines qf position, calibration r~f the Loran-C time difference
to latitude/longitude coordinate converter.
4.8. What are the two categories of Loran-e system accuracy? What do they
mean?
Ans.: Repeatable accuracy to which one can return to a point vis-
ited before; absolute or predictable accuracy of the fix against some
external reference grid such as latitude and longitude.
4.9. Suppose that an Omega receiver processes four I 0.2 kHz signals from
stations with geographic bearing angles (at the receiver) of 31 o, 121 o, 211 c,
and 30 I o using a least-squares algorithm to estimate position change.
(a) What is the GDOP?
Ans.: 3/2Vl.
(b) If the signal phase error associated with each of the four stations is
4 cecs, what is the corresponding position error in kilometers.
Ans.: a,.= 1.25 Km.
(c) If one of the station signals becomes unusable (e.g., due to modal inter-
ference as the path becomes dark), by what factor is the position accu-
racy degraded?
Ans.: a,. becomes larger by a factor of 4/3.
5 Satellite Radio Navigation
5.1 INTRODUCTION
Since the 1960s, the use of satellites was established as an important means of
navigation on Earth. The earliest :;ystems were designed primarily for position
updates of ships, but were also found useful for the navigation of land vehicles.
Beginning in the early 1970s, satellite-navigation systems for aircraft (as well
as other platforms) were under intense development. Those efforts benefited
from the techniques used and the experience gained with the earlier systems.
In the 1980s, systems suitable for aircraft became mature and by 1996 their use
for aircraft navigation was increasing at a widespread and rapid pace.
The satellite-navigation systems described in this chapter are comprised of a
system of satellites that transmit radio signals. Appropriately equipped aircraft
receiving these transmitted signals can derive their three-dimensional position
and velocity and time. Two systems are described in detail, namely the U.S.
Department of Defense's NAVSTAR Global Positioning System (GPS) and the
Russian Federation's Global Orbiting Navigation Satellite System (GLONASS).
The International Civil Aviation Organization (ICAO) and RTCA, Inc. have
defined a more global system that includes these two systems, geostationary
overlay satellites, along with any future satellite navigation systems, in what
has been named the Global Navigation Satellite System (GNSS) [1, 2]. A third
major system, the United States Navy's Transit System, also called the Navy
Navigation Satellite System (NNSS), is a low-altitude Doppler satellite radio
navigation system. In Russia, a similar system was developed, called Tsikada.
Since GPS was fully operational, after 32 years the U.S. Navy ceased operations
of Transit on December 31, 1996 [ 124]. It will not be discussed here further.
(Design details are given in references [3] and [4].)
The systems described in this chapter provide users with a passive means
of navigation; that is, there is no requirement for their equipment to transmit,
only to receive. Both GPS and GLONASS are ranging systems. They provide
both range and range rate (or change in range) measurements. Once initialized,
they provide an instantaneous and continuous navigation solution in a dynamic
environment. Details of these solutions are described later in this chapter.
The advantage of satellite navigation systems is that they provide an accurate
all-weather worldwide navigation capability. The major disadvantages are that
they can be vulnerable to intentional or unintentional interference and tempo-
rary unavailability due to signal masking or lack of visibility coverage. In some
I. Monitor stations that track the satellites' transmitted signals and collect
measurements similar to those that the users collect for their na vigation.
2. A master control station that uses these measurements to determine and
predict the satellites' ephemeris and time history and subsequently to
upload parameters that the satellites modulate on the transmitted signals.
'I
/I
~rc-7.~/
,- - )f 1 .,.r WF~:~{:~:'"~: ,,;'r?~:,~:, '~:tLCo
~lj\;,{\ GROUND
1j! 1 ~ , , ANTENNA
L' ' ~~ USER SEGMENT
MASTER CONTROL
STATION 1-4-___J
CONTROL
SEGMENT
3. Ground antennas that perform the upload and general control of the satel-
lites.
Z. SATELUTE2
X
Figure 5.2 Ranging s;atellite radio-navigation solution.
THE BASICS OF SATELLITE RADIO NAVIGATION 181
filter algorithm. The measurements are not range and range rate (or change
in range), but quantities described as pseudorange and pseudorange rate (or
change in pseudorange). This is because they consist of errors, dominated by
timing errors, that are part of the solution. For example, if only ranging type
measurements are made, the actual measurement is of the form
(5.1)
(5.2)
where X, 1, Y,;, and Z,; are the Earth-centered, Earth-fixed (ECEF) position com-
ponents of the satellite at the time of transmission and Xu. Yu. and Zu are the
ECEF user position components at that time. For the three satellites in Fig-
ure 5.2, Equation 5.2 represents the equations for spheres whose centers are
located at the satellites. The user position is the reasonable intersection of the
three spheres. (There is another solution, but not near the Earth.)
oR;= R,m- R,
182 SATELLITE RADIO NAVIGATION
(5.3)
where oR; is ~he range measurement residual, R;111 is the range measurement
to satellite i, R; is the estimated range to satellite i, lx;, I yi. and I zi are the
components (directional cosines) of the estimated line-of-site (LOS) unit vector
1; between the user and satellite i and oX11 , oY,, and oZ11 are the components of
the vectm;.oX11 of differences between the position solution Xu and the estimated
position X 11 • This vector represents an offset from the intersection of the three
spheres. Figure 5.3 illustrates this linearization in two dimensions, where the
inner circles represent the measured ranges and the outer circles represent the
"computed" ranges based on the estimated position. The shaded areas represent
the range measurement residuals.
Solving three equations representing range measurement residuals from three
satellites give a solution for the position correction vector, provided that the
geometry is sufficient (i.e., the solution exists). If the differences are large, as in
the exaggerated example of Figure 5.3, so that they exceed the range of the lin-
earization, an iterative solution is generally required. This can be accomplished
by using either the same set of measurements or subsequent measurements,
where the user position in the computation of new LOS vectors is propagated
from the previous solution.
(5.4)
where
V; =[X; Y; Z;f is the known satellite i velocity vector
Vu =[Xu Yu Zuf is the unknown user velocity vector
X;= [X; Y; Z;f is the known satellite i position vector
Xu= [Xu Yu Zu]T is the unknown user's position vector
Note that Equation 5.4 is also nonlinear because the LOS vector is also a func-
tion of the user's unknown position. However, linearizing about an estimate of
position and velocity yields
_:__
The second term of the equation can be neglect~d under normal circumstances
because the rate of change of the LOS vector I; is small, which is due to the
large distance to the satellite.
In some precision landing applications, the measurements may be Doppler
count measurements. Then, change in range can be computed by differencing
Equation 5.3 over a time interval [tj _ 1, tj ], resulting in the equation for a change
in range as
184 SATELLITE RADIO NAVIGATION
where the subscript j indicates measurements taken at time ti. Note that over
short time intervals, this change in range measurement can be used to estimate
range rate by dividing by ti - fj _ I.
/::,.fu
-c(/::,.tuj - /::,.tu,J- I) = -c - - (tj - tj _ I) (5.8)
fo
z1
/Earth 's polar axis
'
v
Satellite
x3
i inclination
y
uinol'.
'Jerna\ eC\_ __
X~ ~--- / Yj
n ~~- -------------~
'..Ascending node Equatorial plane
X/
Figure 5.4 T he orhital plan e.
186 SATELLITE RADIO NAVIGATION
APO~EEI
~
·--~~~~----~--~
LINE OF ABSIDES
ASCENDING
NODE
LINE OF NODES
2. Inclination of the orbital plane with the equatorial plane, i (Figure 5.4).
where b is the semiminor axis. Note that for a circular orbit, where the two
axis are equal, the eccentricity is zero. The sixth orbital element is the time
of perigee passage, tp, measured with respect to some arbitrary time scale. tp
establishes the phase of the satellite along the geometric path defined by the
other elements.
Useful Orbital Parameters and Equations The six orbital elements describe
the path of the satellite in its unperturbed orbit. To perform a navigation solu-
tion, however, it is usually better to describe the satellite's position and velocity
in ECEF coordinates. The relationship between the six coordinates and satellite
position is as follows [6]:
The mean motion is the average angular rate of the satellite radius vector r.
It is defined as
(5.1 0)
(5.11)
in radians, which is the basis for the solution of Kepler's equation [6]. Kepler's
equation defines the relationship for the eccentric anomaly, where
E- esin E = M (5.13)
a(l - e 2 )
r =a( I- ecos E ) = - - - - (5.14)
1 + e cos v
where v is the true anomaly in radians. Note from Figure 5.5 that to compute
the cartesian components of the satellite's position (Xp, Yp) in the orbit plane,
the cosine and sine of the true anomaly is required. They can be computed as
~sinE
sin v =' (5.15)
I - e cos E
cos E- e
cos v = (5.16)
I - ecos E
This position in the orbit plane can then be transformed into ECEF coordinates
by performing a Euler transformation through the orientation parameters w, i,
and 0, in that order. Note, however, that 0 is not constant because the Earth
is rotating. It varies from some predefined epoch value 0 0 (right ascension) at
some time to as
(5.17)
in radians, where ~h is the earth's rotational rate, which has a WGS 84 value
of7.2921151467 x 10- 5 radjsec [7, 8].
UV, A,¢) = ~ [I + t, ( A: )" ~ (C,,, cos rnA + S,, ,, sin mA)P,,, (sin ¢) l
(5.18)
where
AE is the WGS-84 semimajor axis of the Earth's ellipsoid =
6378.137km [9]
n,m are degree and order
¢,"A are geocentric latitude and longitude
Pn.mCsin ¢) are Legendre polynomials
Cn,m,Sn,m are geopotential coefficients
Neglecting the effects of longitude (m = 0), which are relatively small compared
to the effects of latitude, the gravity potential due to the second zonal harmonic
is [6]
U2(r,¢)= J-tAioC2 ·
o ( -3s i n
2 I)
¢-- (5.19)
1 r 2 2
3 A 2 C2
Jl E 'O (3 COS 2¢- 1) (5.20)
4r 4
3A~C2.osin 2¢
(5.21)
2r4
Note that the radial force has two components-a constant that adds to the
nominal gravitational force and one that oscillates as function of the satellite's
latitude ¢, where
based on the argument of latitude, the argument of perigee, and the inclination
angle of the satellites orbit. Note that the period of the oscillating force is one-
half the orbit's period. The force in the direction of latitude also oscillates with
latitude.
190 SATELLITE RADIO NAVIGATION
(5.23)
where
is the fractional frequency offset averaged over 7 seconds after the systematic
frequency offset and frequency drift have been removed. i:J.¢(tk) is the change
in clock phase, measured in radians over 7 seconds as illustrated in Figure 5.6.
f 0 is the nominal frequency of the oscillator in Hertz and M is the number of
samples used in the computation.
A typical square root of Allan variance for a good quality crystal oscillator is
shown in Figure 5.7. Three typical stability characteristics are shown, depend-
ing upon averaging time: white frequency noise, flicker frequency noise, and
random walk frequency noise. They are defined using the coefficients hex of
ORBITAL MECHANICS AND CLOCK CHARACTERISTICS 191
Cl)
zc:(
c
~
-
':&:
w
Cl)
c:(
:I:
0..
TIME t- SECONDS
FUCKERFREQUENCYNO~E
j2tn2 h_ 1
1o-12-t----.,...------r------y-----.------,.--------.-
o.o1 0.1 1.0 10.0 100.0 1000.0 10000.0
. h f 2 h h_J h_2
S, (j ) = 2. + 1f + ho + f +f 2
(5.25)
The first two terms define high-frequency phase noise, which can affect signal
tracking, but they are not considered stability terms. The square root of the
Allan variance is a measure of frequency stability. An estimate of clock phase
(or time) stability, in seconds, can be obtained by multiplying the ordinate axis
of Figure 5.7 times the abscissa. However, a more accurate representation in
terms of the coefficients hex is given as a function of time since last measured
or estimated as [12]
(5.26)
At frequencies below about 30 MHz the ionosphere acts almost like a mirror,
bending the path traveled by a radio wave back toward the Earth, thereby allow-
ing long-distance communication. At higher frequencies, such as those used in
satellite radio navigation, radio waves pass through the ionosphere. They are,
nevertheless, affected by it.
The Refractive Index and Phase and Group Velocity The velocity of prop-
agation of a radio wave at some point in the ionosphere is determined by the
density of electrons there. The velocity of a carrier, the pure sinusoidal radio
wave conveying the signal, is actually increased by the presence of the elec-
trons. The greater the density of electrons, the greater the velocity. The net
effect on a radio wave is obtained by integrating the electron density along the
whole path that the signal travels from the satellite to a receiver. The result is
that a particular phase of the carrier arrives at the receiver earlier than it would
have had the signal traveled in complete vacuum. The early arrival is termed a
phase advance.
The increased phase velocity of propagation is related to what is called the
refractive index n by the expression [ 14 J
(5.27)
(5.28)
(5.29)
where Ne is the electron density in electrons per cubic meter, e is the electron
charge, m is the mass of an electron, Eo is the permittivity in free space, andf is
the signal's carrier frequency. Note that the higher this carrier frequency is, the
closer n is to I, and thus the less the ionosphere affects signal propagation. Note
also that the higher the electron density is, the higher the plasma frequency is,
and thus the more the ionosphere affects signal propagation.
On the other hand, the signal that is modulating the carrier (e.g., pseudo-
random noise codes and navigation data) is delayed by the ionosphere. Since
the composite signal can be thought of as being formed by the superposition
194 SATELLITE RADIO NAVIGATION
d I
ng= -(nw)=- (5.30)
dw n
Thus, the relationship of the phase velocity and group velocity (the rate of
energy propagation) satisfies that for a signal passed through a wave guide,
which is
(5.31)
Electron Density and Total Electron Content The electron density is quan-
tified by counting the number of electrons in a vertical column with a cross
section of one square meter [13]. This number is called the total electron con-
tent (TEC). The TEC is a function of the amount of incident solar radiation.
On the night side of the Earth, the free electrons have a tendency to recombine
with the ions, thereby reducing the TEC. As a consequence, the TEC above a
particular spot on the Earth has a strong diurnal variation.
Changes in TEC can also occur on much shorter time scales. One of the phe-
nomena responsible for such changes is the traveling ionospheric disturbance
(TID). TIDs, which have characteristic periods on the order of I 0 minues, are
manifestations of waves in the upper atmosphere believed to be caused in part
by severe weather fronts and volcanic eruptions. There are also seasonal varia-
tions in TEC and variations that follow the sun's 27-day rotational period and
the roughly 11-year cycle of solar activity.
If the carrier frequency is chosen to be high enough so that the fourth-order term
is negligible, then the deviation of n from a free space value of one is inversely
proportional to f 2 . Thus, by making measurements on two widely spaced fre-
quencies and combining them. the electron density Ne can be determined, and
ATMOSPHERIC EFFECTS ON SATELLITE SIGNALS 195
almost all of the ionospheric effect can be removed. This is true whether the
measurements are pseudorange, Doppler, or integrated Doppler measurements,
since they all have an error component that is a function of the refraction
index. However, the correction is applied with a different sign, depending upon
whether the measurement is obtained from the carrier or the modulated signal.
Applications of ionospheric corrections using models or dual frequency mea-
surements are peculiar to satellite radio navigation system, and thus will be
described later in this chapter.
N = I 0\n - I) = --:r-
77.6 (
.P + -4810e)
T- (5.33)
where Pis total pressure in millibars, Tis absolute temperature inK, and e is
partial pressure of water vapor in millibars, where one definition of e is given
as [ 15]
77.6 373256e
N(h) = ------:;- P.fd(h) + T2 .fw(h) (5.35)
I. The path does not follow a straight line. The consequence of this is small
and can be neglected except for very small elevation angles.
2. The wave velocity is slightly lower than it is in a vacuum, producing an
apparent increase in the length of the path given as
11L =c JR (n - I) ds (5.36)
()
where s is the curved abscissa on the path and R is the distance to the
satellite, which can be treated as infinite.
Since the real path does not deviate much from a straight line for all but very
small elevation angles, Equation 5.36 can be made a function of elevation angle
and integrated with respect to altitude. The result is a correction model for
ranging measurements.
In general, the dry and wet terms need to be integrated separately because
they include different functions of altitude and elevation angle. This is the case
when actual surface measurements are used and the ultimate accuracy is desired,
such as in the control segment of the radio-navigation system. However, in the
case of avionics applications, those measurements are not generally available,
and a standard day is used to define the correction model. Then, commensurate
with the accuracy of that standard day, Equation 5.36 becomes
11L = f= (n - I) dh (5.37)
ho sin ¢o
where ho is the altitude of the user and ¢o is the elevation angle of the signal
path to the satellite. An approximation of 11L can be obtained by assuming the
atmosphere is exponential. That is, let
so that
n0 - I
i1L = e~bho (5.39)
b sin ¢o
over the years, the number being tied to budget constraints. The intent, how-
ever, is to provide coverage at all locations on the Earth as nearly to 100% of
the time as possible. Each satellite transmits signals at two frequencies at L-Band
[ 1575.42 (L I) and 1227.6 (L2) MHz] to permit ionospheric refraction corrections
by properly equipped users [ 18]. These signals are modulated with synchronized,
satellite-unique, pseudorandom noise (PRN) codes that provide the instantaneous
ranging capability. Those codes are modulated with satellite position, clock, and
other information, in order to provide the user with that information. Details on
the constellation and signal structure appear in later sections of this chapter.
All equations of Section 5.2 apply to the GPS navigation solution in
that all three measurement capabilities-ranging, Doppler, and integrated
Doppler-exist, and, in general, the solution for the user's clock and clock drift
are required. It is not uncommon for a specific user equipment to use at least
two of the three measurement capabilities to simultaneously solve for position,
velocity, clock offset, and clock drift, in some cases using all satellites in view
[19, 20]. That number can be as high as 12.
System Accuracy GPS provides two positioning services, the Precise Posi-
tioning Service (PPS) and the Standard Positioning Service (SPS) [21 ]. The PPS
can be denied to unauthorized users, but the SPS is available free of change to
any user worldwide. Users that are crypto capable are authorized to use crypto
keys to always have access to the PPS. These users are normally military users,
including NATO and other friendly countries. These keys allow the authorized
user to acquire and track the encrypted precise (P) code on both frequencies
and to correct for intentional degradation of the signal.
Encryption of the precise code provides GPS with an anti-spoofing (A-S)
capability. A-S is not meant to deny the P code to unauthorized users but to pre-
vent the spoofing of the precise code by an unfriendly force. Unfortunately, A-S
denies the P code to unauthorized users. Thus. A-S prevents these users from
correcting for ionospheric refraction, since the L2 signal only carries the P code,
although there are "codeless cross-correlation" techniques that do allow this
measurement [22, 23]. A-S does not prevent the use of the coarse/acquisition
(Cj A) code, which is only carried on the L I signal.
The intentional degradation. on the other hand, is meant to deny accuracy
to an unfriendly force. It is called selective availabilit_v (SA). Unfortunately,
SA also denies accuracy to unauthorized users that are friendly, which is the
entire civil community. The peao~-time policy of the DoD is to provide an
unauthorized accuracy of I 00 meters, 2drms (horizontal accuracy) [24].
Either A-S or SA. or both, may be turned on. If neither is turned on, SPS
accuracy is the same as PPS. In 1996, the U.S. stated that it is its intention to
turn off SA within a decade [135]. More details on GPS system accuracy are
provided later in this chapter.
The GPS Segments GPS has the basic system configuration illustrated in Fig-
ure 5.1. The monitoring and satellite control sites are dispersed around the
L1 AND L2
NAVIGATION
SIGNALS
• MONITOR STATIONS
(5 .40)
where the ith row (for satellite i) of the measurement matrix H (h;) is given as
(5.41)
in terms of Equation 5.3. (The first three elements of Equation 5.41 are the
directional cosines from the user position to the satellite position.) GDOP is
defined as
where trace[·J indicates the sum of the diagonal elements of [·]. If the vector
of pseudorange residuals defined in Equation 5.40 are statistically uncorrelated
with equal !-sigma errors of aPR· then the !-sigma position/time error is
I 1 1 2 , 2
(5.43)
ax. 1 c:c GDOP ·apR=ya~+a~+a:+c-a 11 ,
HLTr = HTEcEF~LTP =
cos El 1 cos Az1
COS El1 COS
-:
Azo
-
cos El 1 sin Az1
cos El2 sin Az2
sin El1
sin El2 -1]
-I
[
cos EIN cos AzN cos EIN sin AzN sin EIN -1
(5.44)
where El; and Az; are the elevation and azimuth angles from the user to satellite
i and where
202 SATELLITE RADIO NAVIGATION
~]
- sin Au sin ¢u COS ¢u
COS Au 0
(5.45)
sin Au cos tPu sin ¢u
0 0
where ¢u and Au are the estimated latitude and longitude of the user. This new
H matrix then replaces the one in Equation 5.42, and fix, fiy, and fiz are replaced
with fiN, fi£, and fih, respectively. Then, the sum of only the first two diagonal
elements in the trace of Equation 5.42 are included for computing HDOP, and
only the third diagonal element is included for computing VDOP. The linear
transformation of Equations 5.44 and 5.45 is only appropriate for the transfor-
mation of location residuals and is not appropriate for total state ECEF coor-
dinates to latitude, longitude, and altitude, which requires a nonlinear transfor-
mation (see References [9] and [23] of Chapter 2).
Adequate coverage is usually defined by the U.S. DoD when PDOP is less
than 6 for elevation angles greater than 5°. However, in some applications of
GPS, such as civil and commercial aviation, this is not adequate coverage. In
these applications, augmentation using pseudolites or geostationary satellites
is necessary. These augmentations are discussed later in this chapter. Because
of these advanced applications, the DOP concept has been replaced with the
concept of availability of accuracy, accounting for satellite failures. This avail-
ability is based upon analytical procedures using the mean-time-between-fail-
ures (MTBF) and mean-time-to-repair (MTTR) characteristics of satellites (and
appropriate augmentations), a concept originally developed by the French space
agency CNES [27, 28] and later extended [29, 30]. With these extensions, cov-
erage is more appropriately defined in terms of probabilities (or availability).
0.98 1( . . ···_
. _1
__ · · · J24 GPSI s. ATEL~ITES ·-
/I: I
~
....1
0.975
-/~
I : : I
L_____
Ll···································
- ___L.__ _
-
..L__----,
co
<! 0.97 t 1E+01 , ,
I
i.
....1
<i:
> 0.965
I:
: ~
m
1 E-01
1 E-o3
r---.....-------r-
-- r ~
I I I I
<! 1
•
I: ! ~ 1 E-05 I - I ···-·
~0.. ~ 'r--..._
0.96
I :' ! 1E-07
• ' I
0.955 I~ 1 ! ~ 1E-09 1'-..1
I: ! ~ 1E-11 I I ..........
0.95 :
I;
I ;;: 1E-13 t-------1- j ~
: 1E-15 ,
0.945 I: ! o 1 2 3 4 5 s 1 a 9 10 11
i I; j NUMBER OF SIMULTANEOUS SV FAILURES
0.94 .
0.5 1.5 2 2.5 .3 3.5 4 4.5 5 5.5 6
N
"ALL-IN-VIEW" DILUTION OF PRECISION OVER THE CONTINENTAL US
0
Vol
:Figure 5.9 Availability of HOOP and VDOP over the continental United States (CONUS).
204 SATELLITE RADIO NAVIGATION
FREQ. --------------------------M'ou-·:
STDS
~~- -------------- -· -- -----------------, I
I
I
I I
I
I
I
I
PROCESSOR lI ·--.
I
-------------------------------------l
hI
MEMORY
110 TX ~An~~~ 1+1 TY ~
I .... , &COMSEC
: I NT!:R-SA TELUTE
:coMMUNICATIONS
DATA I
I
,UPLOAD
I
I
I
I
FORMAT
CONVERTER
I
I
I
I
I
·- -------------- ------------------------ -·
I CTDU :
I
I
I
I
SMCD I
I
I
I
I
I
I
I
POWER I
I
I
I SWITCHING I
I
·---
I
------------· I
EPHEMERIS
--·----·--~I
rt ..L.t..
I I
r-----
DATA
TT&C : :
I
EPS
I
:
'
I
I ADS
I I
I I I I
I
•___________ J
----------J ·-------------~
Figure 5.10 Block diagram of the Block IIR GPS satellite payload electronics (courtesy,
N ITT Aerospace/Communications Division, ITT Defense).
Q
Ul
206 SATELLITE RADIO NAVIGATION
FALCON AFS or
ONIZUKA AFB
• NAVIGATION
PROCESSING
HANDOVER
DATA
EJ
AFSCF
DMA
UPLOAD (TRANSMIT)
GROUND ANTENNAS
•ACTIVE TRACK
•REDUNDANTELECTRON~S
The receiver's channels are periodically calibrated using a dual frequency GPS
signal generator that is slaved to the same frequency reference. This signal gen-
erator is also used to provide signals for fault isolation.
One of the two frequency standards is on hot standby to provide a quick
change-over in case the operational standard exhibits faulty operation. A phase
comparator compares the phases of the two standards to help detect faulty oper-
ation, as well as to provide an initial frequency estimate of the backup standard
for a smooth transition. Both frequency standards are kept alive with a backup
power supply in case of prime power failure.
Meteorological sensors (barometric pressure, temperature, and relative
humidity) are polled by the central processor to permit accurate correction for
tropospheric delay at the MCS. These data, along with the receiver measure-
ments and status information, are forwarded to the MCS over a dedicated secure
communication channel in response to tracking orders received over the same
duplex channels. This channel utilizes a commercially developed SDLC proto-
col to provide error detection and data block re-transmission.
Master Control Station The MCS consists of the processing complex and
controller facility to completely manage and control the operational GPS space
assets and to provide navigation messages. The navigation mission requires an
upload availability of 98%, so redundancy is provided for all mission-critical
equipment. Dual processors with communication controllers and the customary
compliment of peripherals are configured to permit processing of the on-line
navigation processing and satellite control functions with either unit. Personnel
at the MCS control all navigation processing, constellation, and CS assets and
are responsible for the system integrity. This requires established procedures
and efficient access to critical mission data.
The navigation process is illustrated in Figure 5.12. It is based upon a lin-
ear expansion about a reference trajectory, which is obtained by integrating the
equations of motion forward in time from an initial position and velocity state.
The satellite force model used to accomplish this integration includes mathe-
matical expressions for the following effects: ( l) WGS-84 geopotential expan-
sion, (2) sun gravitational attraction, (3) moon gravitational attraction, (4) Earth
gravitational tides, (5) solar flux reaction (including eclipse), and (6) satellite
y-axis acceleration.
Partial derivatives with respect to the epoch states of inertial position and
velocity are generated to reduce the estimation and prediction process to linear
mathematical relationships. Given values for estimated state residuals, evalu-
ation of the first-order expansion provides the position trajectory used to pro-
cess measurements or to generate the navigation message. The inertial-to-ECEF
coordinate transformation matrix is generated to be consistent with externally
supplied Earth orientation data.
The measurement update process consists of data editing, smoothing, mea-
surement model transformation, and estimation steps. The raw MS data are
examined and correlated with receiver fault indicators to ensure track conti-
MS SURVEY DATA INERTIAL PARTIAL DERIVATIVES INITIAL
POSITION
& VELOCITY
~
TRAJECTORY
MEASUREMENT STATE
RESIDUALS RESIDUALS
REFERENCE
PHASE
TRAJECTORY
AIDED
GENERATION
SMOOTHING
DERIVATIVES
A PRIORI POLAR
TIME MOTION
INERTIAL
REFERENCE
DATA
~
A PRIORI INERTIAL POSITION TRAJECTORY
NAVIGATION
UPLOAD GROUND
MESSAGE
PARAMETER
MESSAGE I ... ANTENNA
GENERATION DATA
FITTING
nuity. Both first- and second-difference histories are compared with threshold
values to detect any data inconsistencies. Smoothing is applied to refine the MS
measurements, combining the 1.:5-sec measurements into one low-noise data
point every 15 min. Corrections for the ionosphere and troposphere induced
delays and interchannel receiver delays are made. Then the measurement model
accounts for the Earth's rotation during the signal propagation time, the effects
due to crustal tides, relativistic time distortions and the a priori knowledge
of estimated states to form measurement residuals. Measurement partials with
respect to the estimated states are also evaluated to support the estimation pro-
cess.
A Kalman filter recursive estimator (Chapter 3) is implemented to estimate
the following states: (1) satellite position at epoch time, (2) satellite velocity
at epoch time, (3) satellite clock phase at epoch time, (4) satellite clock fre-
quency at epoch time, (5) satellite clock aging (rubidium frequency standards
only), (6) solar flux, (7) satellite y-axis acceleration bias, (8) MS clock phase
at epoch time, (9) MS clock frequency at epoch time, and (I 0) MS wet tropo-
sphere height. The estimation process is partitioned to conserve computational
resources. A common time scale ensemble of MS clock states is determined and
this common time scale holds the partitions together. The partitioned estima-
tion performance penalty is insignificant when an adequate number of satellites
exist in each partition to maintain a solid MS-time transfer network.
Every eight hours the state of each satellite is used to generate a predic-
tion of the time-scale correction and trajectory in ECEF coordinates. The same
linear expansion and coordinate transformations are used. Navigation message
parameters are obtained as a weighted least-squares fit to these data in accor-
dance with ICD-GPS-2008-PR and the GPS SPS signal specification [7, 8],
and the generated user information is put into a standard format and forwarded
to the GAs for upload. This process is time-phased across satellites to evenly
distribute the MCS work load.
Navigation service integrity is monitored using MS data for four-satellite
position solutions and by comparing the pseudorange and received message
data with MCS expectations. Measurement residuals are evaluated by the MCS
whenever tracking data exist, and performance statistics are maintained as infor-
mation available to system operator personnel and the user community. Each bit
of the received navigation messag,e is compared with the corresponding upload
data-base value to verify proper dissemination.
transmission features. Data from the MCS are buffered on disk prior to satellite
contact as a telemetry data from the satellite prior to transmission to the MCS.
Signal Modulation The PRN codes and data are modulated onto the carriers
using binary phase shift keying (BPSK). This modulation shifts the phase of
the carrier 180° each time there is a change in state of the digitally defined
code and data. First of all, the PRN code is a sequence of 1's and O's, as are
the message data. The message data sequence is modulo-2 added to the code
sequence, which is nothing more than an exclusive or of a code bit and a data
bit, although the data bits transition at a much slower rate than do the codes.
The resulting sequence of I 's and O's are converted to 180" and oc phase shifts
of the carrier, respectively. Since the 180° phase shifts simply change the sign
of the carrier, an equivalent representation is simply an amplitude modulation
of ±1 's. The result is a mathematical representation of the signals as follows:
where A is the Ll P signal amplitude, P(t) and C(t) are the ±I P and C/ A code
PRN sequences, D(t) is the ±I message data bit sequence,.f 1 andh are the Ll
and L2 carrier frequencies, and ¢ 01 and ¢ 02 are the ambiguous Ll and L2 carrier
phases. Note that the two signals are not phase coherent, even though they are
derived from the same frequency reference. There are also slight group (code)
delay differences between the P code and Cj A code modulations and between
the P code modulations at the two frequencies l7l
PRN Code Properties The PRN codes are generated as products (modulo-
2 sums, if expressed as I 's and O's) of two other codes clocked at the same
214 SATELLITE RADIO NAVIGATION
for satellite i, where T, is the C/A code chip width, or the inverse of the C/ A
code chipping rate of 1.023 MHz, Tp is the P code chip width, or the inverse of
the P code chipping rate of 10.23 \1Hz, and 11; is an integer assigned to satellite
i for the C/ A code. In case of the P code, 11; takes on a value between 1 and
37 for 32 satellites and 5 reserved for ground transmitters (GTs) [7, 8]. In the
case of the C/ A code, the 11; are selected values between 1 and I 023 for codes
exhibiting desirable properties [32].
These PRN codes provide the desirable code-division, multiple-access
(CDMA) property that, to an extent, the codes received from the various satel-
lites do not correlate with each other, nor do they correlate with a reference
code in the user's receiver unless the state of the received code matches that of
the reference code. Thus, all satellite signals can be received at the same fre-
quency (except for Doppler differences) and selectively acquired and tracked,
depending upon the selection of the code in the reference code generator. These
PRN codes also exhibit the property that, to an extent, they spread interfer-
ence signals over the signal's bandwidth, providing a degree of interference
rejection. These spread spectrum properties of the GPS PRN codes will be dis-
cussed later in this chapter. Here we will discuss the correlation properties of the
codes.
I. CjA codes. The C/ A codes are Gold codes [ 18], where the G 1 and G2
codes are generated in maximal-length, 10-stage linear feedback shift registers,
each of which generate repeating maximal-length codes of length 1023 chips.
Figure 5.13 represents an implementation of a C/ A coder [32]. In this imple-
mentation the initial state of the G2 register represents the delayed state of that
maximal-length code (the 11; in Equation 5.48) from its initial state of all 1's.
The resulting C/ A (Gold) code is not a maximal-length code. A maximal-
length code x(t) has the autocorrelation property that [33]
2M -I
(5.50)
where k is an integer number and M is the number of stages in the shift register.
This is partly true for the C/ A codes, but they do not have perfect correlation
properties. On the average, a C/ A code has the autocorrelation property that,
for k -J 0 [ 17],
NAVSTAR GLOBAL POSITIONING SYSTEM 215
r------------------------------------------
CODE SET TO
CLOCK ALL 1'S
INITIALIZE
G2 STATE
_..G2(t+niTc)
..._.... .
1
~~~,-L-~~,_~~~.-~~ I
I
I
I
I
I
I
I
I
I
G2CODER :
------------------------------------------·
Figure 5.13 C/ A coder implemented with initial 02 state.
1023
L C(j)C(j + k Tc)
j= I
although the probabilities vary somewhat, depending upon which code is eval-
uated [32]. Also, whereas the maximal-length codes are always balanced, 256
of the I 023 Cj A codes are not. That is [32],
(5 .52)
216 SATELLITE RADIO NAVIGATrON
All of the 32 codes selected for the GPS satellites and the four codes selected
for the ground transmitters are balanced.
The autocorrelation property de:;cribed above carries over to the cross-corre-
lation between the different Cj A codes. This means that the Cj A codes are not
quite orthogonal, and care must be taken during acquisition of the satellite sig-
nals to prevent false acquisitions and false alarms. At zero Doppler difference
between satellite signals, Equation 5.51 above represents a separation of 23.9
dB between signals of equal power. There is also a degree of cross-correlation
between signals at other Doppler differences. This property will be described
later when we discuss the spectral characteristics of the PRN codes.
2. P codes. Whereas the C/ A codes are linear codes, the P codes are non-
linear. That is, the Cj A codes are made up of two maximal-length codes that
are clocked synchronously and are allowed to proceed through all of their I 023
states. This is not true for the P codes. The underlying linear codes of the P
code are short-cycled before creating the product of the codes. This has the
effect of creating an extremely long code. In fact the 37 individual P codes are
simply a one-week piece of a long code that is approximately 38 weeks long.
A typical P code implementation is presented in Figure 5.14 [7]. Note that
there are four shift registers, two each for the X I and X2 code generators. Each
of these 12-stage shift registers that have 2 12 - I = 4095 possible states are
short-cycled, either at 4092 or 4093 states, and reset. Both of the XI shift reg-
isters are reset on X I epochs (every 1.5 seconds), while the X2 shift registers
are reset every 1.5 seconds plus 37 chip clock cycles. All shift registers are reset
at the end of the week. Although the X I and X2 coders repeat every 1.5 and
1.5+ sec, the fact that they are running asynchronously at I 0.23 MHz, their
modulo-2 addition generates an extremely long code. Delaying the X2 code
with respect to the X I code an additional i- 1 chips for the ith satellite pro-
vides codes for each satellite. The count of the X I epochs provides a Z-count
that is used as basic timing for the system to which the data message and the
Cj A coder are synchronized. If th1~ coder were not set at the end of week, the
code would eventually run into the code of the other satellites and return to
the beginning almost 38 weeks later. Unlike the Cj A codes, the P codes have
excellent cross-correlation properties.
3. Spectral characteristics qf" the PRN codes. In the frequency domain the
spectral density of the signal is the spectral density of the PRN codes centered
at ±f;, i = I, 2. At baseband, the spectral density of the P code is
where Tc is the P code chip width, or the inverse of the P code chipping rate
of 10.23 MHz, and Pp; is the appropriate P code carrier power. Actually, this
spectral density is bandlimited in the GPS satellites in order to protect radio
astronomers, so its range is not infinite. This could prevent a user's receiver
:--lAVSTAR GLOBAL POSITIONING SYSTEM 217
COUNT TO
RESET
403,199
RESET (ONE WEEK)
-- ------------------------------'
'' + 3750
X1A EPOCHS
FOR X1
2.5 kHz EPOCH
L----..JX11.5 SEC EPOCHS
15,345,000 CHIPS
DELAY i
15,345,037 CHIPS CHIPS FOR
SATELLITE
'
, X2 CODE GENERATOR :'
~---------------------------------·
Figure 5.14 Typical P coder implementation.
from achieving full correlation. However, all P code receivers also are band
limited. The resulting effect is known as correlation loss due to filtering.
The C/ A code for satellite i is a short one millisecond repeating code whose
spectral density is a line spectrum with components cii• j = -oo to +oo, where j
represents spectral lines I kHz apart and
where P,; is the C/ A code carrier power. The envelope of this line spectrum
takes on the form
218 SATELLITE RADIO NAVIGATION
(5.55)
However, the spectral lines deviate from this envelope significantly. Figure 5.15
illustrates this for the first two lobes of the spectrum. As stated earlier, the C/ A
codes do have a level of cross-correlation at Doppler differences. It has been
shown that this cross-correlation level is approximately equal to the magnitude
of the spectral line component of another Cj A code in the family of l 023 [ 18].
Since the line spectrum shown in Figure 5.15 is typical of all the codes as
far as variation about the envelope is concerned, that spectrum shows typical
levels for this type of cross-correlation-on the order of 2I dB below the carrier
power. Thus, the cross-correlation at Doppler differences can be higher than at
zero Doppler offset. However, its occurrence is quite rare, and a receiver would
never track it because it would disappear as the relative codes move past each
other. These occurrences can cause false alarms during initial signal acquisition.
Frames, Subframes, and TLM and HOW Words The GPS navigation mes-
sage consists of a frame of five 300-bit subframes spanning 30 seconds of time
as illustrated in Figure 5.16 [7, 8]. Each six-second subframe consists of ten
30-bit words, the first two of which repeat in each subframe. These words, also
illustrated in Figure 5 .16, are the Telemetry Word (TLM Word) and the Hand-
Over Word (HOW Word). These two words are generated by the satellite, while
the other eight words are generated by the CS and uploaded to the satellite as
0
-10
N
::I: -20 - - - PRN 2 LINE SPECTRUM
u
co
"C
' -30
>-
!:::
C/l
z -40
w
0
...1
ct -50
et:
t-
(.)
w -60
Cl..
C/l
w
0 -70
0
(.)
~ -80
(.)
-90
-100
0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000
FREQUENCY OFFSET- kHz
Figure 5.15 PRN 2 CI A code spectra l dens ity.
220 SATELLITE RADIO NAVIGATION
SUBFRAME
NUMBER
I
I
10 30-BIT WORDS; 6 SECONDS; ONE SUBFRAME ~I
I I I
rl
L2 FLAGS, URA, WEEK NUMBER AND
TLM HOW CLOCK CORRECTION PARAMETERS
w
2 I I ITLM HOW
EPHEMERIS PARAMETERS
::;:
<(
0:::
u..
w
z
0
(/)
0
l
TIME OF START OF NEXT SUBFRAME LALERT l...LBITS TO CAUSE
FLAG FINAL PARITY
A-S FLAG BITS TO ZEROS
Figure 5.16 GPS navigation message frame and TLM and HOW words.
described in Section 5.5.4. The TLM Word consists of an 8-bit preamble and a
satellite telemetry message. The preamble (1 000 I 011) allows the user's receiver
to synchronize with the subframes to establish time of reception and subsequent
resolution of the Cj A code ambigu:tty. This ambiguity exists because the signal
transit time from the satellite is on the order of 80 msec, while the length of the
repeating Cj A code is only 1 msec. For the most part, the telemetry message
is of no use to the user; it primarily provides real-time telemetry information
from the satellite to the CS.
The HOW Word contains a 17-bit truncated Z-count indicating the time of
NAVSTAR GLOBAL POSITIONING SYSTEM 221
the start of the next subframe in units of six seconds, or four X 1 epochs. It
also contains two flags, an A-S onjoff flag and an alert flag, and a subframe
identification (ID) (1 to 5 for subframes 1-5). The Z-count is used by the user's
receiver to establish GPS time. For P code users the Z-count also provides the
time information required to hand over to P code, and thus the reason for the
name of the word.
The A-S on/off flag alerts non-PPS users as to whether or not they are able
to acquire the P code, and it alerts PPS users as to whether or not they should
encrypt their reference P code [7]. The alert flag indicates to unauthorized users
that the satellite user range accuracy (URA) may be worse than indicated in
subframe I and that they should use the satellite at their own risk [7].
Parity Both the TLM and HOW words contain satellite-generated parity, six
bits per word. The other eight words of each subframe contain MCS-generated
parity. The HOW Word contains two noninformation bearing bits so that the
last two parity bits can be forced to O's, since the parity algorithm overlaps
word boundaries.
The parity algorithm is a (32, 26) Hamming error detection algorithm. It
will detect 1-, 2-, or 3-bit errors. As stated above, the parity overlaps the 30-
bit words, since it is based upon 32 bits, always including the last two bits
of the previous word. Thus, since the parity in the TLM and HOW words are
generated in the satellite, and the parity on the other words are generated in the
MCS, two bits in the HOW Word and Word I 0 of the MCS generated words
are wasted. This is so that both the satellite and the MCS know the state of the
last two bits of the previous word (set to 0) when generating the TLM Word
(satellite) and Word 3 (MCS). The Hamming error detection algorithm used is
specified in [7, 8].
Words 3 through I 0 of each subframe contain the message data. Subframes 1,
2, and 3 repeat every 30 seconds consisting of the same data for nominally one
hour but can change more often or less frequently. Subframes 4 and 5 subcom-
mutate 25 times each, so that a complete data message requires the transmission
of 25 full 1500-bit frames, or 12.5 minutes. Thus every 30 seconds one page of
subframes 4 and 5 are transmitted. These pages consist of less timely data that
change whenever the satellite is uploaded, while sub frames I, 2, and 3 consist
of more timely data that change periodically. This periodicity is either once per
hour, on the hour, once per 4, 6, 12, or 24 (or more) hours on 4-, 6-, 12-, or
24-hour (or more) boundaries. Any period greater than an hour represents a
case when a satellite could not be uploaded within a day of a previous upload,
which would be rare. In addition, the data could change at any time with a
quick upload. Again, this is a rare condition required to change substandard
or erroneous data. One-hour data sets are valid for three additional hours. All
other data sets are valid for two hours past their period of transmission.
Other contents include the GPS week number, URA, health of satellite, various
flags, and data reserved for authorized users. Content details are given in Table
5.3.
t = t, - !::.t,,. (5.56)
where
(5.57)
where afO· an, and af2 are the polynomial coefficients transmitted in subframe
1, toe is the clock data reference time in seconds, and !::.t r is the relativistic term
in seconds computed as
2X,·V, r:.
!::.tr = - - ·-2 - · =FeY A Sll1 Ek (5.58)
c
where the orbit parameters e, A, and E are derived from the data contained in
subframes 2 and 3, and
NAVSTAR GLOBAL POSITIONING SYSTEM 223
-2.j/i
F = = -4.442807633 x 10-IO sec/~ (5.59)
c2
where
is the speed of light. These relativistic effects on the GPS satellite clocks are
described in [35]. It suffices to use the uncorrected value of ts in the evaluation
of the above equations, since, prior to correction, it is within 1 msec of the
corrected value.
Since the clock correction parameters estimated by the MCS are based upon
L I and L2 measurements combined for correction of the delay through the iono-
sphere, and since there is a potential L I /L2 differential bias in the satellite, the
L !-only user must revise his computations of the satellite's clock offset with
[7, 8, 36]
(5.62)
that of a curve fit over the time of validity (fit interval) and not necessarily
an entire orbit. A fit interval flag in subframe 2 indicates whether the interval
is four hours or greater than four hours. lODE has similar meaning for the
ephemeris parameters as IODC has for the subframe I parameters.
The user receiver applies the parameters to a variation on the equations given
in Section 5.3. However, not all the parameters are indicated in that section.
These additional parameters reflect orbit drift and other orbit perturbations such
as the sun and moon gravitational forces and solar radiation pressure. Typical
equations for the computation of the satellite's position in ECEF coordinates
are given in Table 5.5 [7, 8].
The accuracy of curve fit used in generating these parameters and subsequent
truncation is quite good. For the four-hour fit, the user range error (URE) based
on a projection of the curve fit error onto the user range is less than 0.4 meters,
1-sigma. For a six-hour fit, the URE degrades to 1.6 meters, 1-sigma. The equa-
tions provide the satellite's antenna phase center position in the WGS-84 ECEF
reference frame defined in [7, 8].
1. Satellite almanacs. The ephemeris and clock almanac data are a subset
NAVSTAR GLOBAL POSITIONING SYSTEM 225
!1k = !1o + (0- Or:)tk- OEtoc Corrected longitude of ascending node at time
tk accounting for earth's rotation rate !1 e
Xk = x~ cos !1k- y~ sin !1k cos ik ECEF coordinates at time tk
Yk = x~ sin !1 k + y~ cos !1 k cos ik
Z.k = Yk' ..
Sin lk
gitude. Details of the algorithm are given in [7, 8] and [36]. The algorithm
reduces the delay error for the single-frequency users on the order of 50%
to 60% [37].
3. UTC parameters. Page 18 of subframe 4 also provides the parameters
needed to relate GPS time to UTC and notices to the users regarding the
scheduled future or recent past changes due to leap seconds. The param-
eters include the number of integer seconds between UTC and GPS time,
plus first-order polynomial coefficients describing the drift between GPS
time and UTC. However, this drift is kept to a minimum by the MCS by
steering GPS time toward UTC modulo I second. Details of this relation-
ship are described in l7, 8].
dPR; . . . .
-- =
dt R; + cb.t,;
. ~ cb.tu + EPR I (5.64)
dPR; ~ ~ . .
0 -- = 1; · oXu ~ cf1tu + EPR
dt I
(5.65)
Since the L2 delay is greater than the Ll delay, the correction is negative,
decreasing the pseudorange.
3. Tropospheric refraction correction. This correction is made by subtracting
the effects of Equation 5.36. There are many variations to this correction,
depending upon the desired accuracy or complexity.
228 SATELLITE RADIO NAVIGATION
(5.67)
(5.68)
Uncorrected Error Sources After all the corrections stated above, residual
error sources remain. Among the most predominant are SA errors for the unau-
thorized user (SPS mode). The next most predominant error source is the resid-
ual ionospheric refraction error for the Ll-only users. Next is the standard sys-
tem errors in the satellite clock parameters, followed by multipath errors, the
tropospheric refraction, the satellite ephemeris, and finally receiver noise. New
receiver technology has reduced the effects of multipath and receiver noise,
with the exception of operating in a jamming environment.
(5.70)
is a 2Nx8 matrix mapping the position, velocity, and time state residuals into the
2N measurement residuals. The H matrix could be the HLTP matrix of Equation
5.44 if the solution is in the north, east, and up domain. Since the HT is general
not a square matrix, the solution is overly determined. In that case, the solution
for oPVu in a weighted least-squares sense is
(5.71)
(5.72)
· PREAMPLIFIER
RF~ I DOWN-
CONVERTER
I >I I IF
> SIGNAL
PROCESSING
PR,,
"' NAVIGATION
PROCESSING "'
POSI~ ON,
) )~ j' PHASE,
)~
VELOC TY,
'
LO
LOs& CLOCKS
DATA TIME, E TC.
'",
REFERENCE FREQUENCY
OSCILLATOR SYNTHESIZER INTERRUPTS
..
Figure 5.17 Generic GPS avionics receiver functional block diagram.
N
~
......
232 SATELLITE RADIO NAVIGATION
low-noise amplifier (LNA). Its primary function is to set the receiver's noise
figure (see Section 4.2.1) and to reject out-of-band interference. The parameters
that dictate the preamplifier requirements are the unwanted RF environment as
received through the antenna and losses that precede and follow the preamplifier
and desired system noise figure (or noise temperature) as derived from overall
receiver performance requirements. The gain of the preamplifier is not a system
level requirement per se, but a derived requirement that satisfies the system level
requirement.
The unwanted RF environment as received through the antenna affects the
preamplifier in two ways. It can cause damage to the preamplifier electronics,
or it can cause saturation of the preamplifier and circuitry that follows. Nor-
mally, except for damage prevention, one can do nothing to suppress the RF
environment, as passed by the antenna, at frequencies that are in the bandwidth
of the desired GPS signal. That environment is considered to be either jam-
ming or unintentional interference. There are, however, more advanced tem-
poral interference suppression techniques that can be used to suppress narrow-
band interference [121, (Chapter 20)]. Suppression of the RF environment out
of the desired GPS signal band can be accomplished by filtering before, during.
and/or after amplification. When It is accomplished it is based upon a trade-off
between system noise figure requirements, filter insertion loss. and bandwidth
efficiency. Suppression of in-band and out-of-band damaging interference is
usually accomplished with diode-; that provide a ground path for strong sig-
nals. In the case of lightning protection, more complex lightning arrestors may
be used.
The system noise figure is set using a low-noise amplifier (LNA) that pro-
vides enough gain to cause any losses inserted after the LNA to have a negli-
gible effect. An LNA cannot account for losses inserted before its operation or
for its own noise floor.
The reference oscillator provides the time and frequency reference for the
receiver. Since GPS receiver measurements are based on the time of arrival
of PRN code phase and received carrier phase and frequency information, the
reference oscillator is a key function of the receiver. Its output is used by the
frequency synthesizer, which converts the oscillator output to local oscillators
(LOs) and clocks used by the receiver. One or more of those LOs are used by
the downconverter to convert the radio frequency (RF) inputs to IF frequencies.
The signals are easier to process in the IF section of the receiver.
The requirements on reference oscillators for avionics receivers vary depend-
ing upon the avionics application. A high-quality oscillator can be the most sig-
nificant cost item of a modern receiver. Thus, there are compromises made on
oscillator performance. There are some commercial and military applications
where refrence oscillator performance is critical.
Typical requirements parameters applied to reference oscillators are as fol-
lows:
Figure 5.18 OEM GPS receiver card (courtesy, NovAte! Communications, Ltd.)
I 00 mw. The antenna is .1.5 in. in diameter and weighs 5 oz. It meets
the requirements specified in references [42] and [ 122].
'"l Reference oscil/atm: Like most commercial receivers, this receiver uses a
small TCXO as its reference oscillator (on the center of the hoard shown
in Figure 5.18). These small TCXOs have marginal stability and phase
noise characteristics for some GPS applications, hut they are satisfac-
tory when used in conjunction with a multiple-parallel-channel receiver
with tracking loop bandwidths commensurate with dynamic applications.
Note that the reference oscillator output (20.473 MHz) is used directly
as clocks for sampling and signal processing.
3. Synthesi::.erjdmvnconvertcr. A block diagram of the synthesizer and
downconverter is shown in Figure 5.21 ]4.1]. A commercial synthesizer
chip with a programmable divider (prescaler) is used to phase lock a
voltage-controlled oscillator (VCO) to the reference oscillator. The VCO
output frequency is doubled to provide the downconverter LO.
4. Filtering. IF filtering is realized with a 8 MHz bandwidth surface acous-
tic wave (SAW) filter centered at the IF frequency of .15 MHz ]4.1].
Although this SAW filter has a significant group delay, the delay is sta-
ble, and the filter provides excellent rejection of out-of-band noise and
interference.
5. IF sampling and A/D conversion. The IF sampling and A/D conver-
sion process is illustrated in Figure 5.22. This processing includes a
precorrelation automatic gain control (AGC) ]4.1] that controls the level
y t-11CROPROCESSQR
lf:> DOWN~
CONVERTER
LO
FILTERING ANDND
FUNCTION!>
~FILTERING
~ACQUISITION
ALGORITHMS
~TRACKING
LOOPS
LOCK
DETECTION
~DATA RECOVERY
REFERENCE SYNTHESIZER ~MEASUREMENT
OSCILLATOR NCO
PROCESSING
NCO AND
REFERENCE
CLOCKS
Figure 5.19 Functional block diagram of commercial receiver (courtesy, NovAte! commu-
nications, Ltd.)
N
~
......
238 SATELLITE RADIO NAVIGATION
GPS
oo
cNOi PA\N1
PROGRAMMABLE
INPUT
N=9632
SIGNAL PROCESSING
Q__""-PHASE
CLOCK
~DETECTOR
SAMPLE CLOCK
(5.73)
or
l,k, -Q,,k, -l,k, Q,k. l.,k. -Q,k, -l,k. Q,k, ... (5.74)
where Ck. [)" and ¢k are the code, data bit, and signal phase at sample
time tk.
6. Doppler removal (phase rotation). The Doppler removal process of Fig-
ure 5.19 is part of the signal-phase or frequency-tracking function, which
is a complex multiplication between the signal I and Q samples and ref-
erence I and Q samples generated by the carrier number-controlled oscil-
lator (NCO). This NCO is controlled by the microprocessor's portion of
the carrier-tracking loop. Since the Cj A PRN code Doppler [s related
to the carrier Doppler by a factor of 1540, the carrier NCO also outputs
the basic code clock Doppler correction, which is further corrected by
the microprocessor's code tracking loop function with code phase cor-
rections. In some receiver implementations a completely separate NCO
is used to derive the code Doppler and phase.
7. Coder. The implementation of a typical C/ A coder is shown in Figure
5.13. This implementation, in which the G2 state is initialized, allows the
generation of all I 023 codes in the C/ A code family. which is important
for future implementations !32].
8. Correlators. The correlation process for narrow conelator processing is
illustrated in Figure 5.23 [43, 44]. In this process early, punctual, and
late codes are derived in a shift register that shifts the (early) C/ A code
from the coder at a clocking rate defined by the desired early/late cor-
relator spacing (in a fraction of a C/ A code chip). The dual correlation
process is realized by performing a multi-bit exclusive or between the
single-bit PRN codes and the multi-bit I and Q samples. A discrimi-
nator selection process allows the selection of either early and late or
early-minus-late and punctual correlation. Early and late correlation is
used during the signal acquisition process using the maximum (approx-
imately) one chip spacing :N I0 in Figure 5.23) for rapid acquisition.
Early-minus-late and punciUal correlation is used during tracking using
a 0.1 chip spacing (N •·· I) for optimum parallel code (early-minus-late
times punctual) and carrier (punctual) tracking. This dynamic spacing
concept provides fast acquisition and C/ A code-tracking performance
NAVSTAR GLOBAL POSITIONING SYSTEM 241
I SAMPLES
------l~\.6)
K/..~-----------'
.. E-L OR EARLY
""I SAMPLES
PUNCTUAL OR LATE
r-------------------r--~ISAMPLES
Q SAMPLES
~~-----------------r--~E-LOREARLY
Q SAMPLES
PUNCTUAL OR LATE
r-~-PU_N_C_T_U_A_L_O
__R______-r--~QSAMPLES
LATE CODE
DISCRIMINATOR SELECT
~
SAMPLE CLOCK
~ EARLY CODE
(5.77)
where BL is the single-sided tracking loop bandwidth and Bu- is the two-
sided predetection bandwidth. The d for N = I is one-tenth that for the
conventional N = I 0, while the signal-to-noise density for Cj A code
tracking is twice that for P code tracking.
9. Postcorrelation filtering (accumulators). After correlation, the two sets
242 SATELLITE RADIO NAVIGATION
(5.78)
(5. 79)
Military GPS Aviation Receivers In 1996. the standard GPS receiver for mil-
itary avionics was the miniature airborne GPS receiver (MAGRJ produced
by Rockwell International, which also produces another variety intended for
embedded applications. where the entire receiver (the miniature GPS receiver.
MGR). less antenna, is housed inside another avionics assembly such as an
inertial navigation system [51[. The requirements for the MAGR are defined
in [52[. Guidelines for embedded military receivers are specified in [5~ I and
described in [54[. The MAGR receiver is illustrated in Figure 5.24.
""'""'
MICRQ!'flQQE;>l;Qfl
RF F'REAMPLIFIGATION ~-~liONS
IF
AND FILTERING
-FILTERING
-ACQUISITION
ALGORITHMS
1ST
-TRACKING
LO LOOPS
-LOCK
DETECTION
·DATA RECOVERY
REFERENCE
OSCILLATOR
SYNTHESIZER CARRIER
NCO
-MEASUREMENT
PROCESSING
NCO AND
REFERENCE
CLOCKS
Figure 5.27. Dual downconversions are accomplished, one for Ll and one for
L2. The LO for these downconversions are common, converting both RF fre-
quencies to identical IF frequencies. This LO is derived in the synthesizer,
which is driven with the output of an ovenized reference oscillator at a fre-
quency of approximately I 0.95 MHz. The synthesizer also generates in-phase
and quadraphase LOs for later con version to baseband and clocks for the signal-
processing function.
----------- --------------1
154F0
(L1)
)I'V\_---'>•~L1
\61 r ~··· ~
I
I
I TO 17 F0
I
I BPF
I
120 F0 ~
~-:,L2
(L2)
I
I
I
I
I
I
I
I I ---..1~.+.·.'.J...'
i
CONTROL VOLTAGE TO PLL CHIP
FROM PLL CHIP F : 10.23 MHz
0
.tr~~-- .~~~~~~~~~~---~
17.25 F0 (0 DEGREES)
PHASE LOCK
DETECTOR
17 F0 :'
IF IN :
''
''
'''
:__ ----T---------------------------------------------------------------------------- J'
AGC VOLTAGE
• Because of the very high dynamic requirements (9g, lOg/sec [59]), the
MAGR performs carrier frequency tracking, instead of carrier phase
tracking, by minimizing tan- 1(Qk/h)- tan- 1(Qk-J/h- 1) (time differ-
ence of carrier phase error). Data are demodulated differentially by
observing the changes in the sign of h.
• Instead of verification of phase lock, measurement validity is deter-
mined comparing signal-plus-noise power L~= I u;
+ Q~) to a threshold
based upon an estimate of noise power.
• Code-tracking loop aiding processed from corrected external inputs
from an inertial navigation system to achieve a high AJ tracking capa-
bility [52, 60].
• L2 tracking for ionospheric delay corrections is normally performed
sequentially on the fifth channel. However, if it is determined that L I
is jammed, L2 tracking can be performed on some or all channels [52,
60].
6. A-S and SA processing. The MAGR incorporates a precise position-
ing service-security module (PPS-SM) and auxilliary output chip (AOC)
devices to perform A-S and SA processing [52, 59]. The AOC devices are
provided for each channel to allow tracking of the encrypted code when
the receiver is properly authorized. The PPS-SM performs crypto-key pro-
cessing and management for the A-S and SA processing. It operates on
battery power so that keys may be loaded or zeroized without receiver
prime power. A dedicated data path from the PPS-SM to the AOC devices
prevents sensitive data from being handled by the other processors.
1. Selective availability errors. Although these errors are not biases, they
have correlation times that are long enough to be eliminated if the cor-
NAVSTAR CJLOB A L POSITIONI NG SYSTEM 249
DIFFERENTIAL
REFERENCE
STATION
Figure 5.29 Dillerential GPS concept.
rcction update rate is high enough. Typical pseudorange errors are about
JO meters, 1-sigma. but they have the potential to be hi gher.
,.., lonos;Jheric dclurs . These propagation group delay errors can be as high
as 20 to 30 meters during the afternoon hours to I to 6 meters at ni ght.
if not removed usi ng two frequ ency correcti ons. The single freq uency
mode l will reduce this by approximately .'iO'/c . These errors are slowly
varying biases but spatiall y dccorrclatc over larger distances.
J. Tro;}().lplieric dclovs. These propagation delays can be as much as 30
meters to a low-elevation satellite but are quite consistent and can be
modeled. However. variation s in the index of refraction can cause dif-
ferences between the reference station and the use r of I to 3 meter.-; for
low-ele vati on satellites. T hey arc also slow ly varying biases. and they
spatially decorrelate over larger distances.
4. E'p/u' nle ri.l errors. Normally, the difference between the actual satellite
location and the locati on computed from the broadcast e phemeris is small ,
less than I to 3 meters, but this error can be increased signi fica ntly with
selective availability. Ephe meri s errors arc very slowly varying biases but
can spatially decorrelate over large distances .
.'i. So!e!lilc clock errors. The differences between the actual satell ite clock
250 SATELLITE RADIO NAVIGATION
time and that computed from the broadest corrections can become large
if a satellite's clock is misbehaving.
dPRC;(to)
PR;c(t) = PR;(t) + PRC;(to) + (t - to) (5.80)
dt
where PRC;(t 0 ) and dPRC;(t 0 )jdt are the broadcast corrections for satel-
lite i at their time of applicability to.
2. Delta-differential GPS corrections made up of corrections to the broad-
cast corrections applicable to the previous issue of satellite navigation
data (IOD) for a period of time after an IOD change. These delta cor-
NAVSTAR GLOBAL POSITIONING SYSTEM 251
(1) to ensure the higher update rates using existing data links, (2) to provide
additional information required for the precision approach and landing appli-
cation, and (3) to increase the integrity of the broadcast with a much stronger
parity algorithm. The flight test results described above suggest that DGPS,
using differential pseudorange corrections, can meet even Category III preci-
sion approach and landing requirements [68].
Wide Area DGPS RTCA Special Committee SC 159 is also preparing require-
ments for the use of Wide Area DGPS (WADGPS) as part of the FAA's future
Wide Area Augmentation System (WAAS) to achieve a Category I precision
approach and landing capability [69]. WAAS uses a braodcast through a geosta-
tionary satellite to provide corrections over a very wide area, such as the conti-
nental United States (CONUS). Th1~ accuracy of this approach suffers somewhat
because of spatial decorrelation and limited broadcast bandwidth but is expected
to provide the required Category I precision approach accuracy over a region
such as CONUS. To achieve this, the broadcast messages differ significantly
from the RTCM messages because information on ephemeris and ionospheric
errors must be provided to correct for spatial decorrelation of these errors. More
detail on the WAAS is provided in Section 5.7.3.
few techniques available to improve this initialization process. Two of the most
promising are the use of dual-frequency receivers [72] and the use of local near-
Ll transmitting pseudo lites (pseudosatellites) [73 J. The former approach uses
the differential carrier between Ll and L2 to first resolve ambiguities using a
larger beat frequency wavelength, and then transferring the solution to initialize
the Ll ambiguity resolution. The larger wavelength ambiguity takes much less
time to resolve. The second approach takes advantage of the rapidly changing
geometry between the pseudolite and the user, providing more leverage to the
resolution problem. The use of pseudolites for DGPS is discussed further in
Section 5.7.4.
The use of carrier phase techniques for precision approach and landing is still
in development. However, carrier-smoothed-code techniques are much more
robust. These techniques solve for the ambiguity not as an integer but as a
floating-point number. The ultimate accuracy of DGPS for precision approach
still remains unknown.
I. The horizontal CEP (circular error probable) budgeted for the MAGR is
given as [52]
Unauthorized Unauthorized
C/A Code C/A Code DGPS Special
Authorized Ll User Ll User Category I
Error Source Ll/L2 User with SA Without SA User
3. The 2drms budgeted for the commercial avionics user is given as [40]
For the special Category I DGPS user [41], the space/control segment error
budget is reduced to residual SA, residual clock, and spatial decorrelation errors
amounting to 0.5, 0.01, and 0.0 I meters, respectively. The reference station
budget is set at 1.1 meters, resulting in a total spacejcontrol segment/reference
station error budget of 1.21 meters, root-sum-squared. The users' error bud-
get presented in Table 5.7 is an example. The special Category I user has the
choice in allotting his budget between sensor (GPS receiver) error and flight
technical error (FTE), which defines the pilot's or auto-pilot's ability to fly the
prescribed flight path. The total vertical error budget is 9.76 meters, with a prob-
ability of 95%. Thus, the GPS receiver error budget depends upon the assigned
FTE for a given aircraft. Using a 95% probability VDOP (ratio of 95% vertical
navigation error to !-meter rms pseudorange error) of 3.1, the 95% probabil-
ity vertical error for the pseudorange error budget in Table 5.7 is 5.52 meters,
leaving a 95% probability budget for FTE of 8.05 meters. An aircraft with a
good autopilot could use a receiver with larger errors. The HDOP is not speci-
NAVSTAR GLOBAL POSITIONING SYSTEM 255
fied for special Category I DGPS users because the 95% probability horizontal
navigation error budget is so large-33.54 meters, which is easy to achieve
using DGPS techniques.
Time Accuracy The time error budget with respect to universal coordinated
time (UTC) for the MAGR is a GPS system level specification for time trans-
fer [52]. The control segment's budget for maintaining the difference between
GPS time and UTC is 90 nsec, 1-sigma. The 1-sigma GPS time error due to
pseudorange error is
UERE 6.98
a~::. 1 =TDOPx = 1.12x --
c c
= 26.1 ns (5.86)
assuming a TDOP of 1.12 with an elevation mask angle of 5°. The MAGR
is required to output time via a one pulse per second (1- PPS) accurate to the
specified I 00 nsec, !-sigma. The remaining error budget to achieve this is
2
33 33
(
1.12 X ~ )
+ (90 X lQ- 9 )2 = 153.6 ns (5.88)
RELATIVE
10 -
(/)
0:::
w
1-
w
:2
I
u>-
~
:::> 0.1
u
u<(
ABSOLUTE
0.01 1-- -
I I I
0.001 I I I I I
cies are indicative of what might be achieved if differential carrier phase tech-
niques are developed for precision landing applications. Otherwise, accuracies
indicated for differential GPS are more applicable.
than those of GPS that cause the GPS satellites to repeat their ground track
each day. Because of this, the GLONASS ground tracks precess around the
Earth and repeat every 17 orbits lasting 8 whole days plus 32.56 minutes. This
is equivalent to 16 GPS orbits.
That is, when all 24 satellites were to be in the constellation, each would have
a frequency assigned to it with 562.5 kHz separation between satellite signals.
Similarly the L2 P-Code signal is transmitted at a frequency assigned to the
satellite. Each satellite is assigned a unique frequency according to the follow-
ing equation [82]:
That is, when all 24 satellites are in the constellation, each will have a frequency
assigned to it with 437.5 kHz separation between satellite signals. Note that the
ratio of the Ll and L2 frequencies is 9/7, including the frequency separations.
Also note that they are an integer multiple of a common frequency of 62.5 kHz.
GLOBAL ORBITING NAVIGATION SATELLITE SYSTEM (GLONASS) 259
The shift down in frequency is to avoid interference from future Mobile Satellite
Services (MSS) terminals.
Signal Modulation The L I and L2 signals arc both bi-phasc modulated with
the PRN codes and navigation data. The PRN code and navigation data char-
acteristics arc as follows:
5
9
9-BIT SHIFT REGISTER
7
511 kHz
C/A CLOCK
C/ACOOE OUT
1.953 msec of GLONASS system time. The resultant signal spectrum is a line
spectrum centered at the assigned satellite frequency with an envelope equal
to that given in Equation 5.54 with a T, of I /511 ,000 sec, where the lines are
spaced 1 kHz apart, and a spectral null occurs at multiples of 511 kHz. Since
the assigned satellite frequencies are spaced only 562.5 kHz apart, there is an
overlap of signal spectra.
Even though the spectra of the C/ A codes of the different satellites over-
lap, it has very little effect on signal acquisition and tracking, because the user
receivers, when correlating with the code, will track the correct carrier. Spectral
interference will occur, but will be well below the thermal noise level. Adjacent
frequency numbered satellite signals will have a cross-correlation level not to
exceed 48 dB [79]. Because of the separation in frequency, even if full code
correlation between signals occurred for an instant, postcorrelation integration
reduces the effect to that level. The C/ A code only appears on the Ll signals
[77, 80].
2. P code. The CIS has never published the GLONASS P code. However,
it has been determined independently [80]. The GLONASS P code is com-
prised of a 25-stage shift register with tap feedback that would produce a
33,554,431-bit maximal-length sequence, except for the fact that it is short-
cycled to 5, II 0,000 bits and reset to all I 's. It is clocked at a rate of 5.11 MHz
so that it repeats once per second. A functional block diagram of the P code
generator is shown in Figure 5.32. Every satellite generates the same P code.
The !-sec code epochs are synchronized to the 1-msec C/ A code epochs to
ease the handover from one code to the other.
The resultant P code signal spectrum is a line spectrum centered at the
assigned satellite frequency with an envelope equal to that given in Equation
5.53 with aT, of I /5, II 0, 000 sec, where the lines are spaced I Hz apart, which
5.11 MHz
P CLOCK
C/A CODE
HAND OVER
Signal Power The GLONASS ICD for the L I C/ A code signal indicates a
minimum received power of - 161 dEW, which is 1 dB less than specified for
the GPS L I C/ A code [79], although this level has been updated to -160 dBW
[82], which may be the total received C/ A code plus P code power. The L I
P code signal level is not published. The L2 P code received signal power is
-166 dEW [82].
Ephemeris Data The GLONASS ephemeris data are broadcast as ECEF carte-
sian coordinates in position and velocity with lunar/solar acceleration perturba-
tion parameters that are valid over about 0.5 hour [84]. The assumption is that
the user integrates via a fourth-order Runge-Kutta method the motion equations
that include the second zonal geopotential harmonic coefficient. Details of these
equations are given in the GLONASS ICD [79].
Almanacs Even though the ephemeris data differs completely from that of
GPS, the almanac parameters are quite similar as modified Keplerian parame-
ters.
262 SATELLITE RADIO NAVIGATION
Clock Corrections Clock correction parameters are also similar to that of GPS
in terms of clock offsets and clock drift.
l. Multisensor system implies that the GNSS and any augmentations can be
used for navigation, but only after it has been compared for integrity with
another approved navigation system in the aircraft.
2. Supplemental system implies that the GNSS and any augmentations can
be used alone without comparison to another approved navigation sys-
tem. However, another approved navigation system must be on board the
aircraft and usable when the GNSS is not available.
3. Required navigation performance (RNP) system is one that meets all the
requirements without need for any other navigation equipment on board
the aircraft. An RNP system may include one or more navigation sensors
in its definition (e.g., GPS with an inertial reference system, IRS).
GNSS does not add much to the aircraft's navigation system if it is only cer-
tified as a multisensor system service. It can add accuracy as long as the system
it is being compared with meets RNP requirements. This service is also useful
for test purposes. GNSS can, by itself, be certified as a supplemental system
through the use of receiver autonomous integrity monitoring (RAIM) [40] and
possibly oceanic en route [69, 89J. In 1996, based on FAA requirements, as
an RNP system, GNSS requires augmentation, such as combining GPS either
with GLONASS, an independent WAAS-type system, and pseudolites or with
another type of sensor, such as an IRS [69, 891.
to verify their position fixes since antiquity. The thing that is different now is
that computer technology has made it possible to use relatively sophisticated
mathematical methods in performing the consistency checks. Before getting into
the details of one RAIM method, it should be noted that catastrophic failures
are easy to detect with primitive methods, so they are not discussed further
here. It is the more subtle or incipient failures that are treated here; those where
a somewhat out-of-tolerance signal in space causes the user position error to
wander outside some specified limit for the phase of flight in progress. If the
GNSS is a supplemental system, it is sufficient for RAIM to simply detect the
failure and sound an alarm accordingly. If the GNSS is an RNP system, it is
necessary for RAIM to both detect and isolate and exclude the failed source.
This added burden of isolation and exclusion complicates the RAIM (now called
fault detection and exclusion, FDE) problem considerably [90, 91, 135].
RAIM Basics For tutorial purposes it is useful to begin with a simple two-
dimensional example. Suppose that we have three range measurements, each
defining a line of position (LOP) in a plane. Three possible situations are shown
in Figure 5.33. In Figure 5.33a we have the usual situation with good geometry
and consistent measurements (at least within the expected measurement noise).
The result is three intersections (fixes) that are close together. The observer
would then conclude ''no failure" in this situation. In Figure 5.33h we see
another possible situation in which we have favorable geometry but the fixes
are relatively far apart. The observer must conclude that something is wrong
here, and the decision is "failure.'' Note, though, that the information is insuf-
ficient to tell us which measurement is at fault that is, we can do simple error
detection here, but we cannot do fault isolation with just one redundant mea-
surement. Finally, in Figure 5.33c we see an extreme case of poor geometry.
Two of the LOPs are parallel. We can conclude here that measurements I and
2 are consistent (their LOPs are close together), but there is no valid check on
measurement 3. An error in it would go unnoticed. Thus the decision as to a
possible failure in this case is inconclusive. The observer must simply say, "No
(a) Consistent Measurements (b) lnc,:msistent Measurements (c) Poor Geometry Situation
valid integrity check is possible because of poor geometry." All three situations
shown in Figure 5.33 have their counterparts in the more complex GNSS RAIM
setting. Of course the meanings of "close together" and "far apart" need to be
quantified. Also, statistical performance criteria relative to the reliability of the
observer's decision need to be developed. More will be said of these items later.
Work on autonomous means of GNSS failure detection began in earnest dur-
ing the latter part of the 1980s. It was also during this period that the acronym
RAIM (for receiver autonomous integrity monitoring) was coined. We will not
attempt to document all of the technical papers on RAIM that appeared during
this period. One has only to browse through the proceedings of the meetings of
the Institute of Navigation to assess the degree of activity that took place dur-
ing this period and on into the 1990s. A summary of three different methods
is given in [92].
One of many RAIM schemes will now be described; it is easily understood
and can be thought of as a baseline or reference method. While it is a good
scheme, there is no claim that it is the best.
SSE::;; TH (5.94)
SSE> TH (5.95)
decide "failure."
The intuitive rationale for this rule is simply that if the measurements are
consistent, we can expect the residuals to be small; on the other hand. if the
measurements are inconsistent. we can expect SSE to be large because of a poor
least-squares fit. Once the threshold value THis set, the decision rule is quanti-
fied. With this RAIM algorithm, it is easy to set the threshold to yield an alarm
rate that is independent of geometry in the absence of a satellite malfunction.
266 SATELLITE RADIO NAVIGATION
This is usually set at the maximum allowable rate. The RAIM algorithm then
accepts whatever detection probability that results from this threshold setting.
RAIM Specifications There are four key parameters that must be included in
the RAIM specifications:
1. Alarm limit (also called alert limit). Alarm limit refers to the maximum
allowable radial error before the alarm is sounded.
2. Time response of the alarm (i.e., delay to alarm time). Too much delay
can be disastrous in critical situations.
3. Maximum allowable alarm rate in the absence of a satellite malfunction.
There must be a limit to nuisance alarms.
4. Detection probability. This must be close to unity if the RAIM algorithm
is to be effective.
These specifications compete with each other to some extent. For example,
tightening the false alarm rate specification makes it more difficult to meet
the detection probability requirement. The "elastic" in the system that makes it
possible to meet all of the stated four requirements is availability. The RAIM
algorithm can (and indeed must) reject poor detection geometries. RAIM avail-
ability, of course, suffers from such rejections.
RAIM for an RNP System In the RNP application, RAIM will have to be aug-
mented with additional measurement information from outside GPS. Many such
possibilities exist, such as using the combination of GPS and GLONASS, and
the marketplace (and perhaps politics) will determine the mix of sensor infor-
mation to be used in any particular application. Also, it is likely that the ultimate
GPS integrity protection will be provided by a combined W AAS/RAIM sys-
GNSS INTEGRITY AND AVAILABILITY 267
tern. The two systems are complementary, and there is much to be gained by
having the two systems work together [95]. (For a discussion of the W AAS,
see Section 5.7.3.)
GPSjGLONASS RAIM The fact that the two systems operate on different
time scales can be solved by the avionics user by simply adding the time dif-
ference to his solution state vector. This does, however, require an additional
satellite signal source because it adds another unknown to the solution. Further-
more, at least two satellites are required from both systems in the solution. If
there is only one satellite, any error in its pseudorange will simply be assigned
to the solution for the time difference based upon the position solution deter-
mined from the other system. However, this requirement for an additional satel-
lite would not exist continuously, since the time scales of both systems are quite
stable and a reliable time difference solution would remain valid over a long
time. Continuous monitoring may be required for solution integrity, however,
as well as detecting interfrequency errors in the receiver.
The problem with operating with two coordinate systems can be solved over
a period of time and updated as necessary with data-base parameters. For most
phases of flight, the differences appear to be small enough, so they do not mat-
ter [85 J. For precision approach applications, the use of differential corrections
would cancel the differences, including the differences in the time scales.
268 SATELLITE RADIO NAVIGATION
The combined system can also be used in conjunction with the W AAS
described in Section 5. 7 .3, in which case the differences can be broadcast via
the W AAS. If there are a number of failures in either system, the availability
of RAIM and continuity of service could also suffer because, as stated at the
beginning of this section, the two systems' orbits are not synchronized. That is,
if the GLONASS system were to augment the GPS system on one day, it may
not on the next because the ground tracks of the satellites moved with respect
to each other. This could be a problem if the GLONASS system continues to be
unreliable and the number of satellites in orbit do not maintain an operational
status.
W IDE-AREA
REFERENCE
STATION
(WRS)
able for wide area augmentation at the L I frequency were in 1996 [87]. Four
satellites are planned with an edge-of-Earth coverage shown in Figure 5.35 .
Note that many areas have double coverage, while some areas (e.g. , Europe)
have triple coverage. However, at least double coverage is required everywhere
in the service volume to provide the required RNP availability and reliability
[29, 30]. Thus, for CONUS, at least one or two additional GEOs are required
[30] . Unlike the Inmarsat-3 communications satellites, these additional GEOs
may be small single-mission navigation satellites [I 00].
In theW AAS concept, a network of monitoring stations (wide area reference
stations, WRSs) continuously track the GPS (and GEO) satellites and relay the
tracking information to a central processing facility [69, 98]. The central pro-
cessing facility (wide area master station, WMS), in turn, determines the health
and W ADGPS corrections for each signal in space and relays this information,
via the broadcast messages, to the ground Earth stations (GESs) for uplink to
270 SATELLITE RADIO NAVIGATION
Figure 5.35 lnmarsat-3 four ocean-region dep loyment showing 5° elevation contours.
the GEOs. The WMS also determines and relays the GEO ephemeris and clock
state messages to the GEOs. The signal is converted to the L I frequency on
the GEO satellite is then broadcast to the avionics user by the GEO satellite.
Message Format and Content [I 02, 1291 The integrity message contains the
status of each GPS satellite as use/ don '1 use information as well as WADGPS
GNSS INTEGRITY AND AVAILABILITY 271
error corrections and GEO ephemeris and clock data. The messages set into a
format with the capability to include both GPS and GLONASS data, although
GLONASS data are not planned in the FAA system. The magnitude of the
W ADGPS corrections can be also be used as error statistics for the users that
are not applying the corrections in the appropriate phases of flight.
The message data rate differs from that of the standard GPS signal [I 02,
I 03, 104, 129]. It has a symbol rate of 500 symbols per second. A rate 1/2
forward error-correcting (FEC) convolutional code of length seven is used to
reduce the effective data rate to 250 bps, but allowing a 5-dB gain in effective
energy-to-hit ratio over an uncoded 250 bps transmission [105]. Each message
block, shown in Figure 5.36, contains 250 bits, lasting one second. Each frame
contains 24 bits of parity to provide a strong burst error detection capability
as required for high integrity. The higher data rate provides two capabilities.
The first is a required capability to provide an integrity alarm to within 5.2
sec of a signal-in-space fault during Category I precision approach [98]. The
second is to broadcast W ADGPS corrections at a rate commensurate with SA
and ionospheric delay errors.
The various message types are listed in Table 5.9 [129]. There are two types
of correction data-fast and slow. The types 2 through 5 fast corrections are
intended to correct for rapidly changing errors such as GPS SA clock errors,
while the slow corrections are for slower changing errors due to the atmo-
spheric and long-term satellite clock and ephemeris errors. The fast GPS clock
errors are common to all users and will be broadcast as such. Corrections des-
ignated with the maximum positive number indicate not-monitored satellites,
while those designated with the maximum amplitude negative numbers indi-
cate don't-use satellites, which is the integrity indication. Procedures for using
these messages are given in the RTCA MOPS [104, 129].
For the slower corrections, the users are provided with ephemeris and clock
error estimates for each satellite in view (message types 24 and 25). Users are
separately provided with a wide area ionospheric delay model and sufficient
real-time data to evaluate the ionospheric delays for each satellite using that
model (message types 18 and 26). This model is comprised of vertical iono-
spheric delays at a set of grid of points that a user can interpolate to the iono-
spheric pierce points of his pseudorange observations.
~-------------------250BITS-1SECONn-----------------------~l
I I I 212-BIT DATA FIELD
24-BITS
1 PARITY.
I
Type Contents
0 Do not use this GEO for anything (for W AAS testing)
PRN Mask assignments, set up to 52 of 210 bits
2-5 Fast corrections
6 Integrity information
7 UDRE acceleration information
8 Estimated standard deviation message
9 GEO navigation message, (X, Y, Z, time, etc.)
10-11 Reserved for future messages
12 WAAS network/UTC offset parameters
13-16 Reserved for futu,:e messages
17 GEO satellite almanacs
18 Ionospheric grid point masks
19-23 Reserved for future messages
24 Mixed fast corrections/long-term satellite error corrections
25 Long-term satellite error corrections
26 Ionospheric delay corrections
27 WAAS service message
28-62 Reserved for future messages
63 Null message
0.985 -J-------- --+/: ·--·- ············ ·········· --·- -- · · · · ·-·-----· · · · -· · . ····-· 24 GPS + 4 GEO
• I
0.98 / •: .
...........•...... ······r············ ·······-- .. ....... ___ ···············- SATELLITES
~ 0.975 of·····-·-·················------
'
···-··-!········-········-·---·- ·+··········
...J
m 0.97 ---,. ...i... ________ -. l-- ........ ---- ....... --------!- .........___ ..__________ . . J.---------
<(
...J
-----,1: -· ··---···-~----·--- ~
1E+01
~ 0.965 --1~ .... -.. + . . . · · · 1-- ---1 ----t ~ 1E..01 ~
r-..1""---. I· 24GPS
<( :::i 1E..03 1.:'""
0.96
0.955 +···--····· . .
'
··ll···················
-lj- ··-
.... j
' Kt
0:
::> 1E-11
. .J-!J
0 1 2 3 4 5 6 7 8 9 10 11
0.945 .. ·············-····l··--···f···--··········-·····-~---·'··--·---1
I
NUMBER OF SIMULTANEOUS SV FAILURES
11 lr I 1 I II
0.5 1.5 2 2.5 3 3.5 4
~
GNSSSPACESEGMENT
DIFFERENTIAL CORRECTIONS
REFERENCE STATION
Figure 5.38 Pseudolite DGPS concept.
I. They provide a data link for local DGPS corrections and integrity infor-
mation that can be received with a slightly modified GPS avionics receiver
[107, 108].
2. They provide additional ranging signals, just as the W AAS, resulting in
a significant VDOP availability enhancement [109, 110, 111, 112, 113].
This VDOP enhancement is illustrated in Figure 5.39 for two PLs aug-
menting GPS at the FAA Technical Center in Atlantic City, New Jersey.
VDOP is reduced from 2.3 for GPS only to less than 1 in the area of
the runways for GPS augmented with two PLs [113]. Similar enhance-
ments are shown in [ll2] for runways at O'Hare International Airport in
Chicago.
3. They provide a rapid change in geometry that is extremely important for
276 SATELLITE RADIO NAVIGATION
1.7
I
:::.:::
g
(.!)
0
g< -5.--_ __
M
i
1:1::: -10
-10 -5 0 5 10 15
RW13 CROSSTRACK (KM)
Figure 5.39 Illustration of VDOP reduction with 2 PLs at the FAA Technical Center.
kinematic carrier phase ambiguity (both integer and floating point) resolu-
tion techniques [73] and multipath mitigation. This approach is described
in Section 5.5.9.
4. Their signals in space can be more accurate than satellite signals because
of the nonpresence of ephemeris and ionospheric delay errors and reduced
tropospheric delay errors [Ill I].
Along with these advantages there are two significant disadvantages that require
attention: I. The proximity of the PL to the avionics user causes potential inter-
ference to the reception of satellite signals-the well-known near/far problem
[107, 108, 110, 114]. The PL signals become quite strong when the avionics
receiver is near the PL if the PL power is set for reception at a distance. For
example, the received PL signal power increases inversely with the square of
distance from the PL. For a precision approach and landing applications, the
PL power would be set for reception at about 20 nmi. Thus, at 0.1 mile, the
received signal is 40,000 times stronger (46 dB). 2. A PL's location on the
ground could present a problem with the antenna location on the aircraft for
simultaneous reception of the pseudo lite and the satellites [ 113, 115].
prevent the capture of a receiver due to the excessive dynamic range required
for a close-in signal [108]. The pulsing with a relatively low duty cycle reduces
the interference to any signal in the reception band. The strong signal simply
punches holes in the lower powered signals, reducing their received power by
only the duty cycle percentage. For example, if the duty cycle is set at 0.1, the
signal loss of the satellite signal is only I 0 log 10 0.9 = 0.458 dB, although an
additional loss in CjN0 is realized because the pulse power also enters the corre-
lator. This pulse power is reduced significantly by precorrelation pulse clipping
or pulse suppression, techniques that are already implemented in GPS receivers.
Pulsing is also required to prevent multiple PLs interfering with each other. This
adds a requirement of pulse timing so that pulses from two different PLs are
not received simultaneously, in addition to the requirement that received pulses
must be asynchronous with the reception of GPS data bit edges. This type of
timing is possible [113, 115].
Even though pulsing can reduce interference significantly, the strong signal
within the clipped pulses can still cross-correlate with the received satellite sig-
nals, if indeed the PL signals carry GPS-like Cj A codes [ 108, II 0, 113, 114,
115, 116]. This causes another dynamic range problem because the cross-cor-
relation margin between C/ A codes is only on the order of 22 to 24 dB [ 18,
32]. The I 0 dB reduction in average power of the pseudolite signal due to a
10% duty cycle pulsing does little to prevent cross-correlation [ 114 ]. This is
where the frequency offset helps. It has been shown that the cross-correlation
can also occur at frequency offsets but is proportional to the 1 kHz spectral
line component levels of another CjA code in the same family [18, 32]. As it
happens, however, the spectral line components near the null of the spectrum
are down on the order of 70 to 80 dB, as can be observed in Figure 5.15. Thus,
if the PL were to transmit in the null of the GPS satellite C/ A code spectra,
which are all at the same frequency to within 5 kHz due to Doppler differences,
the cross-correlation would be insignificant [ 113, 115]. Cross-correlation peaks
within 4 dB of the 0 offset case can still occur [ 116]. However, if the carrier
frequency jcode frequency ratio of 1540 is maintained, these peaks disappear
very rapidly and simply create an interfering noise that the pulsing mitigates.
Test results back up these interference mitigation theories [ 113 ]. Figure 5.40
20 40 60 80
Average PL to GPS Power (dB)
Figure 5.40 Effects of signal structure on PL interference to GPS satellite signal
reception.
278 SATELLITE RADIO NAVIGATION
PROBLEMS
5.1. A GPS user receiver is tracking four GPS satellites-PRNs 1, 13, 19, and
22. Via the reception of navigation data, the receiver receives the following
ephemeris parameters:
(a) Parameters common to all four satellites:
vA = 5153.619629 meters
e=O
i0 = 0.3055555556 semicircles
wo = 0 semicircles
toe = 345, 600 seconds in the GPS week
280 SATELLITE RADIO NAVIGATION
X 1 = 13,672.46475 km
Y 1 = -6,720.41440 km
Z 1 = 21,755.97535 km
x, 3 = --2,370.46666 km
Y 13 = --23,498.04734 km
zl3 = -12,150.94171 km
X ,9 = - 18,962.99343 km
Y 19 = 6, 971.55345 km
Z 19 = 17,240.21601 km
X 22 = -10,899.89991 km
Y 22 = -14,301.92165 km
z22 = 19,546.60953 km
X= -1, 268.4451896 km
Y = -4,739.4160255 km
Z = 4, 078.0482708 km
R1 = 23,230.69260 km
R 13 = 24,829.03623 km
R19= 24,969.68770 krn
Rn = 20, 578.69009 km
HOOP = 1.241
VOOP = 1.631
POOP = 2.050
TOOP = 0.823
GOOP = 2.208
5.3. After traveling some distance the user then measures the following actual
ranges to the four satellites:
What is the user's new position in ECEF coordinates? What are the
azimuth and elevation angles to the satellites and the HOOP, VOOP, POOP,
TOOP, and GOOP at this new position, which has an approximate position
of N44° 58", W93° 15", at an altitude of 200 meters? Assume a spherical
earth with a radius of 6378.163 kilometers.
Ans.: The user's position solution (near Minneapolis) is as follows:
X = - 255.843602 km
Y = - 4,505.54881 km
Z = 4, 507.55905 km
282 SATELLITE RADIO NAVIGATION
HDOP = 1.304
VDOP = 1.603
PDOP = 2.066
TDOP = 0.8432
GDOP = 2.231
6 Terrestrial Integrated Radio
Communication-Navigation
Systems
6.1 INTRODUCTION
Since the 1970s, many radio communication and navigation systems have
used the same portion of the frequency spectrum and common technology,
such as time synchronous operation, digital modulation, spread spectrum wave
forms, coding and user-borne clock oscillators. Synchronous operation, in con-
junction with signal time-of-arrival measurement, has lead to a direct method
for measuring the range between transmitter and receiver locations in sys-
tems using this technology. For these reasons, integrated relative and absolute
communication-navigation systems, which provide both digital communication
and navigation functions by means of the same wave form, have been widely
developed. These systems typically use the content of digital data and the time
of anival of the messages measured by the receiver, to determine the receiver
platform's position, through some form of multilateration. In general, the posi-
tions are determined in a relative sense within an arbitrary grid, although the
unit positions can be referenced to an absolute, geodetic coordinate system,
such as latitude, longitude, and altitude, through the use of reference stations
whose positions are independently known in the absolute coordinate system. In
addition, the position data may be combined in a Kalman filter with dead-reck-
oning sensor data, such as from an inertial platform, for the purpose of position
extrapolation and calibration of the dead-reckoning sensor enors.
Several types of tenestrial integrated communication-navigation systems
have been developed. One is a decentralized system, in which the operation
is not dependent on any central site or node, and each user in a community of
members determines its own position. Such a system is also called nodeless. A
second type is a centralized system, wherein the operation is dependent on a
central site (node) and may be controlled by it and wherein the determination
of the positions of the users in the community is performed by that central site.
Frequently, it is desired to have the positions of a large number of users known
and tracked at the central site, such as in military or civil command and con-
trol systems. Typically, in such a system, users may obtain their positions by
automatic, periodic, or occasional requests from the central node; hence such
a system is considered nodal. Sys.tems are being developed that exhibit both
nodal and nodeless characteristics and thus become hybrid systems. However,
the fundamental design of these systems is typically based on either the decen-
tralized or centralized concepts.
Typical examples of these systems are represented by the Joint Tactical Infor-
mation Distribution System Relative Navigation (JTIDS RelNav) function and
the position location reporting system (PLRS) and its enhanced versions, whose
principles of operation are described in this chapter. The former is representa-
tive of a decentralized system and the latter is representative of a centralized
system. Applications of these types of systems cover a wide spectrum, includ-
ing the handover of targets between units operating within a common grid,
rendezvous of aircraft or other units, command and control from the viewpoint
of a military commander having knowledge of the position of his forces, and
such specialized purposes as search and rescue and medical evacuation.
The systems described in this chapter were mature and operational in 1996.
For example, by 1996, over 3600 PLRS and enhanced PLRS user units had been
produced and deployed on a variety of U.S. Army and Marine Corps vehicles,
including tanks and helicopters and as manpack units, and about 1500 more
were planned for the future. By 1996, about 500 airborne JTIDS terminals had
been installed on such aircraft as the U.S. E3A, E2C, B-1, F-14D, F-ISC, and
JSTARS, as well as on several aircraft of other NATO countries. About 400
more such terminals had been planned for later installation in various military
aircraft and ships. Also in 1996, a major development was under way by a
consortium of several countries for a JTIDS-like smaller and modular MIDS
terminal that includes the relative navigation function. This reduction in termi-
nal size will make it possible to install it in a large variety of other aircraft.
• II IITfiTI1lffi IITITIITfl1
___~~,~,~,~.~~.~.~~,___ _,.~~,~,~,~I.LI.L__.~~.Willillill~~,~,~,~~--JIWUI.II.I.II.u.u.-
32 chips/pulse encoding five
information bits, continuous phase
shift modulated, pseudo-randomly
encoded
Figure 6.1 JTlDS signal structure.
clock, very precise synchronization of all the individual clocks, although not
necessary for the communications function of JTIDS, allows the receiver to
convert message TOA to an accurate one-way radio range to each of the trans-
mitters and thus support a precise multilateration navigation function.
Interrogator t (o)
TOA(i)
err (i)
Interrogator
Source
Source t (o)
messages received from users of superior quality and establishes its own rank in
the hierarchy on the basis of the qualities of its sources and the measurement
geometry as reflected in the RelNav Kalman filter covariance matrix. Some
units in the community must, of course, have independently known positions
and qualities to get things started.
In the time hierarchy, the net time reference (NTR) transmits the highest
time quality ( 15). Other units transmit time qualities derived from the time vari-
ance developed in their synchronization Kalman filters. Primary users employ
only the RTT technique, while secondary users employ primarily passive syn-
chronization, making recourse to RTT only under certain conditions of poor
geometry. Explicit designation of secondary users is seldom made. Radio-silent
users, unable to participate in RTT message exchanges, automatically assume
secondary user status.
Within the RelNav hierarchy, the equivalent of the NTR is the RelGrid
navigation controller (NC). RelGrid coordinates transmitted by the NC are by
definition perfect as indicated by transmission of the highest relative position
quality of 15. There is no equivalent of the NC in the geodetic frame; that
is, there is no equivalent arbitrarily perfect reference designator. Designation
as a position reference disables geodetic position update; however, this desig-
nation does not connote perfection. The accuracy of the geodetic position in
PPLT messages is characterized by the geodetic position quality which, for all
units, is determined initially by operator entry and subsequently by the RelNav
Kalman filter covariance. The maintenance of the covariance based hierarchy
290 COMMUNICATION-NAVIGATION SYSTEMS
in both time and position is the function of the source selection logic described
in Section 6.2.6.
6.2.6 Mechanization
The overall diagram of the RelNav function in Figure 6.3 shows its three major
subfunctions: source selection, Kalman filter and navigation processing. Each
received PPLI message is processed by the source selection function immedi-
ately following its reception. Host dead-reckoner data are processed in the nav-
igation function to provide an estimated own-unit position at the time of receipt
of each selected PPLI message. Source selection stores the selected PPLI obser-
vations and the associated navigation data to await processing by the RelNav
Kalman filter. The computed range and direction to the source and the mea-
sured range from the TOA of the received PPLI message provide a range error
vector which the filter uses to estimate position, velocity, and other dead-reck-
oner error states. These error estimates are applied to the internal dead-reck-
oner model in the navigation function and corrections are supplied to the host
platform.
Estimated Present Position
~
N
~
292 COMMUNICATION-NAVIGATION SYSTEMS
I Communications
Management Function
I PPLI
IMessage -----4_ Minimum Range Test I
messages with those of the receiving user. In general, both the time and position
qualities of the source must exceed the position quality of the receiving user.
This test is, in a sense, antithetical to the concept of the Kalman filter which
was designed to derive low-variancae estimates from higher-variance observa-
tions; however, it has been found to be essential to the operation of the network
of interactive filters in a RelNav community. It recognizes that the simple qual-
ities transmitted in the PPLI messages do not represent true variances (noises)
in the estimation sense; rather, they represent uncertainties of estimation errors
that are predominantly correlated bias errors. The limited PPLI message bits
available to the RelNav function precluded the more sophisticated approach
of transmitting the separate position terms of the RelNav filter covariance
matrix.
Users enter the network with an operator-assigned initial position quality
indicating the uncertainty of the initial position estimate. In airborne users, posi-
tion quality will degrade with time in accordance with the error signature of the
user's dead reckoner as modeled in the Kalman filter's time propagation of the
covariance matrix. The Kalman filter processes PPLI observations to estimate
the dead-reckoner errors (and reduce their covariance terms) and this will be
reflected in a decreased rate of degradation of position quality. A dynamic bal-
ance supportable by the quality and geometry available from the PPLI message
sources is soon established.
The geometric rank tests of Figure 6.4 recognize that the value of a PPLI
observation is related not only to its position and time qualities but also to the
direction from which it was received. Several sources of equally high quality
all in approximately the same direction provide little more information than
one such source, but each will consume source selection buffer locations and
filter processing time and will crowd out observations of lesser quality but of
greater value by virtue of their directions. The rank test uses the orientation,
eccentricity, and semimjaor axis of the bivariate error ellipse defined by the
horizontal position terms of the RelNav Kalman filter covariance matrix and
the quality-based variance and direction of a received observation to compute a
numerical rank that is stored with each observation. The rank is the approximate
variance of a hypothetical observation lying directly on the extended major axis
that would provide the same benefit as the received observation at its angle off-
axis. In this context, benefit implies the reduction in the major axis of the error
ellipse to be expected of Kalman filter processing of the received observation.
Kalman Filter The JTIDS RelNav Kalman filter is an extended Kalman filter
that estimates linear error states of the navigation process. Table 6.1 presents
a typical state vector for use with an inertial dead reckoner. Two of the state-
vector elements (7, 16) are time and frequency states used only by passive
users. For active users, these two states are carried in a separate synchronization
filter that uses RTT or ETR data. Five states (1, 2, 6, 14, 15) are relative grid
states. Of the remaining nine states, the two horizontal geodetic position error
terms are carried as quaternions, while the third dimension is carried as altitude
294 COMMUNICATION-NAVIGATION SYSTEMS
error and an altimeter instrument scale factor term. The inertial filter models
velocity errors in the north and west directions in the geodetic frame and the
three platform misalignment or tilt states in the local-level frame. A filter for use
with air-data computer/attitude and heading reference system (ADC/ AHRS)
inputs is also included and differs in that it models the dead-reckoner errors
as two wind components, an airspeed instrument scale factor, and azimuth bias
and azimuth gyro drift rate errors.
Figure 6.5 is a flow diagram of the Kalman filter function. Only the mea-
surement geometry, measurement innovation, measurement validity, and filter
characterization features that are peculiar to JTIDS RelNav will be discussed
here. See Chapter 3 for a general discussion of Kalman filters used in multi-
sensor navigation systems.
The source position and the own-unit position stored with it are subtracted
vectorially to obtain a predicted range and three-dimensional direction cosines.
The range error (measurement innovation) is obtained by subtracting the pre-
dicted range from the measured range. The measurement variance derived from
the source qualities and the direction cosines are supplied to the Kalman gain
function, and the range error is supplied to the state and covariance update
function.
The validity tests are intended to detect divergence of the filter and to pro-
vide protection against inconsistent observations. This is particularly impor-
tant to JTIDS RelNav, as compared to other radio-ranging systems because the
ranging sources are typically other, sometimes erroneous, navigating users. The
dilemma facing any validity test is whether the error lies with the local esti-
mate or with the input data. The source selection function makes this decision
more difficult by narrowing the group of sources to the few, typically three to
five, with the best announced qualities (whether true or false). An observation
is rejected if the measurement innovation (the computed range error) is large
compared to the receiving user's filter variances (3-sigma reasonableness test).
A series of observations exhibiting an average innovation exceeding 2-sigma
for this test, triggers a proportional increase of the position terms of the filter
covariance matrix. Eventually, if the recurrence rate of these covariance expan-
sions exceeds a threshold, the process is abandoned and the filter is restarted.
The filter characterization function serves only to convert filter covariance
data to forms more convenient to the source selection and PPLI generation
functions. Position quality in the PPLI message is defined by the JTJDS mes-
sage standard as representing the semimajor axis of the horizontal position error
ellipse and the rank computation requires the semimajor and semiminor axes
and the orientation angle of the error ellipse. Each Kalman cycle, the filter char-
acterization module computes these terms and stores them for use during the
next cycle.
w noise
covariance
gain Geometry ..._ PPLI
II Computation posiUon
Navigation
Correction Careotion & Meaeurement
to t-- Validity T eating
Covariance Predicted
Navigation Extrapolation
model measurement
Error state Measurement Measurement (TOA) ......:: ~ Source
Correctio
state
' transition
matrix
innovation
~
' of
Selection
Function
transition State Transition Updated error state & Error State & Measurement
matrix
- Matrix
Maintenance
covariance Covariance
Update
Innovation Innovation
Computation
N
\C
Ul
N
',c)
C'l
~ Misalignment correction
Platform an_g_ular....::.
rates Misalignment Angle Velocity correction
Update Kalman Filter
14--- Alti~ & au.aternion
J Misalignment
angles
corrections
v
data angular rates
Update
Host 1-<-- Grid
Navigation
System
Interface
IVelocity
'V
AI
V'/
correction
PPLI Message
Altitude
Geodetic Quatemion
Update
1-- 1
Rei Grid
~
Generation
Navigation
1 Quaternions & altitude . .. jUpdate
Geodetic pOSition
Direction Cosine and
Position Geodetic Update Geodetic & Qrid ..:::. Source
Veloaty .& Tilt pOSitiOn &. velOCity Selection
corrections
is selected to correspond with that of the host system. JTIDS initial position
is obtained from the host inertial system, but thereafter JTIDS uses only the
velocity and baro-inertial altitude inputs to compute delta-velocities for use in
the internal inertial mechanization. This model continuously applies the Rel-
Nav filter corrections to velocity and acceleration (via platform misalignment
calibration) to improve upon the solution provided by the host inertial system.
Corrections to position, velocity, and tilts are also returned to the host data
system.
As mentioned earlier, an air-data model is available for use in aircraft having
only an ADC/ AHRS, but JTIDS RelNav has seldom been installed in such
aircraft. Its primary use is as a backup mode to continue RelNav operation and
PPLI transmissions should the host's inertial system fail and force a switch to
air data.
Integration with GPS JTIDS RelNav processes position fixes from an inter-
connected GPS receiver as two uncorrelated, one-dimensional Kalman updates
in the north and east directions. The variance data in the GPS input are used in
RelNav source selection and Kalman gain computations in the same manner as
the qualities in received PPLI messages. If GPS data are of high quality-that
is, if GPS GDOP and signal availability (Chapter 5) are within typical GPS sys-
tem performance criteria-the JTIDS terminal will tend to operate exclusively
on GPS data, and the transmitted PPLI messages will reflect GPS position accu-
racy. A few aircraft with interconnected JTIDS and GPS equipments can, via
JTIDS RelNav, extend the benefits of GPS-based position fixing to an entire
community of JTIDS users.
Dead-Reckoner Errors The inertial filter explicitly models velocity and plat-
form misalignment errors, but higher-order terms are modeled only as dynamic
process noise. The air-data filter models wind and azimuth errors as Markov
processes. To the extent that actual dead-reckoner errors depart from these
assumptions, JTIDS RelNav will experience errors. As one example, winds aloft
typically vary with altitude so a RelNav user coupled with an ADC/ AHRS can
be expected to exhibit temporarily increased errors following a substantial alti-
tude change.
Source Position Errors The best PPLI position quality (15) implies a site sur-
vey of better than 50 ft, !-sigma. More important than accurate survey, however,
are the position qualities assigned at the ground reference sites. These must
reflect a conservative estimate of the position uncertainty. Optimistic position
qualities can result in user validity failures and recurrent filter resets leading to
instability rather than just increased position error.
tor filterng. All positions are available to the cooperating users and to command
centers.
PLRS also provides short message data exchange for both manual and auto-
mated users. All control, measurement reporting, and data exchange are crypto-
graphically secured in a synchronous, anti-jam communications network. Mas-
ter stations (MSs) establish control circuits between radio sets (RSs) and the
MSs via a control network. This control network supports the position location,
navigation aid, and friendly unit identification. From a message flow standpoint,
the control network provides an "order-wire" capability that can be used for
data exchange between users. In addition to the order wire, the control network
also supports user access to a wide range of position location, navigation, and
identification information.
The system has a range of capabilities which support the conduct of coordi-
nated military operations. For the individual tactical user on foot, in a surface
vehicle, or in an aircraft, the system determines and displays to him his accurate
position in real time. It alerts him if he enters a restricted area. It also provides
the user with guidance to predesignated points, to other users, or along corri-
dors in accordance with requests, as well as providing a free text data exchange
capability.
For a tactical commander, the system provides the identification, location,
and movement of all cooperating users within an assigned area of responsibility.
In addition to allowing the commander to monitor the movement of forces,
the system also has the ability to input and modify control measures such as
coordination points, safe corridors, and restricted zones.
For all participants, the system (which operates beyond the line of sight
via integral relays) incorporates electronic counter-counter measures (ECCM)
and provides cryptographically secure digital data communications. Each user
has the capability of sending preassigned short messages to provide data to or
request information from the system and to exchange short free text messages.
A single synchronous community can support over 900 users with a varied
distribution of manpack, surface vehicle, and airborne platforms. System per-
formance is provided within the primary ground operating area, and airborne
users can be located and tracked within a 300-km square extended operating
area. It can interoperate with other communities in adjacent or overlapping geo-
graphical areas. It operates in the UHF band at frequencies from 420 to 450
MHz.
TOA LINK
PORT LINK ••
Figure 6.7
NODE (RADIO SET)
GRID REFERENCE NODE
a set of PORT links that connects RSs (nodes) to the MS either directly or
via one to three relay nodes. Network control and measurement reporting is
transferred over the PORT path. In addition to the bilateral PORT links, one-way
TOA links are utilized to provide the additional multilateration structure needed
for position location and tracking. Since the timing of RS clocks is established
using paired TOA data along the PORT paths (Section 6.3.5), one-way TOAs
can be converted to true range es.timates. In a typical PLRS deployment over
half of the range measurements are based on one-way TOAs.
To initialize the position locarion function and to maintain a relationship
between the internal coordinates and the external military grid reference system
(MGRS) coordinates, the MGRS positions of three or more cooperating RSs are
input to PLRS. These are normally input as three-dimensional fixed reference
positions, and the RSs then become grid reference nodes. The MS may be, but
is not necessarily, one of the grid reference nodes. In addition the system can
operate without any fixed reference RSs as long as the positions of three or
more RSs are regularly input to the position tracking function. In this latter
case the positions may be input and updated by RSs which are moving. This
is termed a dynamic baseline operation. External position sources such as the
global positioning system (GPS) (Chapter 5) can be used to provide the position
reference information to PLRS, but external data sources require the appropriate
coordinate conversion from the respective geoids and datums to MGRS.
and crypto key. Also, ranging may be accomplished with any of the bursts sent,
no matter what their purpose, as far as the message portion of the transmission
is concerned.
PLRS is a synchronous time division multiple access (TDMA) system, which
also employs frequency division multiple access and a srread spectrum wave
form (Figure 6.9). PLRS employs a network that is fully synchronous in three
respects: all RSs maintain timing such that cryptographic resynchronization is
seldom required; all RSs perform actions in a programmed cyclic manner such
that reprogramming RS assignments for relay. ranging, and reporting is seldom
required ; and each RS's time base is maintained with sufficient accuracy such
that one way time of arrival measurements can be translated to ranges by the
MS. Each of these aspects reduces the number of required control transmissions
and makes time available for other syslem functions or for increased system
capacity.
PLRS employs time division multiplexing to permit a large number of users
to utilize the same frequency. Each of the RSs in a network takes turns trans-
mitting its burst while other RSs listen. These timeslots are assigned by the
MS. based on the particular requirements of each user. For example. for a given
----TRANSMISSION BURST
=BOO )lSEC
Figure (l.9 PLRS time division multiple access tTDMA) organization.
304 COMMUNICATION-NAVIGATION SYSTEMS
Timeslot The fundamental time division is the timeslot. The timeslot length is
1.95 msec. The burst transmission accounts for 800 p.,sec, and 600 p.,sec is allo-
cated to RF propagation delay. The remaining time is required for processing
overhead such as message encoding, validation, and guard time.
The Code Resource In addition to the time and frequency separation, PLRS
uses a pseudonoise code resource. These codes provide a different spread spec-
trum pattern for each transmission in the network, thereby eliminating cross talk
and reducing interference.
6.3.5 Measurements
TOA Measurements One-way ranging is made possible by the fully syn-
chronous nature of the PLRS network. Each RS employs a set of time markers
POSITION LOCATION REPORTING SYSTEM 305
RANGE
(PROP)
DELAY
-t------~TRANSMISSSION ~--------1--
p"
-RS#1
I
I~ TOA2 .I RS#2 F"~/
LEGEND
~T= TOA2-TOA1
2 ~TRANSMISSION ~I - TIMING OFFSET
P- PROPAGATION TIME
~-RANGE 1 , 4 - - - - - - - - f
~ ®-RA~E2.4
L-~-T_A_T~_o_N~/~:: ~ ~ ~ : - - - dGE3,<
POSrTIONS
RS#1,
OJ 0 RANGE 1,4
RANGE 2,4
}
..-'\. POSITION OF RS #4
RS#2, RANGE 3,4 L..y"
&RS#3 ELEVATION RS#4
OOE REFERENCE RS
BAROMETRIC PRESSURE - - . - 1 1-4--- ALTITUDE REFERENCES
MEASUREMENT
• .,..__ _ _ _ BAROMETRIC MEAN SEA
OOE NON-REFERENCE R S - - - - 1 LEVEL (MSL) PRESSURE
BAROMETRIC PRESSURE ALTITUDE
MEASUREMENT FILTER UPDATED NON-REFERENCE
- ~ RSALTITUDEANDVERTICAL
......---- VELOCITY ESTIMATES
UPDATED NON-REFERENCE ----1,.,._.
RS ALTITUDE ESTIMATE POSITION REFERENCES
(FIXED)
SINGLE TOA OR [ TRACE
PAIR OF TOA ---------~ FILTER 1-----II•UPDATED 3-D POSITION
MEASUREMENTS AND VELOCITY ESTIMATES
FOR BOTH RSs
frequency to keep all RSs nominally synchronized with the MS, and thereby
with one another.
Finally, the purpose of the TRACE filter is to take each TOA measurement
and partially update the two cooperating RS's position and velocity estimates in
three-dimensional space. Without the supporting processing from the other three
filters, TRACE would not be able to operate successfully. Link value account-
ing is one of the unique byproducts of the TRACE processing. If a particular
TOA link assignment is not aiding the position location accuracy (due to geom-
etry or excess TOAs in a particular direction) of either RS then that TOA link
assignment is replaced by one that may be more beneficial to overall system
accuracy.
These filters provide the partially updated positions necessary to permit the
MS to report a fully updated posJition about half a second after MS receipt of
each user's measurement report.
Control Network Traffic Four types of messages are transmitted over the con-
trol network: user data input, user data output, network control, and measure-
ment report messages. All messages are sent to their appropriate destination,
through relays as necessary, without data content change. Error control is used
in order to insure a low(< 10- 5 ) message error rate. User data input messages
originate at any RS via an inputjoutput (I/0) device. An I/0 device may also
be used to originate queries requiring user data output messages to be sent back
to the RS. This establishes two-way communication between any RS and the
MS. Two way communications between any pair of RSs can also be established
using these message types.
Network control messages and measurement report messages are used by the
MSs to maintain communication with, and exercise control over, each RS. The
network control messages contain commands to the RS, (e.g., link assignments
and timing correction commands). Measurement report messages contain TOA
and altitude measurements and status information. The actual routing of mes-
sages to the proper destination is implicit in the network structure. The MS,
with its knowledge of the connectivity, assigns transmit and receive times to
provide proper relaying of messages to their destinations. RSs performing a
relay function make no distinction between different message types.
User-to-User Traffic The basic PLRS provides for user-to-user data exchange
using two distinct approaches. The first is by way of the MS, using the data
input and data output messages mentioned above. This approach can be used
by any RS operator to send short alerting or coordination messages to any other
user in the community. These messages are stored and forwarded by the MS(s).
The other approach provides local groups of up to eight users the ability to send
short messages to each other, without going through the MSs.
the PLRS control network for monitoring and controlling large communities
of user RSs. In addition to the positioning, position-reporting, navigation aid,
cryptographic key distribution, and status-reporting functions, the control net-
work is also utilized for distributing communications circuit assignments and
monitoring user-to-user communications performance. The EPLRS RSs support
communication network management by implementing the commands, moni-
toring and reporting circuit status, establishing new circuit paths, and controlling
the flow of data packets into and out of the network.
Both duplex (point to point) and group-addressed (broadcast) types of service
are available via the same user RS. Each RS can support up to 30 user circuits
(needlines) simultaneously with a composite (receive plus transmit) information
rate of 4 kbps. Each duplex circuit is capable of supporting acknowledged data
rates of up to 640 bps in each direction. Each group addressed circuit is capable
of supporting nonacknowledged data rates up to 1280 bps. The primary user
interface is via the Army Data Distribution System Interface (ADDSI), which
uses permanent virtual circuit protocols based on CCITT x.25. A !553B data
bus interface is also used for compatibility with existing aircraft and vehicular
systems.
The EPLRS concepts have been proven through live testing with 160 RSs.
In addition, extensive computer modeling and large-scale user community sim-
ulation have been used to confirm extension of performance to operations with
over 500 user RSs. EPLRS is interoperable with the basic PLRS allowing for
mutual support and coordinated operations between basic PLRS- and EPLRS-
equipped users.
PROBLEMS
6.2. The PLRS time of arrival signal processing splits a PN chip into 16 equal
parts. What is the precision of a one-way range measurement in meters:
Ans.: 3.7 meters.
6.3. If the clock offset error between two PLRS RSs involved in a range mea-
surement is 15 nsec, what error in one way range measurement does this
cause?
Ans.: 4.5
6.4. In Problem 6.2, what error would result if this were a two-way range mea-
surement?
Ans.: 19 meters.
6.5. In Problem 6.3, what error would result if this were a two-way range mea-
surement:
Ans.: None.
6.6. In the military grid reference system used for reporting positions to a user
readout in PLRS, a location is reported as 4 decimal digits of casting and
4 decimal digits of northing within a designated 100-km grid square. What
is the maximum error introduced due the precision of this report?
Ans.: 5J2 meters
7 Inertial Navigation
7.1 INTRODUCTION
I. The position and velocity information degrades with time. This is true
whether the vehicle is moving or stationary.
2. The equipment is expensive ($50,000 to $120,000 for the airborne sys-
tems in 1996).
3. Initial alignment is necessary. Alignment is simple on a stationary vehicle
at moderate latitudes, but it degrades at latitudes greater than 75° and on
moving vehicles.
4. The accuracy of navigation information is somewhat dependent on vehicle
maneuvers.
Platform wander a
angle, if not
north pointing
Vehicle azimuth - Vehicle azimuth
relative to platform + from north
---···· ---
Vehicle attitude angles
(roll and pitch)
Roll
I Azimuth Vehicle attitude
I 1 Pitch angle
Pitch gimbal readout
I
I Pitch
I I
L_~
Redundant gimbal Inner roll angle J
drive -------------
(Sec 7.4.2)
ble element. The gimbal angles are measured by transducers, usually resolvers
(Section 7 .4.2), whose outputs indicate the aircraft's roll, pitch, and heading
to the displays, auto-pilot, and sometimes to the computer. In strapdown sys-
tems, attitude angles are mathematically extracted from the analytic platform
transformation matrix (Section 7 .4.1 ).
When the inertial system is turned on, it must be aligned so that the computer
knows the initial position and groundspeed of the vehicle and so that the plat-
form (gimballed or analytic) has the correct initial orientation relative to the
Earth. The platform is typically aligned in such a way that its accelerometer
input axes are horizontal, often with one of them pointed north. As the vehi-
cle accelerates, maneuvers, and cruises, the accelerometers measure changes in
velocity, and the computer faithfully records the position and velocity.
The inertial navigator also contains power supplies for the instruments, a
computer, often a battery to protect against power transients, and interfaces
to a display-and-control unit. The system may be packaged in one or more
modules. Typical gimballed systems in 1968 weighed 50 to 75 lb (excluding
cables), of which 20 lb were for the platform. Steady-state power consumption
was approximately 200 w. First-generation strapdown navigators (early 1980s)
weighed 40 to 50 lb and consumed l 00 to 150 w. In 1996, strapdown systems
weighed 20 to 30 lb and consumed approximately 30 w.
7.3 INSTRUMENTS
7.3.1 Accelerometers
Purpose An accelerometer is a device that measures the force required to
accelerate a proof mass; thus, it measures the acceleration of the vehicle con-
taining the accelerometer. Figure 7.3 shows a black-box accelerometer whose
input axis is indicated. The instrument will supply an electrical output propor-
tional to (or some other determinate function of) the component along its input
axis of the inertial acceleration minus gravitation. If the instrument is mounted
in a vehicle whose inertial acceleration is a and if the vehicle travels in a New-
tonian gravitational field G, (Section 2.2), then the force acting on the proof
mass mp is
F = mpa = F R + mPG + F o
FR Fo
-=a-G--=f (accelerometer output) (7.1)
mp mp
where F R is the force exerted on the proof mass by the restoring spring or
318 INERTIAL NAVIGATION
a-G
-G
Electrical output _ _ 1 -
proportional to --,I 1
xa (a - G)IA I I I
I I I
Inertial acceleration
of accelerometer
G
R (Gravitation) along input axis
~til)(;. z
~
Permanent
magnet
Pendulum and___..--
torquer coil
Capacitive
pick off
plate
Flexure pivot
(flat metal spring)
where
TR is the residual torque applied to pendulum by friction in supports
and connecting wires and by electrical forces, Newton-meters
320 INERTIAL NAVIGATION
e= mb [
k
/y + TR/mb- I¢jmb
I + mbfzlk
l (7.3)
If the stiffness k is high enough, 11 is small, and the instrument measures only
f 1, independent of the presence of a cross-axis acceleration /2. Sensitivity to h
is called cross-coupling and is most serious in a vibration environment when
e andh oscillate in phase and rectify. This rectification is often referred to as
vibropendulous error. The term Trd k is the angle offset due to the presence
of an unwanted torque on the pendulum; it causes an accelerometer bias. I¢/k
is an angular offset of the pendulum due to angular acceleration of the case
around the pivot axis; ¢ is negligible when the accelerometer is mounted on a
mechanical platform, but it is an important source of error in strapdown systems
e
where ¢ and oscillations can rectify. If position calculations are referred to
the center of percussion of the pendulum, the sensitivity to angular acceleration
is reduced (see size effect, Section 7 .4.1 ). The distance from the center of mass
to the center of percussion is I/ mb.
Flexure-pivoted accelerometers are simpler to construct than the older
floated instruments since they do not require adjustment for buoyancy f23, pp.
288-289]. Because they are undamped, they exhibit high-frequency mechanical
resonances. Resonances must therefore be controlled relative to both vibration
inputs and rebalance-loop characteristics. Undamped accelerometers offer the
greatest bandwidth (important for strapdown systems) but must almost always
be supported on a shock-mounted sensor block (Section 7 .4.1) in order to sup-
press high-frequency vibration and to prevent shock damage. Accelerometers
that include fluid damping exhibit reduced bandwidth and additional thermal
sensitivity due to changes in the fluid characteristics.
In navigation-grade accelerometers, a restoring loop maintains the pendu-
lum near null. The restoring servo must be linear and repeatable from l 0 to 25
11-g to 40 g, a range of six to seven orders of magnitude. A digital output can
be obtained by either digitizing the analog output (the current in the torquer
coil) or by pulse-rebalancing with a digital restoring servo. When rebalancing
with pulses of uniform height, pulse width measures incremental velocity Ll V.
INSTRUMENTS 321
In either case, a properly initialized digital counter accumulates the pulses and
stores the velocity change. The rebalance pulse train must not excite accelerom-
eter resonances.
The pivot or flexure supporting the pendulum must provide minimal restraint
for the pendulum in the direction of the input axis while exhibiting high stiffness
in the other two directions. The spring constant of the pivot/flexure generates a
restoring force that reduces the gain of the electronic restoring loop. The spring
constant should be repeatable in order to ensure accuracy, but the high-stiffness
restoring loop dominates. The pivot must not exhibit hysteresis, which may
cause accelerometer biases. Generally, high-quality accelerometers ca