Abstract and Applied Analysis

Volume 2013 (2013), Article ID 535979, 9 pages

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

## Robustness of Operational Matrices of Differentiation for Solving State-Space Analysis and Optimal Control Problems

^{1}Department of Mathematics, Islamic Azad University, Zahedan Branch, Zahedan, Iran^{2}Department of Mathematics, Universiti Putra Malaysia, 43400 Serdang, Malaysia

Received 13 January 2013; Accepted 8 March 2013

Academic Editor: Mustafa Bayram

Copyright © 2013 Emran Tohidi 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

The idea of approximation by monomials together with the collocation technique over a uniform mesh for solving *state-space analysis* and *optimal control* problems (OCPs) has been proposed in this paper. After imposing the Pontryagins maximum principle to the main OCPs, the problems reduce to a linear or nonlinear boundary value problem. In the linear case we propose a monomial collocation matrix approach, while in the nonlinear case, the general collocation method has been applied. We also show the efficiency of the operational matrices of differentiation with respect to the operational matrices of integration in our numerical examples. These matrices of integration are related to the Bessel, Walsh, Triangular, Laguerre, and Hermite functions.

#### 1. Introduction

In the last four decades, numerical methods which are based on the operational matrices of integration (especially for orthogonal polynomials and functions) have received considerable attention for dealing with a huge size of applied mathematics problems such as *state*-*space analysis* and *optimal control*. The key idea of these methods is based on the integral expression
where is an arbitrary basis vector and is a constant matrix, called the operational matrix of integration. The matrix has already been determined for many types of orthogonal (or nonorthogonal) bases such as Walsh functions [1–3], block-pulse functions [4], Laguerre polynomials [5], Chebyshev polynomials [6], Legendre polynomials [7], Hermite polynomials [8], Fourier series [9], Bernstein polynomials [10], and Bessel functions [11]. As a primary research work which was based on the operational matrices of integration, one can refer to the work of Corrington [1]. In [1], the author proposed a method of solving nonlinear differential and integral equations using a set of Walsh functions as the basis. His method is aimed at obtaining piecewise constant solutions of dynamic equations and requires previously prepared tables of coefficients for integrating Walsh functions. To alleviate the need for such tables, Chen and Hsiao [2, 3] introduced an operational matrix to perform integration of Walsh functions. This operational matrix approach has been applied to various problems such as time-domain analysis and synthesis of linear systems, and piecewise constant-feedback-gain determination for optimal control of linear systems and for inverting irrational Laplace transforms.

On the other hand, since the beginning of 1994, the Bernoulli, Chebyshev, Laguerre, Bernstein, Legendre, Taylor, Hermite, and Bessel matrix methods have been used in the works [12–24] to solve high-order linear and nonlinear differential (including hyperbolic partial differential equations) Fredholm Volterra integrodifferential difference delay equations and their systems. The main characteristic of these approaches is based on the operational matrices of differentiation instead of integration. The best advantage of these techniques with respect to the integation methods is that, in the fundamental matrix relations, there is not any approximation symbol, meanwhile in the integration forms such as (1) the approximation symbol could be seen obviously. In other words where is the operational matrix of differentiation for any selected basis such as the previously mentioned polynomials, functions, and truncated series. The readers can see that there is no approximation symbol in (2), meanwhile this can be seen in (1) by using operational matrices of integration. For justifying this expression, one can refer to this subject that after differentiating an th degree polynomial we usually reach to a polynomial which has less than th degree. However, in the integration processes the degree of polynomials would be increased.

In this paper, we generalize a new collocation matrix method that was applied for solving a huge size of applied mathematics models (see for instance [16] and the references therein), to several special classes of systems of ordinary differential equations (ODEs). Two important classes of such systems of ODEs are(i)*State space analysis,*(ii)*Hamiltonian system,*which are the necessary (and also are sufficient in several special cases) conditions for optimality of the solutions of OCPs, originate from the PMP, and have considerable importance in optimal control and calculus of variation.

We again emphasized that the methods that are based on the operational matrices of differentiation are more accurate and effective with regard to the integration ones. We illustrate this fact through several examples for dealing with the previously mentioned systems in the section of numerical examples. It should be noted that one of the best tools for the integration approaches is using high accurate Gauss quadrature rules such as the method of [25, 26]. However, more CPU times are required for using such quadrature rules, and also the matrix coefficient associated to these methods is ill-conditioned usually and should be preconditioned.

