Consider the following MIMO system Go(s) = [1-s/(s+1)² s+3/(s+1)(s+2) 1_s/(s+1)(s+2) s+4/(s+2)²] = Gon(s)[GoD(s)]^-1I (26.7.1) where GoN(s) = [(1_s)(s+2)² (s+1)(s+2)(s+3) (1-s)(s+2)(s+3) (s+1)(s+4)] (26.7.2) GoD(s) = (s + 1)²(s + 2)² (26.7.3) (i) Determine the location of RHP zeros and their

Question

Consider the following MIMO system

[latex]G_{o}(s)=\begin{bmatrix} \frac{1-s}{(s+1)^{2}} & \frac{s+3}{(s+1)(s+2)} \ \frac{1-s}{(s+1)(s+2)} & \frac{s+4}{(s+2)^{2}} \end{bmatrix} =G_{oN}(s)[G_{oD}(s)]^{-1}I[/latex] (26.7.1)

where

[latex]G_{oN}(s)=\begin{bmatrix} (1-s)(s+2)^{2} & (s+1)(s+2)(s+3) \ (1-s)(s+2)(s+3) & (s+1)^{2}(s+4) \end{bmatrix}[/latex] (26.7.2)

[latex]G_{oD}(s)=(s+1)^{2}(s+2)^{2}[/latex] (26.7.3)

(i) Determine the location of RHP zeros and their directions.
(ii) Evaluate the integral constraints on sensitivity that apply without enforcing dynamic decoupling and obtain bounds on the sensitivity peak.
(iii) Evaluate the integral constraints on sensitivity that apply if dynamic decoupling is required and obtain bounds on the sensitivity peak.
(iv) Compare the bounds obtained in parts (ii) and (iii).

Accepted Answer

(i) The zeros of the plant are the roots of [latex]det(G_{oN}(s))[/latex], i.e. the roots of [latex]−s^{6} − 11s^{5} − 43s^{4} − 63s^{3} + 74s + 44[/latex]. Only one of these roots, that located at s = 1, lies in the RHP. Thus [latex]z_{o} = 1[/latex], and

[latex]d\Omega (z_{o},\omega )=\frac{1}{1+\omega ^{2}}d[/latex] (26.7.4)

We then compute [latex]G_{o}(1)[/latex] as

[latex]G_{o}(1)=\begin{bmatrix} 0 & \frac{2}{3} \ 0 & \frac{5}{9} \end{bmatrix}[/latex] (26.7.5)

from where it can be seen that the dimension of the null space is [latex]µ_{z} = 1[/latex] and the (only) associated (left) direction is [latex]h^{T} = [5 − 6][/latex]. Clearly this vector has two nonzero elements, so we could expect that there will be additional design trade-offs arising from decoupling.

(ii) Applying Theorem 24.2 on page 759 part (ii) for r = 1 and r = 2 we obtain respectively

[latex]\frac{1}{\pi }\int_{\infty }^{\infty }{\ln \left|5[S_{o}(j\omega )]_{11}-6[S_{o}(j\omega )]_{21} ight| } \frac{1}{1+\omega ^{2}}d\omega \geq \ln \left|5 ight|[/latex] (26.7.6)

[latex]\frac{1}{\pi }\int_{\infty }^{\infty }{\ln \left|5[S_{o}(j\omega )]_{12}-6[S_{o}(j\omega )]_{22} ight| } \frac{1}{1+\omega ^{2}}d\omega \geq \ln \left|6 ight|[/latex] (26.7.7)

If we impose design requirements, as in Lemma 24.4 on page 760, we have, for the interacting MIMO loop, that

[latex]\left\|[S_{o}]_{11} ight\|_{\infty }+\frac{6}{5}\left\|[S_{o}]_{21} ight\|_{\infty } \geq (\frac{1}{\epsilon _{11}+\frac{6}{5} \epsilon _{21} } ) ^{\frac{\psi (\omega c)}{\pi -\psi (\omega c)} }[/latex] (26.7.8)

[latex]\left\|[S_{o}]_{22} ight\|_{\infty }+\frac{6}{5}\left\|[S_{o}]_{12} ight\|_{\infty } \geq (\frac{1}{\epsilon _{22}+\frac{6}{5}\epsilon _{21} } ) ^{\frac{\psi (\omega c)}{\pi -\psi (\omega c)} }[/latex] (26.7.9)

