Multiobject Holographic Feedback Control of Differential Algebraic System with Application to Power System

A multiobject holographic feedback (MOHF) control method for studying the nonlinear differential algebraic (NDA) system is proposed. In this method, the nonlinear control law is designed in a homeomorphous linear space by means of constructing the multiobject equations (MOEq) which is in accord with Brunovsky normal form. The objective functions of MOEq are considered to be the errors between the output functions and their references. The relative degree for algebraic system is defined that is key to connecting the nonlinear and the linear control laws. Pole assignment method is addressed for the stability domain of this MOHF control. Since there is no any approximation, the MOHF control is effective in governing the dynamic performance stably both to the small andmajor disturbance.The application in singlemachine infinite system (SMIS) shows that this approach is effective in the improvement of stable and transient stability for power system on the disturbance of active power or three-phase short circuit fault.


Introduction
Power system shows strong nonlinear properties, especially in the action of large disturbance such as three-phase short circuit fault.Many effective control theories and methods have been proposed for this nonlinear system in literatures, such as differential geometry theory [1,2], Lyapunov energy function [3,4] and Hamiltonian function [5,6], machine learning, and neural network [7][8][9][10].Feedback linearization based on differential geometry theory is effective in the stability control of the nonlinear power system [2].Different from the approximate linearization, Differential geometry method designs the nonlinear control law by constructing the mapping from the nonlinear space to the homeomorphous linear without approximation.Tan and Wang in [11] designed the adaptive excitation and phase shifter controller based on third-order generator model.Kennedy et al. in [12] designed a nonlinear excitation controller based on a form of state feedback linearization using the geometric approach.With the work of Sastry and Isidori [13], the adaptive control of "minimum-phase" nonlinear system was studied by state feedback exactly linearization and several initial results were derived.
Developing the differential geometry theory, Li et al. in [14] proposed the Multi-index Nonlinear Control (MINC) method.Their research demonstrated that it was not necessary for the system to satisfy the condition of completely exact linearization, which is usually too strict condition to satisfy for most of systems.Li et al. also addressed that the part exact linearization made the selection of output function more flexible to adapt the control requirement of the controlled system, which made MINC method get more satisfied control performance.Li et al. [14] also proposed an effective and practical method of determining output function that was linear combination of state vectors.The precondition of MINC method adopting part exact linearization is that the zero-dynamic must be stable, which is possibly difficult to prove for general system.For solving this problem, Liu et al. in [15] proposed the multiobject holographic feedback (MOHF) control method.This method designed the nonlinear feedback law in linear space by constructing multiobjective equations that satisfy Brunovsky normal form 2 Mathematical Problems in Engineering [16].Furthermore, it is unnecessary to prove the stability of zero-dynamic.These research works are based on the ordinary differential equation (ODE) model.
For most of the complicated power system, it is insufficient that ODE describes wholly the dynamic activity, since the connecting nodes of power system have to be formulated in algebraic equations.The general power system, hence, is usually modeled as differential algebraic equation (DAE).
Several typical methods are described in literatures for studying the differential algebraic system and its application in power systems.Kurina and März in [17] addressed the linear quadratic optimal control problems with constraints described by general linear DAEs with variable coefficients.The optimal feedback control matrix can be obtained by suitably formulating a Riccati DAE system, which is similar to the classical example in which the constraints are described by explicit ODE system.Tsolas et al. in [18] presented the preserving structure model of general power system and suggested that the DA model could be replaced locally by an ODE description.They developed further some geometric characterizations of the region of attraction for a stable equilibrium point.The transient stability was studied in [19,20] based on preserving structure model and the total energy function were shaped.Wang and Chen [21] presented the differential geometric method that used the preserved structure model for DA systems.The Lyapunov function was proposed in [22] for the stability research of multimachine power systems.The energy function presented in [23] as well as [18] showed some properties of Lyapunov function.Kautsky et al. [24] described the numerical method for determining robust solutions of pole assignment by state feedback.The robust pole assignment problem is also discussed in [24].One may see the eigenstructure assignment approach of the state output feedback pole placement for  input  output in [25].The pole assignment problems were proposed by [26][27][28] with projection and deflation methods together with optimization approach, respectively.This paper explores the MOHF control described by the differential algebraic model.As an example of application in power system, the feedback law is designed for single machine infinite system.The structure is organized as follows.Section 2 presents a model of general SIMO differential algebraic system.Section 3 defines the Γ derivative and relative degree for differential algebraic system.Section 4 proposes the MOHF design method and pole assignment approach for this design.Section 5 designs a practical power system based on the preserving structure as an application of the MOHF theory.Section 6 gives the conclusions.

