The system

$$\begin{align*} \dot{x}_1 & = x_4^2 + x_3^3 + u_1 + au_2 \\ \dot{x}_2 & = x_3 \\ \dot{x}_3 & = \sin{x_4}+\cos{x_1}+bu_1 + u_2 \\ \dot{x}_4 & = -x_4 \\ y_1 &= x_1 \\ y_2 &=x_2. \end{align*}$$

The system

$$\begin{align*} \dot{x}_1 & = x_2 \\ \dot{x}_2 & = x_3^2 + x_4 + u_1 + au_2 \\ \dot{x}_3 & = x_4 + bu_1 + u_2 \\ \dot{x}_4 & = -x_4 \\ y_1 &= x_1 \\ y_2 &=x_3 \end{align*}$$

with $ab \neq 1$, has a well-defined vector relative degree (2, 1) and a nonsingular decoupling matrix $\begin{bmatrix} 1 & a \\ b & 1 \end{bmatrix}$.

The model of the MAGLEV system is unstable and nonlinear

$$ m\ddot{x}=mg-\dfrac{K_{c}V^{2}}{x^{2}}, $$

where $x$ is the metal ball position being the system output, $V$ is the system input as the voltage. Other parameters are $m$ as the mass of the metal ball, $K_c$ as constant for magnet circuit, and $g$ is the gravitational acceleration of 9.8 m/s$^2$. A free-body diagram is shown also in the attached image.

The behavior of phytoplankton cells in a continuous reactor is usually described by the Droop model. Cell growth is limited by a nutrient with concentration $S$. The biomass has a concentration $N$ and $Q$ represents the cell quota of assimilated nutrient, expressed as the amount of intracellular nutrient per biomass unit. The dilution rate $D$ corresponds to the flow rate of renewal medium over the volume of the reactor, and $D$ is the input of the system.

We denote $D = D_0 + u$, and the system fits

$$\sum_D \begin{cases} \dot{x}_i = f(x) + ug(x)\\ y=h(x_1) \end{cases}$$

with

$f(x)=\begin{pmatrix} a_2\left(1-\dfrac{1}{x_2}\right)x_1 - D_0x_1\\ a_3\dfrac{x_3}{a_1+x_3} - a_2(x_2 - 1)\\ D_0(1-x_3)-\dfrac{x_1x_3}{a_1+x_3} \end{pmatrix}$

$g(x)=\begin{pmatrix} -x_1\\ 0\\ 1-x_3 \end{pmatrix}$, and $h(x_1)=x_1$, where

$ x_1 = (\rho_m N/S_i);\\ x_2 = (Q/K_Q);\\ x_3 = (S/S_i);\\ a_1 = (K_{\rho}/S_i);\\ a_2 = \mu_m;\\ a_3 = (\rho_m/K_Q). $

The transfer function of the three-mass-system is much more complex than it is for one dominant elasticity (two-mass-system).

$${G_{\rm mech}}(s) = \underbrace{{\dfrac{1} {T_{ \Sigma} \cdot s}}}_{G_{\rm rs}(s)} \cdot \underbrace{ \dfrac{ a_{7} \cdot s^{4} + a_{6} \cdot s^{3}+a_{5} \cdot s^{2} + a_{4} \cdot s + 1}{a_{3} \cdot s^{4} + a_{2} \cdot s^{3}+ a_{1} \cdot s^{2} + a_{4} \cdot s + 1}} _{G_{\rm nrs}(s)}$$

with

$T_{\Sigma} = T_{\rm M} + T_{{\rm L}1} + T_{{\rm L}2}$

and

$\begin{align*} a_{1}&=d_{1}d_{2}T_{{\rm C}1}T_{{\rm C}2}+T_{{\rm L}2}\left(T_{\rm M}+T_{{ \rm L}1}\right) \cdot \frac{T_{C2}}{T_{\Sigma }}+T_{\rm M}\left(T_{{\rm L}1}+T_{{\rm L}2}\right) \cdot \frac{T_{{\rm C}1}}{T_{\Sigma}} \\ a_{2}&=\frac{T_{{\rm C}1}T_{{\rm C}2}}{T_{\Sigma}}\cdot\left(d_{1}T_{{\rm L}2}\left(T_{\rm M}+T_{{\rm L}1}\right)+d_{2}T_{\rm M} \left(T_{{\rm L}1}+T_{{\rm L}2}\right)\right) \\ a_{3}&=\frac{T_{\rm M}T_{{\rm L}1}T_{{\rm L}2}T_{{\rm C}1}T_{{\rm C}2}}{T_{\Sigma}} \\ a_{4}&=d_{1}T_{{\rm C}1}+d_{2}T_{{\rm C}2} \\ a_{5}&=d_{1}d_{2}T_{{\rm C}1}T_{{\rm C}2}+T_{{\rm L}2}T_{{\rm C}2}+\left(T_{{\rm L}1}+T_{{\rm L}2}\right)\cdot T_{{\rm C}1} \\ a_{6}&=\left(\left(d_{1}+d_{2}\right)T_{{\rm L}2}+d_{2}T_{{\rm L}1}\right)\cdot T_{{\rm C}1}T_{{\rm C}2} \\ a_{7}&=T_{{\rm L}1}T_{{\rm L}2}T_{{\rm C}1}T_{{\rm C}2}. \end{align*}$