(iii) If we require dynamic decoupling, expressions (26.7.6) and (26.7.7) simplify, respectively to

[latex]\frac{1}{\pi } \int_{\infty }^{\infty }{\ln \left|[S_{o}(j\omega )]_{11} ight|\frac{1}{1+\omega ^{2}}d\omega } \geq 0[/latex] (26.7.10)

[latex]\frac{1}{\pi } \int_{\infty }^{\infty }{\ln \left|[S_{o}(j\omega )]_{22} ight|\frac{1}{1+\omega ^{2}}d\omega } \geq 0[/latex] (26.7.11)

And, with dynamic decoupling, (26.7.8) and (26.7.9), simplify to

[latex]\left\|[S_{o}]_{11} ight\|_{\infty } \geq (\frac{1}{\epsilon _{11}} )^{\frac{(\omega c)}{\pi -\psi (\omega c)} }[/latex] (26.7.12)

[latex]\left\|[S_{o}]_{22} ight\|_{\infty } \geq (\frac{1}{\epsilon _{22}} )^{\frac{(\omega c)}{\pi -\psi (\omega c)} }[/latex] (26.7.13)

(iv) To quantify the relationship between the magnitude of the bounds, in the coupled and the decoupled situations, we use an indicator [latex]κ_{1d}[/latex], formed as the quotient between the right hand sides of inequalities (26.7.8) and (26.7.12) , i.e.

[latex]κ_{1d}\overset{ riangle }{=} \left(1+\frac{6 }{5}\lambda _{1\epsilon } ight)^{-\frac{\psi (\omega c)}{\pi -\psi (\omega c)} }[/latex] where [latex]\lambda _{1\varepsilon }\overset{ riangle }{=} \frac{\epsilon _{21}}{\epsilon_{11}} [/latex] (26.7.14)

Thus, [latex]\lambda _{1\epsilon }[/latex] is a relative measure of interaction in the direction from channel 1 to channel 2.

The issues discussed above are captured in graphical form in Figure 26.10.

In Figure 26.10 we show a family of curves, each corresponding to a different bandwidth [latex]ω_{c}[/latex]. Each curve represents, for the specified bandwidth, the ratio between the bounds for the sensitivity peaks, as a function of the decoupling indicator, [latex]\lambda _{1\epsilon}[/latex]. We can summarize our main observations, as follows

a) When [latex]\lambda _{1\epsilon}[/latex] is very small, there is virtually no effect of channel 1 into channel 2 (at least in the frequency band [0, [latex]ω_{c}[/latex]]), then the bounds are very close ([latex]κ_{1d}[/latex] ≈ 1).

b) As [latex]\lambda _{1\epsilon}[/latex] increases, we are allowing the off diagonal sensitivity to become larger than the diagonal sensitivity in [0, [latex]ω_{c}[/latex]]. The effect of this manifests itself in [latex]κ_{1d}[/latex] < 1, i.e. in bounds for the sensitivity peak which are smaller than for the decoupled situation.

c) If we keep [latex]\lambda _{1\epsilon}[/latex] fixed, and we increase the bandwidth, then the advantages of using a coupled system also grow.

Also note that the left hand sides of (26.7.8) and (26.7.12) are different. In particular, (26.7.8) can be written as

[latex]\left\|[S_{o}]_{11} ight\|_{\infty } \geq \left(\frac{1}{\epsilon _{11}+\frac{6}{5}\epsilon _{21} } ight)^{\frac{\psi (\omega c)}{\pi -\psi (\omega c)} } -\frac{6}{5} \left\|[S_{o}]_{21} ight\|_{\infty }[/latex]

[latex]\geq \left(\frac{1}{\epsilon _{11}+\frac{6}{5}\epsilon_{21} } ight)^{\frac{\psi (\omega c)}{\pi -\psi (\omega c)} }-\frac{6}{5}\epsilon _{21}[/latex] (26.7.15)

Question 26.6: Consider the following MIMO system Go(s) = [1-s/(s+1)² s+3/(...

Related Answered Questions

Verified Answer:

Verified Answer:

Consider the plant Go(s) = Gon(s)[GoD(s)]^-1 (26.4.2) where GoN(s) = [-5 s² 1 -0.0023]; GoD(s) = [25s + 1 0 0 s(s + 1)²] (26.4.3) (i) Convert to state space form and evaluate the zeros. (ii) Design a pre-stabilizing controller to give static decoupling for reference signals. (iii) Design a ...

Verified Answer: