- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 582879, 22 pages

http://dx.doi.org/10.1155/2013/582879

## Design of Optimal PID Controller with -Routh Stability for Different Processes

^{1}Department of Information Service and Intelligent Control, Shenyang Institute of Automation Chinese Academy of Sciences, Shenyang 110016, China^{2}University of Chinese Academy of Sciences (Graduate School of Chinese Academy of Sciences), Beijing 100039, China

Received 30 November 2012; Revised 27 February 2013; Accepted 11 March 2013

Academic Editor: Shane Xie

Copyright © 2013 XianHong Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper presents a design method of the optimal proportional-integral-derivative (PID) controller with -Routh stability for different processes through Lyapunov approach. The optimal PID controller could be acquired by minimizing an augmented integral squared error (AISE) performance index which contains control error and at least first-order error derivative, or even may contain *n*th-order error derivative. The optimal control problem could be transformed into a nonlinear constraint optimization (NLCO) problem via Lyapunov theorems. Therefore, optimal PID controller could be obtained by solving NLCO problem through interior method or other optimization methods. The proposed method can be applied for different processes, and optimal PID controllers under various control weight matrices and -Routh stability are presented for different processes. Control weight matrix and -Routh stability’s effects on system performances are studied, and different tuning methods’ system performances are also discussed. -Routh stability’s effects on disturbance rejection ability are investigated, and different tuning methods’ disturbances rejection ability is studied. To further illustrate the proposed method, experimental results of coupled water tank system (CWTS) under different set points are presented. Both simulation results and experiment results show the effectiveness and usefulness of the proposed method.

#### 1. Introduction

PID controllers have been widely used in various industrial processes due to simple structure, few parameters, and easy implementation. In [1–43], many tuning methods such as empirical methods, robust methods, frequency methods, numerical methods, and intelligent methods had been presented to design PID controllers. In [3, 4], Ziegler-Nichols (Z-N) methods include continuous cycling (C-C) method and process reaction curve (PRC) method, and they can be called as first Z-N method and second Z-N method, respectively. The PRC method is especialy proposed for processes with time delay. In [5], the Cohen-Coon method is presented for first-order plus time delay (FOPTD) processes. In [6], the integral absolute error (IAE) performance index is used to design the optimal PID controller for FOPTD and second-order plus time delay (SOPTD) plants. In [7–13], the immune algorithm, modified genetic algorithm, modified ant colony algorithm, particle swarm algorithm (PSA), and modified PSA are used to design optimal PID controller for different systems. In [14], multiple tabu search algorithm (MTSA) is applied to design the optimal fuzzy-logic PID controller for load frequency control system. In [15], some formulas are developed to directly compute the optimal PID controller for certain processes. In [16], zero-pole assignment method is directly applied to compute optimal PID controller for continuous-time processes. In [17], direct search algorithm is used to design the optimal PID controller. In [18], the optimal PID controller is presented for FOPTD plants via dimensional analysis and numerical optimization techniques. In [19], a design method is presented to design the optimal PID controller via minimizing load disturbances and solving a constrained optimization problem in frequency domain. In [20, 21], design methods are proposed to develop optimal PI/PID controller via solving a nonconvex optimization problem in frequency domain. In [22], hybrid Taguchi genetic algorithm (HTGA) is used to design optimal PID controller for pulse width modulation (PWM) feedback systems through solving a constrained optimization problem. In [23–25], linear quadratic regulator (LQR) method is proposed to design the optimal PID controller for different systems and processes. In [26], the loop transfer recovery technique is applied to devise multivariable robust optimal PID controller for a nonminimum phase boiler system. In [27], optimal PID controller with norm constraint is presented for the superheated steam temperature system via PSA. In [28], a parameter space approach is proposed to design the minimum variance PID controller for linear time-invariant (LTI) plants. In [29], a design method is proposed to design PID controller by considering a transient performance item. In [30], a design method of multivariable fixed-structure PID controller is proposed for multivariable plants by solving an optimization problem with linear matrix inequality (LMI) constraint. In [31], iterative feedback tuning (IFT) method is used to design the optimal PID controller for achieving fast response to set point variations. In [32], a two-layer online autotuning algorithm is presented to design the optimal PID controller via general predictive control method. In [33], the LQR method is presented to design the optimal PID controller for robot arms. In [39], a design method of PID controller is proposed for SOPTD processes via internal model control (IMC) framework. In [40], analytic rules of the PID controller are presented for FOPTD and SOPTD processes via the IMC method and model reduction method. In [41], a performance assessment method of PID controller is proposed for integral processes via focusing on step responses of the set point change and load disturbance. In [42], a design method of the PID controller with feedforward compensator is proposed for the FOPTD process. In [43], a design method of the multiloop PI/PID controller is presented for square multiple-input and multiple-output (MIMO) process which is decomposed into a series of equivalent effective open-loop process (EOP) with reasonable gain margins and phase margins. In [44], an on-line relay feedback approach is used to identify, assess, and tune the PID controller. In [45], the pole placement method is used to design the PID controller. In [46], the multiobjective genetic algorithm is presented to design the optimal fractional-order PID controller. In [47], a design method is presented for optimal PID controller with dynamic performance constrained.

Design methods of the PID controller in [1–43] can be roughly divided into the following cases: the empirical methods, intelligent algorithm methods, numerical parametric methods, classic LQR/LQG methods, robust design methods, IMC methods, and pole placement methods. Many design methods of the PID controller are proposed for first-order and second-order processes. In this paper, a novel systematic design method is presented to design optimal PID controller with *ε*-Routh stability for different processes. The optimal PID controller with *ε*-Routh stability is proposed by minimizing augmented integral squared error (AISE) performance index. The proposed method is unlike classic optimal control design methods which need to solve the Riccati equation (RE). The proposed optimal control problem is equivalently transformed into the NLCO problem via Lyapunov theorems. Therefore, optimal controller parameters can be obtained by solving a NLCO problem. The proposed method is used to design optimal PID controller for a third-order, SOPTD, and FOPTD processes, respectively. The proposed method is used to control the liquid level of the CWTS under different set points. The simulation results and experiment results show effectiveness and usefulness of the proposed method. The proposed method is convenient to design the optimal PID controller which can provide high performance for control system.

