Hand Calculations logo Hand Calculations All help pages ▾

Gps Dilution Of Precision⚠ unverified

Aerospace / Navigation · Compute the geometric dilution of precision (GDOP)

Parameters

InputSymbolUnitDefaultDescription
hdophdop1.0Horizontal dilution of precision, dimensionless
vdopvdop1.0Vertical dilution of precision, dimensionless
OutputSymbolUnitDescription
resultGDOPGeometric dilution of precision, dimensionless

The science & history

Understanding the Parameters

Derivation (Approaching a Proof)

DOP falls out of the least-squares GPS solution. Linearising the pseudorange equations about an estimated position gives an observation matrix $\mathbf{G}$ (each row a unit line-of-sight vector to a satellite, plus a 1 for the clock term). The covariance of the position/time estimate is

$$\operatorname{cov} = \sigma_{\text{range}}^2\,(\mathbf{G}^\top\mathbf{G})^{-1},$$

so the geometry alone determines the error shape through $(\mathbf{G}^\top\mathbf{G})^{-1}$. The diagonal elements of this matrix are the squared DOPs for each coordinate:

$$HDOP^2 = q_{xx}+q_{yy}, \quad VDOP^2 = q_{zz}, \quad TDOP^2 = q_{tt},$$

where $q_{ii}$ are the diagonal entries of $(\mathbf{G}^\top\mathbf{G})^{-1}$. The composite DOPs are then just sub-sums of these variances (because variances of independent coordinate errors add):

$$PDOP = \sqrt{HDOP^2 + VDOP^2} = \sqrt{q_{xx}+q_{yy}+q_{zz}},$$ $$GDOP = \sqrt{PDOP^2 + TDOP^2} = \sqrt{q_{xx}+q_{yy}+q_{zz}+q_{tt}}.$$

The calculator computes the first line — PDOP — the root-sum-square of horizontal and vertical DOP. To reach true GDOP one would add the $TDOP^2$ (clock) term under the root, which the registry formula omits. The root-sum-square structure is simply the Pythagorean addition of independent error variances. $\qquad\blacksquare$

Dimensional check. $HDOP$ and $VDOP$ are dimensionless geometry factors, so $\sqrt{HDOP^2+VDOP^2}$ is dimensionless ✓ — DOP is a pure multiplier.

History and Development

Related Concepts: Barometric Altitude, Great-Circle Distance, Cross Track Error, Magnetic Heading, Circular Orbit Velocity

Notes: This is PDOP, not GDOP — $\sqrt{HDOP^2+VDOP^2}$ is the position DOP; true GDOP = $\sqrt{PDOP^2+TDOP^2}$ adds the clock term the registry omits (misnamed). DOP is a dimensionless geometry amplifier: $\sigma_{\text{pos}}\approx DOP\times$ ranging error. Small = satellites well spread; large = clustered. VDOP > HDOP typically (satellites only above the horizon), so GPS altitude is the weak axis. From $\operatorname{cov}=\sigma^2(\mathbf G^\top\mathbf G)^{-1}$. Defaults ⇒ $\sqrt2$.

← Back to the workspace  ·  All help pages  ·  Getting started