The remainder of this paper is organized as follows. In Section 2, the considered problems such as *state-space analysis* and *Hamiltonian system* are introduced. In Section 3, the fundamental matrix relations together with the method of obtaining approximate solutions are described. In Section 4, several numerical examples are provided for confirming high accuracy of the proposed method. The last Section is devoted to the conclusions.

#### 2. Problems Statement

In this section two types of problems are considered. In the first subsection, we show that how the *Hamiltonian systems* can be obtained in both linear and nonlinear forms. In the second subsection, we introduce a general form of *state-space analysis* problems.

##### 2.1. Hamiltonian Systems

###### 2.1.1. Linear Quadratic Optimal Control Problems

In this part, we consider the following linear optimal control problem (OCP): where , , and . The control is an admissible control if it is piecewise continuous in for . Its values belong to a given closed subset of . The input is derived by minimizing the quadratic performance index , where is positive semidefinite matrix and is a positive definite matrix. We consider Hamiltonian for system (3) as where is the costate vector.

According to the Pontryagin’s maximum principle, we have [27] The optimal control is computed by [27] where and are the solution of the Hamiltonian system:

###### 2.1.2. Nonlinear Quadratic Optimal Control Problems

Consider the nonlinear dynamical system with denoting the state variable, the control variable, and is the given initial state at . Moreover, and are two continuously differentiable functions in all arguments. Our aim is to minimize the quadratic objective functional subject to the nonlinear system (8), for , positive semidefinite and positive definite matrices, respectively. Since the performance index (9) is convex, the following extreme necessary conditions are also sufficient for optimality [28]: where is referred to the Hamiltonian. Equivalently, (10) can be written in the form of where is the costate vector with the th component , and with , .

Also the optimal control law is obtained by

For solving such a two-point boundary value problem (TPBVP) in (11), we apply a similar collocation method that was proposed in [29].

##### 2.2. State Space Analysis Problems

In this part, we consider the following *state space analysis* problem:
where , , and are known, meanwhile is unknown. The goal is to obtain the approximation of in (13). The previously mentioned system (13) is similar to *Hamiltonian system* (7) and the scheme of their solutions is the same.

*Remark 1. *We recall that the main goal of this paper is to approximate the solution of the systems (7), (11), and (13) by applying a new matrix method which is based on the operational matrix of differentiation and also the uniform collocation scheme in the parts of *Hamiltonian systems* and *state space analysis* problems.

#### 3. Fundamental Matrix Relations and Method of the Solution

In this section, by using the collocation points and the matrix relations between the monomials and their derivatives, we will find the approximate solution of the system (7) expressed in the truncated monomial series form (assuming that and and also together with is independent of time , that is, and ) so that and ; are the unknown coefficients.

Let us consider the desired solutions and , of (7) defined by the truncated monomial series (14). We can write the approximate solutions, which are given in relation (14) in the matrix form where and , .

The matrix form of the relation between the matrix and its th derivative is so that By using the relations (15) and (16), we have the following relations: Thus, we can express the matrices and as follows: where

Now, we can restate the system (7) in the matrix form where Applying the following collocation points in (21): yields to equations as follows: All of the these equations can be written in the following matrix form: where where denotes the Kronecker product and is the identity matrix of dimension .

With the aid of relation (19) and the collocation points (23), we gain which can be written as where If the relation (28) is substituted into (25), the fundamental matrix equation is obtained as Thus, the fundamental matrix equation (30) corresponding to (7) can be written in the form which corresponds to a linear system of algebraic equations in the unknown monomial coefficients so that By the aid of the relation (19), the matrix form for the boundary conditions which are given in (7) can be written as

Finally, by replacing the rows of the matrices by the last rows of the matrices , we obtain the new augmented matrix The unknown monomials coefficients which exist in are determined by solving this linear system, and hence , are substituted in (14). Therefore, we find the approximated solutions

We can easily check the accuracy of the method. Since the truncated monomial series (14) are the approximate solutions of (7), when the functions , and their derivatives are substituted in (7), the resulting equation must be satisfied approximately; that is, for , and , ( positive integer).

