Incremental holistic controller

The trace of inverse coefficient matrix method in state space is used to implement and simplify an iterative design of a control system with main interest in increasing dynamical performance quality or to check if stability danger exists. The τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\tau $\end{document}-locus is a helpful design element. Various applications and examples are presented for illustration.

The term of holistic controller can be renamed in more detail by Trace of inverse matrix method, i.e., as specified in [3,5], For the closed-loop control system with the output controller K y one has tr[(A + BK y C) − where the controller is utilized to approach a change of the control system dynamics. Executing the first derivative with respect to K and K y , for the state controller and for the output controller, one has, Eq. (17) and Eq. (18), respectively, and Holistic controller yields a proposal for the controller design improving the dynamic performance of the control system as far as stability and oscillation properties are concerned. An increment change in the state controller K = −αB T A −2T can be performed.
A norm-oriented τ -(root)-locus is developed in what follows where the sum of eigenvalues of the A −1 matrix (not λ i [A]) are balanced versus the controller norm. A norm-limited output controller and the τ -assignment is described as well.
The main idea is to use tr Eq. (6) which is equated with zero to obtain the extremum, usually the minimum.
Consider two poles at −8 and +0.1, one of them unstable. The resulting τ yields a command to shift to the left in the same way as if the pole were at −0.1. This corresponds to τ → min. Consider the shift result to −0.08, τ goes from negative high to a negative bigger value but τ continues reducing its size. Crossing the origin then τ would jump to positive high, the direction would keep all right.
The question is of interest when τ is a correct indicator. In the range of correct indicator, it is a simple scalar number, easily to be expressed with syms parameter for the state-space data. Eigenvalues, on the contrary, are always difficult to derive and to manipulate. Trend towards stability or danger is useful and possible even if it is not precise. Anyhow, in engineering systems in many cases one has to investigate additionally the numerical effort or computer-assisted simulation [2]. Well-established stability criteria are useful for checking stability of a given system. However, modifying to use for design facilities, the task of system stability often requires some additional efforts. A similar situation results in this connection.
Consider, e.g., a three-dimensional system 1/[(a + jb)(a − jb)c]. Then, τ is 2a/(a 2 + b 2 ) + 1/c and 2a/b 2 + 1/c for a small. The stability border is a = 0, of course. If parameters of the system are changed in such a way that all the poles move towards the imaginary axis, the stability border is yet reached at a = 0 and one has τ = 2a/b 2 + 1/c = 1/c. Smaller values of a show instability yet. The τ -stability border is yet crossed at τ > 1, τ = 0 is not the stability border. For τ < 1/c the system is unstable all along, one eigenvalue is located in the right half-plane.
Trace of inverse matrix method is a good help to move the design activity towards stability increase or decrease. The correct ap-

Fig. 1. Example 1 numerically solved. Diamonds indicate the eigenvalues of A for the initial setup when starting the ten steps
proach has to be decided with additional considerations of physical or engineering domain. Hence, the result can be outlined as "extremum"{tr[A −1 ]} either max or min in a formula.

Norm-oriented and stability-oriented τ -locus
} as an equality condition. Selection of a weighting parameter μ and the τ -locus of corresponding eigenvalues yields an opportunity for valuating stability and control effort. Since τ is a replacement only, the advantages and disadvantages have to be rigorously investigated. Intending for − i 1/λ i [A + BK] → min K including the norm as a condition, one has Differentiation ∂ ∂K providing the minimum, using Eq. (13) and (17), yields Analytical solution is only possible for very low order, with respect to commonly having A and A T in the same relation. With respect to the power −2, solution is only possible in a numerical way.
The parameter μ is a factor in the τ -locus, balancing the norm of K versus τ using their gradients. No single eigenvalues are considered.
The τ -locus can also be used to represent the closed-loop eigen- The τ -characteristics as displayed in Fig. 1 are very different from the classical root loci. The characteristic starts at small μ.
Example 1 Eq. (6) is solved numerically for an unstable system of order n = 6, without executing a neglection and with ten μ steps λ[A(K, μ)]. Figure 1 shows the result while increasing μ following the arrows, including a subplot norm of K versus μ; or some alternative measure valuating the control effort for the current application.
To develop the example, having selected the initials for μ, K, and A, find K ∝ −B T A −2T /(2μ), augment K with K, get A := A + BK and plot λ[A] for illustration, yielding Fig. 1.
Repeating the derivation leading to the previous section, now for the special case of output state controller, one finds K y = −B T (A + BK y C) −2T C T /(2μ).  Fig. 2). Their minimum is b GK and provides the optimum μ.