Relative Degree of the Differential Algebraic System
The definition of the relative degree for ordinary differential equation system was given by Sastry and Isidori in [13].
We redefine it here in order to design differential algebraic system.Firstly, a special matrix is given as where   is an unit matrix of  × .This matrix Λ(, ) exists in the condition of Assumption 3 holding.The definition of Γ derivative is formulated as follows.
Definition 4. The product of the Jacobi matrix of output function ℎ(, ) with respect to  and , the matrix Λ(, ), and the vector field (, ) are defined as Γ derivative of ℎ(, ) with respect to (, ) in differential algebraic system; that is The higher Γ derivatives are then defined as Correspondingly, the Γ derivative of the Γ   ℎ(, ) ( = 1, 2, . . ., ) to (, ) is defined as An immediate definition of the relative degree for the differential algebraic systems is yielded as the following.
Definition 5.The DA system of ( 1) is said to have relative degree  to the output function  = ℎ(, ), if and only if the following conditions holds: With these definitions of Γ derivative and relative degree, we can establish the homeomorphous mapping between nonlinear space and linear in the MOHF control design of next section.

MOHF Control Design
In this section, we describe the MOHF design method in detail.We establish, firstly, the adaptive model and then discuss the rules of choosing objective function.In the following section we propose the structure of multiobjective equation and design feedback control law.The problem of pole assignment is solved in Section 4.4.

Adaptive Model. Consider the reference model:
The adaptive model is derived from ( 9) and ( 1): where The reference adaptive output function of plant ( 1) is followed as where   = ℎ(  ,   ) is the output of ( 9) and also the reference output of plant (1).

Objective Functions.
Objective functions are the variables that need to be constrained in control.These functions are considered as the form where and   and   are the th component of output function in (12) and its reference trajectory.In this MOHF control, the output functions ( 12) are determined for the purpose of getting satisfied control object.They include not only the practical outputs of the system but also other variables related to control properties.We consider them as the linear combination of state vectors and algebraic variables, which means where  ∈   , and are the constant matrices decided by the dynamic and stability properties of plant.
Remark 6.The choice of objective function acts importantly in MOHF control.It decides the effect of control and the complexity of feedback law.Two factors fall into main considerations.The first, also the most important, is the relations between objective functions and control performance.The second is the practicability of the designed control law, since the excessively complicated law is probably difficult to carry out and maybe declines the robustness.Based on both of the considerations, it is of high priority to let the relative degree hold in determining the objective functions.Compared with higher , this can make feedback law more simple and hence easy to be used in practical system.For most of system, ( 17) is not a hard condition if the objective functions are arranged appropriately.Our following design is based on this condition satisfied.

