Abstract

This paper proposes a new method to design the observer for state delay systems such that (i) state estimation errors converge to zero quickly and (ii), at the same time, a quadratic performance measurement of the deviation of estimates from the actual states is minimized for reducing large error during the transient period of observation. The proposed new approach fuses the merits of both the orthogonal functions approach and evolutionary optimization. One illustrative example is given to verify the effectiveness and efficiency of the proposed new optimization method on performance improvement of state estimations. From the illustrative example, in addition to the asymptotical convergence of the estimated state errors, the performance index for the proposed optimal design approach is clearly much lower than that of the nonoptimal design method.

1. Introduction

Control systems are commonly used in a variety of machines, products, and processes to regulate system output. Most control system design approaches that exploit state feedback control theory are based on the assumption of having sensors to measure the output of a system under control. In short, all the state variables are assumed to be available outputs. However, for many control systems, state variables are not accessible to direct measurement, or the number of measurement devices is limited, possibly due to cost considerations. In order to effectively use state feedback, developing an approach to estimating the system states is inevitable. As a result, an important problem regarding linear multivariable control systems is designing an observer for a given dynamic system. Therefore, many researchers have proposed various methods for designing observers ([121] and references therein). On the other hand, time-delay systems are commonly seen in industrial processes, electrical and mechanical systems, economic and population growth, biological systems, etc. [2229]. Hence, due to their importance and the many predicaments encountered in analyzing and designing them, many interesting research topics have emerged, such as observer design for time-delay systems [30]. In order for the estimated states to be close and converge fast to the actual states of the observed system, many researchers have proposed various observer design methods for time-delay systems ([16, 3036] and references therein). These existing observer design methods guarantee that the state estimation errors can asymptotically converge to zero. Although the existing approaches can be applied to improve the convergence of the estimated states, they may lead to unsatisfactory estimation results, since a large estimation error may occur during the transient period of observation [19].

For delay-free systems, some research has been devoted to considering the observer design issue of reducing the large error during the transient period of observation [14, 6, 19]. However, to the authors’ best knowledge, for the observer design of time-delay systems, there is no extant literature investigating the issue of minimizing estimation errors during the transient period of observation. Therefore, the purpose of this paper is to study the observer design problem for state delay systems by concurrently considering the following two design issues: (i) the state estimation errors quickly converge to zero, and (ii) a quadratic performance measurement of the deviation of the estimates from the actual states is minimized so as to reduce the large error occurring during the transient period of observation.

The remainder of the paper is organized as follows. Section 2 describes the problem statement. A new method of designing an optimal observer for state delay systems is proposed in Section 3. An illustrative example is provided in Section 4. Finally, Section 5 gives some conclusions.

2. Problem Statement

Consider the following time-delay system:with , where is a known constant time delay, k is a constant vector, is the n-dimensional state vector, is the r-dimensional control vector, is the m-dimensional output vector, and is a given time function vector. In addition, , , , and are constant matrices of appropriate dimensions.

For estimating state vector , the state observer for the system (1a), (1b), and (1c) is designed as [37]:where is the estimate of is the estimate of , and is the observer gain matrix to be designed.

The vector of state estimation error is defined aswhere with denoting the state estimation error vector when belongs to , with denoting the state estimation error vector when belongs to , and is the final time which is long enough for state estimation error to decrease to almost zero. Consequently, for

It can be shown thatin which That is,

A quadratic performance measurement of estimation error is defined aswhere the matrices and are specified positive definite weighting matrices. In addition, , and denote the transpose vectors of , and , respectively. Performance measurement is exploited to design a proper observer that can avoid large estimation error during the transient period of observation.

The optimal observer design problem for the time-delay system is to directly find the observer gain matrix such that the state estimation error vector converges to zero, while at the same time the quadratic performance measurement defined by (7) is minimized.

3. Optimal Design of the Observer Gain Matrix

For the state observer in (2a), (2b), and (2c), the problem of guaranteeing the state estimation error vector converges to zero, under the situation that the observer gain matrix has been given beforehand, to derive an asymptotic stability criterion for checking whether the state estimation error vector can asymptotically converge to zero or not. For delay-free systems, some researchers [14, 6, 19] have considered the following two issues concerning the design of linear observers: (i) the eigenvalues of the linear observers are specified to satisfy desired asymptotic convergence performance, and (ii) a quadratic performance measurement of state estimation error is minimized so as to reduce the transient estimation error. But, for time-delay systems, the eigenvalues of the linear observers cannot be specified to satisfy desired convergence performance. Hence, in what follows, an LMI-based (linear-matrix-inequality-based) asymptotic stability criterion is presented to analyze whether the state estimation error vector can asymptotically converge to zero or not, where the observer gain matrix has been specified beforehand.

Theorem 1. The state estimation error vector in (4) asymptotically converges to zero when the observer gain matrix has been specified beforehand if, for the specified observer gain matrix , an symmetric positive definite matrix exists such that the following inequality is satisfied:where

Proof. Let the Lyapunov function candidate for the error dynamic system described by (4) be [38]in whichThen, the time derivative of givesBy using (8), it can be shown that Therefore, the state estimation error vector converges to zero asymptotically.

