Nonlinear rational systems/models, also known as total nonlinear dynamic systems/models, in an expression of a ratio of two polynomials, have roots in describing general engineering plants and chemical reaction processes. The major challenge issue in the control of such a system is the control input embedded in its denominator polynomials. With extensive searching, it could not find any systematic approach in designing this class of control systems directly from its model structure. This study expands the U-model-based approach to establish a platform for the first layer of feedback control and the second layer of adaptive control of the nonlinear rational systems, which, in principle, separates control system design (without involving a plant model) and controller output determination (with solving inversion of the plant U-model). This procedure makes it possible to achieve closed-loop control of nonlinear systems with linear performance (transient response and steady-state accuracy). For the conditions using the approach, this study presents the associated stability and convergence analyses. Simulation studies are performed to show off the characteristics of the developed procedure in numerical tests and to give the general guidelines for applications.
1. Introduction
This section justifies the reasons for designing controllers for rational models by introducing model expression and representations, achieved results in model identification, and a critical review of controller-designing approaches.
Assign a triplet X,f,h, where X is an irreducible real affine variety and f,h are mapping functions. A system Σ, with input U∈ℝm and output Y∈ℝr, is defined as polynomial/rational, while the functions f=fα∣α∈U and h:X→ℝr both on X are mappings from input space to state space and from state space to output space polynomial/rational, respectively. That is, for polynomial systems, hi∈A for all i=1,…,r where A is the algebra of all polynomials on the variety X, and for rational systems, hi∈Q for all i=1,…,r where Q is the algebra of all rational functions on the variety X.
For a single-input and single-output (SISO) nonlinear dynamic rational system, it can be generally modelled with a ratio of two polynomials [1, 2].
(1)yk=NpkDpk+ek=NpYk−1,Uk−1,Ek−1DpYk−1,Uk−1,Ek−1+ek=∑j=1numpnjkθnj∑j=1denpdjkθdj+ek,where yk∈ℝ, uk∈ℝ, and ek∈ℝ denote measured output, input, and model error/noise/uncertainties, respectively, at time instant k=1,2,…. Npk and Dpk are real valued and smooth numerator and denominator polynomials, respectively. Yk−1∈ℝn⊃yk−1,…,yk−n, Uk−1∈ℝn⊃uk−1,…,uk−n, and Ek−1∈ℝn⊃ek−1,…,ek−n denote the delayed outputs, inputs, and model noises, respectively. pnjk∈ℝ and pdjk∈ℝ for regression terms θnj∈ℝ and θdj∈ℝ, respectively, are the coefficients and num and den for numbers of total regression terms of the polynomials. The major properties of the rational model (1) are summarised below:
It is also defined as a total nonlinear model [2] as it covers many different linear and nonlinear models as its subsets (such as NARMAX (nonlinear autoregressive moving average with exogenous input) models [3] and intelligent models for neurofuzzy systems [4]). Rational systems have been observed in general engineering, chemical processes, physics, biological reactions, and econometrics; for example, rational models are a class of mechanistic models in describing catalytic reactions in chemical kinetics [5, 6]; metabolic, signal, and genetic networks in systems biology [1]; and movement of satellites in Earth orbit [1]. There have also been reports of rational modelling applications [7–9].
This is more concise in structure than a polynomial; the example below uses a Taylor series expansion to approximate a simple rational model below.
(2)yk=sinuk−11+y2k−3=sinuk−11−y2k−3+y4k−3.
The other characteristic of the rational models is the power to quickly change the model output while input has small variations. Consider a simple system output below
(3)yk=11+uk−2.
Clearly the model output will be dramatically increased, as the input uk−2 approaches −1. This comes from the function of the denominator.
Introducing a denominator polynomial makes the model concise in describing complexity and adds more functions in describing nonlinearities. On the other side, in contrast to polynomial systems, this makes identification and control system design noticeably different and more difficult with the inherent nonlinear parameters and control inputs [2]. Therefore, comprehensive studies of this class of systems in theoretical and application aspects are required. This study takes the pioneer step towards the control of rational systems.
1.2. Model Identification
Model identification has been relatively mutual to some extent. So far, the identification aspect has gone through data-driven model structure detection, parameter estimation, and model validation from noise-contaminated input and output data. The major work on rational model identification is summarised in the following categories: linear least squares (LLS) algorithms for parameter estimation—extended LLS estimator [10], recursive LLS estimator [11], orthogonal LLS structure detector and estimator [12], fast orthogonal algorithm [13], and implicit least squares algorithm [14], and nonlinear least squares algorithms—prediction error estimator [15] and globally consistent nonlinear least squares estimator [16]. Other algorithms include the following categories: back propagation (BP) algorithm [17] and enhanced linear Kalman filter (EnLKF) [18].
There are two model validation methods: higher order correlation tests [19] and omnidirectional cross-correlation tests [20].
A summary of the representative publications till 2015 can be found in a survey of rational model identification [2].
1.3. Controller Design
As surveyed above, rational models have been increasingly used to represent nonlinear dynamic plants. Consequently, the control system design should have been considered on the agenda in the follow-up studies. However, up to now, there is no reference found for designing such controllers directly referred to the model analytical knowledge. The paramount difficulty is that part of the controller output is embedded in the denominator polynomial Dpk. For example, yk=0.5yk−1−yk−2uk−1+0.1u5k−1/1+y2k−2+0.2u2k−1. With extensive investigations through major academic publication searching engines, it can be concluded that this study is the first trial with analytical approaches to design a controller for rational systems.
Regarding controller design approaches possibly referred to the rational systems, these could be the reduction of rational model structure complexity, which are neural network models, linear approximation models, linearization, and iterative learning control and U-model enhanced control. A brief critical review of the approaches is presented.
Reference [21] on neural controllers is probably the first publication relating to control of rational models. However, the design approach has merely used rational models as extreme nonlinear examples; it has not designed controllers by taking the model structure into consideration (even if known in advance), except for taking the models as the representatives of complex nonlinear dynamic systems.
Piecewise linearization [22, 23] around operating points has been widely studied to simplify controller-designed procedures when plants are subject to mild nonlinear dynamics. It should be mentioned that a group of piecewise linear models can be admitted as a linear model, with varying order and parameters in different operating intervals. The promising property is using linear control design strategies directly. However, it could induce inaccuracy and dynamic uncertainty because of ignoring some inherent nonlinearities from their original nonlinear representations. Further, this method may also increase computational burden/complexity while overborrowing piecewise linear intervals to match severe nonlinearities.
Pointwise linearization has been claimed by neural network-based control and/or adaptive control, which uses linear models to approximate predominant dynamics around an operating point or every input-output dynamic gain at each time instance and then employs a neural network to determine the error induced by the linearization [24, 25]. Once again, it uses linear control system design to construct nonlinear control systems. However, this involves online neural network learning or online model iden parameter estimation, and therefore, the constructed nonlinear control system is operated under adaptive principles (the controller parameters are updated with the neural network output), even for deterministic nonlinear plants. The other related issue is the selection of neural network topology, which has no systematic procedure available to find the best fitted neural network representative.
Feedback linearization is a well-developed subject [26]. A general SISO nonlinear system is described as
(4)ẋ=fx=gxuy=hxwhere x is the state vector and u and y are the input and output, respectively. f⋅, g⋅, and h⋅ are real valued and smooth mapping functions. With this model structure, a series of analogies with some fundamental features of linear control systems have been established, which provides a very useful concept in the design of nonlinear control systems using linear design methodologies. Obviously, the model has u in an explicit position. The studied nonlinear rational model has no such explicit expression for input u to be designed, and this immediately reveals that the methodologies rooted in the approach, although useful references, are not directly applicable in designing control of nonlinear rational systems. The other input-output linearization techniques [27] have had similar requirements for an explicit u expression and special skills for state variable transformation.
Iterative learning/data-driven control/model-free control is another possible control system design methodology in avoiding model structure complexity. The approaches do not require a clear plant model structure but still need plants with some mild conditions in control [28, 29]. Again, if a rational model is available, it is wasteful without using the model information in the control system design. It is believed, particularly for man-made engineering systems/products (built up by rules/models), that any repetitive process and motion has a model existing in operation even though the model is yet to be identified.
U-Model-based control has claimed to radically relieve the dependence of plant model-oriented design foundation. The use of the plant model is effectively reduced as a reference for converting to U-model and accordingly to work out the control output [30]. U-Model-based control assumes the feasibility of using linear system design procedures to design the control of nonlinear dynamic plants with assigned response performances. The U-model control platform is illustrated in Figure 1.
U-Model-based control system design.
The U-model systematically converts smooth (polynomial and extended including transcendental functions) models, derived from principles or identified from measurements, into a type of U-based model to equivalently describe plant input-output relationship, so that it establishes a general platform to facilitate control system design and dynamic inversion. It should be mentioned that there is nothing lost with the derived U-models from the original nonlinear models. The difference between the two types of model expressions is that those original nonlinear models could be obtained from principles, such as Newton’s law, or identified directly from measured data; the U-models are derived from the original models in control design-oriented expressions. Regarding the U-control (U-model-based control) research status, Zhu and Guo [31] have brought forward a fundamental framework in terms of pole placement control for nonlinear systems. More recently, U-control has been expanded to general predictive control [32] and sliding mode control [30]. To accommodate the U-control of state space models, a backstepping algorithm is being expanded to extract the controller output within multiloop U-models. With the nature of separating control system design (specifying closed-loop performance) and controller output calculation (by resolving plant dynamic inversion through U-model), it can be forecast that the other classical control issues could be similarly formulated within a general and concise framework.
1.4. Organisation of the Study
The remaining study is organised into five major sections. Section 2 is used to define a generic framework of control-oriented U-model for representing smooth nonlinear dynamic plants. It is then expanded by including a rational model and transcendental functions as its subsets to lay a basis for applying linear control system design techniques. Section 3 proposes a general pole placement controller for nonlinear rational systems within the U-model framework. Section 4 shows design of an adaptive UPPC for the control of stochastic nonlinear rational systems. Section 5 tests a number of typical rational systems with the developed procedures and shows the exemplary procedures for potential users.
2. U-Model: A Generic Framework of Control-Oriented Nonlinear Plant Models2.1. U-Model Foundation: Polynomial [30]
Consider a general polynomial description of
(5)yk=NpYk−1Uk−1=∑i=0Lpikθi,where yk∈ℝ and uk∈ℝ denote the plant output and input, respectively, at time instant k=1,2,…. Np⋅∈ℝ is a real valued and smooth polynomial function and Yk−1∈ℝn⊃yk−1,…,yk−n and Uk−1∈ℝn⊃uk−1,…,uk−n denote the delayed outputs and inputs, respectively. pik∈ℝ denotes the model structure variables, e.g., uk−2y3k−1, uk−1u2k−3, yk−2yk−3, and θi∈ℝ denote the coefficients. To convert the above polynomial into a U-model, which is a polynomial with an argument of control input uk−1 (also called controller output while talking about control system design), it gives [30]
(6)yk=∑j=0Mλjkujk−1,where degree M is of controller output uk−1 and λk=λ0k…λMk∈ℝM+1 is the time-varying parameter vector, a function of absorbing past inputs Uk−2, outputs Yt−1, and parameters θnj in the original polynomial. An example illustrating the conversion to U-model from an ordinary polynomial is shown here. Consider a polynomial,
(7)yk=0.2yk−1yk−3+0.5uk−1uk−3−0.9yk−2u2k−1.
Rearrange polynomial (7) with
(8)yk=λ0k+λ1kuk−1+λ2tu2k−1,where λ0k=0.2yk−1yk−3, λ1k=0.5uk−3, and λ2k=−0.9yk−2.
Clearly, the time-varying λjk is absorbing the past inputs/outputs and parameters of the original polynomial, associated with ujk−1.
Property 1.
Assign φ:ℝL+1→ℝM+1 a U-mapping to convert the classical polynomial expression of (5) to its U-expression of (6) and the inverse be φ−1, that is
(9)fpi,θi→φfuj,λj
Thus, it has good mapping properties [30].
2.2. U-Mode: Rational
With reference to (1), its deterministic parametric rational expression is given below:
(10)yk=Np∗Dp∗=∑j=1numpnjkθnj∑j=1denpdjkθdj.
Its U-model realisation can be determined by removing the denominator to the left-hand side of (10); it gives
(11)ykDp∗=Np∗.
Convert (11) into U-model form to yield
(12)yk∑i=0Lγikuik−1=∑j=0Mλjkujk−1,where λjk∈ℝ is a function of past inputs Uk−2 and outputs Yk−1 and parameters θnj in the numerator polynomial. Similarly, γik∈ℝ is a function of past inputs Uk−2 and outputs Yk−1 and parameters θdj in the denominator polynomial. M and L are the degrees of the model input uk−1 in the numerator and denominator, respectively. Here is a simple example to show the conversion of
(13)yk=yk−11uk−1.
Inspection of (12) gives
(14)ykγ1kut−1=λ0k,where γ1k=1 and λ0k=yk−1.
In the following sections of the controller design, it is required to make a dynamic inversion of (12) to solve for roots.
There are many standard root-solving algorithms for such polynomial equations [30].
Remark 1.
Compared with polynomial U-realisation, it can be noted that rational model U-realisation is an implicit expression of yk due to the multiplicative item ykDpk.
2.3. U-Model: Extended
To describe more general nonlinear terms including those transcendental functions, define the extended U-model below:
(15)ykfbuk−1=fauk−1,where fbuk−1∈ℝ and fauk−1∈ℝ are smooth functions. In general, these can be expressed as
(16)fbuk−1=∑jfbjuk−1,fauk−1=∑jfajuk−1.
Here is a simple example to show its U-model representation; consider
(17)yk=yk−1sinuk−11+cos2uk−1.
For its U-model of (15), it gives
(18)fbuk−1=γ0k+γ1kcos2uk−1,fauk−1=λ1ksinuk−1,where γ0k=1, γ1k=1, and λ1k=yk−1.
3. Pole Placement Controller: A Show Case of the Design Procedure
The control objective is, for a desired trajectory vk, to determine a control input ut to drive the underlying system output yk to follow the desired trajectory vk with an acceptable performance (such as transient response and steady-state error), while all the inputs and outputs of the control system are bounded within the permitted ranges.
3.1. U-Control System Design
In general, there are three steps in the U-control system design routine:
Form a proper linear feedback control system structure, as shown in Figure 2. The controller, in the dashed line block, consists of two functions, the invariant controller Gc1 and the dynamic inverter Gp−1. The plant model is Gp.
U-Model control system.
Design the invariant controller Gc1 by linear control system approach. By letting Gp=1, therefore Gp−1=1, and specifying the desired closed-loop transfer function G, it gives Gc1=G/1−G and the invariant controller output vt is the desired output while the plant model is a unit constant.
Determine dynamic inverter Gp−1 to work out the controller output uk−1. Assuming the plant model is bounded-input/bounded-output (BIBO) stable and the inverse of Gp exists, expressing the plant model Gp in forms of U-model, lettingyk=vk in the U-model, gives model (15) in expression of vkfbuk−1=fauk−1. To determine control input ut−1 is to find the inverse by resolving the equation of vkfbuk−1−fauk−1=0.
It should be noted that the arrow line from the plant to the dashed line block represents the U-model update from the plant model at each time instance.
Proposition 1.
Generality: the U-model-based control allows a once-off design for all linear and nonlinear polynomial models. This is due to the controller Gc1 design being independent of model Gp.
Proposition 2.
Simplicity: the U-model-based control requires no repeated computation if a plant model is changed. Again, this is due to the controller Gc1 design being once-off and independent of model Gp, and changes to plant model Gp only change the U-model to resolve different roots. In comparison, almost all classical and modern control approaches are plant model-based designs; that is, the controller design is a function of both system performance and plant; accordingly, if the plant model is changed, the controller must be redesigned.
Proposition 3.
Feasibility for controller design of rational systems: this can be proved directly from Proposition 1 and U-realisation of the rational model in (12).
In formality, the U-adaptive control is very similar to deterministic U-model control. The difference is that the plant model is required to be estimated or updated online in the adaptive control.
For simplicity, but without losing generality, in formulation of the U-model (polynomial), once the invariant controller is designed, the real controller output can be determined by letting
(19)vk=∑j=0Mλjkujk−1.
Then resolving one of the roots from
(20)vk−∑j=0Mλjkujk−1=0.
3.2. Stability and Robust Analysis of U-Model Control Systems
There are two typical situations: ideal case—deterministic systems without modelling error and disturbance, and nonideal case—deterministic systems with modelling errors and/or disturbance.
Theorem 1. (bounded-input, bounded-output (BIBO) stability of deterministic U-model control systems).
Regarding the U-model control system shown in Figure 2, it is BIBO stable and tracks the bounded reference signal r properly while the following conditions are satisfied:
Invariant controller Gc1 is closed-loop stable; that is, all poles of the closed loop are located with the unit circle.
Plant model Gp is a bounded-input/bounded-output (BIBO).
The inverse of the plant model Gp−1 exits.
Proof.
With reference to Figure 2, it has Gp−1Gp=1 from the conditions (ii) and (iii). Accordingly, the closed-loop transfer function is given in terms of Gc1/1+Gc1, which is stable from (i), and thus, the tracking performance is given by rGc1/1+Gc1.
Remark 2.
This establishes a framework for designing control for both linear and nonlinear dynamic plants. It is feasible, simple, general, and with no repetition of controller design on changes to the plant model, except the computation of the inversion of the changed plant U-model polynomial. In other words, this is a new methodology for minimising the complexity induced by the plant model in control system design, which is particularly important for nonlinear plants. U-model, as a universal dynamic inverter, is the key to achieve the goals.
Theorem 2. (BIBO stability of uncertain U-model control systems).
Regarding the U-model control system structured in Figure 2, modelling error and/or disturbance εUt can be treated as an external disturbance as shown in Figure 3. It is BIBO stable and tracking the reference signal with a bounded error while the following conditions are satisfied:
Invariant controller Gc1 is closed-loop stable.
Plant model Gp is a bounded-input/bounded-output (BIBO).
The inverse of the plant model Gp−1 exits.
The upper bound of modelling error and/or disturbance εUt is satisfied with the conditions of small gain robust stability [33].
Proof.
In Figure 3, Gp−1Gp=1; this gives y=rGc1/1+Gc1+εU/1+Gc1.
Then the stability of Figure 3 is the same as in Figure 2 while the upper bound εUt is satisfied with the small gain robust stability criterion.
Remark 3.
It should be noted that the tracking error is determined by εU/1+Gc1; therefore, a properly designed Gc1 will have a degree of robustness against uncertainties/disturbance.
Uncertain U-model control system.
4. Design of Pole Placement Controller
A classical approach [34] has been selected to formulate the U-model-enhanced pole placement controller (UPPC) [30, 31]. Here, a further refined version of UPPC is presented. Within the U-model framework, closed-loop control system performance is independently specified without involving the plant model. Therefore, the classical version involving plant model can be simplified as below.
(21)Rvk=Trk−Syk,and
(22)R=zn+r1zn−1+⋯+rn,T=t0zm+t1zm−1+⋯+tm,S=s0zl+s1zl−1+⋯+sl,with rk for reference, vk for invariant controller Gc1 output, and yk for plant output. The polynomials R, S, and T, with backward shift operator z−1 and proper orders (n, m, and l), are used to specify closed-loop control system performance.
To guarantee that the control system is realistically implementable, specify
(23)OS<OR⇔l<n,OT≤OR⇔m≤n,where the operator O∗=order∗ denotes the order of the concerned linear polynomial.
With reference to (19), two control roles can be assigned with negative feedback −R/S for stabilising the closed-loop system with requested dynamics and feedforward T/R for reducing steady-state errors. The structured control system is shown in Figure 4.
Structured UPPC control system.
For designing an invariant controller, let vt=yt in (19); thus, it gives the closed-loop transfer function
(24)yk=TR+Srk=TAcrk.
Accordingly, the required design task is to assign the closed-loop denominator polynomial Ac and the numerator polynomial T.
It should be noted that after Ac is specified (by customers and/or designers), a routine for resolving Diophantine is needed to work out the parameters of polynomials R and S from the following relationship:
(25)R+S=Ac.
To achieve zero steady state, T can be designed with
(26)T=Ac1.
The detailed design procedure and examples can be refereed to [31].
Remark 4.
Compared with classical pole placement control design procedures [34], the UPPC is more concise and independent of the plant model, which results in the UPPC being generalised to any plant model structure and once-off designed. For each different plant model, this task is merely the resolving of the U-model to obtain one of the roots as the operational controller output. The relevant comparison details can be referred in [30].
5. U-Model-Based Pole Placement Control with Adaptive Parameter Estimation
The U-model-based adaptive control schematic diagram is shown in Figure 5. Again, this U-model adaptive control is different from those classical adaptive/self-tuning control approaches in terms of control structure. The feedback controller parameters are not tuned and thereafter are fixed: the only adaptation is to update U-model parameters to accommodate the plant model parameter variation and/or external disturbance, which is consistent with Propositions 1, 2, and 3.
Adaptive U-model control system.
In general, an adaptive control system can be considered as a two-layer system, that is:
Layer 1: conventional feedback control
Layer 2: adaptation loop
In this study, the UPPC presented in Section 3 is selected to form the conventional feedback control. Thus, this section mainly develops this adaptation loop formulation.
In recursive formulation, there are two ways to estimate the U-model parameters in the adaptation loop.
Indirect parameter estimation: estimate the original rational model parameters (θnjk,θdjk) first and then convert into U-model parameters λjk. The challenging issue is that classical recursive least squares estimation algorithms give biased estimates and recursive rational model estimators need noise variance information in advance [11, 18].
Direct parameter estimation: estimate the U-model parameters λjk directly. The challenging issue is that the parameters λjk, while converted from a rational model, are time varying at every sampling time. It has been proved [35] that for time-varying stochastic models, the parameter estimation errors (PEE) with the well-known forgetting factor least squares (FFLS) algorithm are bounded and the FFLS is capable of reducing the squared measurement error (the difference between measured output and model-predicted output); even the time-varying parameter estimates are not converged to their real values.
In this study, a FFLS estimator [36] is selected with the following formulations:
(27)εUk=yk−ΨTkλ^k−1,Kk=Pk−1Ψkρ+ΨTkPk−1Ψk,λ^k=λ^k−1+KkεUk,Pk=Ι−KkΨTkPk−1,where vector λ^k=λ^0kλ^1k…λ^MkT∈ℝM is the estimate of λk; εk∈ℝ is the error, that is, the difference between the measured output and the model-predicted output; Kt∈M+1×1 is the weighting factor vector indicating the effect of εt to change the parameter vector; Ψk=1uk−1…uMk−1T∈ℝM is the input vector at time k − 1; ρ is the forgetting factor (a number less than 1, e.g., 0.99 or 0.95, represents a trade-off between fast tracking and noisy estimate), the smaller the value of ρ, the quicker the information in previous data will be forgotten; and Pk∈ℝM+1×M+1 is the covariance matrix.
In presenting the stability of the proposed adaptive U-control, expand the virtual equivalent system (VES) concept and methodology [37] for the analysis, which is an alternative insight and judgement of the stability/convergence for adaptive control systems. Following the similar arguments as shown before, we assume Gp−1Gp=1, and the invariant controller Gc1 is well defined to stabilise conventional feedback control systems and track the bounded reference signal in terms of mean squares. Then for a slow time-varying parameter model (because it is converted from its original time-invariant parameter model referred to in (5) and (6)), the U-model parameter estimation errors εUt are bounded with FFLS or the other recursive algorithms [35, 38]. In this case, using Figure 3 again, knowing εUt includes U-model parameter estimation errors. Hence, in terms of VES, the adaptive control system can be treated as a summation of two subsystems of
(28)y=y1+y2=rGc11+Gc1+εU1+Gc1.
As εUt is bounded, the adaptive control system is stable and the tracking control error will converge to a bounded compact set around zero, whose size depends on the ultimate bounds of estimation error εU.
Remark 5.
The U-model provides a platform for simplifying control system design, and VES provides a platform for simplifying the analysis of stability and convergence of general adaptive control systems.
6. Simulation Studies
Four case studies have been conducted to initially validate the new design procedure. It should be made clear that there is no comparison result that can be provided as this is the first study in the control of such nonlinear rational systems.
As described before, the design is split into two stages, design invariant control Gc1 (thus, vk by pole placement) and determination of the controller output uk−1 by resolving plant U-model equation.
To design the pole placement controller, assign the characteristic equation
(29)Ac=z2−1.3205z+0.4966.
Factorisation of (29) gives the closed-loop poles as 0.6603±0.2463i; this gives a decayed oscillatory response (ωn=1, ζ=0.7), which is a commonly used dynamic response index. For steady-state error performance, making its error zero gives
(30)T=Ac1=1−1.3205+0.4966=0.1761.
From the causality condition, specify the structures of R and S with
(31)R=z2+r1z+r2,S=s0z+s1.
Form a Diophantine equation with polynomials Ac, R, and S [30] to yield
(32)r2+s1=0.4966,r1+s0=−1.3205.
To make polynomial R stable and having the requested response, assign r1=−0.06 and r2=0.0005, which give two poles z−0.05andz−0.01. Then the coefficients of polynomial S are resolved in the Diophantine equation of (32) as follows.
(33)s0=−1.2605,s1=0.4961.
Consequently, controller (19) can be recursively implemented to calculate the virtual controller output vt:
(34)vk+1=0.06vk−0.0005vk−1+0.1761rk−1+1.2605yk−0.4961yk−1.
Case 1 (feasibility test of U-control of nonlinear rational systems).
Consider a rational system modelled by
(35)yk=0.5yk−1uk−1+u3k−11+y2k−1+u2k−1,
where yk is the plant output and uk is the input of the model or controller output. This is used to test deterministic feedback control. The model structure has been typically investigated in system identification. Accordingly, its U-realisation can be expressed as
(36)yk1+y2k−1+u2k−1=0.5yk−1uk−1+u3k−1.
To obtain the dynamic inverter Gp−1 output, that is, the controller output ut, let yk=vk; then it gives rise to
(37)vk1+y2k−1+u2k−1=0.5yk−1uk−1+u3k−1.
To determine the control input uk−1, form a U-model equation from (37) as
(38)λ0k−λ1kuk−1+λ2ku2k−1−λ3ku3k−1=0,where
(39)λ0k=vk1+y2k−1,λ1k=0.5yk−1,λ2k=vt,λ3k=1.
In this simulation, the operation time length was configured with 400 sampling points and the reference was a sequence of multiamplitude steps. The achieved output response and controller output are shown in Figures 6(a) and 6(b), respectively.
Case 2 (test of external disturbance).
Consider a stochastic rational system modelled by
(40)yk=0.5yk−1uk−1+u3k−11+y2k−1+u2k−1+ek,
where yk is the plant output, uk is the input of the model or controller output, and ek is Gaussian noise representing an unknown disturbance acting on the controlled plant output.
This case study was used to test adaptive feedback control. The feedback control loop has been designed as in Case 1; that is, all configurations for feedback control were kept as those used in Case 1. For the adaptation loop, the disturbance was configured with ek~N0,0.01, the initial covariance matrix with Pk=106I4, and the forgetting factor with ρ=0.95 to deal with fast time-varying parameter estimation; the initial parameter vector was randomly assigned with λ^0=λ^00λ^10λ^30λ^40T=0.30.20.10.1T; and the input vector was specified with Ψk=1uk−1u2k−1u3k−1T. The achieved output response and controller output are shown in Figures 7(a) and 7(b), respectively.
Case 3 (test of internal parameter variation).
The same model structure as Case 1 is used, but the parameter associated with yk−1 and uk−1 is time varying representing internal parameter disturbances, such as worn parts in mechanical and electrical systems.
(41)yk=akyk−1uk−1+u3k−11+y2k−1+u2k−1.
In simulation, all the setups were the same as those used in Case 1. The parameter variation was configured as
(42)ak=0.9,120≤k≤250,0.5,otherwise.
The adaption loop, specified as in Case 2, was used to follow the plant model internal structure variation. The achieved output response and controller output are shown in Figures 8(a) and 8(b), respectively. Inspecting the simulation results, the output of the systems are seen to track the reference signals after a short transient phase. U-model parameter estimation is shown in Figure 9. It should be noted that this estimated parameter vector is to achieve smaller squared error between the measured output and model-predicted output. Therefore, the estimates are not converged to those real time-varying parameters in the U-model. In the future, studies to deal with time-varying parameter estimation will be conducted in terms of reducing both squared measurement errors and squared dynamic errors [39].
Case 4 (feasibility test of U-control of extended nonlinear rational systems).
This study is used to test the U-control of extended rational systems with transcendental input and delayed output.
(43)yk=0.5yk−1+sinuk−1+uk−11+exp−y2k−1,
where yk is the plant output and uk is the input of the model or controller output. Accordingly, the extended U-model can be expressed as
(44)yk1+exp−y2k−1=0.5yk−1+sinuk−1+uk−1.
With the same controller designed in (44) above, assigning the output yk of (44) with the desired output vk of (34) gives
(45)vk1+exp−y2k−1=0.5yk−1+sinuk−1+uk−1.
Therefore, the control input uk−1 can be solved by
(46)vk1+exp−y2k−1−0.5yk−1−sinuk−1−uk−1=0.
The achieved output response and controller output are shown in Figures 10(a) and 10(b), respectively. Once again, the computational experiment confirms the feasibility of U-control.
Plant output and control input.
Plant output response
Control input
Plant output and control input.
Plant output
Control input
Plant output and control input.
Plant output
Control input
U-Model parameter estimates.
Plant output and control input.
Plant output
Control input
7. Conclusions
A fundamental question is raised in this study and those for the other U-model-enhanced controls: after two generations of plant model- (polynomial and state space) centered control system design research/applications, what is the next generation of development? Should the research for new model structures continue, or should control systems be designed without such plant model requirements (possibly implying separation of control system design and controller output determination)?
One of the feasible choices in the future progression could be the U-control design methodology, which radically reduces the complexity of plant model-oriented design methods. The proposed U-control method provides a platform (1) with a universal control-oriented structure to represent existing models, (2) separating closed control system design from plant model structure (no matter whether linear or nonlinear or polynomial or state space), (3) where all well-developed linear control system design methods can be expanded in parallel to nonlinear plant models, (4) with a supplementary to all existing control design methods. Accordingly, this study is a show case using the U-model framework to design the control of the nonlinear rational systems with classical linear design approaches. Further study on the rational model control could derive concise algorithms for robust and adaptive control with reference to the recent research development [40, 41].
This foundation work has put an emphasis on formulation of structure in a systematic approach. Rigorous mathematical considerations should be followed to establish a comprehensive description and explanation.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors acknowledge Dr. Steve Wright for English proof reading. Finally, the authors are grateful to the anonymous reviewers for their constructive comments and suggestions with regard to the revision of the paper.
NemcovaJ.2009AmsterdamCentrum Wiskunde & Informatica (CWI)ZhuQ. M.WangY. J.ZhaoD. Y.LiS. Y.BillingsS. A.Review of rational (total) nonlinear dynamic system modelling, identification, and control201546122122213310.1080/00207721.2013.8497742-s2.0-84928755361BillingsS. A.2013Chichester, West SussexWiley, John & Sons10.1002/97811185355612-s2.0-84939520302WangL. X.1994Englewood Cliffs, NJ, USAPrentice HallDimitrovS. D.KamenskiD. I.A parameter estimation method for rational functions199115965766210.1016/0098-1354(91)87027-72-s2.0-0026226398KamenskiD. I.DimitrovS. D.Parameter estimation in differential equations by application of rational functions199317764365110.1016/0098-1354(93)80052-o2-s2.0-0027622717FordI.TitteringtonD. M.KitsosC. P.Recent advances in nonlinear experimental design19893114960x10.1080/00401706.1989.104884752-s2.0-0024606752PontonJ. W.The use of multivariable rational functions for nonlinear data representation and classification199317101047105210.1016/0098-1354(93)80086-32-s2.0-0027675811KambhampatiC.MasonJ. D.WarwickK.A stable one-step-ahead predictive control of nonlinear systems200036448549510.1016/S0005-1098(99)00173-92-s2.0-0033881654BillingsS. A.ZhuQ. M.Rational model identification using an extended least-squares algorithm199154352954610.1080/002071791089341742-s2.0-0026218753ZhuQ. M.BillingsS. A.Recursive parameter estimation for nonlinear rational models199116367BillingsS. A.ZhuQ. M.A structure detection algorithm for nonlinear dynamic rational models1994591439146310.1080/002071794089231402-s2.0-0000719723ZhuQ. M.BillingsS. A.Fast orthogonal identification of nonlinear stochastic models and radial basis function neural networks199664587188610.1080/002071796089216622-s2.0-0030192825ZhuQ. M.An implicit least squares algorithm for nonlinear rational model parameter estimation200529767368910.1016/j.apm.2004.10.0082-s2.0-18744412956BillingsS. A.ChenS.Identification of non-linear rational systems using a prediction-error estimation algorithm198920346749410.1080/002077289089101432-s2.0-0024628610MuB. Q.BaiE. W.ZhengW. X.ZhuQ. M.A globally consistent nonlinear least squares estimator for identification of nonlinear rational systems20177732233510.1016/j.automatica.2016.11.0092-s2.0-85009466002ZhuQ. M.A back propagation algorithm to estimate the parameters of nonlinear dynamic rational models200327316918710.1016/S0307-904X(02)00097-52-s2.0-0037367056ZhuQ. M.YuD. L.ZhaoD. Y.An enhanced linear Kalman filter (EnLKF) algorithm for parameter estimation of nonlinear rational models201748345146110.1080/00207721.2016.11862432-s2.0-84969800031BillingsS. A.ZhuQ. M.Nonlinear model validation using correlation tests19946061107112010.1080/002071794089215132-s2.0-0001658853ZhuQ. M.ZhangL. F.LongdenA.Development of omni-directional correlation functions for nonlinear model validation20074391519153110.1016/j.automatica.2007.02.0102-s2.0-34547600511NarendraK. S.ParthasapathyK.Identification and control of dynamical systems using neural networks19901142710.1109/72.802022-s2.0-002539956718282820RomanchukB. G.SmithM. C.Incremental gain analysis of piecewise linear systems and application to the antiwindup problem19993571275128310.1016/S0005-1098(99)00023-02-s2.0-0033165186OzkanL.KothareM. V.GeorgakisC.Model predictive control of nonlinear systems using piecewise linear models2000242–779379910.1016/S0098-1354(00)00376-82-s2.0-0034661317TsujiT.XuB. H.KanekoM.Adaptive control and identification using one neural network for a class of plants with uncertainties199828449650510.1109/3468.6867112-s2.0-0032124215ZhuQ. M.MaZ.WarwickK.Neural network enhanced generalised minimum variance self-tuning controller for nonlinear discrete-time systems1999146431932610.1049/ip-cta:199903642-s2.0-0032598248IsidoriA.MarconiL.SerraniA.New results on semiglobal output regulation of nonminimum-phase nonlinear systemsProceedings of the 41st IEEE Conference on Decision and Control, 2002December 2002Las Vegas, NV, USA1467147210.1109/CDC.2002.1184726SlotineJ. J. E.LiW.1991LondonPrentice-HallLiY. Q.HouZ. S.FengY. J.ChiR. H.Data-driven approximate value iteration with optimality error bound analysis201778798710.1016/j.automatica.2016.12.0192-s2.0-85010380397FliessM.JoinC.Model-free control201386122228225210.1080/00207179.2013.8103452-s2.0-84888882608ZhuQ. M.ZhaoD. Y.ZhangJ. H.A general U-block model-based design procedure for nonlinear polynomial control systems201647143465347510.1080/00207721.2015.10869302-s2.0-84941800860ZhuQ. M.GuoL. Z.A pole placement controller for nonlinear dynamic plant2002216646747610.1177/095965180221600603DuW. X.WuX. L.ZhuQ. M.Direct design of a U-model-based generalized predictive controller for a class of non-linear (polynomial) dynamic plants20122261274210.1177/09596518114096552-s2.0-84863048057KravarisC.WrightR. A.Deadtime compensation for nonlinear processes19893591535154210.1002/aic.6903509142-s2.0-0024732521AstromK. J.WittenmarkB.19952ndReading, MA, USAAddison-WesleyDingF.ChenT.Performance bounds of forgetting factor least-squares algorithms for time-varying systems with finite measurement data200552355556610.1109/TCSI.2004.8428742-s2.0-16344391042SoderstromT.StoicaP.1989Hemel HempsteadPrentice Hall InternationalZhangW. C.On the stability and convergence of self-tuning control–virtual equivalent system approach201083587989610.1080/002071709034874212-s2.0-77951564730DingF.XuL.ZhuQ. M.Performance analysis of the generalised projection identification for time-varying systems201610182506251410.1049/iet-cta.2016.02022-s2.0-85006797844KalabaR.TesfatsionL.Time-varying linear regression via flexible least squares1989178-91215124510.1016/0898-1221(89)90091-62-s2.0-0024861578NaJ.HerrmannG.ZhangK. Q.Improving transient performance of adaptive control via a modified reference model and novel adaptation20172781351137210.1002/rnc.36362-s2.0-84987679594NaJ.MahyuddinM. N.HerrmannG.RenX. M.BarberP.Robust adaptive finite-time parameter estimation and control for robotic systems201525163045307110.1002/rnc.32472-s2.0-84943200482