Design of MOHF Control Law.
In this section, we use the objective functions in (13) to construct the multiobjective equation that is in agreement with Brunovsky normal form [16] and design the nonlinear feedback law in linear space.The multiobjective equation is formed as where are Brunovsky constant matrices.This system ( 18) is called the extra system.The feedback law V is called extra system control law.It is designed as the following form according to linear optimal theory The feedback matrices K = ⌊− 1 − 2 ⋅ ⋅ ⋅ −  ⌋ can be yielded by the next pole assignment.
The control law  of system (1) is accordingly called as the inner system control law, which is bridged with V by ( 5)- (7).From (18) we have and from (20), the inner feedback control law Δ is obtained: Remark 7.This MOHF method uses the extra control law V to govern the objective functions to their equilibriums.The homeomorphous transformations of ( 21)- (22) transfer this law to equivalent input Δ of nonlinear system to constrain the nonlinear dynamic oscillation.Since this transformation from nonlinear linear space to linear holds all the state information without approximation, this approach may enhance the dynamic stability when V is designed appropriately.The feedback and output constant matrix are designed by the following pole assignment.
For the stability control of power system, the aim of control is to make the system reinstate to its original equilibrium or achieve a new equilibrium after the disturbance deleted.So we consider that the reference trajectories of ( 23)-( 25) are equilibrium; that is Substituting ( 24) into ( 23) and ( 25), the linear adaptive model is as follows: which is expressed concisely as Selecting the nonsingular linear transformation matrix   , let Δ x = x − x =   Δ x; then (29) are transformed to the following controllable normal form where From ( 29), (30), and (31) we have that Assume that ( 30)-(31) hold condition (17), we get the feedback law (34) Substituting (34) into (30) we get where The characteristic polynomial of (31) is Remark 8. Equation ( 27) describes the linear asymptotical system of plant (1) in the neighborhood of hyperbolic singular point.According to Hartman-Grobman theory, the above pole assignment to the linear system can achieve locally stability control of (1).
Remark 9.This MOHF design method supplies an inner adaptive reference track that is the equilibrium state.In case that the feedback law constrains the objective functions to its reference tracks, the stability of the system is held.

Application
As an example of application, we designed MOHF law for single machine infinite system (SMIS) based on a DA model.
In Figure 1, the solid lines show the simulink results of MOHF control; the dotted lines are the results of linear optimal control (LOC), and the dash-dotted lines depict the dynamics of adjustment when the power system stabilizer (PSS) is implemented in the power system.As the overshoot of rotor speed Δ of PSS is much larger than both of the MOHF control and LOC, additional plot (Figure 1(c)) is used to illustrate the dynamics precisely.
Figure 1(a) shows the response of the electromagnetic power   of the generator when the disturbance of mechanical power imposed on the SMIS.From these simulink curves one can see that when the system is exerted by a constant disturbance from the active power, the MOHF method can adjust the rotor angle by supporting an equivalent input Δ to improve the electromagnetic power and attain rapidly new power balance with good stable and dynamic properties.This dynamic process is accompanied with smaller oscillation and overshoot than LOC and PSS method, which is beneficial for protecting the generator against the mechanical shock.Figure 1(b) is the response of rotor speed whose precise comparison of MOHF and LOC is shown in Figure 1(c).Oscillations of small magnitude and low frequency may present limitations on power transfer capability in some cases and affect the stability [29,30].Figures 1(b) and 1(c) demonstrate that the MOHF control is more effective in damping frequency oscillations of system.Figure 1(d) is the plot of the generator terminal voltage deviation.It gives the result that the MOHF control method adjusts the output voltage to stable state in very short time with zero error.Although the overshoot is bigger than LOC, it contents the power system criterion.Figures 1(b), 1(c), and 1(d) show that with the control of MOHF, the system attains better electric qualities in regards to voltage fluctuation and low frequency oscillation than the LOC and PSS when this system is subjected to an active power disturbance.
Figure 2 shows the responses when the three-phase short circuit fault occurs at one second in the infinite bus of SIMS and is solved after 0.1 second.
The dynamics of   , Δ, and Δ  are showed in Figures 2(a), 2(b), 2(c), and 2(d), respectively, with the solid blue lines for MOHF, red dotted lines for LOC, and orange dash-dotted lines for PSS.From Figure 2(a) one can see that, when the serious fault acts on the system, the dynamic feedback law Δ of MOHF can effectively restrain the oscillation of power and make the active power   rapidly reinstate its balance point.From these figures one may find that when the control constant matrices are adjusted, the MOHF still contains good transient response with zero error while the LOC with error.This result indicates that the MOHF method holds better robustness than LOC.