If ( positive integer) is prescribed, then the truncation limit is increased until the difference at each of the points becomes smaller than the prescribed , see [24].

*Remark 2. *We must recall that a similar approach can be applied for the *state space analysis* problem (13). Moreover, as we say before, for solving a general nonlinear system of ODEs such as (11), we apply a generalization of the collocation method that was proposed in [29].

#### 4. Numerical Examples

In this section, several numerical examples are given to illustrate the accuracy and effectiveness of the proposed method. All calculations are designed in MAPLE 13 and run on a Pentium 4 PC Laptop with 2 GHz of CPU and 2 GB of RAM. In this regard, in tables and figures, we report the absolute error functions associated to the trajectory and control variables and also the approximated values of performance index. In the first example, we provide an OCP that was recently considered by a new method (which is based on the operational matrix of integration of Triangular functions) [30] and reach to more accurate results. Also, in the second example, we consider another OCP (or *Hamiltonian system*) with time variant dynamical system, in which our results have more accuracy and credit with regard to methods [30, 31]. Moreover, we consider a nonlinear OCP as our third numerical illustration. In the fourth example, we provide a *state space analysis* problem together with a full comparison with the methods that are based on the operational matrices of integration such as Bessel [11] and Laguerre [32].

*Example 3 (see [30] linear Hamiltonian system). *Consider the problem of minimizing
subject to
The purpose is to find the optimal control which minimizes (37) subject to (38). The Optimal value of performance index for this problem is and also exact solutions have been given in [30] as
Since the objective function of this OCP is convex, therefore the following necessary conditions (i.e., *linear Hamiltonian system*) for optimality are also sufficient:
Hence, we need to solve the previous system of differential equations such that the obtained numerical solution is the optimal solution of problem (37)-(38). It should be noted that according to (6) the optimal control is computed by , where is the solution of the previous system.

We solve this problem by using our proposed method in the cases of , and 8. The approximated solutions corresponding to these values of are provided below
The associated performance indexes for the selected values of are , , , , and . We provide the , associated to our proposed method (PM) and a new method that is based on the operational matrix of integration of Triangular functions [30] for different values of in Table 1. It can be seen from this table that our obtained results for such considered values of (i.e., 4, 5, 6, 7, and 8) are the same and equal to the obtained results of [30] for higher values of such as , and 64 in computation of . Moreover, our results corresponding to the are more accurate with regard to the method of [30] even by choosing lower values of .

*Example 4 (see [30, 31] linear Hamiltonian system). *Consider the linear time-varying system
with the cost functional
The problem is to obtain the optimal control which minimizes (43) subject to (42). The optimal control is
where is the feedback controller gain matrix and the solution of the Riccati equation [30]
According to the optimality conditions (5) and (6) we have
We first solve the previous system and obtain the numerical solutions and for , and then solve (45) by ODE solver commands which exist in MAPLE 13. Since, then our numerical results of are equal to , that is, . The numerical results of system (46), which are obtained by the proposed method could be deduced as
Also, the exact solution of the Riccati equation (45) at the uniform mesh in the interval are , , , , , , , , , , and . In Table 2, we provide the absolute values of errors at the selected points for the previously considered values of together with the same errors associated with other methods [30, 31]. Again, we can see the accuracy of method with regard to the methods that are based on operational matrices of integration.

*Example 5 (nonlinear Hamiltonian system). *As our third illustration, consider the following nonlinear optimal control problem:
Trivially , , , , , and . As mentioned in Section 2.2, we solve the following system of ordinary differential equations:
Also the optimal control law is given by

Similar to the linear cases, we suppose that the state and costate variables could be written in terms of linear combination of monic polynomials which are defined in Section 3, with the unknown monomial coefficients. These coefficients will be determined after imposing the previous system of differential equations at the uniform mesh in the interval . In other words, applying these collocation points to the main system together with the considered boundary conditions on and transforms the basic problem to the corresponding system of nonlinear algebraic equations. By assuming different values of such as 5, 7, and 9, we solve the previously mentioned system. In Table 3, we provide the approximated performance index , which is obtained by our proposed method and also the difference between the approximated for the considered values of .

*Example 6 (see [11, 32] state-space analysis). *We consider a linear-time invariant state equation
where
We are given that the input is the unit step function in the interval and the initial state is
The exact solution for (51) is
We solve this problem by using our proposed method in the cases of , and 20. The approximated solutions corresponding to the , and 10 are provided later

