PDF (letter size)

## Using Mason rule to obtain the transfer function from state space description

September 21, 2015   Compiled on October 15, 2023 at 4:40am

Given the observable form \begin {align*} x^{\prime } & =\begin {pmatrix} 0 & 0 & 0 & -\alpha _{0}\\ 1 & 0 & 0 & -\alpha _{1}\\ 0 & 1 & 0 & -\alpha _{2}\\ 0 & 0 & 1 & -\alpha _{3}\end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}\\ x_{3}\\ x_{4}\end {pmatrix} +\begin {pmatrix} \beta _{0}\\ \beta _{1}\\ \beta _{2}\\ \beta _{3}\end {pmatrix} u\\ y & =\begin {pmatrix} 0 & 0 & 0 & 1 \end {pmatrix}\begin {pmatrix} x_{1}\\ x_{2}\\ x_{3}\\ x_{4}\end {pmatrix} +\left [ \gamma \right ] u \end {align*}

Show that the transfer function is $$\frac {Y}{U}=\gamma +\frac {\beta _{3}s^{3}+\beta _{2}s^{2}+\beta _{1}s+\beta _{0}}{s^{4}+\alpha _{3}s^{3}+\alpha _{2}s^{2}+\alpha _{1}s+\alpha _{0}}$$ using Mason rule.

The transfer function can ofcourse be found using $$\frac {Y}{S}=C\left ( sI-A\right ) ^{-1}B+D$$ which gives\begin {align*} \frac {Y}{U} & =\begin {pmatrix} 0 & 0 & 0 & 1 \end {pmatrix} \left ( \begin {pmatrix} s & 0 & 0 & 0\\ 0 & s & 0 & 0\\ 0 & 1 & s & 0\\ 0 & 0 & 1 & s \end {pmatrix} -\begin {pmatrix} 0 & 0 & 0 & -\alpha _{0}\\ 1 & 0 & 0 & -\alpha _{1}\\ 0 & 1 & 0 & -\alpha _{2}\\ 0 & 0 & 1 & -\alpha _{3}\end {pmatrix} \right ) ^{-1}\begin {pmatrix} \beta _{0}\\ \beta _{1}\\ \beta _{2}\\ \beta _{3}\end {pmatrix} +\gamma \\ & =\gamma +\frac {\beta _{3}s^{3}+\beta _{2}s^{2}+\beta _{1}s+\beta _{0}}{s^{4}+\alpha _{3}s^{3}+\alpha _{2}s^{2}+\alpha _{1}s+\alpha _{0}} \end {align*}

But we want to use Mason rule here. The ﬁrst step is to write the equations so that the nodes variable are on the left side of the equation. The node variables are the states $$x_{1},x_{2},x_{3},x_{4}$$. From the matrix equations we obtain\begin {align*} x_{1}^{\prime } & =-\alpha _{0}x_{4}+\beta _{0}u\\ x_{2}^{\prime } & =x_{1}-\alpha _{1}x_{4}+\beta _{1}u\\ x_{3}^{\prime } & =x_{2}-\alpha _{2}x_{4}+\beta _{2}u\\ x_{4}^{\prime } & =x_{3}-\alpha _{3}x_{4}+\beta _{3}u \end {align*}

Solving for $$x$$ now, and in the process we change $$x^{\prime }$$ to $$sx$$ by taking Laplace transforms of all variables, and we add the output equation as well\begin {align*} X_{4} & =-\frac {s}{\alpha _{0}}X_{1}+\frac {\beta _{0}}{\alpha _{0}}U\\ X_{1} & =sX_{2}+\alpha _{1}X_{4}-\beta _{1}U\\ X_{2} & =sX_{3}+\alpha _{2}X_{4}-\beta _{2}U\\ X_{3} & =\left ( s+\alpha _{3}\right ) X_{4}-\beta _{3}U\\ Y & =X_{4}+\gamma U \end {align*}

We now draw the Mason diagram, putting $$U$$ on the left most node and $$Y$$ on the right most node (the input and output). Here is the result

