The dynamic engine model, with parameters for a 1.6 liter, 4-cylinder fuel injected engine, is a two-input/two-output system, given by the following differential equations: $\dot{P}=k_P(\dot{m}_{ai}-\dot{m}_{ao}),$ where $k_p=42.40$

$\dot{N}=k_N(T_i-T_L),$ where $k_N=54.26$

$\dot{m}_{ai}=(1+0.907\theta+0.0998\theta^2)g(P)$

$g{P}= \begin{cases} 1, & P<50.6625 \\ 0.0197(101.325P - P^2)^{\frac{1}{2}}, &P \geq 50.6625 \end{cases}$

$\dot{m}_{ao} = -0.0005968N-0.1336P+0.0005341NP+0.000001757NP^2$

$m_{ao} = \dot{m}_{ao}(t-\tau)/(120N), \ \tau=45/N$

$T_i = -39.22+325024m_{ao}-0.0112\delta^2+0.635\delta+(0.0216+0.000675\delta)N(2\pi/60)-0.000102N^2(2\pi/60)^2$

$T_L = (N/263.17)^2+T_d$.

Coupled nonlinear differential equations describing a process involving a continuous flow stirred tank reactor are given by

$$\begin{align*} \dot{C}_1 &= -C_1u + C_1(1-C_2)e^{C_2/\Gamma} \\ \dot{C}_2 &= -C_2u + C_1(1-C_2)e^{C_2/\Gamma}\dfrac{1+\beta}{1+\beta-C_2}. \end{align*}$$

In these equations, the state variables $C_1$ and $C_2$ represent dimensionless forms of cell mass and amount of nutrients in a constant volume tank, bounded between zero and unity. The control $u$ is the flow rate of nutrients into the tank, and is the same rate at which contents are removed from the tank. The constant parameters $\Gamma$ and $\beta$ determine the rates of cell formation and nutrient consumption; these parameters are set to $\Gamma$= 0.48 and $\beta$ = 0.02 for the nominal benchmark specification.

The Euler-Lagrange equations of motion of the system are given as follows

$$\begin{pmatrix} (M+m) & mL\cos{q_1} & 0\\ mL\cos{q_1} & mL^2+\Theta & \dfrac{\Theta}{2}\\ 0 & \dfrac{\Theta}{2} & \Theta \end{pmatrix} \begin{pmatrix} \ddot{x}\\ \ddot{q_1}\\ \ddot{q_2} \end{pmatrix} + \begin{pmatrix} -mL\sin{q_1\dot{q_1}^2}\\ -mLg_g\sin{q_1}\\ 0 \end{pmatrix} = \begin{pmatrix} Q_x\\ Q_1\\ Q_2 \end{pmatrix},$$

where $Q_x (N)$ is the generalized force pushing the cart in the horizontal “$x$” direction, $Q_1$ and $Q_2$ are torques in $(N · m)$ rotating the beam of the crane around a horizontal axis orthogonal to “$x$” and counter-rotating the hamper at the free end of the beam to avoid turning out the worker from the hamper. $L (m)$ denotes the lenght of the crane’s beam, $g_g$ ($m/s^2$) is the gravitational acceleration, $m$ ($kg$) and $\Theta$ $(kg · m^2)$ denote the momentum (with respect to its own center of mass that was supposed to be on the rotational axle) and the mass of the hamper.

Image below shows the block diagram of a discrete-time system.

$$\begin{align*} H_1(z) &=\dfrac{0.2z^{-1}}{z^{-1}-0.21z^{-2}} \\ H_2(z) &=\dfrac{0.1z^{-1}}{1-1.1z^{-1}+0.3z^{-2}} \\ H_3(z) &=\dfrac{0.3z^{-1}}{1-0.4z^{-1}} \end{align*}$$

The block diagram of a third-order nonlinear discrete time system adopted by Fakhouri for identification evaluation is shown below.

$$\begin{align*} H_1(z) &=\dfrac{0.1z^{-1}}{1-0.5z^{-1}} \\ H_2(z) &=\dfrac{0.1z^{-1}}{1-1.3z^{-1}+0.42z^{-2}} \\ H_3(z) &=\dfrac{1.0z^{-1}}{1-0.7z^{-1}} \end{align*}$$