digital-brain/content/zettels/norms.md

8.3 KiB

+++ title = "Systems and Signals Norms" author = ["Thomas Dehaeze"] draft = false +++

Tags :

\[ \SI{1}{\meter\per\second} \]

Resources:

Definition

A norm of \(e\) (which may be a vector, matrix, signal of system) is a real number, denoted \(\|e\|\), that satisfies the following properties:

  1. Non-negative: \(\|e\| \ge 0\)
  2. Positive: \(\|e\| = 0 \Longleftrightarrow e = 0\)
  3. Homogeneous: \(\|\alpha \cdot e\| = |\alpha| \cdot \|e\|\) for all complex scalars \(\alpha\)
  4. Triangle inequality: \(\|e_1 + e_2\| \le \|e_1\| + \|e_2\|\)

Vector Norms

  • Vector 1-norm (Sum Norm): \[ \|a\|_1 \triangleq \sum_i |a_i| \]
  • Vector 2-norm (Euclidean Norm): \[ \|a\|_2 \triangleq \sqrt{\sum_i |a_i|^2} \]
  • Vector p-norm: \[ \|a\|_p \triangleq \left( \sum_i |a_i|^p \right)^{1/p} \]
  • Vector \(\infty\text{-norm}\) (Max Norm): () \[ \|a\|_\infty \triangleq \max_i |a_i| \]

Matrix Norms

A norm on a matrix \(\|A\|\) is a matrix norm if, in addition to the four norm properties, it also satisfies the multiplicative property: \[ \|AB\| \le \|A\| \cdot \|B\| \]

  • Sum matrix norm: \[ \|A\|_\text{sum} \triangleq \sum_{i,j} |a_{ij}| \]
  • Frobenius matrix norm (Euclidean Norm): \[ \|A\|_F \triangleq \sqrt{\sum_{i,j} |a_{ij}|^2} = \sqrt{\text{tr}(A^H A)} \]
  • Max element norm: (which is not a matrix norm) \[ \|A\|_\text{max} \triangleq \max_{i,j} |a_{ij}| \]

Induced Matrix Norms

Induced matrix norms are important because of their close relationship to signal amplification in systems.

Consider the figure below where \(w\) is the input vector, \(z\) the output vector and where the "amplification" or "gain" of the matrix \(A\) is defined by the ration \(\|z\|/\|w\|\).

{{< figure src="/ox-hugo/induced_matrix_norm.png" >}}

The maximum gain for all possible input directions is given by the induced norm: \[ \|A\|_{ip} \triangleq \max_{w \neq 0} \frac{\|Aw\|_p}{\|w\|_p} \]

Thus, the induced norm gives the largest possible "amplification" of the matrix. The following equivalent definition is also used: \[ \|A\|_{ip} = \max_{\|w\|_p \le 1} \|Aw\|_p \]

Signal Norms

For signals, we may compute the norm in two steps:

  1. "Sum up" the channels at a given time using a vector norm. For a scalar, we simply take the absolute value.
  2. "Sum up" in time using a temporal norm.

We normally use the same p-norm both for the vector and the signal.

  • 1-norm in time (Integral Absolute Error): \[ \|e(t)\|_1 = \int_{-\infty}^{\infty} \sum_i |e_i(\tau)| d\tau \]
  • 2-norm in time (Quadratic Norm): \[ \|e(t)\|_2 = \sqrt{\int_{-\infty}^{\infty} \sum_i |e_i(\tau)|^2 d\tau} \]
  • \(\infty\text{-norm}\) in time (Peak value in time): \[ \|e(t)\|_\infty = \max_\tau \left( \max_i |e_i(\tau)| \right) \]
  • Power-Norm or RMS-Norm: \[ \|e(t)\|_\text{pow} = \lim_{T\to \infty} \sqrt{\frac{1}{2T} \int_{-T}^T \sum_i |e_i(\tau)|^2 d\tau} \]

Signal Interpretation of Various System Norms

Consider a system \(G\) with input \(d\) and output \(e\), such that: \[ e = G d \]

