Model Based Adaptive Piecewise Linear Controller for Complicated Control Systems

A model based adaptive piecewise linear control scheme for industry processes with specifications on peak overshoots and rise times is proposed. It is a gain stabilized control technique. Large gain is used for large tracking error to get fast response. Small gain is used between large and small tracking error for good performance. Large gain is used again for small tracking error to cope with large disturbance. Parameters of the three-segment piecewise linear controller are found by an automatic regulating time series which is function of output characteristics of the plant and reference model.The time series will be converged to steady values after the time response of the considered system matching that of the reference model. The proposed control scheme is applied to four numerical examples which have been compensated by PID controllers. Parameters of PID controllers are found by optimization method. It gives an almost command independent response and gives significant improvements for response time and performance.


Introduction
Gain and phase stabilized are two conventional design methods for feedback control systems.They can be analyzed and designed in gain-phase plots to get wanted gain and phase margins or gain and phase crossover frequencies [1,2].The gain crossover frequency is closely related to the system bandwidth (or response time).The phase margin is closely related to performance (or peak overshoot).In general, fast response time and good performance cannot be obtained simultaneously for some feedback control systems.For example, the altitude control system of the airframe needs large gain for fast response time and low gain for good performance.It is in conflict with another.Simple and effective way to solve this problem and provides better results than those of compensated by linear controllers are generally expected.This is the motivation of this paper.
Variable structure control is a switching control method for feedback control systems [3][4][5][6][7].It gives good performance and robustness for coping with system uncertainty.But it suffers from chattering problem and state measurements.In this paper, a fast response system and a good performance system are selected for switching.An adaptive switching algorithm is used.There is no discontinuous connection between two systems.Therefore, there is no chattering problem.Gain scheduling has been used successfully to control nonlinear systems for many decades and in many different applications, such as autopilots and chemical processes [8][9][10].It consisted of many linear controllers for operating points to cope with large parameter variations.This concept will be expanded for response time and performance.Operating points are replaced by fast response and good performance conditions and interpolation for gain evaluation is replaced by an adaptive switching point.It is determined by the filtered command tracking errors.Nonlinear controller's syntheses using inverse describing function for use with hard nonlinear system have developed for several researchers [11][12][13].They are complicated but effective for nonlinear systems.
In this paper, a simple three-segment piecewise linear controller is proposed.It is symmetry and has two switching points.It is equivalent to two quasilinear gains; that is, there are two equivalent systems.One is the fast but bad performance system, and the other is the slow but good performance system.Switching points are used to select optimal gain.They are functions of the dynamic of tracking error.It is easy to analyse and design.Based on the information of the reference model, parameters of three-segment piecewise linear controller can be tuned in the online manner.An automatic regulating time series [14] is used for tuning the parameter.It can be retuned again if the system is needed.Furthermore, it gives an almost reference input independent response.It implies that the controlled system is similar to a linear system.
The proposed control scheme will be illustrated by a 2nd order altitude control system and applied to four numerical examples those have been compensated by optimized PID controllers.Time responses show that the proposed method gives significant improvements on response time and performance than those compensated by PID controllers.

