4.1.7 Sensor Error Sources
This chapter is intentionally cross-cutting. The individual sensor pages explain the physics and error mechanisms of gyroscopes, accelerometers, magnetometers, star trackers, Sun sensors, and GNSS. Here, those effects are placed into one common engineering framework so they can be used consistently in simulation, calibration, covariance design, Kalman filtering, Monte Carlo analysis, FDIR, and mission-level GNC verification.
The central engineering question is not simply, “What is the sensor accuracy?” It is: which errors are deterministic, which are stochastic, which vary with time or environment, which can be calibrated, which must be estimated online, which should appear in the process model, which should appear in the measurement covariance, and which must instead be handled by validity or fault logic?
A spacecraft sensor never provides the exact physical quantity that exists in reality. Every measurement is affected by imperfections associated with the sensing element, electronics, calibration, installation, spacecraft environment, timing system, signal processing, and stochastic processes.
The general idea is
\[ \boxed{ \text{True Physical Quantity} \rightarrow \text{Sensor Physics} \rightarrow \text{Deterministic Errors} \rightarrow \text{Stochastic Errors} \rightarrow \text{Sampling / Processing} \rightarrow \text{Measured Quantity} } \]A useful generic sensor model is
\[ \boxed{ \mathbf z_m = h(\mathbf x) + \mathbf e_d + \mathbf e_s } \]where
- \(\mathbf z_m\) = measured quantity,
- \(h(\mathbf x)\) = ideal measurement predicted from the true state,
- \(\mathbf e_d\) = deterministic or calibratable error,
- \(\mathbf e_s\) = stochastic error.
For a three-axis vector sensor, a more detailed model is
\[ \boxed{ \mathbf z_m = (\mathbf I+\mathbf S+\mathbf M) \mathbf z_{true} + \mathbf b + \mathbf n } \]where
\[ \mathbf S= \begin{bmatrix} s_x&0&0\\ 0&s_y&0\\ 0&0&s_z \end{bmatrix} \]represents scale-factor errors,
\[ \mathbf M \]represents misalignment, non-orthogonality and cross-axis coupling,
\[ \mathbf b \]is sensor bias, and
\[ \mathbf n \]represents stochastic measurement noise.
For real spacecraft sensors, this can be extended further:
\[ \boxed{ \mathbf z_m(t_k) = \mathcal Q \left[ \mathcal S \left( h(\mathbf x(t_k-\tau)) + \mathbf e_{cal} + \mathbf e_{env} + \mathbf e_{dyn} + \mathbf n \right) \right] } \]where \(\mathcal S\) represents sampling, \(\mathcal Q\) quantisation, and \(\tau\) measurement latency.
4.1.7.1 Why Sensor Errors Matter in GNC
Sensor errors propagate directly into the quantities used by guidance, navigation and control.
The complete chain is
\[ \boxed{ \text{Sensor Error} \rightarrow \text{Measurement Error} \rightarrow \text{State-Estimation Error} \rightarrow \text{Guidance Error} \rightarrow \text{Control Error} \rightarrow \text{Spacecraft Performance} } \]For example, a gyroscope bias produces attitude drift:
\[ b_g \rightarrow \delta\omega \rightarrow \delta\theta. \]An accelerometer bias produces velocity and position drift:
\[ b_a \rightarrow \delta v \rightarrow \delta r. \]A magnetometer bias distorts the measured geomagnetic vector and therefore affects attitude estimation and magnetic control.
Star-tracker centroiding, alignment or thermal errors directly influence spacecraft attitude knowledge.
GNSS range, clock, geometry and antenna errors propagate into position, velocity and timing estimates.
Consequently, sensor modelling is not merely an instrumentation problem. It is fundamental to the entire GNC error budget.
4.1.7.2 Sensor Error Classification
Sensor errors can be classified in several ways.
A useful top-level classification is
\[ \boxed{ \text{Sensor Errors} = \text{Deterministic} + \text{Stochastic} + \text{Environmental} + \text{Dynamic} + \text{Digital} + \text{Installation} + \text{Operational} } \]Deterministic errors
These can often be characterized and calibrated:
- constant bias,
- scale-factor error,
- axis misalignment,
- non-orthogonality,
- cross-axis sensitivity,
- nonlinearity,
- known temperature dependence,
- mounting alignment error,
- fixed latency.
Stochastic errors
These cannot be predicted exactly for an individual measurement:
- white noise,
- random walk,
- bias instability,
- correlated noise,
- coloured noise,
- flicker noise,
- shot noise,
- clock noise.
Environmental errors
Examples include:
- temperature,
- radiation,
- magnetic contamination,
- vibration,
- mechanical stress,
- electromagnetic interference,
- stray light,
- Earth albedo,
- thermal gradients.
Dynamic errors
Examples include:
- finite bandwidth,
- response lag,
- motion blur,
- vibration sensitivity,
- aliasing,
- slew-rate limitations.
Digital errors
Examples include:
- quantisation,
- finite numerical precision,
- timestamp discretisation,
- packet loss,
- sample jitter.
Installation errors
Examples include:
- sensor-to-body misalignment,
- lever-arm displacement,
- antenna phase-centre offset,
- structural deformation.
4.1.7.3 Accuracy, Precision, Resolution and Repeatability
These terms should not be used interchangeably.
Accuracy
Accuracy describes closeness to the true value.
\[ \boxed{ e=z_m-z_{true} } \]Small measurement error generally corresponds to high accuracy.
Precision
Precision describes how tightly repeated measurements cluster.
A sensor may be highly precise but biased.
Resolution
Resolution is the smallest input change that can be represented or detected.
For an \(N\)-bit sensor over range
\[ x_{min}\le x\le x_{max}, \]an approximate digital resolution is
\[ \boxed{ \Delta x = \frac{x_{max}-x_{min}}{2^N} } \]depending on encoding convention.
Repeatability
Repeatability describes the ability to reproduce measurements under nominally identical conditions.
Stability
Stability describes how sensor characteristics change with time.
These distinctions are essential when reading sensor datasheets.
4.1.7.4 Bias
Bias is an additive offset between the sensor output and the ideal measurement.
\[ \boxed{ z_m=z_{true}+b } \]For a vector sensor,
\[ \boxed{ \mathbf z_m = \mathbf z_{true} + \mathbf b. } \]Bias may arise from:
- electronics,
- manufacturing tolerances,
- mechanical stress,
- calibration residuals,
- temperature,
- ageing,
- radiation,
- magnetic contamination.
4.1.7.5 Constant Bias
The simplest model assumes
\[ \boxed{ b(t)=b_0. } \]Thus
\[ z_m=z_{true}+b_0+n. \]Constant bias can often be estimated during calibration.
However, even a small constant bias can become extremely important when the measurement is integrated.
For a gyro,
\[ \delta\theta(t)\approx b_g t. \]For an accelerometer,
\[ \delta v(t)\approx b_a t \]and
\[ \boxed{ \delta r(t)\approx\frac12b_a t^2. } \]Therefore an apparently small inertial-sensor bias may produce large navigation errors.
4.1.7.6 Bias Repeatability
A sensor may not return to exactly the same bias after:
- power cycling,
- thermal cycling,
- launch,
- mechanical shock,
- long shutdown periods.
This is called bias repeatability or turn-on-to-turn-on bias variation.
A useful model is
\[ \boxed{ b_0\sim\mathcal N(\bar b,\sigma_{b0}^2). } \]Each Monte Carlo realization can therefore begin with a different initial bias.
4.1.7.7 Bias Drift
Real bias usually changes slowly with time:
\[ \boxed{ b=b(t). } \]The measurement becomes
\[ z_m=z_{true}+b(t)+n(t). \]Bias drift may be caused by:
- temperature,
- electronics ageing,
- mechanical relaxation,
- radiation,
- changing operating conditions.
Bias drift is especially important for inertial sensors.
4.1.7.8 Bias Random Walk
A common stochastic model is
\[ \boxed{ \dot b=w_b } \]where
\[ w_b(t)\sim\mathcal N(0,q_b). \]In discrete time,
\[ \boxed{ b_{k+1}=b_k+w_{b,k}. } \]More carefully, for a continuous white driving process with spectral density \(q_b\),
\[ \boxed{ b_{k+1} = b_k+\sqrt{q_b\Delta t}\,\xi_k } \]where
\[ \xi_k\sim\mathcal N(0,1). \]This model is widely used in Kalman filters.
4.1.7.9 First-Order Gauss-Markov Bias
Not all bias processes should be represented as unbounded random walks.
A first-order Gauss-Markov model is
\[ \boxed{ \dot b = -\frac{1}{\tau_b}b+w_b } \]where \(\tau_b\) is the correlation time.
Its discrete form is
\[ \boxed{ b_{k+1} = e^{-\Delta t/\tau_b}b_k+w_k. } \]This gives a correlated but bounded stochastic process.
4.1.7.10 Bias Instability
Bias instability describes slowly varying fluctuations that are not adequately represented by simple white noise.
It is commonly identified using Allan deviation.
Conceptually,
\[ \boxed{ \text{Bias Instability} \rightarrow \text{Slowly Varying Measurement Offset}. } \]It is particularly important for gyroscopes and accelerometers.
Bias instability should not automatically be treated as identical to bias random walk. They represent different stochastic behaviours.
4.1.7.11 Scale-Factor Error
Scale-factor error causes the sensor response to differ from the ideal proportionality.
For a scalar measurement,
\[ \boxed{ z_m=(1+s)z_{true}. } \]If
\[ s=0.01, \]the scale-factor error is approximately \(1\%\).
For a three-axis sensor,
\[ \boxed{ \mathbf z_m = (\mathbf I+\mathbf S)\mathbf z_{true}. } \]Scale-factor error becomes particularly important when the measured quantity is large.
4.1.7.12 Scale-Factor Asymmetry
Some sensors have different scale factors for positive and negative inputs:
\[ s_+\neq s_-. \]Therefore,
\[ z_m= \begin{cases} (1+s_+)z, & z\ge0\\ (1+s_-)z, & z<0. \end{cases} \]A single scale-factor coefficient may therefore be insufficient for high-fidelity modelling.
4.1.7.13 Scale-Factor Temperature Dependence
Scale factor can vary with temperature:
\[ \boxed{ s(T) = s_0+k_{s1}(T-T_0)+k_{s2}(T-T_0)^2+\cdots } \]so that
\[ z_m=[1+s(T)]z_{true}. \]This is one reason calibration tables are sometimes functions of temperature rather than fixed constants.
4.1.7.14 Axis Misalignment
A sensor's physical axes are not perfectly aligned with the assumed coordinate frame.
If
\[ C_B^S \]maps a body-frame quantity into the sensor frame,
\[ \boxed{ \mathbf z^S=C_B^S\mathbf z^B. } \]An error in \(C_B^S\) produces a mounting or alignment error.
For small angular misalignment,
\[ C\approx I-[\delta\boldsymbol\theta\times] \]for one common convention.
Therefore,
\[ \boxed{ \delta\mathbf z \approx -[\delta\boldsymbol\theta\times]\mathbf z. } \]The sign depends on the adopted frame/transformation convention.
4.1.7.15 Non-Orthogonality
Ideally, the three sensor axes are mutually perpendicular:
\[ \hat x\cdot\hat y = \hat y\cdot\hat z = \hat z\cdot\hat x = 0. \]Manufacturing imperfections can cause small deviations.
A non-orthogonality matrix can be represented as
\[ \boxed{ M= \begin{bmatrix} 0&m_{xy}&m_{xz}\\ m_{yx}&0&m_{yz}\\ m_{zx}&m_{zy}&0 \end{bmatrix}. } \]The sensor model becomes
\[ \mathbf z_m=(I+M)\mathbf z_{true}+\cdots. \]4.1.7.16 Cross-Axis Sensitivity
An input along one physical axis may generate an output on another axis.
For example,
\[ z_{m,x} = z_x + c_{xy}z_y + c_{xz}z_z. \]The complete model is
\[ \boxed{ \mathbf z_m=C_{cross}\mathbf z_{true}. } \]Cross-axis sensitivity can result from sensing-element construction, electronics, mounting, or imperfect calibration.
4.1.7.17 Combined Calibration Matrix
Scale factor, misalignment and cross-axis effects can be grouped into one calibration matrix:
\[ \boxed{ \mathbf z_m = A\mathbf z_{true} + \mathbf b + \mathbf n. } \]where
\[ A = \begin{bmatrix} a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33} \end{bmatrix}. \]Calibration seeks estimates
\[ \hat A,\qquad\hat{\mathbf b}. \]The corrected measurement becomes
\[ \boxed{ \hat{\mathbf z}_{true} = \hat A^{-1} (\mathbf z_m-\hat{\mathbf b}). } \]This formulation is useful for accelerometers, gyroscopes and magnetometers.
4.1.7.18 Nonlinearity
A real sensor may not maintain a perfectly linear input-output relationship.
Instead of
\[ z_m=Kz, \]the response may be
\[ \boxed{ z_m = K_1z + K_2z^2 + K_3z^3 +\cdots. } \]Nonlinearity is often specified relative to full scale.
It becomes particularly important near the sensor's operating-range limits.
4.1.7.19 Hysteresis
Hysteresis means the output depends not only on the current input but also on its previous history.
Conceptually,
\[ \boxed{ z_m(t) = h[z(t),\text{history}]. } \]The output for increasing input may differ from the output for decreasing input.
Hysteresis may appear in:
- magnetic sensors,
- mechanical sensing elements,
- thermal systems.
4.1.7.20 Dead Zone
A sensor may have a region in which small input changes produce little or no observable output.
\[ z_m= \begin{cases} 0,& |z|4.1.7.21 Saturation
Every physical sensor has a finite operating range.
If
\[ z_{min}\le z\le z_{max}, \]the measured value may be represented as
\[ \boxed{ z_m = \operatorname{sat} (z,z_{min},z_{max}). } \]When the input exceeds the measurement range, information is lost.
Examples include:
- gyro rate exceeding full scale,
- accelerometer acceleration exceeding range,
- magnetometer saturation near strong magnetic fields,
- photodetector saturation,
- GNSS RF front-end overload.
A realistic simulation should generate a saturation flag when appropriate.
4.1.7.22 Quantisation
Digital sensors represent continuous measurements using discrete levels.
For quantisation step
\[ \Delta, \] \[ \boxed{ z_q = \Delta \operatorname{round} \left( \frac{z}{\Delta} \right). } \]Quantisation error is
\[ e_q=z_q-z. \]For an ideal uniform quantiser,
\[ -\frac{\Delta}{2} \le e_q < \frac{\Delta}{2}. \]Under suitable assumptions,
\[ \boxed{ \sigma_q^2\approx\frac{\Delta^2}{12}. } \]This approximation should not be blindly applied when the signal is highly correlated with the quantisation process.
4.1.7.23 Resolution vs Quantisation
Resolution and quantisation are related but conceptually different.
Resolution describes the smallest distinguishable input increment.
Quantisation describes the process of mapping a continuous value onto discrete output levels.
A sensor may also have effective resolution poorer than its digital bit depth because of analog noise.
Thus,
\[ \boxed{ \text{Digital Bits} \neq \text{Effective Measurement Resolution}. } \]4.1.7.24 White Noise
A common sensor-noise model is zero-mean Gaussian white noise:
\[ \boxed{ n(t)\sim\mathcal N(0,\sigma^2). } \]In discrete form,
\[ z_k=z_{true,k}+n_k. \]with
\[ n_k\sim\mathcal N(0,R). \]For vector measurements,
\[ \boxed{ \mathbf n_k\sim\mathcal N(\mathbf0,R). } \]The covariance matrix need not be diagonal.
4.1.7.25 Measurement Noise Covariance
For a three-axis sensor,
\[ R= \begin{bmatrix} \sigma_x^2&R_{xy}&R_{xz}\\ R_{yx}&\sigma_y^2&R_{yz}\\ R_{zx}&R_{zy}&\sigma_z^2 \end{bmatrix}. \]If axis errors are independent,
\[ R= \operatorname{diag} (\sigma_x^2,\sigma_y^2,\sigma_z^2). \]Assuming diagonal covariance without justification can underestimate filter uncertainty when measurement errors are correlated.
4.1.7.26 Noise Density
Sensor noise is often specified as an amplitude spectral density.
Typical units include
\[ \frac{^\circ/s}{\sqrt{\mathrm{Hz}}}, \qquad \frac{m/s^2}{\sqrt{\mathrm{Hz}}}, \qquad \frac{nT}{\sqrt{\mathrm{Hz}}}. \]A generic relationship is
\[ \boxed{ \sigma^2 = \int S_n(f)\,df. } \]For approximately white noise over effective bandwidth \(B\),
\[ \boxed{ \sigma \propto N_d\sqrt{B}. } \]The exact numerical conversion depends on whether the specification uses one-sided or two-sided PSD conventions and on the actual filter bandwidth.
4.1.7.27 Power Spectral Density
The power spectral density describes how noise power is distributed with frequency.
\[ \boxed{ S_x(f). } \]White noise approximately has
\[ S_x(f)=\text{constant}. \]Coloured noise has frequency-dependent PSD.
PSD analysis helps distinguish different stochastic processes that cannot be understood from RMS noise alone.
4.1.7.28 Coloured Noise
Real sensor noise is frequently temporally correlated.
Unlike white noise,
\[ E[n_k n_{k-j}] \neq0 \]for some \(j\neq0\).
Examples include:
- low-frequency drift,
- \(1/f\) noise,
- thermally correlated noise,
- structural vibration,
- electronics filtering.
A filter that assumes white measurement noise may become inconsistent if the real noise is strongly correlated.
4.1.7.29 Flicker Noise
Flicker noise is often associated with approximately
\[ \boxed{ S(f)\propto\frac1f. } \]It dominates many sensors and electronic systems at low frequencies.
It is related conceptually to the slowly varying behaviour observed as bias instability.
4.1.7.30 Random Walk
Random walk arises when white noise is integrated.
\[ \boxed{ x_{k+1}=x_k+w_k. } \]Its variance grows with time.
This behaviour is fundamental in inertial navigation.
4.1.7.31 Angle Random Walk
Gyroscope white rate noise integrates into attitude uncertainty.
Conceptually,
\[ \boxed{ \text{Gyro White Noise} \xrightarrow{\int} \text{Angle Random Walk}. } \]Attitude uncertainty therefore increases with integration time.
4.1.7.32 Velocity Random Walk
Accelerometer white noise integrates into velocity uncertainty:
\[ \boxed{ \text{Acceleration Noise} \xrightarrow{\int} \text{Velocity Random Walk}. } \]Further integration propagates the error into position.
4.1.7.33 Rate Random Walk
Rate random walk describes a lower-frequency stochastic process in gyroscope output, distinct from ordinary white rate noise.
It is one of several noise mechanisms that can appear in Allan-deviation analysis.
4.1.7.34 Allan Variance
Allan variance is widely used to characterize inertial-sensor stochastic errors.
For averaging interval \(\tau\),
\[ \boxed{ \sigma_A^2(\tau) = \frac12 E [ (\bar y_{k+1}-\bar y_k)^2 ]. } \]The Allan deviation is
\[ \boxed{ \sigma_A(\tau) = \sqrt{\sigma_A^2(\tau)}. } \]Plotting
\[ \log\sigma_A \quad\text{against}\quad \log\tau \]reveals different stochastic regimes.
4.1.7.35 Interpreting an Allan-Deviation Plot
Different slope regions may correspond to different processes such as:
- quantisation noise,
- white noise,
- bias instability,
- random walk,
- rate ramp.
Conceptually,
\[ \boxed{ \text{Time-Series Data} \rightarrow \text{Allan Deviation} \rightarrow \text{Noise Identification} \rightarrow \text{Stochastic Sensor Model}. } \]Coefficient extraction depends on the exact Allan-variance convention and units being used, so numerical conversion factors should be applied carefully.
4.1.7.36 Temperature Effects
Temperature affects nearly every spacecraft sensor.
A bias model can be written as
\[ \boxed{ b(T) = b_0 + k_{b1}(T-T_0) + k_{b2}(T-T_0)^2 +\cdots. } \]Similarly,
\[ s=s(T). \]A complete measurement may therefore be
\[ \boxed{ z_m = [1+s(T)]z_{true} + b(T) + n. } \]Temperature effects may arise from:
- electronics,
- mechanical expansion,
- refractive changes,
- detector sensitivity,
- oscillator drift,
- structural alignment.
4.1.7.37 Thermal Gradient Effects
Uniform temperature and thermal gradients are different problems.
A sensor may have acceptable bias stability at uniform temperature but experience alignment errors when different parts of the spacecraft structure are at different temperatures.
For precision attitude determination,
\[ \boxed{ \text{Thermal Gradient} \rightarrow \text{Structural Deformation} \rightarrow \text{Sensor-to-Payload Misalignment}. } \]This is especially important for star trackers and precision payloads.
4.1.7.38 Ageing
Sensor characteristics may change over mission life:
\[ b=b(t_{life}), \qquad s=s(t_{life}). \]Ageing can affect:
- electronics,
- oscillators,
- detector sensitivity,
- magnetic materials,
- optical components,
- mechanical structures.
Long-duration missions may therefore require in-orbit recalibration.
4.1.7.39 Radiation Effects
The space radiation environment can produce:
- total ionizing dose effects,
- displacement damage,
- single-event effects,
- detector hot pixels,
- electronics drift,
- memory corruption,
- temporary sensor upset.
Conceptually,
\[ \boxed{ \text{Radiation} \rightarrow \text{Sensor Parameter Drift / Transient Fault / Permanent Degradation}. } \]Radiation should therefore be considered both as an error source and as a reliability issue.
4.1.7.40 Vibration and Mechanical Effects
Launch and on-orbit vibration can influence sensor measurements.
Sources include:
- reaction wheels,
- cryocoolers,
- pumps,
- mechanisms,
- thruster firing,
- structural modes.
A vibration signal may enter the sensor bandwidth and appear as real measurement content.
Thus,
\[ \boxed{ \text{Mechanical Disturbance} \rightarrow \text{Sensor Motion} \rightarrow \text{Measured Disturbance}. } \]This is not always sensor noise—the sensor may be correctly measuring actual structural motion.
4.1.7.41 g-Sensitivity of Gyroscopes
Some gyroscopes exhibit sensitivity to linear acceleration:
\[ \boxed{ \boldsymbol{\omega}_m = \boldsymbol{\omega} + K_g\mathbf a + \mathbf b_g + \mathbf n_g. } \]This may matter during:
- launch,
- thruster firing,
- vibration,
- high-acceleration maneuvers.
4.1.7.42 Lever-Arm Effects
A sensor mounted away from the spacecraft centre of mass experiences rotational acceleration.
For sensor offset \(\mathbf r\),
\[ \boxed{ \mathbf a_s = \mathbf a_O + \dot{\boldsymbol\omega}\times\mathbf r + \boldsymbol\omega\times (\boldsymbol\omega\times\mathbf r) } \]for a rigid body-fixed sensor.
Therefore an accelerometer away from the reference point measures both translational and rotational effects.
GNSS antennas also require proper treatment of antenna location relative to the spacecraft reference point.
4.1.7.43 Electromagnetic Interference
Spacecraft electronics can interfere with sensors.
Potential sources include:
- switching converters,
- transmitters,
- reaction-wheel electronics,
- high-current cables,
- magnetorquers,
- digital buses.
Magnetometers are particularly sensitive.
Conceptually,
\[ \boxed{ \mathbf B_m = \mathbf B_{Earth} + \mathbf B_{spacecraft} + \mathbf n. } \]4.1.7.44 Optical Contamination and Stray Light
Optical sensors are affected by environmental optical errors.
Examples include:
- direct Sun,
- Earth limb,
- Moon,
- Earth albedo,
- spacecraft reflections,
- contamination,
- radiation-induced detector defects.
These effects apply particularly to:
- star trackers,
- Sun sensors,
- cameras.
The measurement may become biased or completely invalid.
Therefore a realistic model requires a validity state, not merely increased Gaussian noise.
4.1.7.45 GNSS-Specific Environmental Errors
GNSS contains several error sources that differ from conventional inertial sensors:
- satellite orbit error,
- satellite clock error,
- receiver clock error,
- propagation effects,
- multipath,
- antenna phase-centre variation,
- RF interference,
- geometry,
- tracking error,
- cycle slips.
Thus GNSS error cannot generally be represented only by
\[ \mathbf r_m=\mathbf r+\mathbf n. \]That simple model is appropriate only at solution-level simulation fidelity.
4.1.7.46 Sampling
Digital sensors provide measurements at discrete times:
\[ \boxed{ t_k=kT_s. } \]The sampled measurement is
\[ z[k]=z(t_k). \]Sampling introduces important GNC considerations:
- update rate,
- synchronization,
- aliasing,
- latency,
- missed samples.
4.1.7.47 Aliasing
If the sampling frequency is insufficient, high-frequency signal content can appear as false lower-frequency content.
The Nyquist condition for a band-limited signal is
\[ \boxed{ f_s>2f_{max}. } \]In practice, anti-aliasing filters are normally required before sampling.
Aliasing is especially important when spacecraft vibration frequencies overlap sensor bandwidth.
4.1.7.48 Sensor Bandwidth
A real sensor cannot respond instantaneously.
A simple first-order model is
\[ \boxed{ G(s)=\frac{1}{\tau s+1}. } \]The cutoff frequency is approximately
\[ \boxed{ f_c=\frac{1}{2\pi\tau}. } \]Bandwidth affects:
- dynamic tracking,
- phase lag,
- noise,
- control-loop interaction.
A higher bandwidth is not automatically better because it can also admit more noise.
4.1.7.49 Latency
A sensor measurement may correspond to an earlier physical state.
\[ \boxed{ z_m(t)=h[x(t-\tau)]+n(t). } \]For a rotating spacecraft, attitude error caused by timing error is approximately
\[ \boxed{ \delta\boldsymbol\theta \approx \boldsymbol\omega\,\delta t. } \]Latency therefore becomes increasingly important during fast maneuvers.
4.1.7.50 Latency vs Timestamp Error
These should be distinguished.
Known latency
If
\[ \tau \]is known accurately, the estimator can often compensate for it.
Timestamp error
If the measurement epoch itself is uncertain,
\[ \delta t\neq0, \]the measurement is associated with the wrong spacecraft state.
Therefore,
\[ \boxed{ \text{Known Latency} \neq \text{Timing Uncertainty}. } \]4.1.7.51 Sampling Jitter
The actual sample time may differ from its nominal value:
\[ \boxed{ t_k=kT_s+\delta t_k. } \]For rapidly changing signals,
\[ \delta z \approx \dot z\,\delta t. \]Sampling jitter can therefore become a measurable sensor error.
4.1.7.52 Synchronization Error Between Sensors
Sensor fusion assumes measurements are associated with consistent epochs.
Suppose
\[ z_1=z_1(t) \]but
\[ z_2=z_2(t-\Delta t). \]Directly combining them as simultaneous measurements introduces inconsistency.
This matters for:
- GNSS + IMU,
- gyro + star tracker,
- camera + LiDAR,
- differential GNSS between spacecraft.
4.1.7.53 Dropout
Sensors do not always provide measurements.
Define
\[ \boxed{ valid_k\in\{0,1\}. } \]Then
\[ valid_k=0 \]means the measurement should not be used.
Dropout can arise from:
- star-tracker blinding,
- Sun-sensor eclipse,
- GNSS outage,
- communication failure,
- sensor reset,
- invalid signal quality.
4.1.7.54 Outliers
Not every erroneous measurement follows a Gaussian distribution.
An outlier may be represented conceptually as
\[ \boxed{ z_m=z_{true}+n+o } \]where \(o\) is a large occasional error.
Sources include:
- false star identification,
- GNSS multipath,
- cycle slip,
- detector artifact,
- packet corruption.
Outliers require robust validation rather than merely increasing \(R\).
4.1.7.55 Stuck Measurements
A sensor can fail while continuing to output a plausible-looking constant value:
\[ \boxed{ z_k=z_{k-1}. } \]A simple range check may not detect such a fault.
Detection may require:
- rate-of-change checks,
- comparison with predicted dynamics,
- redundant sensors,
- innovation monitoring.
4.1.7.56 Spike Errors
A transient sensor fault may produce
\[ z_k=z_{true,k}+A_{spike}. \]Spike detection can use:
- thresholding,
- innovation gates,
- temporal consistency,
- median filtering where appropriate.
Care is needed not to reject genuine rapid spacecraft motion.
4.1.7.57 Noise vs Fault
A fundamental distinction is
\[ \boxed{ \text{Expected Random Variation} \neq \text{Sensor Fault}. } \]Noise belongs within the expected statistical measurement model.
Faults correspond to abnormal behaviour outside the assumed model.
The GNC architecture should distinguish these cases because the response is different.
4.1.7.58 Measurement Validity Logic
A realistic sensor output should include more than a numerical value.
A measurement packet can contain
\[ \boxed{ \mathcal Z = \{ \mathbf z,\, t,\, R,\, valid,\, quality,\, status \}. } \]This allows the estimator to know:
- what was measured,
- when it was measured,
- its uncertainty,
- whether it should be trusted.
4.1.7.59 Innovation Monitoring
For an estimator,
\[ \boxed{ \mathbf y_k = \mathbf z_k - h(\hat{\mathbf x}_k^-) } \]is the innovation.
Its covariance is
\[ \boxed{ S_k = H_kP_k^-H_k^T+R_k. } \]Large innovations may indicate:
- bad measurements,
- incorrect sensor covariance,
- model error,
- maneuver not represented in dynamics,
- sensor fault.
4.1.7.60 Mahalanobis Gating
A normalized innovation statistic is
\[ \boxed{ d^2 = \mathbf y^TS^{-1}\mathbf y. } \]A measurement can be rejected when
\[ \boxed{ d^2>\gamma } \]where \(\gamma\) is selected from an appropriate statistical threshold.
This is preferable to applying the same raw-error threshold to measurements with different covariance.
4.1.7.61 Residual Monitoring
With redundant measurements, residuals can be formed between sensors.
For two sensors measuring the same quantity,
\[ \boxed{ r=z_1-z_2. } \]A growing residual may indicate sensor degradation.
However, the expected residual covariance must include uncertainty from both sensors.
4.1.7.62 Redundancy
Spacecraft often carry redundant sensors.
Possible configurations include:
- cold redundancy,
- warm redundancy,
- hot redundancy,
- multiple sensing heads,
- dissimilar sensor redundancy.
Redundancy provides:
\[ \boxed{ \text{Fault Detection} + \text{Fault Isolation} + \text{Continued Operation}. } \]4.1.7.63 Common-Mode Errors
Redundant sensors do not guarantee independent errors.
If both sensors experience the same disturbance,
\[ e_1\approx e_2, \]comparison may not reveal the problem.
Examples include:
- common temperature environment,
- common structural deformation,
- shared clock error,
- shared GNSS propagation errors,
- spacecraft magnetic contamination.
This is a major consideration in redundancy architecture.
4.1.7.64 Correlated Sensor Errors
For two measurements,
\[ E[e_1e_2^T] = P_{12}. \]If
\[ P_{12}\neq0, \]their errors are correlated.
The combined covariance is therefore not simply diagonal.
Ignoring correlation can cause an estimator to become overconfident.
4.1.7.65 Calibration
Calibration estimates parameters relating sensor output to the physical truth.
Conceptually,
\[ \boxed{ \text{Known Reference} + \text{Sensor Measurement} \rightarrow \text{Parameter Estimation} \rightarrow \text{Calibration Model}. } \]Typical calibration parameters include:
- bias,
- scale factor,
- axis alignment,
- non-orthogonality,
- temperature coefficients,
- optical distortion,
- antenna offset.
4.1.7.66 One-Point Calibration
A one-point calibration can estimate an offset.
If the reference value is \(z_{ref}\),
\[ \boxed{ \hat b=z_m-z_{ref}. } \]This does not independently determine scale factor.
4.1.7.67 Two-Point Calibration
With two known inputs,
\[ z_{ref,1},\qquad z_{ref,2}, \]and outputs
\[ z_{m,1},\qquad z_{m,2}, \]a linear model
\[ z_m=az+b \]can estimate both scale factor and offset.
4.1.7.68 Multi-Position Calibration
Three-axis inertial and magnetic sensors frequently require measurements in many orientations.
The calibration problem becomes
\[ \boxed{ \min_{A,b} \sum_k \| \mathbf z_{m,k} - A\mathbf z_{ref,k} - \mathbf b \|^2. } \]This can estimate:
- scale factors,
- cross-axis terms,
- bias,
- alignment parameters.
4.1.7.69 Temperature Calibration
Sensor parameters can be measured over a thermal range:
\[ \boxed{ \{T_i,b_i,s_i\} \rightarrow b(T),s(T). } \]The flight software can then apply temperature compensation.
This requires an accurate sensor temperature measurement and an appropriate calibration model.
4.1.7.70 Ground Calibration vs In-Orbit Calibration
Ground calibration
Allows controlled reference inputs and environmental testing.
In-orbit calibration
Can account for:
- launch shifts,
- zero-g environment,
- thermal conditions,
- ageing,
- structural alignment changes.
Therefore,
\[ \boxed{ \text{Ground Calibration} + \text{In-Orbit Estimation} } \]is often more powerful than either alone.
4.1.7.71 Calibration Residual Error
Calibration does not eliminate uncertainty completely.
If
\[ \hat b=b+\delta b, \]then after correction,
\[ z_{corr}=z_m-\hat b \]contains
\[ -\delta b. \]Therefore calibration parameters themselves require uncertainty models.
4.1.7.72 Sensor Error States in a Kalman Filter
Slowly varying errors can be estimated as states.
For an INS,
\[ \boxed{ \delta\mathbf x = \begin{bmatrix} \delta\mathbf r\\ \delta\mathbf v\\ \delta\boldsymbol\theta\\ \mathbf b_a\\ \mathbf b_g \end{bmatrix}. } \]Additional states may include:
- scale factors,
- clock bias,
- clock drift,
- alignment errors.
The appropriate state vector depends on mission requirements and observability.
4.1.7.73 Process Noise vs Measurement Noise
These are fundamentally different.
Process noise
Represents uncertainty in state evolution:
\[ \boxed{ \mathbf x_{k+1} = f(\mathbf x_k)+\mathbf w_k } \]with
\[ \mathbf w_k\sim\mathcal N(0,Q). \]Measurement noise
Represents uncertainty in sensor observations:
\[ \boxed{ \mathbf z_k = h(\mathbf x_k)+\mathbf v_k } \]with
\[ \mathbf v_k\sim\mathcal N(0,R). \]Thus,
\[ \boxed{ Q\neq R. } \]Confusing the two produces incorrect filter tuning.
4.1.7.74 Observability of Sensor Errors
Including a bias state in an EKF does not guarantee that the filter can estimate it.
A parameter must influence measurements in a distinguishable way.
Conceptually,
\[ \boxed{ \text{State Included in Filter} \not\Rightarrow \text{State Observable}. } \]Observability depends on:
- sensor combination,
- spacecraft motion,
- reference geometry,
- measurement availability.
4.1.7.75 Error Propagation Through Integration
Some sensor errors become increasingly severe because GNC algorithms integrate their measurements.
Gyro:
\[ \boxed{ \delta\theta = \int\delta\omega\,dt. } \]Accelerometer:
\[ \boxed{ \delta v = \int\delta a\,dt } \]and
\[ \boxed{ \delta r = \int\delta v\,dt. } \]This explains why inertial navigation requires external aiding.
4.1.7.76 Attitude Error Coupling into Acceleration
Accelerometer measurements must be transformed using the estimated attitude.
If the attitude estimate contains
\[ \delta\boldsymbol\theta, \]then the transformed acceleration also contains error.
For a small attitude error,
\[ \boxed{ \delta\mathbf a \approx -[\mathbf f\times] \delta\boldsymbol\theta } \]up to the adopted frame/error convention.
Thus gyro errors can indirectly produce position errors.
4.1.7.77 Frame-Transformation Errors
Sensor measurements often require transformations:
\[ \boxed{ \mathbf z^B = C_S^B\mathbf z^S. } \]If
\[ C_S^B \]is uncertain, the transformed measurement contains additional error.
This applies to:
- IMUs,
- magnetometers,
- Sun sensors,
- star trackers,
- GNSS antenna reference points.
4.1.7.78 Error Covariance Transformation
Suppose measurement covariance is known in frame \(A\):
\[ P^A. \]After rotation into frame \(B\),
\[ \boxed{ P^B = C_A^B P^A (C_A^B)^T. } \]Covariance must be transformed with the measurement.
Rotating the vector but leaving its covariance in the original frame is inconsistent.
4.1.7.79 Linear Error Propagation
For
\[ \mathbf y=f(\mathbf x), \]small perturbations satisfy
\[ \boxed{ \delta\mathbf y \approx J\delta\mathbf x } \]where
\[ J=\frac{\partial f}{\partial\mathbf x}. \]Therefore covariance propagates as
\[ \boxed{ P_y \approx JP_xJ^T. } \]This is one of the fundamental tools of GNC uncertainty analysis.
4.1.7.80 Root-Sum-Square Error Budget
If scalar error contributors are independent,
\[ \boxed{ \sigma_{total} = \sqrt{ \sigma_1^2+ \sigma_2^2+ \cdots+ \sigma_n^2 }. } \]This is commonly used in engineering error budgets.
However, RSS is not valid without modification when errors are correlated.
4.1.7.81 Covariance-Based Error Budget
For correlated vector errors,
\[ \boxed{ P_{total} = \sum_iP_i + \sum_{i\neq j}P_{ij}. } \]This is more general than simply adding independent \(1\sigma\) values by RSS.
4.1.7.82 1σ, 2σ and 3σ
For a one-dimensional Gaussian distribution,
\[ x\sim\mathcal N(\mu,\sigma^2), \]the intervals
\[ \mu\pm\sigma, \qquad \mu\pm2\sigma, \qquad \mu\pm3\sigma \]represent progressively larger probability regions.
However, the corresponding probabilities for multidimensional error ellipsoids are not identical to the one-dimensional interpretation.
4.1.7.83 RMS Error
Root-mean-square error is
\[ \boxed{ RMSE = \sqrt{ \frac1N \sum_{k=1}^{N} (z_k-z_{true,k})^2 }. } \]RMSE combines bias and random variation.
For a scalar error with mean bias \(b\) and variance \(\sigma^2\),
\[ \boxed{ RMSE^2=b^2+\sigma^2. } \]4.1.7.84 Monte Carlo Sensor Error Modelling
Monte Carlo simulation is essential for evaluating combined sensor uncertainties.
For run \(j\), sample uncertain parameters such as
\[ \boxed{ b^{(j)},\; s^{(j)},\; \delta\theta^{(j)},\; \sigma^{(j)},\; \tau^{(j)}. } \]Then simulate
\[ \boxed{ \text{Truth} \rightarrow \text{Sensor Errors} \rightarrow \text{Estimator} \rightarrow \text{Guidance} \rightarrow \text{Controller} \rightarrow \text{Mission Outcome}. } \]This captures nonlinear interaction between errors that cannot always be understood from an analytical RSS budget.
4.1.7.85 Fixed vs Run-to-Run vs Time-Varying Errors
This distinction is critical in Monte Carlo simulation.
Fixed parameter
Constant for every simulation:
\[ p=p_0. \]Run-to-run uncertainty
Randomly sampled once per Monte Carlo run:
\[ p^{(j)}\sim\mathcal D_p. \]Examples:
- mounting error,
- initial bias,
- calibration residual.
Time-varying uncertainty
Changes throughout a simulation:
\[ p=p(t). \]Examples:
- white noise,
- bias drift,
- dropout,
- temperature.
These should not be modelled identically.
4.1.7.86 Worst-Case Analysis vs Monte Carlo
Worst-case analysis asks whether requirements remain satisfied under extreme parameter combinations.
Monte Carlo estimates statistical performance.
Therefore,
\[ \boxed{ \text{Worst Case} \neq \text{Monte Carlo}. } \]Both may be required for safety-critical spacecraft functions.
4.1.7.87 Sensor Error Budget
A complete sensor error budget can contain:
| Error Source | Type | Typical Model | Calibration Possible? | Time Dependence |
|---|---|---|---|---|
| Bias | Additive | \(b\) | Yes | Constant / varying |
| Scale factor | Multiplicative | \(1+s\) | Yes | Slow |
| Misalignment | Geometric | \(C\) | Yes | Slow |
| White noise | Stochastic | \(n\) | No | Fast |
| Bias drift | Stochastic | \(b(t)\) | Partly | Slow |
| Quantisation | Digital | \(Q(z)\) | No | Sample dependent |
| Saturation | Nonlinear | \(\operatorname{sat}(z)\) | No | Input dependent |
| Temperature | Environmental | \(b(T),s(T)\) | Often | Slow |
| Latency | Timing | \(z(t-\tau)\) | Often | Fixed / varying |
| Dropout | Availability | \(valid=0\) | No | Event dependent |
4.1.7.88 Error Sources by Sensor Type
| Error Source | Gyro | Accelerometer | Magnetometer | Star Tracker | Sun Sensor | GNSS |
|---|---|---|---|---|---|---|
| Bias | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ |
| Scale factor | ✓ | ✓ | ✓ | —/limited | ✓ | — |
| Misalignment | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ |
| White noise | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ |
| Random walk | ✓ | ✓ | possible | — | — | clock/process dependent |
| Temperature | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ |
| Quantisation | ✓ | ✓ | ✓ | ✓ | ✓ | ✓ |
| Saturation | ✓ | ✓ | ✓ | detector | detector | RF/tracking |
| Environmental interference | vibration | vibration | magnetic | stray light | albedo/light | RF/multipath |
| Dropout | possible | possible | possible | common | eclipse | visibility |
| Timing error | ✓ | ✓ | ✓ | critical | relevant | critical |
The exact relevance and magnitude are sensor- and implementation-dependent.
4.1.7.89 Gyroscope Error Model
A useful complete model is
\[ \boxed{ \boldsymbol\omega_m = (I+S_g+M_g)\boldsymbol\omega + \mathbf b_g + K_{ga}\mathbf a + \mathbf n_g. } \]with
\[ \dot{\mathbf b}_g = -\frac1{\tau_g}\mathbf b_g + \mathbf w_{bg} \]or a random-walk model where appropriate.
Important gyro errors include:
- bias,
- angle random walk,
- bias instability,
- rate random walk,
- scale factor,
- non-orthogonality,
- g-sensitivity,
- temperature,
- saturation.
4.1.7.90 Accelerometer Error Model
A representative model is
\[ \boxed{ \mathbf f_m = (I+S_a+M_a)\mathbf f + \mathbf b_a + \mathbf n_a. } \]Additional effects may include:
- bias drift,
- vibration,
- temperature,
- nonlinearity,
- cross-axis sensitivity,
- saturation,
- lever arm.
4.1.7.91 Magnetometer Error Model
A representative model is
\[ \boxed{ \mathbf B_m = A_m\mathbf B_{true} + \mathbf b_h + \mathbf B_{spacecraft} + \mathbf n_m. } \]where \(A_m\) can represent soft-iron, scale-factor and alignment effects, while \(\mathbf b_h\) represents hard-iron offset.
Magnetic contamination is often as important as the intrinsic sensor error.
4.1.7.92 Star-Tracker Error Model
A high-level model is
\[ \boxed{ q_m = \delta q_{err} \otimes q_{true}(t-\tau). } \]The error quaternion can contain contributions from:
- centroiding,
- catalogue uncertainty,
- optics,
- alignment,
- thermal deformation,
- timing,
- motion blur,
- processing.
A validity flag should represent blinding and loss of tracking.
4.1.7.93 Sun-Sensor Error Model
For a vector-level Sun sensor,
\[ \boxed{ \mathbf s_m = C_{err}\mathbf s_{true} + \mathbf b_s + \mathbf n_s. } \]For photodiode-based coarse Sun sensors, errors can additionally arise from:
- cosine-response deviation,
- scale-factor mismatch,
- albedo,
- shadowing,
- saturation,
- temperature,
- dark current.
During eclipse,
\[ \boxed{ valid_{Sun}=0. } \]4.1.7.94 GNSS Error Model
At solution level,
\[ \boxed{ \mathbf r_m = \mathbf r + \mathbf b_r + \mathbf n_r } \]and
\[ \boxed{ \mathbf v_m = \mathbf v + \mathbf b_v + \mathbf n_v. } \]At raw-observable level,
\[ \boxed{ \rho_i = R_i + c(\delta t_r-\delta t_i) + I_i+T_i + \Delta\rho_{orbit} + \Delta\rho_{mp} + b_i + n_i. } \]For spacecraft applications, geometry, antenna orientation, timing and signal availability are essential components of the error model.
4.1.7.95 Sensor Error and Control-System Performance
Measurement noise enters the feedback loop.
Conceptually,
\[ \boxed{ \text{Sensor Noise} \rightarrow \text{Estimated State Noise} \rightarrow \text{Control Command Variation} \rightarrow \text{Actuator Activity}. } \]Excessive sensor noise may therefore increase:
- reaction-wheel activity,
- thruster firing,
- propellant consumption,
- jitter.
Filtering must balance noise rejection against latency and dynamic response.
4.1.7.96 Sensor Error and Guidance Performance
Guidance operates on estimated states rather than truth.
\[ \boxed{ \hat{\mathbf x} = \mathbf x + \delta\mathbf x. } \]Therefore,
\[ \mathbf u_G = G(\hat{\mathbf x}) \]differs from the ideal guidance command.
For RPO, relative-navigation uncertainty affects:
- approach corridor,
- braking,
- hold-point detection,
- collision avoidance,
- capture readiness.
4.1.7.97 Sensor Error and Autonomous Mode Logic
Sensor errors also influence discrete spacecraft decisions.
Examples include:
\[ \boxed{ \text{Sun Detected?} } \] \[ \boxed{ \text{Star Tracker Valid?} } \] \[ \boxed{ \text{GNSS Fix Valid?} } \] \[ \boxed{ \text{Capture Conditions Satisfied?} } \]Thus sensor uncertainty affects not only continuous estimation but also spacecraft state-machine transitions.
4.1.7.98 Requirement Allocation
Mission-level performance requirements should be flowed down to sensor requirements.
Conceptually,
\[ \boxed{ \text{Mission Requirement} \rightarrow \text{GNC Requirement} \rightarrow \text{Navigation Error Budget} \rightarrow \text{Sensor Requirement}. } \]For example,
\[ \sigma_{position,mission} \]may be allocated among:
- GNSS,
- IMU,
- timing,
- alignment,
- estimator,
- modelling errors.
Sensor selection should therefore begin from system requirements rather than datasheet comparison alone.
4.1.7.99 Calibration vs Estimation vs Filtering
These terms solve different problems.
Calibration
Determines systematic sensor parameters.
\[ \boxed{ \text{Calibration} \rightarrow b,S,M,\ldots } \]Estimation
Determines unknown spacecraft states and sometimes sensor errors.
\[ \boxed{ \text{Estimator} \rightarrow \hat x,\hat b,\ldots } \]Filtering
Reduces unwanted measurement variation and estimates state using models and measurements.
They complement one another.
4.1.7.100 Error Mitigation Strategies
Sensor errors can be addressed using several mechanisms:
- better sensor hardware,
- mechanical/thermal design,
- careful placement,
- electromagnetic cleanliness,
- preflight calibration,
- temperature compensation,
- in-orbit calibration,
- filtering,
- sensor fusion,
- bias estimation,
- redundancy,
- FDIR,
- measurement gating,
- adaptive covariance,
- operational constraints.
No single technique removes all sensor errors.
4.1.7.101 Sensor Error Simulation Architecture
A useful generic simulation chain is
Truth State
|
v
Ideal Measurement h(x)
|
v
Frame Transformation
|
v
Scale Factor
|
v
Misalignment / Cross-Axis Error
|
v
Nonlinearity
|
v
Bias
|
v
Bias Drift
|
v
Environmental Effects
|
v
White / Coloured Noise
|
v
Bandwidth / Sensor Dynamics
|
v
Saturation
|
v
Quantisation
|
v
Sampling
|
v
Latency / Timestamp
|
v
Dropout / Fault Logic
|
v
Measurement + Covariance + Validity
|
v
Navigation Filter
The exact ordering should reflect the physics of the particular sensor; this chain is a modelling framework rather than a universal physical sequence.
4.1.7.102 Simulink Implementation
A reusable Simulink sensor-error subsystem can contain:
Truth Signal
|
v
Coordinate Transform
|
v
Calibration Matrix
|
+------ Bias
|
+------ Bias Drift
|
+------ Temperature Model
|
+------ Noise
|
v
Sensor Dynamics
|
v
Saturation
|
v
Quantizer
|
v
Zero-Order Hold
|
v
Transport Delay
|
v
Validity / Fault Logic
|
v
Sensor Output
Outputs should ideally include
\[ \boxed{ \mathbf z_m,\quad R,\quad t_m,\quad valid,\quad status. } \]4.1.7.103 Low-, Medium- and High-Fidelity Models
Level 1 — Simple Measurement Model
\[ \boxed{ z_m=z_{true}+b+n. } \]Suitable for early controller development.
Level 2 — Engineering Sensor Model
\[ \boxed{ \mathbf z_m = (I+S+M)\mathbf z_{true} + \mathbf b(t,T) + \mathbf n. } \]Add:
- saturation,
- quantisation,
- sampling,
- latency,
- dropout.
Suitable for GNC design and Monte Carlo analysis.
Level 3 — Physics-Based Sensor Model
Model:
- sensing physics,
- environment,
- spacecraft geometry,
- electronics,
- signal processing,
- acquisition/tracking,
- detailed failure mechanisms.
Suitable for high-fidelity SIL/HIL and verification.
4.1.7.104 Sensor Datasheet Parameters
When reading a datasheet, do not look only at “accuracy.”
Important parameters include:
- measurement range,
- resolution,
- bias,
- bias repeatability,
- bias stability,
- noise density,
- random walk,
- scale-factor error,
- scale-factor stability,
- nonlinearity,
- cross-axis sensitivity,
- alignment accuracy,
- bandwidth,
- update rate,
- latency,
- temperature range,
- thermal coefficients,
- radiation tolerance,
- vibration/shock qualification,
- power,
- mass,
- interface,
- calibration information.
Sensor-specific pages then add parameters such as star-tracker exclusion angles or GNSS signal/tracking capabilities.
4.1.7.105 Datasheet Value vs Simulation Parameter
A datasheet specification is not always directly usable as a simulation standard deviation.
For example,
\[ \boxed{ \text{Noise Density} \neq \text{Per-Sample }\sigma } \]without knowing bandwidth and convention.
Similarly,
\[ \boxed{ \text{Maximum Error} \neq 1\sigma. } \]And
\[ \boxed{ \text{Bias Stability} \neq \text{Constant Bias}. } \]The statistical meaning of each specification must be understood before using it in Monte Carlo simulation.
4.1.7.106 Units and Conversions
Sensor error modelling is especially vulnerable to unit mistakes.
Examples include:
\[ ^\circ/s \leftrightarrow rad/s \] \[ \mu g \leftrightarrow m/s^2 \] \[ nT \leftrightarrow T \] \[ arcsec \leftrightarrow rad. \]For angular conversion,
\[ \boxed{ 1^\circ=\frac{\pi}{180}\;rad } \]and
\[ \boxed{ 1\;arcsec = \frac{\pi}{180\times3600}\;rad. } \]Every simulation parameter should therefore carry an explicit unit.
4.1.7.107 Error Correlation Time
The magnitude of an error is not enough to characterize it.
Two noise processes may have the same RMS but very different time behaviour.
Thus a stochastic model may require both
\[ \boxed{ \sigma } \]and
\[ \boxed{ \tau_c } \]where \(\tau_c\) is correlation time.
This strongly affects estimator performance.
4.1.7.108 Sensor Error State Machine
Some sensor errors are mode-dependent.
A sensor can transition through
\[ \boxed{ \text{OFF} \rightarrow \text{INITIALIZING} \rightarrow \text{ACQUIRING} \rightarrow \text{VALID} \rightarrow \text{DEGRADED} \rightarrow \text{FAILED}. } \]Examples:
- star tracker acquiring stars,
- GNSS waiting for fix,
- Sun sensor entering eclipse,
- sensor recovering after reset.
This is often more realistic than assuming continuous availability.
4.1.7.109 Adaptive Measurement Covariance
Sensor quality can change with operating conditions.
Instead of constant
\[ R=R_0, \]use
\[ \boxed{ R=R(\text{quality},\text{geometry},\text{environment}). } \]Examples:
- GNSS covariance increases with poor geometry,
- star-tracker covariance may increase during high slew,
- Sun-sensor accuracy may degrade near FOV limits.
This allows the estimator to weight measurements according to their actual quality.
4.1.7.110 Sensor Fault Injection
Verification should deliberately inject faults.
Examples include:
\[ b\rightarrow b+\Delta b \] \[ s\rightarrow s+\Delta s \] \[ valid\rightarrow0 \] \[ z\rightarrow z_{stuck} \] \[ z\rightarrow z+A_{spike} \] \[ \tau\rightarrow\tau+\Delta\tau. \]This allows verification of:
- estimator robustness,
- FDIR,
- safe mode,
- abort logic,
- redundancy switching.
4.1.7.111 Verification Against Truth
Simulation provides access to truth.
Therefore measurement error can be calculated as
\[ \boxed{ \mathbf e_k = \mathbf z_{m,k} - \mathbf z_{true,k}. } \]Useful statistics include:
\[ \bar{\mathbf e}, \] \[ \sigma_e, \] \[ RMSE, \] \[ e_{max}, \]and error PSD/Allan deviation where appropriate.
Estimator consistency can additionally be evaluated using covariance-normalized metrics.
4.1.7.112 Sensor Error Budget to Navigation Error Budget
Individual sensor specifications are not the final engineering objective.
The important chain is
\[ \boxed{ \text{Sensor Errors} } \] \[ \downarrow \] \[ \boxed{ \text{Measurement Covariance} } \] \[ \downarrow \] \[ \boxed{ \text{Navigation State Covariance} } \] \[ \downarrow \] \[ \boxed{ \text{Guidance / Control Uncertainty} } \] \[ \downarrow \] \[ \boxed{ \text{Mission Performance}. } \]This is the connection between sensor engineering and spacecraft GNC.
4.1.7.113 Complete Generalized Sensor Model
A useful high-level model combining the major error classes is
\[ \boxed{ \begin{aligned} \mathbf z_m(t_k) ={}& \mathcal Q \Big[ \mathcal D \big( A(T,t)\, C_{B}^{S}(t)\, h[\mathbf x(t_k-\tau)] \\ &+ \mathbf b(T,t) + \mathbf e_{env}(t) + \mathbf e_{dyn}(t) + \mathbf n(t) \big) \Big] \end{aligned} } \]subject to
\[ \boxed{ z_{min}\le z_m\le z_{max} } \]and
\[ \boxed{ valid_k\in\{0,1\}. } \]Here:
- \(h(\mathbf x)\) generates the ideal physical measurement,
- \(C_B^S\) represents sensor installation/frame transformation,
- \(A\) represents scale factor, non-orthogonality and cross-axis effects,
- \(\mathbf b\) represents bias,
- \(\mathbf e_{env}\) represents environmental errors,
- \(\mathbf e_{dyn}\) represents dynamic response effects,
- \(\mathbf n\) represents stochastic noise,
- \(\mathcal D\) represents sensor dynamics/sampling,
- \(\mathcal Q\) represents quantisation.
Not every sensor requires every term.
4.1.7.114 Complete Error-Management Chain
At spacecraft-system level:
\[ \boxed{ \text{Sensor Physics} \rightarrow \text{Error Characterization} \rightarrow \text{Calibration} \rightarrow \text{Measurement Model} } \] \[ \downarrow \] \[ \boxed{ \text{Covariance Modelling} \rightarrow \text{Sensor Fusion} \rightarrow \text{Innovation Monitoring} } \] \[ \downarrow \] \[ \boxed{ \text{FDIR} \rightarrow \text{Navigation Estimate} \rightarrow \text{Guidance} \rightarrow \text{Control}. } \]4.1.7.115 Final Engineering Perspective
The most important point of this chapter is:
\[ \boxed{ \text{Sensor Error} \neq \text{White Noise Only}. } \]A realistic spacecraft sensor can contain
\[ \boxed{ \text{Bias} + \text{Drift} + \text{Scale Factor} + \text{Misalignment} + \text{Nonlinearity} + \text{Noise} + \text{Environment} + \text{Timing} + \text{Digital Effects} + \text{Availability}. } \]These errors also behave differently:
\[ \boxed{ \text{Constant} \neq \text{Run-to-Run} \neq \text{Time-Varying} \neq \text{Fault}. } \]The complete spacecraft engineering philosophy should therefore be
\[ \boxed{ \text{Characterize} \rightarrow \text{Calibrate} \rightarrow \text{Model} \rightarrow \text{Estimate} \rightarrow \text{Monitor} \rightarrow \text{Mitigate}. } \]And ultimately,
\[ \boxed{ \text{Sensor Error Model} \rightarrow \text{Navigation Uncertainty} \rightarrow \text{GNC Performance} \rightarrow \text{Mission Safety and Success}. } \]This gives 4.1.7 a different and important role from the individual sensor pages: the gyro, accelerometer, magnetometer, star-tracker, Sun-sensor, and GNSS pages explain the errors of each instrument, while Sensor Error Sources explains the common mathematical and engineering framework connecting all of those errors to calibration, stochastic modelling, Kalman filtering, Monte Carlo analysis, FDIR, and spacecraft-level GNC performance.