An interesting comparison between our presented method (PM) and the method of [11] in the absolute value of the errors at the uniform mesh for , and 20 is considered in Table 4. Moreover, the error histories in the computational interval associated with our method for and 5 together with the error history of the method [11] (Bessel Integration) for and also the error history of the method [32] (Laguerre Integration) for are depicted in Figure 1. From this figure one can see the applicability and high accuracy of the presented method with respect to the methods which are based on operational matrices of integration such as [11, 32].

#### 5. Conclusions

The aim of this paper is to present an indirect approach for solving optimal control problems using truncated monomial series together with the collocation method on a uniform mesh. Our method applies an operational matrix of differentiation which has more efficiency with respect to the integration ones. Operational matrices of differentiation have several specific properties that other integration ones do not have them. One of them is that through using differentiation ones, we solve our problem directly and do not need to integrate the dynamical system. Another property is that the fundamental relations, which are based on differentiation matrices, are the exact relations, meanwhile the methods which are based on integration matrices impose the approximation to the main problem. These properties are shown through several numerical examples such as state-space analysis and specially optimal control problems.

#### Conflict of Interests

The authors declare that they do not have any conflict of interests in their submitted paper.

#### Acknowledgment

The authors are grateful to the comments of the referees, which have led to several improvements in the presentation of this paper.

#### References

