---
title: Polar RMS Error (PE) in Instrument Design
url: https://www.emergentmind.com/topics/polar-rms-error-pe
type: topic
---

# Polar RMS Error (PE) in Instrument Design

Polar RMS Error (PE) quantifies the root-mean-square deviation in the Stokes vector measured by a polarimetric system due to the cumulative effect of physical parameter uncertainties in its optical components. This scalar metric provides a comprehensive single-number figure of merit for polarimetric accuracy as propagated through the full modulation–demodulation pipeline, allowing direct comparison with system-level requirements and facilitating tolerance allocation and error budgeting within instrument design [1207.4241].

## 1. Mathematical Definition of Polar RMS Error

The polarimetric RMS error (PE) is rigorously defined with respect to a reference or “design” input Stokes vector, $S_{\text{in,def}}$. The measured vector, $S_{\text{in}}^*$, is related to the ideal through a total system response

$$
S_{\text{in}}^* = (I_4 + \Delta X) S_{\text{in,def}}
$$

where $\Delta X$ is a $4 \times 4$ matrix collecting all first-order perturbations from imperfect optics and systematics. The error vector is

$$
\epsilon = S_{\text{in}}^* - S_{\text{in,def}} = \Delta X S_{\text{in,def}}
$$

The squared norm is then

$$
\|\epsilon\|^2 = \sum_{k=1}^4 [\Delta X S_{\text{in,def}}]_k^2
$$

and the polarimetric RMS error is

$$
\text{PE} = \|\epsilon\| = \sqrt{\sum_{k=1}^4 [\Delta X S_{\text{in,def}}]_k^2}
$$

as introduced in Equation 20 of [1207.4241].

## 2. Modulation–Demodulation Framework and the Response Matrix

A polarimeter’s signal chain is expressed in terms of a modulation matrix $O$ (dimension $m \times 4$ for $m$ modulation states) and a corresponding demodulation matrix $D$ ($4 \times m$). Their product yields the system’s response matrix:

$$
X = D O
$$

Ideally, $X = I_4$, but in the presence of non-idealities, perturbations are captured in $\Delta X$:

$$
X = I_4 + \Delta X
$$

The measured Stokes vector after demodulation becomes $S_{\text{in}}^* = X S_{\text{in}}$. Perturbations impacting the response, including all optical-element-induced first-order variations and their propagation through the analysis chain, are integrated within $\Delta X$.

## 3. Error Propagation: From Physical Parameter Errors to PE

Each optical element $z$ is characterized by physical parameters $p_{j,z}$ with associated uncertainties $\delta p_{j,z}$. The Mueller matrix for an element is decomposed as

$$
M_z(p_{1,z}+\delta p_{1,z}, \ldots) = M_{z,\text{main}} + \sum_{j=1}^{j_{\max}(z)} \delta p_{j,z} \; m_{j,z}
$$

Assuming small, independent errors, the total instrument Mueller matrix is approximated by expanding to first order in each $\delta p_{j,z}$. The modulation matrix itself is likewise decomposed, and the demodulation matrix is computed as the Moore–Penrose pseudo-inverse of $O_{\text{main}}$ (Eq. 11):

$$
D^* = (O_{\text{main}}^T O_{\text{main}})^{-1} O_{\text{main}}^T
$$

Each perturbation term propagates to the response matrix via

$$
\Delta X = \sum_{z,j} \delta p_{j,z} \Delta X_{j,z}, \quad \text{with} \quad \Delta X_{j,z} = D^* o_{j,z}
$$

The overall error in the output Stokes vector is the sum of errors from individual physical parameter perturbations acting on the design input Stokes vector:

$$
\epsilon = \sum_{z,j} \delta p_{j,z} (\Delta X_{j,z} S_{\text{in,def}})
$$

Assuming mutual independence of the $\delta p_{j,z}$, the total RMS error is

$$
\text{PE}^2 = \sum_{k=1}^4 \left[ \sum_{z,j} \delta p_{j,z} (\Delta X_{j,z} S_{\text{in,def}})_k \right]^2
$$