This paper is arranged as follows: in Section 2, preliminary and optimal control problem statements are described. In Section 3, CWTS model is established, and systematic design procedures of the optimal PID controller are presented in detail. Characteristics and robustness of - stability are analyzed, and the interior method is recommended to solve the NLCO problem of the proposed optimal PID controller. In Section 4, optimal PID controllers and simulation results are presented for different processes. The control weight matrix and *ε*-Routh stability’s effects on system performances are studied. The system performances and disturbance rejection ability of different tuning methods are investigated. -Routh stability’s effects on disturbance rejection ability are discussed as well. To further validate the proposed method, the proposed method is used to control the liquid level of the CWTS under different set points. In Section 5, study contents are reviewed, and some conclusions are made.

#### 2. Preliminary and Optimal Control Problem Statements

*Preliminary*. Based on Lyapunov theorems, a systematic tuning method is proposed to design the optimal PID controller with *ε*-Routh stability. The following theorems and *ε*-Routh stability definition are utilized to design the optimal PID controller.

Theorem 1 (see [1, 2]). *Sufficient and necessary conditions for linear time-invariant system asymptotically stable in a large scope are that: for any given matrix , there is a matrix that satisfies the Lyapunov algebraic equation (LAE): .*

Theorem 2 (see [1, 2]). *The linear time-invariant (LTI) system , and is asymptotically stable, for any given matrix , then the performance index will have equivalent relationships:
*

*Definition of **-Routh stability*. is characteristic polynomial of the linear time-invariant control system. For any nonnegative real , if all elements of the Routh array’s first column of characteristic polynomial can satisfy the following inequality relationships, then the linear time-invariant control system is considered to possess the -Routh stability:

where , and are the elements of the Routh array’s first column. It calls Routh array’s first column as the Routh column, and the general computation formulas of the Routh column can be founded in [48, 49]. It could be known that the *ε*-Routh stability degenerates into the Routh-Hurwitz stable when the is zero . The motivation for the definition of -Routh stability comes from stability analysis of linear systems. As known, if any one element in Routh column of polynomial is zero or negative real, then linear control systems may bring oscillation behaviors or unstable behaviors. Each element in the Routh column of the polynomial is much closer to the zero , at least it is more possible to bring oscillation behaviors. Consequently, the definition of -Routh stability is motivated to describe how far away each Routh column from the zero . It can be known that -Routh stability will imply system robustness in some degree.

*Optimal Control Problem Statements.* We assume that the processes can be described by general transfer function :
where , and are system parameters, and and are the denominator and numerator, respectively. The denominator and numerator’s orders are the and , respectively. For most practical processes, the number of poles is greater than the number of zeros , and and are the positive integers. It could be known that for the all-pole plants , additional for the type I plants , and for the type II plants , and so on. The diagram of the control system is shown in Figure 1, and the control system can be described by the following equation:
where is control error, and are reference command and system output, respectively, is PID controller, and , , and are controller parameters. Equation (4) can be equally expressed as
Therefore, error transfer function can be obtained as
where is the control system’s characteristic polynomial. The error differential equation can be obtained by inversing Laplace transform of (6):
where is the differential operator. When reference command is constant signal or zero signal , or even the piecewise constant signal, the error differential equation will yield
Considering the polynomial differential operator and , error differential equation (8) could be expressed as
Therefore, (9) can be equally expressed as
For the symbol simplicity, the error differential equation (10) is equally expressed as
where is the decision parameter. Parameter is obtained as
Thus, the corresponding state-space model of the differential equation (11) can be obtained:
where is a state vector, is a state matrix, is a parameter vector, is a zero vector, and is an unit matrix. They are obtained as
To obtain the optimal PID controller, it needs to select suitable optimized performance indices. The performance indices such as the integral squared error (ISE), integral absolute error (IAE), integral time absolute error (ITAE), and integral time squared error (ITSE) are often used for criteria when designing optimal PID controller. The precision and steady-state property of control system are directly reflected by control error. In this paper, augmented integral squared error (AISE) in [34, 36] is used to design an optimal PID controller. Hence, optimal PID controller not only stabilizes the control system with -Routh stability, but also minimizes the performance index. Therefore, the proposed optimal control problem can be stated as (i) the optimal PID controller stabilizes the control system with -Routh stability; (ii) the optimal PID controller minimizes the given performance index. Eventually, the proposed optimal control problem can be depicted by the following equation:
where is a control weight matrix, and at least a positive semidefinite real symmetric matrix, and is the linear function of the controller parameters , , and , and , , and are optimal controller parameters. Constraint is the control system’s state constraint, constraint guarantees the control system stable in a large scope, and constraint guarantees the control system with the -Routh stability. However, there are three independent decision variables among decision variables ; the rest decision variables are interrelated variables. Constraint reflects decision variables’ inner relationships. Constraint can assure the existence of feasible solutions and can reduce the searching space. Constraint can also assure PID controller to be physically achieved. In some actual control system, PID controller parameters are limited in certain range. Thus, the controller parameters cannot be too large. The parameters , , and are the scope of controller parameters.

#### 3. CWTS Model and Proposed Design Procedures

The proposed design method is used to design optimal PID controller which is employed to control the liquid level of the CWTS. The experiment setup is shown in Figure 2. It can use the experimental setup to construct different physical processes which are shown in Figure 3.

The physical process in Figure 3 contains the interacting behavior because the liquid level of tank-2 depends on the liquid level of tank-1. In this paper, it does not consider nonlinearities of the coupled water tank system. Thus, the water tank system is considered as the linear process. Based on the volume balance principle, the liquid level equation of tank-1 can be obtained: where is the flow of control valve, is the flow of valve-1, is the flow of valve-3, is the hydraulic capacity of tank-1, and is the liquid level of tank-1. Likewise, the liquid level equation of tank-2 can be obtained: where is the flow of valve-2, is the hydraulic capacity of tank-2, and is the liquid level of tank-2. The water flows equation of valve-1, valve-2, and valve-3 can be obtained: where and are the liquid resistance of valve-1 and valve-2, respectively, and is the liquid resistance of valve-3. Thus, the liquid level model of water tank-1 is acquired: where is the liquid level transfer function of tank-1. The transfer function shows that physical process is the second-order process with a negative zero. System parameters are obtained as