The Model Based Control Scheme
Figure 1 shows the proposed model based control scheme in which reference model represents the wanted system response, adaptive piecewise linear controller is a gain stabilized compensation, and autotuning process is used to find the switching point ( 1 ) of the adaptive piecewise linear controller.The piecewise linear controller can get fast response and performance simultaneously.It is a threesegment piecewise linear controller.The switching point ( 1 ) is dependent on the filtered dynamic of the tracking error ().Parameters  1 ,  2 ,   , and   are found by the autotuning process for matching outputs of the reference model.Several identifying processes will be processed to find them.In general, the controlled system can be described by peak overshoot and rise time of step response.The concept of the proposed control scheme will be discussed by an illustrating example in next subsections.
2.1.The Adaptive Piecewise Linear Controller. Figure 1 shows the three-segment piecewise linear nonlinearity.Now, the problem is to find values of   and   for switching points + 1 and − 1 and gains ( 1 ,  2 ) for two selected systems; that is, one is the fast but bad performance system and the other is the slow but good performance system.Note that switching points + 1 and − 1 are not fixed and will be determined by the absolute value of the command tracking error.
Consider a 2nd altitude control example described by The transfer function of the closed-loop system with loop gain  is Poles locations and natural frequency (  ) for two loop gains ( 1 ,  2 ) are They are overdamped and underdamped systems.Time responses of them are shown in Figure 2 for  =  1 (smalldot-line) and  =  2 (large-dot-line) in which  represents the reference input and  represents the output.The strategy for gain switching is (1) large gain ( 2 ) for large tracking error to get fast response and (2) small gain ( 1 ) for small tracking error () to get good performance.It is a variable structure system and can be achieved by selecting a proper switching point  1 of the piecewise linear controller.For example, the optimal switching point  1 is selected as 0.525 for  = 1 to get fast response and good performance.Large gain ( 2 ) is used for || >  1 and small gain ( 1 ) is used for || ≤  1 .
Step response is shown in Figure 2 (solid-line) also for  = 1.It shows that the controlled system using adaptive gain can give fast response and good performance simultaneously.However, it is not true that  is equal to 5, 10, and 50, respectively.Those step responses are shown in Figure 3. Naturally, another optimal switching point  1 for  = 5, 10, and 50 can be selected for getting good performance.They are 2.625, 5.250, and 26.250 for  = 5, 10, and 50, respectively.They are true for step responses from zeros to 5, 10, and 50 only.Another possible way for the switching point can be dependent on the tracking error ().A possible switching rule for  1 is found as  1 = 0.925|| for good performance.Figure 4 shows time responses for  = 1, 5, 10, and 50, respectively.It can be seen that the switching rule gives an input command () independent of results.However, they are slower than those results shown in Figures 2 and 3.
One possible way to speed up the time response is to enlarge the large gain phase in the beginning.A low-pass filter () =   /(   + 1) is used for the absolute tracking error () to get  1 .Figure 5 shows that faster response than that shown in Figure 4 is got for   = 1.0465 and   = 1/  .The switching point  1 is shaped for speeding up the responses while keeping performance unchanged.Figure 5 shows that the rise time of the controlled system is quite close to that of the controlled system using  =  2 ; that is, the rise time of the system is dependent on  2 .Note that the natural frequency (  ) for  =  2 is used to find   also.Therefore, it is needed to determine   only.
Figure 6 shows time responses for  = 1, 5, 10, and 50, respectively.It shows that almost input independent responses are got.It implies that the controlled system is similar to the system controlled by linear controller.This is the major merit of the proposed method.

The Automatic Tuning Process.
In general, one can use rise time and peak overshoot of step response to describe the characteristics of the controlled system.They will be used as specifications to find parameters of the adaptive piecewise linear controller, that is,  1 ,  2 ,   , and   .

Automatic Tuning Process for Matching Peak Point.
Parameter   = 1.0465 of the low-pass filter

Adaptive piecewise linear controller
Autotuning process

Reference model
Linear compensator   /(   + 1) of the adaptive piecewise linear controller used in Section 2.1 was found by use of the following equations: where   () is an automatic regulating time series [14];  and  are controlling parameters of the regulating property;  is the specification of the peak value;  is the peak point found by using   =   ();  is simulation period of one step response; and  is the th step responses.Equation (4) gives that   () will be converged to a steady-state value after  being matched by .
for three initial guesses   (0).They show that   () are converged to 1.0465 quickly.Figure 7 shows also that the final system matches the peak specification .

Automatic Tuning Process for Matching Peak Point and
Rise Time.In the previous subsection, specification  is used to find parameter   .Another specification for rise time   can be also used to tune high gain  2 and   .Automatic tuning equations are described as follows: where  is regulating parameter and  2 ,  1 , and   will be replaced by  2  2 ,  1 / 2 , and   / 2 , respectively.Figure 8 shows the online automatic tuning process for   4) and ( 7) are  = 0.6,  = 0.5,  = 5,  = 30 sec,   (0) = 0.7, and   (0) = 1.00.Internal values of each tuning interval are given in Table 1.Table 1 gives that system response is converged into wanted specification quickly.Note that system response shown in Figure 8 is faster than that shown in Figure 7.  Design procedures for the proposed method using the proposed adaptive piecewise linear controller can be deduced as follows.
Step 1. Select two loop gains for fast response and good performance, respectively.In general, high loop gain ( =  2 ) is for fast responses and low gain ( =  1 ) is for good performance.The rise time of the system with high gain meets the design specification.The peak overshoot of the system with low gain meets the design specification.
Step 2. Determine parameters of low-pass filter () =   /(   + 1) to find the optimal switching point  1 .The natural frequency (  ) for the high gain system ( =  2 ) is used to find   .The natural frequency (  ) is closely related to the rise time.Equations ( 4)-( 8) are used to find parameters  1 ,  2 ,   , and   for given specifications  and   .Another method for finding parameter   can be used by the optimization method using performance index formulated by integration of the absolute error (IAE) and integration of the square error (ISE) oronline parameterized method [15,16].
The proposed control scheme will be applied to four numerical examples which have been compensated by optimized PID controllers.