For single-parameter contributions,

$$
\text{PE}_{j,z} = |\delta p_{j,z}| \cdot \| \Delta X_{j,z} S_{\text{in,def}} \|
$$

and the total combined error is

$$
\text{PE}_{\text{total}} = \sqrt{ \sum_{z,j} \text{PE}_{j,z}^2 }
$$

## 4. Underlying Assumptions and Approximations

The formalism underlying PE calculation is predicated on several explicit assumptions (as detailed in [1207.4241]):

- All physical-parameter errors $\delta p_{j,z}$ are sufficiently small for a first-order Taylor expansion to be accurate.
- The parameter uncertainties are statistically independent, justifying root-sum-square (RSS) combination.
- The demodulation step uses only the unperturbed modulation matrix $O_{\text{main}}$ in its pseudo-inverse; inclusion of first-order terms in the inversion would yield non-linear $\delta p$ dependence and is not treated.
- The linearized Mueller calculus is assumed adequate for the full error budget.

This framework thus delivers a tractable, analytically transparent mapping from low-level physical tolerances to the high-level system error metric.

## 5. Workflow Implementation in Polarimetric Systems

The M&m's code operationalizes this formalism in the following stepwise workflow:

1. The optical train is specified as an ordered list of $n$ elements, with associated physical parameters and uncertainties.
2. Element-wise, $M_{z,\text{main}}$ and the set of first-order weight matrices $m_{j,z}$ are provided.
3. The total system's “static” Mueller matrix is built via ordered multiplication of $M_{z,\text{main}}$ for each modulation state.
4. In parallel, each $m_{j,z}$ term is propagated as a first-order deviation, accounting for its correct position in the optical train.
5. First rows of each $M_{\text{tot},i}$ yield $O_{\text{main}}$ and the corresponding error-propagation matrices $o_{j,z}$.
6. Pseudo-inverse demodulation matrix $D^*$ is computed.
7. Each error contribution $\Delta X_{j,z} = D^* o_{j,z}$ is formed, and all $\delta p_{j,z} \Delta X_{j,z}$ are summed to yield $\Delta X$.
8. For the design input Stokes vector, the error vector $\epsilon$ and PE are computed as

   $$
   \epsilon = \Delta X S_{\text{in,def}}, \quad \text{PE} = \sqrt{\sum_k \epsilon_k^2}
   $$

This modular architecture enables both single-parameter sensitivity studies and holistic error budgeting for complex systems.

## 6. Role and Interpretation in System Performance Analysis

PE serves as a scalar measure of polarimetric accuracy for a given input Stokes state. It can be directly juxtaposed with mission or science-driven requirements (e.g., stipulations such as "we require $\le 10^{-3}$ of full-scale") [1207.4241]. Its utility extends to trade studies aimed at tolerance allocation: tighter requirements on $\delta p_{j,z}$ are imposed where their associated $\text{PE}_{j,z}$ contributions are dominant. Alternatively, analysis can be conducted over a suite of input states to determine the worst-case $\|\epsilon\|$ or by comparison of $\Delta X$ to a requirements matrix $\Delta X_{\text{req}}$.

This approach formalizes traceability from physical parameter tolerances through to system-level polarimetric fidelity, supporting both design iteration and requirements verification.

## 7. Comparison with Alternative Metrics and Broader Implications

PE offers a single-number summary of system error, suitable for both specification compliance and optimization. While equally valid to directly analyze the full $\Delta X$ matrix or to explore state-dependent error structure, practice as documented in [1207.4241] deems PE the preferred figure of merit for most trade studies. Its analytical structure—clear partitioning of sources, explicit error propagation, and the ability to decompose by element or mechanism—provides actionable insight into instrument design. A plausible implication is that this metric is well-suited for both initial design phases and ongoing assessment as tolerances or requirements evolve.

Source: https://www.emergentmind.com/topics/polar-rms-error-pe