##### 3.1. Proposed Design Procedures

Without losing generalities, a third-order controlled process is considered in the unity feedback system. The process and PID controller have forms, respectively,
Thus, open-looped transfer function of control system yields
Then, the error equation and of the control system are obtained:
where is the control system’s characteristic polynomial. The error differential equation can be obtained by inversing Laplace transform of the error equation:
where * δ* is a differential operator. When reference command is the constant signal, or zero signal , or even piecewise constant signal, then the error differential equation will yield
For symbol simplicity, error differential equation is equally expressed as
where is the decision parameter. The decision parameter yields
Thus, the corresponding state-space model of error differential equation (27) can be obtained:
where is state vector, is state matrix, is the middle decision variables, and the . Therefore, the optimal control problem is formulated as
If optimal control problem (30) is solved, then the proposed optimal PID controller could be obtained. The constraint assures that the control systems asymptotically stable in a large scope. The determinant is . Based on LTI stablility theory, constraint is equivalently transformed into the following inequality:
The constraint assures that the control system possesses the -Routh stability. The characteristic polynomial is . According to -Routh stability definition, the constraint is formulated by the following inequalities:
The inequalities (31) and (32) assure control system asymptotic stability and -Routh stability, respectively. Therefore, the control system is asymptotically stable in a large scope; based on the Lyapunov Theorems 1 and 2, performance index is equivalently transformed into the performance index
where is a positive definite real symmetric matrix, and matrix meets the LAE: . The performance index is determined by the matrix and initial state. Therefore, the proposed optimal control problem (30) can be equivalently transformed into the following NLCO problem:
The proposed optimal control problem (30) is equivalently transformed into a NLCO problem via applying the Lyapunov theorems. The optimal PID controller with -Routh stability can be obtained by solving NLCO problem (34). The design problem of optimal PID controller with -Routh stability is deduced into the issue that: for the given control weight matrix and initial states, we pursue a suitable matrix to minimize performance index

*J*. We assume that the control system is static at the beginning, then a unit step command is inputted into the system, and the initial states can be obtained: , and . The control weight matrix is a diagonal matrix with . Apparently, NLCO problem (34) has five middle decision variables. In fact, it only has three independent variables. Different optimization methods such as the Newton method, quasi-Newton method, Lagrange method, conjugate gradient method, interior point method, linear programming method, particle swarm algorithm, genetic algorithm, evolution algorithm, and other intelligent algorithms in [50–63] have been well established to solve the NLCO problem. Therefore, optimal parameters and proposed PID controller are gotten:

It is noted that (i) the proposed method is also suitable for first-order, second-order, and high-order processes; (ii) the proposed method is also suitable for nonunit feedback system. The non-unit feedback system can be transformed into unit feedback system through transfer function equivalent principle; (iii) the proposed method is also suitable for processes with time delay which can be approximated by transfer function; (iv) the proposed method can be used for different processes.

##### 3.2. Characteristic Analysis of -Routh Stability

The -Routh stability constraint will imply control system robustness in some degree. The -Routh stability constraint could avoid too conservative system performances and may help to improve the robustness of control systems. The closed-loop characteristic polynomial is . The polynomial coefficient is a linear mapping function of the PID controller parameters , , and . The parameters , , and are constrained in intervals , , and respectively. Therefore, the PID controller parameters can be seen as the uncertain bounded parameters in some degree. In this sense, polynomial can be seen as interval polynomial with three uncertain parameters , , and . To analyze the robustness of -Routh stability, without losing generality, interval polynomial has the general form:
where and are the uncertain coefficients of the interval polynomial . The uncertain coefficients , and are continuous mapping functions of uncertain bounded parameters , and :
where is the uncertain parameter space, and are the lower bound and upper bound of uncertain parameter , and is continuous mapping function of polynomial for uncertain parameters. Based on -Routh stability definition for interval polynomial , the following inequalities can be obtained:
where is the element of first column of Routh array. According to the interval polynomial stable criterion in [36, 64], if the disjoint set of feasible set and uncertain parameters space is not a empty set, then the interval polynomial will satisfy the interval polynomial stable criterion:
Therefore, the interval polynomial is Routh-Hurwitz stable in the feasible set , and the **-**Routh stability could guarantee the robustness of linear systems in a certain degree. For any positive real , two inequalities for -Routh stability are obtained:
The inequality (I) and inequality (II) have the feasible spaces and , and the first inequality of inequalities (I) and (II) has the feasible set
Then, any point in the feasible set will have the relationship
Likewise, any point in the feasible set will have the relationship
For the feasible set , there at least a point which will satisfy the relationship
Therefore, the feasible sets and will have the following relationship:
Then, the overall feasible spaces and will have following relationships:
The overall feasible set will be gradually condensed with the *ε* increasing. With real *ε* sequence , one will get following relationships:
Therefore, the feasible set will has the following relationship:
When the feasible set is strictly condensed with , then the limit of feasible space will yield the following relationships:
If the tends to infinity, then the limit of the feasible set will be condensed into a point or an empty set . If the feasible set is condensed into a point , then the interval polynomial degenerates into a usual polynomial. It reflects that the point is a fixed point in the stable region of linear control systems. If the feasible set tends to empty set , therefore it has a real to satisfy the following relationships:
when the feasible set has the relationship . However, this is a special event which can be depicted by the probability event
It assumes that all the probability events are independent and identically distributed. Hence, the probability can be written as
Considering probability event first, the probability will have
The probability for the probability event will yield
Therefore, the probability of the probability event yields

##### 3.3. Interior Method for NLCO Problem