By using the Schur complete formula [39], (8) becomes the following LMI form:where is the identity matrix. In this form of (13), the powerful LMI approach can be adopted to find the solution of (8) when the observer gain matrix has been previously specified.

By using the asymptotic stability criterion of (8), when the observer gain matrix has been previously given, only the property of asymptotic convergence of the state estimation error vector can be guaranteed. However, this observer design result is not satisfactory, as a large estimation error may occur during the transient period of observation [19]. Hence, the problem considered here is how to specify the observer gain matrix in (2a), (2b), and (2c) so that the constraint of the LMI-based condition of the state estimation error vector asymptotically converging to zero in (13) can be satisfied, while the quadratic performance measurement of estimation error in (7) is simultaneously minimized. The procedure for designing the observer gain matrix can be stated as follows.

Step 1. Check the LMI-based constraint condition of the state estimation error vector asymptotically converging to zero in (13).

Step 2. Minimize the quadratic performance measurement of estimation error in (7) for the error dynamic system described by (4).

In order to help design the observer gain matrix, the orthogonal functions approach (OFA) is employed in this paper. The OFA has been successfully applied to investigate various problems of systems and control [22, 4042]. The key characteristic of OFA is that it converts integral or differential equations into algebraic equations. As a result, the OFA has become very popular in computation because the dynamic equations of the system can be transformed into a set of algebraic equations whose solutions lead to the solution of the original problem [22, 4042]. Suppose that the vectors , and can be developed approximately in terms of truncated orthogonal functions (OF) representation as described by (14)-(16), respectively,where is the number of the terms needed for the OF, represents the OF basis vector, while are the orthogonal functions, and , and denote the coefficient vectors. In addition, , and denote the coefficient matrices.

For , using the OF representations of and in (14) and (16), and applying the following integral property of the OF [22, 40] (5) can be cast into the formwhere , and the matrix is equivalent to the matrix in (17) with is the integration-operational matrix of the OF, in which the entries of matrix depend on the definite choice of the OF basis vector [22, 40].

Using the Kronecker product, the solution of (18) can be obtained aswhere denotes the Kronecker product [43], , , , and is the identity matrix.

For , applying the OF representation of in (15) and using the integral property of the OF in (17), integrating (6) giveswhere , matrix denotes matrix in (17) with , and matrix is the time-delay-operational matrix of orthogonal functions, in which the elements of the delay matrix depend on the specific choice of the OF basis vector [22, 40].

Applying the Kronecker product, the solution of Eq. (20) can be obtained aswhere .

Now, substituting the truncated OF representations of the estimation error vectors and in (14) and (15) into (7), the quadratic performance measurement of (7) becomesin which , and the constant matrices and denote, respectively, the product-integration-matrix of two OF basis vectors having different time intervals [22]. From (19), (21), and (22), it can be seen that the OFA-based computational approach does not limit the sizes of both and , where is the given known constant time delay and is the final time which is given by the control engineer for desiring state estimation error decreased to almost zero.

From the aforementioned procedures for solving and of (18) and (20), respectively, it is clear that if the observer gain matrix is given, then and can be ascertained, and thus the quadratic performance measurement value in (22) relevant to the observer gain matrix can be calculated. Given another matrix of observer gain, another performance measurement value in (22) can be obtained. That is, the performance measurement value of the algebraic form in (22) is in fact dependent on the observer gain matrix , indicatingwhere represent the elements of the observer gain matrix . Thus, the optimal observer design problem of the time-delay system in (1a), (1b), and (1c) is to find the optimal for a symmetric positive definite matrix making the LMI in (13) holds (i.e., such that the property of the state estimation error vector asymptotically converging to zero in the theorem is satisfied). At the same time, the performance measurement of the algebraic form in (22) is minimized. This is equivalent to the constrained static optimization problem:subject to and subject to the constraint of there being a symmetric positive definite matrix to let the LMI in (13) holds for and , in which are the given positive real numbers. This means that, by using the OFA, the optimal observer design problem for the time-delay system described by (1a), (1b), and (1c) can be substituted by constrained static optimization problems with constrained algebraic equations. It greatly simplifies the optimal observer design problem, where the state estimation error vector converges to zero, and at the same time the quadratic performance measurement of (7) is minimized. Both genetic algorithm (GA) and particle swarm optimization (PSO) have acquired considerable attention concerning their ability as effective optimization techniques solving complex problems and have been successfully utilized in many areas [29, 4446]. The GA is adopted in this paper to deal with the complex problem of searching for the optimal observer gain matrix in the algebraic form of (24) subject to the constraint of there being a matrix such that the LMI in (13) holds, where (24) is a complex nonlinear function having continuous variables. In addition, for the parameter optimization problems of complex nonlinear functions having continuous variables, it has been shown that the Hybrid Taguchi-Genetic Algorithm (HTGA) method can obtain better results than existing approaches [47]. The HTGA method is therefore employed in this paper to design the optimal observer for a time-delay system.

Optimal Observer Design Procedures

Step 1 (parameter setting). Input is as follows: population size , number of generations, mutation rate , and crossover rate . Output is as follows: the value of in (22), and the observer gain matrix of dimension .

