Abstract
For model-free time-delay systems, an analytical method is proposed to characterize the stabilizing PID region based on the frequency response data. Such characterization uses linear programming which is computationally efficient. The characteristic parameters of the controller are first extracted from the frequency response data. Subsequently, by employing an extended Hermite-Biehler theorem on quasipolynomials, the stabilizing polygon region with respect to the integral and derivative gains and is described for a given proportional gain in term of the frequency response data. Simultaneously, the allowable stabilizing range of is derived such that the complete stabilizing set of the PID controller can be obtained easily. The proposed method avoids the complexity and inaccuracy of the model identification and thus provides a convenient approach for the design and tuning of the PID controller in practice. The advantage of the proposed algorithm lies in that the boundaries of the stabilizing region consist of several simple straight lines, the complete stabilizing set can be obtained efficiently, and it can be implemented automatically in computers.
1. Introduction
Time delay is ubiquitous in many control systems [1]. It may cause the degradation of the control performance and may even render the unstable controlled structure [2, 3]. Hence, the problem of time delay has been studied quite extensively in the literature. The closed-loop characteristic equation of time-delayed systems has an infinite number of roots, which makes the analysis and design extremely challenging.
In industry control, the majority of control systems continue to be operated by PID controller despite the recent advances in control theory and implementation. This popularity stems from their structural simplicity and robust performance in a wide range of operating conditions [4]. Accordingly, the design of PID controllers for the time delay processes is of great importance. Conventional approaches to design a PID controller often involve the development of a mathematical model with time delay, that is, transfer function or state space model, which is obtained by using some approximation or simplification. It inevitably causes the identification error between the mathematical model and the corresponding controlled plant. Thus, the control result may not be reliable. Moreover, some variables or parameters are unavailable or difficult to be obtained in practical applications. Another kind of control technique for the systems with time delay is to design a control system by directly using the original input and output data of the controlled plant which can be easily measured experimentally, that is, model-free control design. Most of the model-free approaches are based on fuzzy logic and neural networks [5, 6]. On the other hand, model-free approaches received special attention in the area of PID autotuning due to many advantages for tuning PID parameters within a closed-loop setting. Many popular tuning rules such as the Ziegler-Nichols [7] procedure and the current practice in PID design for process control are based on the frequency response of the plant. Most of the model-free control methods are the single tuning ones; that is to say, they can only satisfy a specified criterion, but not several required design criteria simultaneously. However, in the practical application, most of the industrial controllers are required to satisfy different performance criteria simultaneously, either in time domain (e.g., overshoot and settling time) or frequency domain (e.g., phase margin and gain margin). In order to reach this purpose, a natural idea is to first present the stabilizing regions of the PID controller and then determine the PID controller by finding the intersection of the areas of the control parameters meeting each required criterion under the constraint of the resultant stabilizing regions. In the recent years, some effective approaches on determining the stabilizing region of the PID controllers have been reported for the system with time delay based on the transfer function model [8–11]. The approaches in [12] and [13] were developed to determine the entire set of stabilizing first-order/PID controller only based on the frequency response of a given LTI plant. In [14], the D-composition approach is employed to determine the stabilizing PID set for the time-delay plant without the detailed model. Although the graphical method based on the D-composition technique is feasible for the determination of the stabilizing PID region, the characterization of the stability boundaries is very complex and cannot be described by using a simple function, which causes difficulty in the practical application. In addition, the procedure for deriving the stabilizing PID region has to be finished manually and cannot be implemented automatically by the computer.
In this paper, a computationally efficient strategy is proposed to generate the linear programming characterization of the stabilizing PID controllers for the time-delay plant based on frequency response data. Based on the extension of the Hermite-Biehler Theorem, the stabilizing polygon region with respect to is derived for a fixed value. Moreover, in order to present the entire stabilizing set of the PID controller efficiently, the allowable stabilizing range of is also derived on the basis of the frequency response data. The results provide the important approach for the PID controller design satisfying different performance criterion simultaneously only based on the frequency response data of the plant.
The paper is organized as follows. The problem statement is given in Section 2. The method to determine the characteristic parameters necessary for the determination of the stabilizing PID controller is addressed in Section 3. In Section 4, the stabilization of the PID controller for the model-free linear time-delay system is analyzed, and the algorithm of determining the set of all stabilizing sets is proposed. Numerical examples are provided to illustrate the main results in Section 5. Finally, the conclusion is given in Section 6.
2. Problem Statement
Consider the feedback configuration with an LTI plant and a PID controller as shown in Figure 1, where is the reference input, is the system output, and is a PID controller with the following form: where is a single-input single output (SISO) linear time-invariant (LTI) plant with the frequency response data, and its detailed model is unknown.
The objective of this paper is to directly determine the complete set of the PID controller that can guarantee the system stability based on the frequency response data of the plant, without constructing a state space or transfer function model.
3. Determination of the Characteristic Parameters
In this section, the characteristic parameters of the plant that are required to determine the stabilizing set of the PID controller will be presented. Assume that the plant has the following transfer function: where is the time delay and are polynomials with real coefficients and high degrees and , respectively. Let , , and denote the numbers of right-half-plane(RHP), left half-plane(LHP), and imaginary-axis zeros of , respectively, and let , , and be the numbers of its right half-plane(RHP), left half-plane(LHP), and imaginary-axis poles of , respectively.
The time-delay can be easily determined by imposing the step signal in the input. The lag time between the step signal and the output response signal is .
The frequency response of only is characterized as where and are the real and imaginary parts of the frequency response of the plant. For the stable plant, the frequency response data can be obtained in terms of multiple point relay test method in [15]. For the unstable plant, the frequency response data cannot be presented directly. However, we can find a controller that can make the closed-loop system stable and then obtain the frequency response data of the plant by measuring the frequency response data of the stable closed-loop system.
The Nyquist and Bode diagrams can be easily sketched based on the frequency response of . From the high frequency slope of the Bode magnitude plot, the relative degree can be presented as follows: where
Then, we compute and . If there exists abrupt increasing change for times at () in the Bode magnitude plot, must have the imaginary-axis poles . To determine the number of multiple imaginary-axis poles, the plot of is drawn for increasing from 1 till the abrupt increasing change at disappears. Assume that when the abrupt increasing change disappears. Thus, we have Conversely, If there exists abrupt decreasing change for times at () in the Bode magnitude plot, must have the imaginary-axis zeros . By the similar lines as the computation of , can also be obtained from the Bode magnitude plot.
Take to be a sufficiently large integer. The net change of the phase angle of as increases from 0 to , denoted as , can be easily obtained from the phase plot. Moreover, when is sufficiently large, the following equation holds (1)When the plant is stable, we have . Then, from (7), can be obtained. (2)When the plant is unstable, the controller that can make the closed-loop system stable, is introduced. The closed-loop transfer function is given by
Transform into the form where and are the numerator and denominator of the transfer function of , respectively.
Take the denominator of (9) as Considering the argument of in , where can be chosen in the interval . Due to the stability of the system, all zeros of are in the left half-plane. From the extension of the extended Hermite theorem [11], the net change in phase of in is , where is the highest order of . Thus, the net change of the phase angle of as increases from to is where and denote the numbers of RHP and LHP of the controller . According to the frequency response data of , can be easily derived. The values of , , , and have been known. Hence, the value of can be obtained from (11).
From the frequency response data of , we have From (12), the argument is obtained. Furthermore, in terms of (3), Hence, the value of can be derived from (13).
In the case ,
4. Stabilization Analysis of the PID Controller
The closed-loop characteristic function of unit feedback system in Figure 1 is
Multiplying (15) by , we have where Multiplying (16) by , we have
Substituting into (18) yields where and are the real and imaginary parts of . The numerator of in (16) is given by Substituting into yields where and represent the real and imaginary parts of . It is seen that and have the same zeros. Substituting into and multiplying it by yields where Actually, in (23) is the numerator of (18).
Lemma 1. If all zeros of lie in the left half-plane, the net change of the phase angle of for increasing from to satisfies the following:
Proof. In terms of the results in [11], it is known that if all zeros of and lie in the left half-plane, and have exactly real zeros in the interval , and the zeros of and interlace. Denote the minimal and maximal zeros of in as and . It can be derived that
As , we have
where is the coefficient of the highest term of . In the case of even, from (31), we have
Thus, we have
In the case of odd, from (31), we have
Thus, (28) also holds for odd.
When , we have
where is the coefficient of the highest term of .
Definition 2. Take to be the real and distinct zeros of in . These roots are also the roots of in . Define as follows.(1)If the plant has the imaginary-axis zero at , then define .(2)If the plant has a zero at the origin, then define , where is given in (20).(3)For all other .
Definition 3. Let . The signature is denoted by
Theorem 4. The necessary and sufficient conditions for the stability of closed-loop system are that there exists such that where , and let be the real and distinct zeros of in . For a fixed , if there exists one string satisfying (38), the stabilizing set of is the intersection of the following inequalities: In addition, for a given LTI time-delay plant with , the stabilizing values of have to satisfy another condition .
Proof. If , are the roots of , then the net change from to can be described as where Take as the roots of in . It is obtained that Since we have It can be obtained that In terms of Lemma 1, if all the zeros of () lie in the left-half plane, then It follows that Substituting (44) into (47) yields From (24), is the odd function, and it has the same roots in and . Furthermore, from (35), has no real zero in when , even and it has one real zero in when odd. Hence, by considering the roots of in the interval . It is easily known that Thus, the expression (48) can be rewritten as where . By taking , the expression (37) can be obtained. This completes the proof.
The following theorem is given to determine the allowable stabilizing range of .
Theorem 5. The necessary condition for the value leading to the existence of the stabilizing region is that has at least distinct roots in , where
By Combining Theorem 4 and Definition 3, the proof of Theorem 5 can be easily proved. Thus, its proof procedure is omitted here.
Based on the results presented in Theorems 4 and 5, the algorithm to determine the stabilizing sets of PID controller based on the frequency response data of the plant is given as follows.
Step 1. Determine the characteristic parameters , , , , , , and .
Step 2. Present the curve of , and with respect to based on the frequency response data of , where .
Step 3. Choose the proper values of and .
Step 4. Compute the allowable stabilizing range of based on Theorem 5.
Step 5. For a grid point of , compute the real and distinct zeros of in (20) in the interval and denote them by in the increasing order.
Step 6. According to Theorem 4 and Definition 3, find satisfying (38).
Step 7. Present the stabilizing set of by computing the intersection of the inequalities determined by (39) with respect to .
Step 8. Go to Step 5 with another grid point of till all the grid points are considered.
5. Numerical Study
Example 6. First, we consider the plant in [11]. Collecting the plant frequency response data, the Nyquist plot and the Bode diagrams of the plant are shown in Figures 2 and 3, respectively.
According to Section 3, the characteristic parameters of the plant are firstly obtained. From the Bode magnitude plot in Figure 3, we have Let , and the net change of the phase is depicted in bode plot as follows: Since the plant is stable, from the following equation: It follows that .
Choose and to find the allowable stabilizing range of . It is seen from Theorem 5 that the following equation must have at least 11 zeros in as follows:
By the plot in Figure 4, it is easy to verify that such a condition holds within the interval . Thus, the allowable stabilizing range of can be obtained.
Choose within . The zeros of are
From Theorem 4, we have Thus, the corresponding stings is
By using (39), it is concluded that the stabilizing corresponding to have to simultaneously satisfy the following: From these inequalities, the stabilizing region of is sketched in Figure 5. The complete stabilizing set of the PID controller can also be derived by sweeping over the allowable stabilizing range of .
Example 7. Let us consider the plant only with the frequency response data, whose Nyquist plot and Bode diagrams are shown in Figures 6 and 7, respectively.
Following the same lines as Example 6, it can be obtained, that and . Let , and let . From Theorem 5, it is derived that the feasible range of is . For within the resultant range, the stabilizing region of is shown in Figure 8. By sweeping over the values in the interval , the complete stabilizing set of the PID controller is presented, which is shown in Figure 9.
6. Conclusion
In this paper, we have proposed an analytical method to characterize the stabilizing PID region for a linear time-delay plant without parametric model using linear programming. Several characteristic parameters, which are required for the determination of the stabilizing PID controller, were firstly provided by using the frequency response information. By employing an extended Hermite-Biehler theorem applicable to quasipolynomials, the 2D parametric region in the space was analytically developed, and its boundaries consist of only several straight lines. Meanwhile, the allowable stabilizing range of is developed. Thus, the 3D visualization of all the stabilizing PID controllers can be conveniently presented only based on the frequency response data of the controlled plant. The results are applicable to arbitrary SISO linear time-delay system, including stable and unstable plants, even the plants with imaginary zeros or poles. The proposed algorithm can be extended to carry out the optimal design of the PID controllers and the design of the PID controllers that can satisfy various performance indices without detailed models.
Acknowledgments
This paper is supported by the National Science Foundation of China (61273116), National Science Funds of Zhejiang Province (Y1111012), and the National Nature Science Funds for Distinguished Young Scholar under Grant (61025016).