Without losing generality, NLCO problem (34) of the optimal PID controller can be written in the following general forms:
where , and are the smooth and differentiable functions, is the decision variable, and , , and denote the number of the decision variable, equality constraint, and inequality constraint, respectively. The various optimization methods in [50–63] are developed to solve optimization problems. The optimization methods such as the Newton methods, conjugate gradient methods, steepest descent methods, interior point methods, genetic algorithms, and particle swarm algorithms in [50–63], are well established to solve constraint optimization problems. Based on interior point methods in [56–61], the interior method is recommended to solve the NLCO problem. Then, the NLCO problem (56) could be transformed into the following form:
where is the barrier parameter, the slack vector is set to be positive, and is an expanded constraint inequality. It introduces the Lagrange multipliers and for the barrier problem (57):
where is Lagrange function, and are Lagrange multipliers for constraints and , respectively. Based on the Karush-Kuhn-Tucker (KKT) optimality conditions [57–59], the optimality conditions for NLCO problem (56) can be expressed as
where is the diagonal matrix and its elements are the components of the vector **,** and is a vector of all ones, and are the Jacobian matrices of the vectors and , respectively, is the grand of function , and is a diagonal matrix and its elements are the components of vector . The system (59) is the KKT condition of the NLCO problem (56). When in the search approach, it should remain , . To obtain the iteration direction, it can make the point satisfying the KKT conditions(59); then the following system will be obtained:
System (60) is obtained by ignoring higher-order incremental and replaces the nonlinear terms with linear approximation in system (59). The system (60) is written in the matrix form:
where is the Hessian matrix in system. Finally, the new iterate direction is obtained via solving the system (61), which is the essential process in the interior point method. Thus, the new iteration point can be gotten in the next iteration:
where is the step size. Choosing the step size holds the , in search process. In this paper, the interior point method is recommended to solve the NLCO problem. The interior algorithm framework is depicted as follows.

*Interior Algorithm Framework for NLCO Problem. One has the following.*

*Step 1. *Choose an initial iteration point in the feasible region set , and the .

*Step 2. *Construct current iterate; we have the current iterate values , , , and of the primal variable , the of the slack variable , the multipliers and , respectively.

*Step 3. *Calculate the Hessian matrix of the Lagrange system , and the Jacobian matrices and are of the vectors and in the current iterate .

*Step 4. *Solve the linear system (61) and construct the iterate direction . Solving the linear matrix equation (61), we obtain the primal solution , multipliers solution , , and also the slack variable solution .

*Step 5. *Choosing the step size holds the , in search process, . Update the iterate values: .

*Step 6. *Check the ending conditions in region . If it does not satisfied, go to Step 3, else minimum of NLCO problem is obtained in feasible region .

*Step 7. *End.

#### 4. Simulation and Experiment Results

The proposed method is used to design the optimal PID controller with -Routh stability for processes [7], [15], and [6], respectively. is a third-order all-pole process, and and are the SOPTD and FOPTD processes, respectively,

In order to use the proposed method, time delay of plant.2 and plant.3 is approximated as and , respectively. Thus, transfer functions of and are the and , respectively. The proposed method is used to design the optimal PID controller with -Routh stability for different processes, and optimal PID controllers under different control weight matrices and -Routh stability are shown in Table 1.

For the proposed method, the plant.1’s step response and system performances under various control weight matrices are shown in Figure 4 and Table 2, respectively. The control weight matrix has effects on system performances. Figure 4 and Table 2 obviously show that system performances are affected by control weight factors , , , and . The results show that the response velocity will be heightened by increasing weight factor or reducing weight factors , , and . When control weight factor (with , , and fixed) is increased, the overshoot is increased, and the rise time, delay time, peak time, and setting time are reduced. Likewise, the overshoot will be increased, and the rise time, delay time, and peak time will be reduced by reducing weight factors , , and (with fixed). Response results show that system performances are affected by control weight matrix.

For th-order processes, it can be inferred that overshoot will be heightened, and the rise time, delay time, and peak time will be reduced by increasing control weight factor (with other factors fixed) or reducing other weight factors (with fixed). On the contrary, the overshoot will be reduced, and the rise time, delay time, and peak time will be increased by reducing weight factor (with other factors fixed) or increasing other factors (with fixed). Therefore, system performances are affected by control weight factors . In all, control weight matrix has effects on system performances.

For plant.2, step response results and system performances under different -Routh stability are shown in Figure 5 and Table 3, respectively. The response results show that system performances are affected by -Routh stability (with control weight matrix fixed). The results show that response velocity will be heightened by increasing . When the (with control weight matrix fixed) is increased, the overshoot is increased, and the rise time, delay time, and peak time are reduced. However, setting time is gradually increased and then is decreased by increasing . Therefore, the rise time, delay time, and peak time will be increased, and the overshoots will be reduced by reducing the . In all, system performances are affected by the -Routh stability too.

##### 4.1. System Performances of Different Tuning Methods

The step response and system performances of different tuning methods for plant.1, plant.2, and plant.3 are shown in Figures 6, 7, and 8 and Table 4, respectively. For plant.1, the results obviously show that the Z-N method, IA-ISE method [7], and IA-ITSE method [7] show oscillatory behaviors, but the proposed method shows smoothly response without any oscillation. The Z-N method has small rise time, delay time, and peak time compared with the proposed method, IA-ISE method, and IA-ITSE method. However, the Z-N method has the largest overshoot and setting time among these tuning methods. The overshoot and setting time of the Z-N method is about 11.08% and 6.41 sec, and the overshoot and setting time of the proposed method, IA-ISE method, and IA-ITSE method are about 4.08% and 3.54 sec, 4.24% and 3.13 sec, and 5.27% and 4.88 sec, respectively. Thus, these tuning methods own advantages for plant.1.