For performance, we may want the output signal \(e\) to be "small" for any allowed input signals \(d\). We therefore need to specify:

  1. What \(d\) are allowed. (Which set does \(d\) belong to?) Some possible inputs signal sets are:
    • \(d(t)\) consists of impulses \(\delta(t)\).
    • These generate step changes in the states.
    • \(d(t) = \sin(\omega t)\) with fixed frequency
    • \(d(t)\) is bounded in energy \(\|d(t)\|_2 \le 1\)
    • \(d(t)\) is bounded in power \(\|d(t)\|_\text{pow} \le 1\)
    • \(d(t)\) is bounded in magnitude \(\|d(t)\|_\infty \le 1\)
  2. What we mean by "small". (Which norm should be use for \(e\)?) To measure the output signal, we may consider the following norms:
    • 2-norm (energy): \(\|e(t)\|_2\)
    • \(\infty\text{-norm}\) (peak magnitude): \(\|e(t)\|_\infty\)
    • Power: \(\|e(t)\|_\text{pow}\)

We now consider which system norms result from the definition of input classes and output norms (Table 1).

Table 1: System norms for sets of inputs signals and three different output norms
\(\delta(t)\) \(\sin(\omega t)\) \(\vert\vert d \vert\vert_2\) \(\vert\vert d \vert\vert_\infty\) \(\vert\vert d \vert\vert_\text{pow}\)
\(\vert\vert e \vert\vert_2\) \(\vert\vert G(s) \vert\vert_2\) \(\infty\) \(\vert\vert G(s) \vert\vert_\infty\) \(\infty\) \(\infty\)
\(\vert\vert e \vert\vert_\infty\) \(\vert\vert g(t) \vert\vert_\infty\) \(\overline{\sigma}(G(j\omega))\) \(\vert\vert G(s) \vert\vert_2\) \(\vert\vert g(t) \vert\vert_1\) \(\infty\)
\(\vert\vert e \vert\vert_\text{pow}\) 0 \(\frac{1}{\sqrt{2}} \overline{\sigma}(G(j\omega))\) 0 \(\le \vert\vert G(s) \vert\vert_\infty\) \(\vert\vert G(s) \vert\vert_\infty\)

System Norms

\(\mathcal{H}_\infty\) Norm

Consider a proper linear stable system \(G(s)\). The \(\mathcal{H}_\infty\) norm is the peak value of its maximum singular value: \[ \|G(s)\|_\infty \triangleq \max_{\omega} \overline{\sigma}(G(j\omega)) \]

In terms of signals, the \(\mathcal{H}_\infty\) norm can be interpreted as follows:

  • it is the worst case steady-state gain for sinusoidal inputs at any frequency
  • it is equal to the 2-norm in the time domain: \[ \|G(s)\|_\infty = \max_{d(t)} \frac{\|e(t)\|_2 \neq 0}{\|d(t)\|_2} = \max_{\|d(t)\|_2 = 1} \|e(t)\|_2 \]

\(\mathcal{H}_2\) Norm

Consider a strictly proper system \(G(s)\). The \(\mathcal{H}_2\) norm is:

\begin{align*} \|G(s)\|_2 &\triangleq \sqrt{\frac{1}{2\pi} \int_{-\infty}^{\infty} \text{tr}\left(G(j\omega)^HG(j\omega)\right) d\omega} \\\ &= \sqrt{\frac{1}{2\pi} \int_{-\infty}^{\infty} \sum_i {\sigma_i}^2(G(j\omega)) d\omega} \end{align*}

In terms of signals, the \(\mathcal{H}_\infty\) norm can be interpreted as follows:

  • it is a measure of the expected RMS value of the output to white noise excitation

The \(\mathcal{H}_2\) is very useful when combined to [Dynamic Error Budgeting]({{< relref "dynamic_error_budgeting" >}}).

As explained in (Monkhorst 2004), the \(\mathcal{H}_2\) norm has a stochastic interpretation:

The squared \(\mathcal{H}_2\) norm can be interpreted as the output variance of a system with zero mean white noise input.

Bibliography

Monkhorst, Wouter. 2004. “Dynamic Error Budgeting, a Design Approach.” Delft University.

Skogestad, Sigurd, and Ian Postlethwaite. 2007. Multivariable Feedback Control: Analysis and Design. John Wiley.

Toivonen, Hannu T. 2002. “Robust Control Methods.” Abo Akademi University.

Zhang, Weidong. 2011. Quantitative Process Control Theory. CRC Press.