Synthetic parametric model

Introduction

On this page you will find a purely synthetic parametric model. The goal is to have a simple parametric model which one can use to experiment with different system orders, parameter values etc.

System description

The parameter $\varepsilon$ scales the real part of the system poles, that is, $p_i=\varepsilon a_i+jb_i$. For a system in pole-residue form

$H(s,\varepsilon) = \sum_{i=1}^{n}\frac{r_i}{s-p_i} = \sum_{i=1}^{n}\frac{r_i}{s-(\varepsilon a_i+jb_i)}$

we can write down the state-space realisation $H(s,\varepsilon) = \widehat{C}\Big(sI-\varepsilon \widehat{A}_\varepsilon - \widehat{A}_0\Big)^{-1}\widehat{B}+D$ with

$\varepsilon \widehat{A}_\varepsilon + \widehat{A}_0 = \varepsilon \left[\begin{array}{ccc} a_1 & & \\ & \ddots & \\ & & a_n\end{array}\right] +\left[\begin{array}{ccc} jb_1 & & \\ & \ddots & \\ & & jb_n\end{array}\right] ,$

$\widehat{B} = [1,\ldots,1]^T,\quad \widehat{C} = [r_1,\ldots,r_n],\quad D = 0.$

Notice that the system matrices have complex entries.

For simplicity, assume that $n$ is even, $n=2k$, and that all system poles are complex and ordered in complex conjugate pairs, i.e.

$p_1 = \varepsilon a_1+jb_1, p_2 = \varepsilon a_1-jb_1, \ldots, p_{n-1} = \varepsilon a_k+jb_k, p_n = \varepsilon a_k-jb_k,$

and the residues also form complex conjugate pairs

$r_1 = c_1+jd_1, r_2 = c_1-jd_1, \ldots, r_{n-1} = c_k+jd_k, r_n = c_k-jd_k.$

Then a realization with matrices having real entries is given by

$A_\varepsilon = \left[\begin{array}{ccc} A_\varepsilon^1 & & \\ & \ddots & \\ & & A_\varepsilon^k\end{array}\right], \quad A_0 = \left[\begin{array}{ccc} A_0^1 & & \\ & \ddots & \\ & & A_0^k\end{array}\right], \quad B = T\widehat{B}, \quad C = \widehat{C}T^*, \quad D = 0,$

with $T = \left[\begin{array}{ccc} T_0 & & \\ & \ddots & \\ & & T_0 \end{array}\right]$ and $T_0 = \frac{1}{\sqrt{2}}\left[\begin{array}{cc} 1 & -j\\ 1 & j \end{array}\right]$.

Numerical values

We construct a system of order $n = 100$. The numerical values for the different variables are

• $r_i = 1$,

• $a_i$ equally spaced in $[10, 10^3]$,

• $b_i$ equally spaced in $[10, 10^3]$,

• $\varepsilon$$\in [1/50,1]$.

In MATLAB the system matrices are easily formed as follows

n = 100;

a = -linspace(1e1,1e3,n/2);

b = linspace(1e1,1e3,n/2);