The plant.2 is a SOPTD plant with dead time ratio . The results in Figure 7 show that dynamic performances of the Z-N method, A-H method [19], gain and phase method [15], and proposed method have different features. The Z-N method, A-H method [19], gain and phase method [15], and the proposed method have small rise time, delay time and peak time, setting time, and overshoot, respectively. The proposed method has the largest rise time, delay time, and peak time, and the A-H method and Z-N method have the largest setting time and overshoot, respectively. The longest rise time, delay time, peak time, and setting time are about 4.778 sec, 1.841 sec, 6.622 sec, and 4.43 sec, respectively, and the largest overshoot is about 43.09%. The gap between the largest and the smallest items of the rise time, delay time, peak time, and setting time is small, but only the gap between the largest and the smallest overshoot is about 13 times of the smallest ones. Thus, the comprehensive performances of Z-N method are poorer than rest methods.

The plant.3 is a FOPTD plant with large dead time ratio . The results in Figure 8 show that all the tuning methods show oscillation response. The overshoot of the Z-N method (C-C) and Z-N method (PRC) are about 18.18% and 21.17%, respectively. The results demonstrate that the proposed method has the smallest peak time 20.87 sec and the smallest overshoot 6.67%, and the Cohen-Coon method has the smallest rise time 17.252 sec and the smallest delay time 14.306 sec, and the Z-N method (C-C) has about the smallest setting time 51.60 sec. However, the Z-N method (PRC) has the largest rise time 20.286 sec, delay time 14.635 sec, and peak time 26.43 sec, and the Cohen-Coon method has the largest setting time 67.05 sec and largest overshoot 38.06%. For plant.3, the Cohen-Coon method shows poorer system performances than other tuning methods. By studying different tuning methods’ dynamic performances, it can be known that the proposed method could provide good system performances.

##### 4.2. Disturbance Rejection Ability of Different Methods

Disturbance response results for plant.1, plant.2, and plant.3 are shown in Figures 9, 10, and 11, respectively. For plant.1, the results show that the proposed method, Z-N method, IA-ISE and IA-ITSE methods have good ability to reject the disturbances, and they show roughly the same disturbance response amplitude (DRA) and disturbance attenuation time (DAT). The proposed method has the best disturbance rejection ability among the methods. The proposed method does not show oscillatory behaviors, but Z-N method, IA-ISE method, and IA-ITSE method show oscillatory behaviors. For plant.2, the disturbance response results of the proposed method, Z-N method, A-H method, and gain and phase method are shown in Figure 10. The results show that these methods have good ability to reject disturbances and they have roughly the same DAT. The DRAs of the proposed method, Z-N method, A-H method, and gain and phase method are about 0.13, 0.10, 0.10, and 0.11, respectively. The proposed method has larger DRA than the rest methods, but the proposed method shows smooth disturbance response. For plant.3, disturbance response results in Figure 11 show that these tuning methods have roughly the same ability to reject disturbances. The proposed method, Z-N method (C-C), Z-N method (PRC), and Cohen-Coon method have roughly the same DRA and DAT. Among these tuning methods, the DRA of the Z-N method (PRC) is about 0.275, the proposed method, Z-N method (C-C), and Cohen-Coon method have the same DRA which is about 0.25. Overall, as known from disturbance response results, it could be that the proposed method, A-H method, and gain and phase method roughly have roughly the same ability to reject disturbances.

To study -Routh stability’s effects on disturbance rejection ability, the disturbance response results under various are shown in Figure 12. The disturbance response results in Figure 12 clearly show that the *ε* has effects on the DRA and DAT. The DAT will be continuously reduced when the increases from 0.5 to 5.0. However, the DAR will be gradually reduced and then be increased when the increases from 0.5 to 5.0. The DAT reflects how fast the disturbances will be eliminated, and DRA reveals how much disturbance amplitude of the control system. Thus, -Routh stability has effects on the disturbance rejection ability.

##### 4.3. Experiment Results for the CWTS

The proposed PID controller is utilized to control the liquid level of a CWTS with interactive behavior. The proposed PID controllers’ experimental results of the CWTS under various control weight matrix and different set points are shown in Figures 14, 15, 16, 17, 18, 19, 20, and 21. The theoretical model of the CWTS is already obtained in previous section; therefore the identification model can be acquired via open-loop experiment. The liquid level’s real-time data is stored in the disk of monitoring computer, and the date can be accessed on-line for process control (OPC) database. The identification model of CWTS is identified as via least square method (sample time is 1 second) when open ratio of control valve, output valve-1, output valve-2, and connected valve-3 is set at 75%, 100%, 30%, and 100%, respectively. The steady-state liquid level is about 294.00 mm, and it will take about 12 minutes to reach this liquid level. The actual liquid level and the identification model’s liquid level are shown in Figure 13. The identification model is consistent with the actual model of the CWTS; thus, the identification model describes the actual model accurately.

In the closed-loop experiment, the proposed method is used to control the liquid level of the CWTS under different set points. Proposed PID controllers under different control weight matrices and -Routh stability are used in the experiments. The proposed PID controllers are used to control liquid level of the CWTS with interactive behaviors. The first group experiments’ setting point of the liquid level is set at 300.00 mm and second group experiments’ setting point is set at 160.00 mm. The corresponding experiment results of the first group and second group experiments are shown in Figures 14–21, respectively. In Figures 14–21, the solid red line is the dynamic response of the liquid level , the solid magenta line is the dynamic response of the liquid level , and the solid cyan line is the control signal. In the first group experiments, the dynamic response results of the proposed method shown in Figures 14–17 obviously reveal that the proposed method provides good system performances for the CWTS with interactive behaviors. In first group experiments, the setting time for the CWTS is less than 7 minutes, and the overshoot of the proposed method is less than 4% under various control weight matrices of the proposed method.

In the second group experiments, the dynamic response results of the proposed method shown in Figures 18–21 clearly reveal that the proposed method provides good system performances and robustness. The proposed PID controller is designed for the liquid level of set point at 294 mm. The actual model of the CWTS has been changed in the second group experiments because the liquid level set point has made a huge change. In the second group experiments, the setting time of the proposed method for the CWTS is less than 3.5 minutes, and the overshoot of the proposed method is less than 12%. Figures 18 and 20 show that the liquid level has a small fluctuation; however, overall system performances are still accepted. The system parameters of liquid model will show perturbations when the set point of liquid level makes a huge change. The proposed method can still provide good system performances when the set point of the liquid level changes from 300 mm to 160 mm. Furthermore, the dynamic response results in Figures 18–21 also reflect that the proposed method shows good robustness to rejection model perturbation, system parameters perturbation, and set point variation.