Numerical Examples
Example 1.Consider a stable plant that has the transfer function [17] It is a second order dynamic plus a pure time delay (SOPDT).
In this example, parameters of the PID controller are It is designed by optimization method first and closed with it.And then low loop gain  1 = 0.50 in cascade is selected and high loop gain  2 = 2.587 in cascade is selected for the system which is just in the sustaining oscillating condition.Time responses using low gain ( 1 = 0.50) and high gain ( 2 = 2.587) are shown in Figure 9.The oscillation frequency is   = 1.5708 rad/s.It shows an overdamped system and a zero-damped system.An optimal parameter   for switching point  1 will be selected to coordinate two systems for matching system specifications.Now, apply the proposed control scheme to the system using  1 = 0.5000,  2 = 2.587,   = 0.6366.
() is converged to be 1.1385 within three period simulations.The online autotuning process is shown in Figure 10.The time response is shown in Figure 9 also.It can be seen that the proposed method can give fast response and good performance simultaneously.It is the combination of overdamped and zero-damped systems with  1 .Zerodamped system is used for fast responses and overdamped system is used for good performance.Naturally, it is almost input command () independent also.
Simulation results of the proposed method and four other methods are presented for comparisons.They are methods of Ziegler-Nichols [18,19] for finding PI and PID compensators, Tan et al. [20,21] for finding PID compensator, and Majhi [17] for finding PI compensator.Parameters of five found compensators are given as follows.
(1) Proposed method:  2. From Table 2 and Figure 11, one can see that the proposed method gives faster response and better performance than those of other methods presented.
Example 2. Consider the sixth-order plant [17]  () = 1 ( + 1) 6 . ( Three parameters of the optimized PID controller are The performance index 50IAE + ISE is used for finding them.The dominate pole for high gain  2 = 1.500 is −0.0848 ± 0.7898.The natural frequency is   = 0.7943 rad/s.It is an underdamped system.The low gain  1 = 0.05 is selected.It is an overdamped system.Parameters for the proposed control scheme found are  Figure 12 shows time responses using optimization method and the proposed method for  = 1, 5, 10, and 50.It can be seen that the proposed method gives better result than that of the optimization method.Simulation results of the proposed method and four other methods are presented for comparisons also.They are rules of Ziegler-Nichols [18,19] for finding PI and PID compensators, Ho et al. [22] for finding PID compensator, and Majhi [17] for finding PI compensator.Parameters of five found compensators are given as follows. (1) The proposed method:  Time responses are shown in Figure 13.Integral of the square error (ISE) and integral of the absolute error (IAE) are given in Table 3. From Table 3 and Figure 13, one can see that the proposed method gives faster response and better performance than those of other methods.
Example 3. Consider the very high order plant [17]  () = 1 ( + 1) 20 . ( Three parameters of the optimized PID controller are The performance index 50IAE + ISE is used for finding them. The dominate pole for high gain  2 = 2.000 is +0.0089 ± 0.1965 (  = 0.1967).Note that it is an unstable system.The low gain  1 = 0.200 is selected.It is an overdamped system.Parameters found for the proposed control scheme are Figure 14 shows time responses using optimization method and the proposed method for  = 1, 5, 10, and 50.It can be seen that the proposed method gives better result than that of the optimization method.Note that the compensated system is a combination of an overdamped system and an unstable system.That is, the proposed control scheme can stabilize the system and gives fast response and good performance.Simulation results of the proposed method and four other methods are presented for comparisons.They are methods of Ziegler-Nichols for finding PI and PID compensators [18,19], Zhuang and Atherton [23] for finding PI compensator, and Majhi [17] for finding PI compensator.Parameters of five found compensators are given as follows.
Integral of the square error (ISE) and integral of the absolute error (IAE) are given in Table 4. Time responses are shown in Figure 15.From Table 4 and Figure 15, one can see that the proposed method gives faster time responses and better performance than those of other mentioned methods.Example 4. Consider an electrohydraulic velocity/position servo control system [24] shown in Figure 16.The relation between the servo spool position  V and the input voltage  is in the form of where  V is the valve gain,  V is the damping ratio of the servo valve, and  V is the natural frequency of the servo valve.In general, (23) can be approximated by  V =  V  for large  V .The relation between the valve displacement   and the load flow rate   is governed by the well-known orifice law [24]: where   is the volumetric displacement;   is the total leakage coefficient;   is the total volume of the oil;   is the bulk modulus of the oil; and  is the velocity of the motor shaft.The torque balance equation for the motor is in the form of where   is the viscous damping coefficient and   is the external load disturbance which is assumed to be dependent upon the velocity of the shaft.The mathematical model of the considered system is shown in Figure 17.System parameters are given as follows: The control configuration for velocity and position servo control of the considered system is shown in Figure 18, in which inner loop and outer loop adaptive nonlinear controllers are included.Design results of the velocity control loop are discussed as follows.Time responses of the controlled system using the found PI controller are shown in Figure 19.

(b) Parameters of Inner Loop Adaptive Nonlinear Controller.
Low gain ( 1 = 1) and high gain ( 2 = 9.223) are selected.The low gain case is the optimized result and the high gain case is the controlled system in the sustaining condition (  = 312.71rad/s).The   = 312.71rad/s gives   = 0.003182.  = 2.2889 is found by (12) using  = 1.001.Time responses for low gain, high gain, and adaptive gain are shown in Figure 19.Rise times of the optimization method and the proposed method are 0.0202 sec and 0.0124 sec, respectively.It shows that the proposed method can give faster response than that controlled by the optimized method.The gain/phase margins, phase/gain crossover frequencies, and rise times are given also in Table 5.It gives that controlled system using two methods has the same robustness while Figure 19 shows that the proposed method gives faster response.
Design results of the position control loop are discussed as follows.
(a) Outer Loop PI Controller.The PI controller is first found by the optimizations toolbox of MATLAB for minimizing the integration of absolute errors (IAE), integration of square errors (ISE), and zero peak overshoot.Parameters of the PI controller are   = 18.506 and   = 0.3666.Time responses of the controlled system using the found PI controller are shown in Figure 20.
(b) Parameters of Outer Adaptive Nonlinear Controller.Low gain ( 1 = 1) and high gain ( 2 = 7.877) are selected.The low gain case is the optimized result and the high gain case is the controlled system in the sustaining condition (  = 91.95rad/s).The   = 91.95rad/s gives   = 0.0010875.  = 13.5 is found by (12) using  = 1.001.Time responses for low gain, high gain, and adaptive gain are shown in Figure 20.Rise times of the optimization method and the proposed method are 0.0513 sec and 0.0334 sec, respectively.It shows that the proposed method can give faster response than that controlled by the optimized method.Gain/phase margins, phase/gain crossover frequencies, and rise times are given also in Table 6.It gives that controlled system using two methods has same robustness while Figure 20 shows that the proposed method gives faster response.

Conclusions
The proposed adaptive piecewise linear controller has shown that controlled systems were reference input independent and both good performance and fast response were obtained simultaneously.Three-segment piecewise linear controller provided a switching algorithm for low gain and high systems, that is, low gain for performance and high gain for  Four numerical examples were designed and comparisons were made with six famous online computing and control methods.They have illustrated better performance and faster response provided by the proposed method than those of other mentioned methods.

Figure 7 :
Figure 7: Online automatic tuning process of the illustrating example using different initial   (0).

Figure 8 :
Figure 8: Online automatic tuning process of the illustrating example for matching peak point and rise time specification.

Figure 9 :
Figure 9: Step responses of Example 1 for constant gains ( = 0.5 and 2.587) and adaptive gain with  1 .

Figure 15 :
Figure 15: Comparisons with other methods for Example 3.

Figure 16 :
Figure 16: Block diagram of the electrohydraulic system.

Figure 17 :
Figure 17: Mathematical model of the electrohydraulic system.

Figure 18 :Figure 19 :
Figure 18: Control configuration of velocity and position servo control system.

Table 1 :
Internal values of each tuning interval.

Table 2 :
The ISE and IAE of five methods for Example 1.

Table 3 :
The ISE and IAE of five methods for Example 2.

Table 4 :
The ISE and IAE of five methods for Example 3.   is a constant for specific hydraulic motor;   is the supply pressure;   is the load pressure; and   is the valve flow gain which varies at different operating points.The following continuity property of the servo valve and motor chamber yields

Table 5 :
Gain/phase margins, phase/gain crossover frequencies, and rise times.

Table 6 :
Gain/phase margins, phase/gain crossover frequencies, and rise times.There are zero-damped, underdamped, and unstable systems used in Examples 1, 2, and 3 individually to get fast responses in large tracking phases.