Abstract
This article refers to practical and application-oriented design aspects arising in the field of chordal automation where stereographic projection is used. Several actions are discussed when high-order systems are involved. For average-order systems there exist elegant analytically concise methods to present and solve the problems. The chordal stability margin design is extended to pole assignment conditions, gradient analysis and actuation signal limits, also adequate for high-order systems.
Zusammenfassung
Dieser Beitrag befasst sich mit praktischen und anwendungsorientierten Aspekten der Auslegung im Bereich der chordalen Automation beim Einsatz von stereographischer Projektion. Es werden mehrere Prozesse behandelt, an denen Systeme höherer Ordnung beteiligt sind. Für Systeme durchschnittlicher Ordnung gibt es elegante, analytisch prägnante Methoden, die Probleme darzustellen und zu lösen. Die Auslegung des chordalen Stabilitätsrandes wird ausgeweitet auf die Bedingungen der Eigenwertvorgabe, die Analyse mit Gradienten und die Begrenzung der Stellgröße; Methoden, welche auch für Systeme höherer Ordnung anwendbar sind.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
An automatic control setup, given as a single-loop linear control system with a two-valued input and output, is considered in detail, see Fig. 1. The dynamic property results from the closed-loop poles \(s\) of the transfer functions, i.e., of \(1 + G(s)K(s) =0\) or \(G(s) = -K(s)^{-1}\). If \(s =j\omega\) was fulfilled by any solution, the system would stay in an unwanted steady-state oscillation (if there was no otherwise solution \(\Re e \,s>0\)). That is, a necessary condition, to simplify matters, reads as \(G( j\omega)\neq -K^{-1}(j\omega)\), \(G(j\omega)\) and \(-K( j\omega)^{-1}\) must be distinctly separated, \(\min_{\omega} |G(j\omega) + K^{-1}( j\omega)|>0\). This makes clear that an interrelation of \(G\) and the inverse \(K\) is responsible for the closed-loop performance. To consider inverses of proper transfer functions is not common with reference to frequency plots routing to infinity. But the stereographic projections keeps it in the usual range.
From Fig. 1 one has
Now, the Riemann sphere is considered. In spite of the analytical effort, relations on a sphere are well suited for requirements in automatic control [2, 3]. The Riemann sphere is situated at the complex plane of frequency plots with its south pole at the origin. Using stereographic projection onto the Riemann sphere, two points \(G_{1}\) and \(G_{2}\) of a frequency plot in the complex plane are mapped to \(\mathbf{g}_{1}\) and \(\mathbf{g}_{2}\) on the Riemann sphere via an arrow pointing to the north pole, i.e. map \(G\mapsto \mathbf{g}\) corresponds with \(\mathbf{g} = [ \Re e\, G, \Im m\,G, |G|^{2} /(1 +| G|^{2})]\).
The straight distance inside the sphere, Fig. 2, results from basic algebra
For automatic control purposes, the chordal distance \(\rho[G, K]\) is the inverse of the maximum singular value of M (conditionally weighted by \(w_{ik}\)).
Comparing Eq. (7) and Eq. (45), the extension of Eq. (8), and Eq. (36) is important, especially the correspondences \(G_{1}\), \(G_{2}\) and \(-K^{-1}\), \(G\), respectively. They are used for purpose in automatic control
The chordal stability margin is defined by the minimum versus frequency
Presupposing stability of the nominal system of Fig. 1, a given family of uncertain (perturbed) plants \(G_{p}\) is robustly stable if [1–3]
This relation has been slightly simplified. The chordal stability margin must be bigger than the chordal uncertainty radius \(r_{G}\). Using the definitions Eq. (46) and Eq. (38) of the Appendix A,
For illustration, an example for \(K(s) = (1 + 2s)/(2 + s)\), \(G(s) = 1/s^{2}\) is given in Figs. 3 and 4.
Now, the main part of this article aims at application-oriented problems characterized especially by high-order systems. The chordal stability margin is extended by a gradient method. An optimization uses eigenvalue assignment and actuation effort conditions. The unstructured chordal stability margin is also combined with the unstructured design method using the inverse trace of system matrix in state space.
2 Gradient of the chordal stability margin \(b_{GK}\)
Define the polynomial setup of the controller \(K(j \omega)=\frac{Z(j\omega)}{N(j\omega)}\), \(N(j\omega)=\mathbf{b}^{T}\boldsymbol{\omega}\), \(Z(j\omega)= \mathbf{a}^{T}\boldsymbol{\omega}\), \(\boldsymbol{\omega}\triangleq (1\ j\omega\ -\omega^{2}\ - j\omega^{3}\ \ldots) ^{T} \in \mathcal{C}^{n}\), \(\omega\in \mathcal{R}^{1}\), \(\boldsymbol{\omega}^{H}\equiv\boldsymbol{\omega}^{\star T}\)
\(\boldsymbol{\Omega}_{r} \triangleq\Re e\,\boldsymbol{\Omega}\), \(\boldsymbol{\Omega}_{i}\triangleq \Im m\,\boldsymbol{\Omega}\), \(|N|^{2} =\mathbf{b}^{T}\boldsymbol{\Omega} \mathbf{b}\), \(|Z|^{2}=\mathbf{a}^{T}\boldsymbol{\Omega} \mathbf{a}\).
where
With \(b_{i}\) as one of the coefficients of \(\mathbf{b}\) in the denominator polynomial, responsible for the poles of the controller, and \(\mathbf{e}_{i}\) the unit vector
Changing the order of operations one has
An extension to vector-valued gradients \(\frac{\partial}{\partial \mathbf{b}}\) or \(\frac{\partial}{\partial \mathbf{a}}\) for uncertain coefficients of \(G\) or \(K\) can easily be stated, irrespective of the system order.
3 Eigenvalue assignment as a condition
Consider a controller \(K\) which order is selected intentionally higher, using a surplus of e.g. \(c_{1}\) and \(c_{2}\) of \(K(s)\). Assume a first-order plant and a second-order controller
and assigned poles of the control system
Then the system is underdetermined referring to the parameters \(c_{1}\) and \(c_{2}\). They can be arbitrarily chosen without affecting the eigenvalue assignment. The task is
The optimum controller results as
4 Valuating constraints
An important issue is valuating constraints, e.g. the maximum actuation signal \(u_{m}\) arising from a reference step function. What is the design process? In Fig. 7 lines of given constant values \(b_{GK}\) are repeated in the basis plane \(c_{1}\), \({c}_{2}\). A plane of resulting \(u_{m}\) is depicted. Above one of these lines of constant \(b_{GK}\), a single line of values of \(u_{m}\) is plotted as a thick line. This line demonstrates \(u_{m}(c_{1},c_{2})\) with no influence on the implemented eigenvalue assignment and chosen chordal stability margin \(b_{GK}\).
5 Structured rectangular uncertainty versus chordal uncertainty
The question arises what is the chordal uncertainty measure \(\delta[G_{p}, G]\) for a given structured rectangular uncertainty in \(G_{p} = \frac{2s+b_{z}}{s^{2}+3s+a_{N}}\); \(a_{N}=5\), \(b_{z} = 4\), \(K(s) = \frac{s +2}{20s+1}\), with parameters \(a_{N}\), \(b_{Z}\) uncertain less than a certain amount, see Fig. 8.
Another question is: Which rectangular uncertainty is equivalent to a chordal one. Which part of the rectangular parameters is included if e.g. \(\delta [G, G_{p}]< 0.03\) is given. The result is depicted in Fig. 9. For illustration purpose the set of eigenvalues of the closed-loop system are included in this figure when a so-called uncertainty walk is executed along the border of the chordal uncertainty.
6 Graphically supported design
Figure 10 is derived from Eq. (58). The function \(N\) is depicted versus frequency before having selected the minimum. The areas between the function and the abscissa are \(\int_{0}^{\infty}(\cdot)^{*}(\cdot)d\omega\). In the figure, 6.3 stands for \(|GS|^{2}\) and 2.5 for \(| GKS|^{2}\), for the outputs caused by disturbance and reference in the steady state, respectively. These values represent the root mean square energy content of the signals. According to Parseval Theorem they correspond to the integral of squared impulse response. Usual requirements are increasing \(b_{GK}\), observing all the actuations and selected limits. Besides, also the squared initial values can be observed if the frequency axis is scaled up to higher values; the initial actuation signal, 4 and 1, from KS and S, caused by reference and disturbance, respectively. When changing parameters during the design process, special limits, e.g. \(\|wKS\|_{\infty}\), can be taken into consideration.
The influence of a varying design parameter \(v\) is shown in Fig. 11. A subplot demonstrates the frequency-dependent \(\rho(G, K)(\omega; v) \) and the local minimum in \(b_{GK}\).
An alternative could be: Varying \(k_{i}\), elements of K, from \(1/b^{2}_{GK} \to\min_{k_{i}}\) from Eq. (48) one finds
Then, \(b_{GK}\) and \(\|KS\|_{\infty}\) are displayed having an eye on the limits.
7 Chordal stability margin for holistic eigenvalue assignment
The trace of the negative inverse closed-loop coefficient matrix in state-space system design is the sum of inverse poles [4]. Shifting the poles to the left as much as possible
is required, conditionally augmented by the norm of \(\mathbf{K}\). Now, when \(b_{GK}\) is included
has to be executed. Figure 12 depicts the result with an additional condition \(c_{2} > 2\).
8 Conclusion
The presented methods are intended to demonstrate high-order control systems and graphically supported design tools, adequate for applications. The methods are selected of simple structure in order to simplify the access.
Change history
24 August 2017
An erratum to this article has been published.
References
El-Sakkary, A. K. (1989): Estimating robustness on the Riemann sphere. Int. J. Control, 49, 561–567.
Qiu, L., Zhou, K. (2013): Preclassical tools in postmodern control. IEEE Control Syst. Mag., 33(4), 26–38.
Vinnicombe, G. (2001): Uncertainty and feedback: \(H_{\infty}\) loop shaping and the \(\nu\)-gap metric. London: Imperial College Press.
Weinmann, A. (2001): Control system design based on holistic eigenvalue allocation. E&I, Elektrotech. Inf.tech., 118, 167–173.
Acknowledgements
Open access funding provided by TU Wien (TUW).
Author information
Authors and Affiliations
Corresponding author
Additional information
An erratum to this article is available at https://doi.org/10.1007/s00502-017-0511-5.
Appendix A
Appendix A
This appendix is considered as a clearly arranged presentation of interrelations which might be confound with each other, and to get familiar with the intermediate derivations, even though they are known in the literature [1–3],
1.1 A.1 Singular values
The singular value obeys \(\sigma [\mathbf{G}] \triangleq+\sqrt{\lambda [\mathbf{G}^{H}\mathbf{G}]}\) and its maximum \(\sigma_{\max}[\mathbf{G}] = +\sqrt{\lambda _{\max}[ \mathbf{G}^{H}\mathbf{G}]}\equiv \|\mathbf{G}\|_{s}\). The \(H_{\infty}\) norm is \(\|\mathbf{G}(s)\|_{\infty}\triangleq \sup_{\omega}\sigma_{\max}[\mathbf{G}(j\omega)]\). The maximum singular value is a generalization of the absolute value of scalars for matrices.
For any matrix \(\mathbf{X}\) or a special vector \(\mathbf{X} = (a, b)^{T}\) one has \(\mathbf{X}^{H}\mathbf{X}= (a, b)^{*} (a, b)^{T} = (a^{*} a + b^{*}b)\) and
1.2 A.2 Chordal uncertainty measurement \(\delta[G_{p}, G]\)
Note that \(p^{*}(j\omega) = p(-s)|_{s= j\omega}\) and \(a(s) \cdot a(-s)|_{s= j_{\omega}} = \mathbf{a}^{T}\boldsymbol{\omega}\cdot \mathbf{a}^{T}\boldsymbol{\omega}^{*} \stackrel{\text{(13)}}{=} \mathbf{a}^{T}\boldsymbol{\Omega} \mathbf{a}\). Note also \(G^{*}\triangleq G(-s)|_{s= j\omega}\) (for functions rational real \(P^{*}(j\omega ) = \text{conj} (P(j\omega) = P(-s)|_{s= j\omega})\).
There is the rather strange definition of the distance of two transfer functions where each one is a function of \(s = j\omega\), but it turns out similar Eq. (7)
For the case of uncertain plant \(G_{1}= G_{p}\) and \(G_{2} = G\).
For \(G_{1}\) close to \(G_{2}\) and for small uncertainty, the 1 in Eq. (33) can be neglected
For the nominal \(G\), the family of perturbed \(G_{p}\) and \(\delta [G(s),G_{p}(s)]\) one has
1.3 A.3 Chordal distance \(\rho[G,K]\) and stability margin
Further correspondencies referring to \(b_{GK}\) are given as follows. For \(G_{1} = -K^{-1}\) and \(G_{2} = G\) in Eq. (36), \(b_{GK}\) is defined as the smallest chordal distance between \(G\) und \(-K^{-1}\)
where \(\min_{\omega} |G^{-1}| = \max_{\omega}|G| = \|G\|_{\infty}\) was used.
We start a new derivation with Eq. (45) and we multiply numerator and denominator with \(|K|\)
Alternatively, from Eqs. (12), (46) and (47) it follows \(b_{GK} \triangleq \min_{\omega}\frac{1}{\sqrt{N}}\) and
The chordal distance \(p[G,K]\) of Eq. (52) equals \(1/\sqrt{N}\) of Eq. (56).
Take notice of the different interpretations of \(b_{GK}\). The origin resulting of automatic control Eq. (39); the chordal distance in Eq. (45); the minimum value in frequency domain in Eq. (12); an expression using \(H_{\infty}\)-norm Eq. (48) and the graphical correspondence using all the components of \(\mathbf{M}\) in Eq. (58).
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Weinmann, A. Chordal automation towards high-order systems. Elektrotech. Inftech. 134, 276–281 (2017). https://doi.org/10.1007/s00502-017-0507-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00502-017-0507-1