Simulation results and experiment results reveal that the proposed method can provide good system performances for different processes. The proposed method also demonstrates good disturbance rejection ability and robustness for different processes. Therefore, the proposed method is useful to design the optimal PID controller. In short, the proposed method is suitable for different processes and can provide good system performances. The simulation results and experiment results show the effectiveness and usefulness of the proposed method.

#### 5. Conclusions

In this paper, a systematic method is presented to design the optimal PID controller with -Routh stability for different processes via Lyapunov approach. The optimal PID controller is obtained by minimizing AISE performance index which contains the control error and at least first-order error derivative, or even may contain the th*-*order error derivative. The proposed optimal control problem could be equivalently transformed into a NLCO problem via Lyapunov theorems. Therefore, the optimal PID controller can be obtained by solving the NLCO problem via the interior method or other optimization methods. Optimal PID controllers under different control weight matrices and -Routh stability are presented for different processes. For the proposed method, control weight matrix and -Routh stability’s effects on system performances are analyzed. Different tuning methods’ system performances are discussed.

The control weight matrix and -Routh stability have large effects on system performances. For third-order processes, system performances and response velocity are affected by the control weight factors , , , and . The peak overshoot and response velocity will be heightened by increasing weight factor (with other factors fixed). Likewise, the overshoot and response velocity will be heightened by reducing weight factors , , and (with fixed). On the contrary, the overshoot and response velocity will be reduced by reducing weight factor (with other factors fixed) or increasing weight factors , , and (with fixed). Through the study, for th-order processes, the above conclusions still hold; it can be inferred that system performances are affected by the control weight factors . The -Routh stability (with control weight matrix fixed) has effects on system performances. The results show that the overshoots and response velocity will be heightened by increasing the . Overshoots will be increased, and the rise time, delay time, and peak time will be reduced by increasing the . On the contrary, overshoots will be reduced, and the rise time, delay time, and peak time will be increased by reducing the .

Disturbance rejection ability of different tuning methods is also investigated, and simulation results reveal that the proposed method shows good disturbance rejection ability. The -Routh stability’s effects on disturbance rejection ability are investigated. The disturbance response results clearly show that the has effects on the DRA and DAT. The DAT will be continuously reduced by increasing the . However, the DAR will be gradually reduced and then be increased by increasing the .

To further validate the proposed method, the proposed method is utilized to control the liquid level of CWTS with interactive behaviors. The experimental results under different set points and control weight matrices are presented in detail. Both simulation results and experimental results show the effectiveness and usefulness of the proposed method.

#### Acknowledgments

A partial content of this paper, as the previous work and earlier version paper, was presented at the 2009 International Conference on Information Engineering and Computer Science (ICIECS-2009), Wuhan, China December 19-20, 2009. This work is supported by the Research Projects no. 2007BAF09B01 and 2009A A04Z155. The authors would like to thank the reviewers for suggestions.

#### References