Step 2 (initialization). First, randomly generate the initial population with the chromosomes of form for executing the HTGA. Next, compute the solutions of and by applying (19) and (21), and compute by using (22) which is defined as the fitness function of the HTGA. Then, the fitness values of the initial population feasible for the LMI-based constraint in (13) are calculated.

Step 3 (HTGA method execution [47]). By integrating (19), (21), and (22), the HTGA method is utilized to search for the optimal observer gain matrix , where a penalty on the fitness value is given for the chromosome violating the LMI-based constraint in (13).

Step 4 (whether or not the stopping criteria are met). If yes, go to Step 5. Otherwise, go back to Step 3.

Step 5. The optimal gain matrix of the observer and the optimal value of J are found.

Remark 2. For studying the stability of time-delay systems, both delay-independent and delay-dependent criteria have been proposed in the literature. The purpose of both delay-independent and delay-dependent criteria is to guarantee that the steady-state errors are reduced. But the purpose of the proposed observer design method with the quadratic performance measurement of (7) is to yield the improvement in observer transient error performance. Thus, the key point of this paper is to propose a new observer design method for reducing the large error occurring during the transient period of observation, not to study the stability of time-delay systems. Therefore, this proposed design approach can be expected to give an adequately damped response such that a small overshoot can be obtained; thus, this provides better physical performance characteristics in the time domain. Though the proposed theorem is a delay-independent stability criterion, the proposed observer design method can also combine delay-dependent stability criterion to reduce the large error occurring during the transient period of observation.

Remark 3. To conquer the inherent difficulty of giving a priori proper bound of parameters in which the optimal solution is placed, it is desirable to allow the HTGA to dynamically expand the search space. To this end, the search-space expansion schemes developed in [48] are also integrated into the HTGA approach to design the optimal observer gain matrix In addition, in the HTGA approach, the bound of parameters of the observer gain matrix can be specified to avoid the case of a high gain value from practical considerations.

4. Illustrative Example

One example is given in this section to verify the efficiency of the proposed method for designing the optimal observer gain matrix such that (i) the state estimation error vector converges to zero and (ii), at the same time, the quadratic performance measurement of algebraic (7) is minimized for reducing the large error occurring during the transient period of observation. The following evolutionary environments of the HTGA are adopted in this example; the maximum number of generation is 100, the crossover rate is 0.95, the population size is 30, and the mutation rate is 0.01. The type of orthogonal function used in this example is the shifted Chebyshev series with [22, 40].

Consider the time-delay system described by (1a), (1b), and (1c), whereIn this example, , and By using the proposed optimal observer design method, and applying the LMI Toolbox [39], the optimal observer gain matrix can be obtained asand the symmetric positive definite matrix can be derived asFor the optimal observer, the value of is

To the authors’ best knowledge, for the observer design of time-delay systems, there is no extant literature investigating the issue of minimizing estimation errors during the transient period of observation. That is, for state delay systems, the extant approaches do not consider to minimizing estimation errors during the transient period of observation. So, no existing methods can be used to compare the proposed optimal design approach. Thus, only the comparisons of the state estimation errors of optimal estimation and nonoptimal estimation are given in the simulation. Letting and solving the LMI in (13) to find the matrices and , then the nonoptimal observer gain matrix can be obtained aswhereand the value of is The comparisons of the state estimation errors and for the to-be-estimated states and of optimal estimation and nonoptimal estimation are given in Figures 1 and 2, respectively. It can be seen that, by using the proposed optimal design method, the estimated state errors can quickly converge to zero, and at the same time the transient estimation errors have been significantly reduced during the observation period. Furthermore, the performance measurement value for the proposed optimal design method is 93.64% lower than that for the nonoptimal design approach only using the LMI method. That is, the proposed optimal observer design method can obtain a smaller estimation error during the transient period of observation.

5. Conclusions

A new approach has been proposed in this paper to design the observer gain matrix for state delay systems such that (i) the constraint of the LMI-based condition of the state estimation error vector asymptotically converging to zero in (13) can be satisfied and (ii) at the same time the quadratic performance measurement of (7) is minimized for reducing the large error occurring during the transient period of observation. One illustrative example has been given to verify the effectiveness and efficiency of the proposed optimization method on the improvement of transient performance of state estimation errors. Based on the simulation results, it is evident that not only do the estimated state errors quickly converge to zero, but also the performance measurement values for the proposed optimal design approach are significantly lower than those of the nonoptimal design approach only using the LMI method. The constraint of the LMI-based condition makes the state estimation error vector asymptotically converge to zero so that the steady-state errors are reduced. Besides, the proposed observer design method with the quadratic performance measurement gives a penalty for the transient error. It has been shown in this paper that there exists a physically realizable, optimally designed observer gain that can yield the improvement in observer transient error performance. Therefore, this proposed design approach can be expected to give an adequately damped response such that a small overshoot can be obtained; thus, this provides better physical performance characteristics in the time domain. Last but not least, to the authors’ best knowledge, this paper is a pioneering work in dealing with the observer design issue of transient estimation performance improvement for state delay systems.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.