- M. S. Corrington, “Solution of differential and integral equations with Walsh functions,”
*IEEE Transactions Circuit Theory*, vol. CT-20, no. 5, pp. 470–476, 1973. View at Google Scholar · View at Scopus - C. F. Chen and C. H. Hsiao, “Time-domain synthesis via Walsh functions,”
*Proceedings of the Institution of Electrical Engineers*, vol. 122, no. 5, pp. 565–570, 1975. View at Google Scholar · View at Scopus - C. F. Chen and C. H. Hsiao, “WALSH series analysis in optimal control,”
*International Journal of Control*, vol. 21, no. 6, pp. 881–897, 1975. View at Google Scholar · View at Scopus - N. S. Hsu and B. Cheng, “Analysis and optimal control of time-varying linear systems via block-pulse functions,”
*International Journal of Control*, vol. 33, no. 6, pp. 1107–1122, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. Hwang and Y. P. Shih, “Laguerre series direct method for variational problems,”
*Journal of Optimization Theory and Applications*, vol. 39, no. 1, pp. 143–149, 1983. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. N. Elnagar, “State-control spectral Chebyshev parameterization for linearly constrained quadratic optimal control problems,”
*Journal of Computational and Applied Mathematics*, vol. 79, no. 1, pp. 19–40, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. N. Paraskevopoulos, “Legendre series approach to identification and analysis of linear systems,”
*Institute of Electrical and Electronics Engineers*, vol. 30, no. 6, pp. 585–589, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Th. Kekkeris and P. N. Paraskevopoulos, “Hermite series approach to optimal control,”
*International Journal of Control*, vol. 47, no. 2, pp. 557–567, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. N. Paraskevopoulos, P. D. Sparis, and S. G. Mouroutsos, “The Fourier series operational matrix of integration,”
*International Journal of Systems Science*, vol. 16, no. 2, pp. 171–176, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. H. Doha, A. H. Bhrawy, and M. A. Saker, “Integrals of Bernstein polynomials: an application for the solution of high even-order differential equations,”
*Applied Mathematics Letters*, vol. 24, no. 4, pp. 559–565, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. N. Paraskevopoulos, P. G. Sklavounos, and G. Ch. Georgiou, “The operational matrix of integration for Bessel functions,”
*Journal of the Franklin Institute*, vol. 327, no. 2, pp. 329–341, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. H. Bhrawy, E. Tohidi, and F. Soleymani, “A new Bernoulli matrix method for solving high-order linear and nonlinear Fredholm integro-differential equations with piecewise intervals,”
*Applied Mathematics and Computation*, vol. 219, no. 2, pp. 482–497, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - A. H. Bhrawy and M. A. Al-Shomrani, “A shifted Legendre spectral method for fractional order multi point boundary value problems,”
*Advances in Differerence Equations*, vol. 2012, article 8, 2012. View at Google Scholar - E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A Chebyshev spectral method based on operational matrix for initial and boundary value problems of fractional order,”
*Computers & Mathematics with Applications*, vol. 62, no. 5, pp. 2364–2373, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Gülsu, B. Gürbüz, Y. Öztürk, and M. Sezer, “Laguerre polynomial approach for solving linear delay difference equations,”
*Applied Mathematics and Computation*, vol. 217, no. 15, pp. 6765–6776, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Sezer, S. yalçinbaş, and N. Şahin, “Approximate solution of multi-pantograph equation with variable coefficients,”
*Journal of Computational and Applied Mathematics*, vol. 214, no. 2, pp. 406–416, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. Tohidi, “Legendre approximation for solving linear HPDEs and comparison with Taylor and Bernoulli matrix methods,”
*Applied Mathematics*, vol. 3, no. 5, pp. 410–416, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - E. Tohidi, A. H. Bhrawy, and Kh. Erfani, “A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation,”
*Applied Mathematical Modelling*, vol. 37, no. 6, pp. 4283–4294, 2013. View at Publisher · View at Google Scholar - E. Tohidi, “Bernoulli matrix approach for solving two dimensional linear hyperbolic partial differential equations with constant coefficients,”
*American Journal of Computational and Applied Mathematics*, vol. 2, no. 4, pp. 136–139, 2012. View at Publisher · View at Google Scholar - F. Toutounian, E. Tohidi, and S. Shateyi, “A collocation method based on Bernoulli operational matrix for solving high order linear complex differential equations in a rectangular domain,”
*Abstract and Applied Analysis*, Article ID 823098, 2013. View at Google Scholar - F. Toutounian, E. Tohidi, and A. Kiliçman, “Fourier operational matrices of differentiation and transmission: introduction and applications,”
*Abstract and Applied Analysis*, Article ID 198926, 2013. View at Google Scholar - S. Yalçinbaş, M. Aynigül, and M. Sezer, “A collocation method using Hermite polynomials for approximate solution of pantograph equations,”
*Journal of the Franklin Institute*, vol. 348, no. 6, pp. 1128–1139, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. A. Yousefi and M. Behroozifar, “Operational matrices of Bernstein polynomials and their applications,”
*International Journal of Systems Science*, vol. 41, no. 6, pp. 709–716, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Yuzbasi,
*Bessel polynomial solutions of linear differential, integral and integro-differential equations [M.S. thesis]*, Graduate School of Natural and Applied Sciences, Mugla University, 2009. - O. R. N. Samadi and E. Tohidi, “The spectral method for solving systems of Volterra integral equations,”
*Journal of Applied Mathematics and Computing*, vol. 40, no. 1-2, pp. 477–497, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - E. Tohidi and O. R. N. Samadi, “Optimal control of nonlinear Volterra integral equations via Legendre polynomials,”
*IMA Journal of Mathematical Control and Information*, vol. 30, no. 1, pp. 67–83, 2013. View at Google Scholar - L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze, and E. F. Mishchenko,
*The Mathematical Theory of Optimal Processes*, Wiley Interscience, 1962. View at MathSciNet - H. P. Geering,
*Optimal Control with Engineering Applications*, Springer, Berlin, Germany, 2007. - E. Tohidi, O. R. N. Samadi, and M. H. Farahi, “Legendre approximation for solving a class of nonlinear optimal control problems,”
*Journal of Mathematical Finance*, vol. 1, pp. 8–13, 2011. View at Publisher · View at Google Scholar - K. Maleknejad and H. Almasieh, “Optimal control of Volterra integral equations via triangular functions,”
*Mathematical and Computer Modelling*, vol. 53, no. 9-10, pp. 1902–1909, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. L. Wang, R. Y. Chang, and S. Y. Yang, “Analysis and optimal control of time-varying systems via generalized orthogonal polynomials,”
*International Journal of Control*, vol. 44, no. 4, pp. 895–910, 1986. View at Google Scholar · View at Scopus - F. C. Kung and H. Lee, “Solution of linear state-space equations and parameter estimation in feedback systems using laguerre polynomial expansion,”
*Journal of the Franklin Institute*, vol. 314, no. 6, pp. 393–403, 1982. View at Google Scholar · View at Scopus