- J. J. D’Azzo and C. H. Houpis,
*Linear Control System Analysis and Design*, McGraw-Hill, New York, NY, USA, 2nd edition, 1981. - M. D. Tong,
*Linear System Theory and Design*, University of Science and Technology of China press, Hefei, China, (Chinese), 2004. - J. G. Ziegler and N. B. Nichols, “Optimum settings for automatic controllers,”
*ASME Transactions*, vol. 64, pp. 759–768, 1942. - J. G. Ziegler and N. B. Nichols, “Process lags in automatic control circuits,”
*Transactions of the ASME*, vol. 65, pp. 433–444, 1943. - G. H. Cohen and G. A. Coon, “Theoretical consideration of related control,”
*Transactions of the ASME*, vol. 75, pp. 827–834, 1953. - C. R. Madhuranthakam, A. Elkamel, and H. Budman, “Optimal tuning of PID controllers for FOPTD, SOPTD and SOPTD with lead processes,”
*Chemical Engineering and Processing*, vol. 47, no. 2, pp. 251–264, 2008. View at Publisher · View at Google Scholar · View at Scopus - D. H. Kim and J. H. Cho, “Intelligent tuning of PID controller with disturbance function using immune algorithm,” in
*Proceedings of the IEEE Annual Meeting Fuzzy Information Processing*, vol. 1, pp. 286–291, June 2004. View at Scopus - Y. P. Wang, D. R. Hur, H. H. Chung, N. R. Watson, J. Arrillaga, and S. S. Matair, “Design of an optimal PID controller in AC-DC power system using modified genetic algorithm,” in
*Proceedings of the 2000 International Conference on Power System Technology*, vol. 3, pp. 1437–1442, December 2000. - Q. Zeng and G. Tan, “Optimal design of PID controller using modified ant colony system algorithm,” in
*Proceedings of the 3rd International Conference on Natural Computation (ICNC '07)*, vol. 5, pp. 436–440, August 2007. View at Publisher · View at Google Scholar · View at Scopus - Z. L. Gaing, “A particle swarm optimization approach for optimum design of PID controller in AVR system,”
*IEEE Transactions on Energy Conversion*, vol. 19, no. 2, pp. 384–391, 2004. View at Publisher · View at Google Scholar · View at Scopus - H. Hu, Q. Hu, Z. Lu, and D. Xu, “Optimal PID controller design in PMSM servo system via particle swarm optimization,” in
*Proceedings of the 31st Annual Conference of IEEE Industrial Electronics Society (IECON '05)*, pp. 79–83, November 2005. View at Publisher · View at Google Scholar · View at Scopus - T. H. Kim, I. Maruta, and T. Sugie, “Particle swarm optimization based robust PID controller tuning scheme,” in
*Proceedings of the 46th IEEE Conference on Decision and Control (CDC '07)*, pp. 200–205, New Orleans, La, USA, December 2007. View at Publisher · View at Google Scholar · View at Scopus - M. I. Solihin, Wahyudi, M. A. S. Kamal, and A. Legowo, “Optimal PID controller tuning of automatic gantry crane using PSO algorithm,” in
*Proceedings of the 5th International Symposium on Mechatronics and its Applications (ISMA '08)*, Amman, Jordan, May 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Pothiya and I. Ngamroo, “Optimal fuzzy logic-based PID controller for load-frequency control including superconducting magnetic energy storage units,”
*Energy Conversion and Management*, vol. 49, no. 10, pp. 2833–2838, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. Zhuang and D. P. Atherton, “Automatic tuning of optimum PID controllers,”
*IEE Proceedings D*, vol. 140, no. 3, pp. 216–224, 1993. View at Scopus - G. Herjdlfsson and A. S. Hauksddttir, “Direct computation of optimal PID controllers,” in
*Proceedings of the 42nd IEEE Conference on Decision and Control*, pp. 1120–1125, Maui, Hawaii USA, December 2003. - S. Daley and G. P. Liu, “Optimal PID tuning using direct search algorithms,”
*Computing and Control Engineering Journal*, vol. 10, no. 2, pp. 51–56, 1999. View at Scopus - S. Tavakoli and M. Tavakoli, “Optimal tuning of PID controllers for first order plus time delay models using dimensional analysis,” in
*Proceedings of the 4th International Conference on Control and Automation*, pp. 942–946, Montreal, Canada, June 2003. View at Scopus - H. Panagopoulos, K. J. Astrom, and T. Hagglund, “Design of PID controllers based on constrained optimization,”
*IEE Proceedings*, vol. 149, no. 1, pp. 32–40, 2002. - C. Hwang and C. Y. Hsiao, “Solution of a non-convex optimization arising in PI/PID control design,”
*Automatica*, vol. 38, no. 11, pp. 1895–1904, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. L. Xue, D. H. Li, and C. D. Lv, “Optimal design of PID controller based on sensitivity constraint,” in
*Proceedings of the 5th World Congress on Intelligent Control and Automation*, Hangzhou, China, June 2004. - C. H. Hsieh and J. H. Chou, “Design of optimal PID controllers for PWM feedback systems with bilinear plants,”
*IEEE Transactions on Control Systems Technology*, vol. 15, no. 6, pp. 1075–1079, 2007. View at Publisher · View at Google Scholar · View at Scopus - R. T. J. O’Brien and J. M. Howe, “Optimal PID controller design using standard optimal control techniques,” in
*Proceedings of the 2008 American Control Conference*, pp. 4733–4738, Seattle, Wash, USA, June 2008. - J. B. He, Q. G. Wang, and T. H. Lee, “PI/PID controller tuning via LQR approach,” in
*Proceedings of the 37th IEEE Conference on Decision and Control*, pp. 1177–1182, Tampa, Fla, USA, December 1998. - G. R. Yu and R. C. Hwang, “Optimal PID speed control of brush less DC motors using LQR approach,” in
*Proceedings of the IEEE International Conference on Systems, Man and Cybernetics (SMC '04)*, pp. 473–478, October 2004. View at Scopus - H. Zargarzadeh, M. R. Jahed Motlagh, and M. M. Arefi, “Multivariable robust optimal PID controller design for a non-minimum phase boiler system using loop transfer recovery technique,” in
*Proceedings of the 6th Mediterranean Conference on Control and Automation (MED '08)*, pp. 1520–1525, Ajaccio, France, June 2008. View at Publisher · View at Google Scholar · View at Scopus - P. Han, Y. Huang, Z. Z. Jia, D. F. Wang, and Y. L. Li, “Mixed
*H*_{2}/*H*_{∞}optimal pid control for superheated steam temperature system based on PSO optimization,” in*Proceedings of the International Conference on Machine Learning and Cybernetics (ICMLC '05)*, pp. 960–964, Guangzhou, China, August 2005. View at Scopus - P. B. Dickinson and A. T. Shenton, “A parameter space approach to constrained variance PID controller design,”
*Automatica*, vol. 45, no. 3, pp. 830–835, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. C. Basilio and S. R. Matos, “Design of PI and PID controllers with transient performance specification,”
*IEEE Transactions on Education*, vol. 45, no. 4, pp. 364–370, 2002. View at Publisher · View at Google Scholar · View at Scopus - M. Saeki, “Fixed structure PID controller design for standard
*H*_{∞}control problem,”*Automatica*, vol. 42, no. 1, pp. 93–100, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. Lequin, M. Gevers, M. Mossberg, E. Bosmans, and L. Triest, “Iterative feedback tuning of PID parameters: comparison with classical tuning rules,”
*Control Engineering Practice*, vol. 11, no. 9, pp. 1023–1033, 2003. View at Publisher · View at Google Scholar · View at Scopus - M. Xu, S. Li, C. Qi, and W. Cai, “Auto-tuning of PID controller parameters with supervised receding horizon optimization,”
*ISA Transactions*, vol. 44, no. 4, pp. 491–500, 2005. View at Scopus - G. L. Luo and G. N. Saridis, “LQ design of PID controllers for robot arms,”
*IEEE Journal of Robotics and Automation*, vol. RA-1, no. 3, pp. 152–159, 1985. View at Scopus - X. H. Li, H. B. Yu, and M. Z. Yuan, “Design of an optimal PID controller based on Lyapunov approach,” in
*Proceedings of the International Conference on Information Engineering and Computer Science (ICIECS '09)*, Wuhan, China, December 2009. View at Publisher · View at Google Scholar · View at Scopus - K. J. Astrom and T. Hagglund, “The future of PID control,”
*Control Engineering Practice*, vol. 9, pp. 1163–1175, 2001. - X. H. Li, H. B. Yu, M. Z. Yuan, and J. Wang, “Design of robust optimal proportional-integral-derivative controller based on new interval polynomial stability criterion and Lyapunov theorem in the multiple parameters' perturbations circumstance,”
*IET Control Theory and Applications*, vol. 4, no. 11, pp. 2427–2440, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - H. Biao and R. Kadali, “System identification: conventional approach,” in
*Dynamic Modeling, Predictive Control and Performance Monitoring*, vol. 374 of*Lecture Notes in Control and Information Sciences*, pp. 9–29, Springer Press, Berlin, Germany, 2008. View at Publisher · View at Google Scholar - F. P. Ti,
*System Identification*, Zhe Jiang University Press, Hangzhou, China, (Chinese), 2004. - D. E. Rivera, M. Morari, and S. Skogestad, “Internal model control. 4. PID controller design,”
*Industrial and Engineering Chemistry Process Design and Development*, vol. 25, no. 1, pp. 252–265, 1986. View at Scopus - S. Skogestad, “Simple analytic rules for model reduction and PID controller tuning,”
*Journal of Process Control*, vol. 13, no. 4, pp. 291–309, 2003. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - M. Veronesi and A. Visioli, “Performance assessment and retuning of PID controllers for integral processes,”
*Journal of Process Control*, vol. 20, no. 3, pp. 261–269, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. Visioli, “A new design for a PID plus feedforward controller,”
*Journal of Process Control*, vol. 14, no. 4, pp. 457–463, 2004. View at Publisher · View at Google Scholar · View at Scopus - H. P. Huang, J. C. Jeng, C. H. Chiang, and W. Pan, “A direct method for multi-loop PI/PID controller design,”
*Journal of Process Control*, vol. 13, no. 8, pp. 769–786, 2003. View at Publisher · View at Google Scholar · View at Scopus - K. K. Tan, T. H. Lee, and X. Jiang, “On-line relay identification, assessment and tuning of PID controller,”
*Journal of Process Control*, vol. 11, no. 5, pp. 483–496, 2001. View at Publisher · View at Google Scholar · View at Scopus - Q. G. Wang, Z. Zhang, K. J. Astrom, and L. S. Chek, “Guaranteed dominant pole placement with PID controllers,”
*Journal of Process Control*, vol. 19, no. 2, pp. 349–352, 2009. View at Publisher · View at Google Scholar · View at Scopus - L. Meng and D. Xue, “Design of an optimal fractional-order PID controller using multi-objective GA optimization,” in
*Proceedings of the Chinese Control and Decision Conference (CCDC '09)*, pp. 3849–3853, June 2009. View at Publisher · View at Google Scholar · View at Scopus - X. H. Li, H. B. Yu, and M. Z. Yuan, “A new design method of optimal PID controller with dynamic performances constrained,”
*Control and Intelligent Systems*, vol. 40, no. 4, pp. 1–21, 2012. - I. J. Nagrath and M. Gopal,
*Control Systems Engineering*, Wiley Eastern, New Delhi, India, 1975. - G. H. Hostetter, C. J. Savant Jr., and R. T. Stefani,
*Design of Feedback Control System*, CBS college, 1982. - T. F. Coleman and Y. Li, “On the convergence of interior-reflective Newton methods for nonlinear minimization subject to bounds,”
*Mathematical Programming*, vol. 67, no. 2, pp. 189–224, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. F. Shanno, “Conditioning of quasi-Newton methods for function minimization,”
*Mathematics of Computation*, vol. 24, pp. 647–656, 1970. View at Publisher · View at Google Scholar · View at MathSciNet - W. W. Hager and H. Zhang, “A new conjugate gradient method with guaranteed descent and an efficient line search,”
*SIAM Journal on Optimization*, vol. 16, no. 1, pp. 170–192, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. E. Dennis, M. El-Alem, and K. Williamson, “A trust-region approach to nonlinear systems of equalities and inequalities,”
*SIAM Journal on Optimization*, vol. 9, no. 2, pp. 291–315, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. W. Hager and D. T. Phan, “An ellipsoidal branch and boun d algorithm for global optimizati on,”
*SIAM Journal on Optimization*, vol. 20, no. 2, pp. 740–758, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. I. M. Gould, “The generalized steepest-edge for linear programming,”
*Combinatorics and Optimization, Technical Report*83-2, University of Waterloo, Waterloo, Canada, 1983. - Y. Zhang, “Solving large-scale linear programs by interior-point methods under the MATLAB environment,” Tech. Rep. TR96-01, Department of Mathematics and Statistics, University of Maryland, Baltimore, Md, USA, 1995.
- R. J. Vanderbei and D. F. Shanno, “An interior-point algorithm for non-convex nonlinear programming,”
*Computational Optimization and Applications*, vol. 13, no. 1–3, pp. 231–252, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Y. Benson and D. F. Shanno, “Interior-point methods for nonconvex nonlinear programming: regularization and warmstarts,”
*Computational Optimization and Applications*, vol. 40, no. 2, pp. 143–189, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. A. Waltz, J. L. Morales, J. Nocedal, and D. Orban, “An interior algorithm for nonlinear optimization that combines line search and trust region steps,”
*Mathematical Programming*, vol. 107, no. 3, pp. 391–408, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. B. Richard, M. E. Hribar, and J. Nocedal, “An interior point algorithm for large-scale nonlinear programming,”
*SIAM Journal on Optimization*, vol. 9, no. 4, pp. 877–900, 1997. View at Zentralblatt MATH · View at MathSciNet - A. Forsgren, P. E. Gill, and M. H. Wright, “Interior methods for nonlinear optimization,”
*SIAM Review*, vol. 44, no. 4, pp. 525–597, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. E. Goldberg,
*Genetic Algorithms in Search, Optimization and Machine Learning*, Addison-Wesley, Reading, Mass, USA, 1989. - H. H. Holger and T. Stutzle,
*Stochastic Local Search: Foundations and Applications*, Morgan Kaufmann, San Francisco, Calif, USA, 2004. - X. H. Li, H. B. Yu, and M. Z. Yuan, “Robust stability of interval polynomials and matrices for linear systems,” in
*Proceedings of the IASTED International Conference Modelling, Identification, and Control*, pp. 138–147, Phuket, Thailand, November 2010.