The forward paths from $$U$$ to $$Y$$ are\begin {align*} F_{1} & =\frac {\beta _{0}}{\alpha _{0}}\\ F_{2} & =\gamma \\ F_{3} & =\left ( -\beta _{2}\right ) \left ( s\right ) \left ( -\frac {s}{\alpha _{0}}\right ) \\ F_{4} & =\left ( -\beta _{1}\right ) \left ( -\frac {s}{\alpha _{0}}\right ) \\ F_{5} & =\left ( -\beta _{3}\right ) \left ( s\right ) \left ( s\right ) \left ( -\frac {s}{\alpha _{0}}\right ) \end {align*}
We now need to calculate the associated $$\Delta _{k}$$ for each of the above forward loops. $$\Delta _{k}$$ is found by removing $$F_{k}$$ from the graph and then calculating the main mason $$\Delta$$ of what is left in the graph. When remove $$F_{1}$$ no loops remain, hence $$\Delta _{1}=1$$. When removing $$F_{2}$$ all the loops remain, hence $$\Delta _{2}=1-\left ( L_{1}+L_{2}+L_{3}\right )$$, and when removing $$F_{3},F_{4},F_{5}$$ no loops remain, hence $$\Delta _{3}=\Delta _{4}=\Delta _{5}=1$$. Therefore$\frac {Y}{U}=\frac {F_{1}\Delta _{1}+F_{2}\Delta _{2}+F_{3}\Delta _{3}+F_{4}\Delta _{4}+F_{5}\Delta _{5}}{1-\left ( L_{1}+L_{2}+L_{3}\right ) }$ There is no other combinations of loops. The above becomes\begin {align*} \frac {Y}{U} & =\frac {F_{1}+\gamma \left ( 1-\left ( L_{1}+L_{2}+L_{3}\right ) \right ) +F_{3}+F_{4}+F_{5}}{1-\left ( L_{1}+L_{2}+L_{3}\right ) }\\ & =\gamma +\frac {F_{1}+F_{3}+F_{4}+F_{5}}{1-\left ( L_{1}+L_{2}+L_{3}\right ) }\\ & =\gamma +\frac {\frac {\beta _{0}}{\alpha _{0}}+\left ( -\beta _{2}\right ) \left ( s\right ) \left ( -\frac {s}{\alpha _{0}}\right ) +\left ( -\beta _{1}\right ) \left ( -\frac {s}{\alpha _{0}}\right ) +\left ( -\beta _{3}\right ) \left ( s\right ) \left ( s\right ) \left ( -\frac {s}{\alpha _{0}}\right ) }{1-\left ( \left ( -\frac {s}{\alpha _{0}}\right ) \left ( \alpha _{2}\right ) \left ( s\right ) +\left ( -\frac {s}{\alpha _{0}}\right ) \left ( s+\alpha _{3}\right ) \left ( s\right ) \left ( s\right ) +\left ( -\frac {s}{\alpha _{0}}\right ) \left ( \alpha _{1}\right ) \right ) }\\ & =\gamma +\frac {\frac {\beta _{0}}{\alpha _{0}}+\beta _{2}\frac {s^{2}}{\alpha _{0}}+\beta _{1}\frac {s}{\alpha _{0}}+\beta _{3}\frac {s^{3}}{\alpha _{0}}}{1-\left ( -\frac {s^{2}}{\alpha _{0}}\alpha _{2}-\frac {s^{2}}{\alpha _{0}}\left ( s+\alpha _{3}\right ) -\frac {s}{\alpha _{0}}\alpha _{1}\right ) }\\ & =\gamma +\frac {\frac {\beta _{3}s^{3}+\beta _{2}s^{2}+\beta _{1}s+\beta _{0}}{\alpha _{0}}}{1+\frac {s^{3}+s^{2}\alpha _{3}+s^{2}\alpha _{2}+s\alpha _{1}}{\alpha _{0}}} \end {align*}
Hence$G\left ( s\right ) =\frac {Y}{U}=\gamma +\frac {\beta _3s^3+\beta _2s^2+\beta _1s+\beta _0}{s^3+s^2\alpha _3+s^2\alpha _2+s\alpha _1+\alpha _0}$