Chordal projection
Example 2 In Fig. 2 the straight line distances between G and K are the basis for yielding b GK , where min ω |G −1 | = max ω |G| = G ∞ was used, One tries to guarantee b GK > r G , the chordal uncertainty radius [1,4]. The chordal projection does not provide a design proposal for the controller. But chordal projection is used to improve the design approximation by τ in such a way that μ is selected observing b GK and/or satisfying the chordal uncertainty radius.

Advantageous diagnosis with uncertainty walk
Useful diagnosis can be obtained by using the τ -method. Consider, e.g., a continuous-time system of n-th order with n parameters and distinct eigenvalues in the nominal case. Including uncertainties of the parameters at a specific closed cyclic (spherical) boundary, each eigenvalue moves along a closed route when displayed in the complex plane. This route is termed uncertainty walk. The uncertainty walks are different as far as the approximation to the instability border is concerned, or where the best opportunity exists for improving the performance quality.
Example 4 Using the τ -method, the n eigenvalues can be replaced roughly by τ -values. The question arises and can be affirmed in many cases, if the τ -value can be used to roughly replace the n uncertainty routes by a single route of the τ -value, or if at least a first diagnosis can be stated. In this example, these properties are displayed and illustrated in a 3D-plot (Fig. 3) or in an unwrapped display (Fig. 4).

Numerical coefficients and pseudo inverse
As a final application, consider a desired change in some dominant eigenvalues in magnitude d, which has to be performed by changes K of the state controller K, see Eq. (9) and Eq. (11). Using multiple differentiation with respect to s yields the coefficients of powers of s to compare the given and desired increments of some dominant eigenvalues with the unknown changes in the controller parameters.
If you were handicapped by numbers of 20 and more, you can start with comparison of replacing powers of s by any numerical numbers of adequate height. Selecting small changes of eigenvalue increments and controller parameter increments you can use the Taylor expansion of the determinant. In Eq. (16) you simply use more comparisons of coefficients than necessary and you directly obtain the coefficients for the left pseudo-inverse matrix M L .
The incremental eigenvalue assignment adaptation is explained with a simple example. This, of course, is classical eigenvalue assignment, only Taylor expansion is unusual. Consider the system, e.g., with closed-loop equation You want to change the conjugate complex pole pair at e.g. −1 ± j3, such that it is moved to the left and towards the negative real axis with a quantity d = −0.05 ± j 0.1. This change increment has to be provided by Eq. (16)
Use is made first by deriving the coefficients of the powers of s by successive differentiation with respect to s, and second by stating the equations resulting from Eq. If the comparison of the factors is performed by replacing s by numerical values, there should be used more than three numerical values inserted. This is useful for applying pseudo inverse matrix.
In Line 10 substituting s by a numerical value produces a vector r ∈ R q×1 . In Line 12 the underbraced expression [aa] delivers a row m ∈ R 1×3 and q times inserted a matrix M ∈ R q×3 . [bb] delivers a vector v ∈ R 3×1 . Thus we have r = Mv which is usually overdetermined with respect to numerical uncertainties. Solving for the given data r and M, one finds the unknown v = (M T M) −1 M T r = M L r.

Conclusion
Changing the set of inverse eigenvalues λ i [A −1 ] to τ , the trace of the inverse, is a strong simplification with disadvantages and advantages. The latter is the opportunity of algebraic representation with syms variables, with an indicator for a simple setup. It has been discussed if this change is worth while. The item is how τ can substitute the eigenvalues and which advantages and problems arise.
The incremental small step eigenvalue assignment and its determinant Taylor expansion, not only for powers of s but for numerical values, directly provides coefficients for the pseudo inverse matrix. The τ -locus with the parameter μ valuates the magnitude of the controller norm for the purpose of direct control system design. The design result is also valuated by the stability margin of the chordal projection method.
In the case of uncertainty walk, the τ -method provides a replacement of λ i by τ .

Funding Note Open access funding provided by TU Wien (TUW).
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.