Conclusion
This paper has explored the principles and design method of multiobject holographic feedback control for nonlinear differential algebraic system and applied them to design the feedback law of single machine infinite bus power system.The concept of relative degree and Γ derivative were defined which are used to connect the output function and equivalent input.Based on these works, the homeomorphous mapping from nonlinear space to linear was established.By constructing the multiobject equation which is in accordance with Brunovsky normal form, the control objects of nonlinear system were transformed to linear space to implement control.Pole assignment is used to design the output constant matrix to attain the stability.
Since this transformation from nonlinear space to linear employed no approximation and held all of the nonlinear information, this MOHF control may inhibit the oscillation more effectively than LOC and PSS approaches while the small or large disturbance imposes on the system.The practical application in the SMIS demonstrates the good Mathematical Problems in Engineering effect of this inhibition.Unnecessary proving the stability of zero-dynamics, the MOHF control based on differential algebraic model may be helpful in the research of more general electrical power system and other similar nonlinear differential algebraic systems.
of output function (2) is The curve of Δ for MOHF control compared to LOCResponse of voltage error to power disturbance The curve of Δ

Figure 1 :
Figure 1: System responses to the disturbance of the mechanical power stepping up 20% regulation.
Rotor speedResponse of rotor speed to three-phase short circuit fault ×10 The curve of Δ for MOHF control compared to LOC of voltage error to three-phase short circuit fault 0The curve of Δ

Figure 2 :
Figure 2: System responses to the disturbance of a three-phase short circuit fault.

Figure 3 :
Figure 3: System responses to adjustment of the feedback matrix.

Figure 2 (
Figure2shows the responses when the three-phase short circuit fault occurs at one second in the infinite bus of SIMS and is solved after 0.1 second.The dynamics of   , Δ, and Δ  are showed in Figures 2(a), 2(b), 2(c), and 2(d), respectively, with the solid blue lines for MOHF, red dotted lines for LOC, and orange dash-dotted lines for PSS.From Figure2(a) one can see that, when the serious fault acts on the system, the dynamic feedback law Δ of MOHF can effectively restrain the oscillation of power and make the active power   rapidly reinstate its balance point.Figure 2(b) and Figure 2(c) show that the dynamic response of rotor speed for MOHF holds smaller maximum overshoot, shorter settling time, and fewer frequency oscillation than LOC and PSS.From Figure 2(d) one can see that the dynamic reactive power compensation of the MOHF ensures the output voltage of the generator with good dynamic and stable performance.It is clear from Figure 2 that the MOHF controls the dynamic behavior of system in serious disturbance more effectively than LOC and PSS.Robust Analysis.Figures 3(a) and 3(b) are the responses while the control constants are adjusted.The solid lines show the results of the MOHF control and the dotted lines are results of the LOC.This figure (Figure 3(a)) shows the simulation result with the output constant  3 of MOHF changed from  3 = −15 of Figure 2 to this figure's  3 = −10, other constants keep no change, and the feedback matrix K of LOC is adjusted from  = [75 10 −2000] to  = [70 10 −2000].Figure 3(b) shows the result of constants  1 and  from  1 = 500 to  1 = 505 and  from  = [70 10 −2000] to  = [70 12 −2000].From these figures one may find that when the control constant matrices are adjusted, the MOHF still contains good transient response with zero error while the LOC with error.This result indicates that the MOHF method holds better robustness than LOC.
The second algebraic equation shows the relations of the voltage, active power, and (or) reactive power flows among the nodes. ∈   is state vector. ∈   is algebraic variable. ∈  is the input of the plant. ∈   is the output function which is decided ]  is the state vector of the generator and  = [    ]  the algebraic constraint variable. is the rotor angle,  is the rotor speed, and    is the transient internal voltage of the generator.  and   are direct-axis and quadrature-axis synchronous reactance of the generator, respectively.   is direct-axis transient reactance.  ∑ =   +   +   ,    ∑ =    +   +   , and   ∑ =   +   +   , where,   is the transformer reactance and   the transmission line reactance.  is output terminal voltage of the generator and  is the voltage at infinite bus.   is the time constant of the field-winding with open-circuit armature winding. is the constant of damper winding and   the inertia time constant of the rotor.