Research Article  Open Access
Optimal Control Systems by TimeDependent Coefficients Using CAS Wavelets
Abstract
This paper considers the problem of controlling the solution of an initial boundaryvalue problem for a wave equation with timedependent sound speed. The control problem is to determine the optimal sound speed function which damps the vibration of the system by minimizing a given energy performance measure. The minimization of the energy performance measure over sound speed is subjected to the equation of motion of the system with imposed initial and boundary conditions. Using the modal space technique, the optimal control of distributed parameter systems is simplified into the optimal control of bilinear timeinvariant lumpedparameter systems. A waveletbased method for evaluating the modal optimal control and trajectory of the bilinear system is proposed. The method employs finite CAS wavelets to approximate modal control and state variables. Numerical examples are presented to demonstrate the effectiveness of the method in reducing the energy of the system.
1. Introduction
Dynamic stability related to parametric resonance is a very important factor in structural dynamics. For example, instability caused by parametric resonance is believed to be the reason for the famous Tacoma bridge collapse in 1940 [1]. A suitable control of the coefficients may provide an effective protection against this phenomenon.
Control in the coefficients is known to be a very effective method in structures governed by elliptic equations [1]. However, not much information is known about the effect produced by control in coefficient for hyperbolic equations [2, 3]. In this paper, we study a control problem for a structure dynamic system governed by a hyperbolic equation where the control is a time dependent coefficient.
The model considered in this work is motivated by recent developments in the area of smart materials [4]. The properties of these materials can be changed by applying external fields, such as electrical, magnetic, or temperature; this is referred to as a phase transformation.
A structure made with such a material is considered where control consists of eliminating a transient disturbance in the structure by varying the material properties in the response to the deformation. The modal dynamic of a structure is governed by a scalar wave equation, where the control variable is the sound speed in the medium. The basic bilinear optimal control problem becomes the minimization of the energy function of the system in a given period of time with a minimum sound speed. Using modal expansion, the optimal control of the distributed parameter system is reduced to the optimal control of a bilinear timevarying lumped parameter system. The parameterization approach is used to approximate the statevariable and each component of the control variable using finiteterm wavelets with unknown coefficients. Therefore, the quadratic problem is transformed into a mathematical programming problem with the objective of minimizing the unknown coefficients to give suboptimal solution of the problem. A necessary condition for the optimality of the unknown coefficients is derived as a system of linear algebraic equations for which the solution is used to obtain the optimal control sound speed and optimal state function.
The bilinear system is a kind of nonlinear system where some related problems such as optimal control are much more difficult to solve than those of linear systems. In literature, many authors [5–9] have tried various methods to overcome the difficulties of solving bilinear systems. In this paper, the focus will be on obtaining the optimal state solution of a wave equation governed by a bilinear system using CAS wavelets taking advantage of some needed properties of this type of wavelets [10, 11]. Compared to conventual method such as Fourier series or finite elements, CAS wavelets with their local properties enable arbitrary functions (even with discontinuity) to be approximated more efficient. To demonstrate the effectiveness of the proposed approach, numerical results will show confirm that the proposed method significantly minimizes the energy of the system.
2. Optimal Control Problem Setting
Let be an open, bounded, and simply connected subset of dimensional Euclidean space . Let denote a given time interval with finite terminal time . Consider the wave equation, defined on : where is the Laplacian operator, and is the disturbance of position and time . The wave speed is assumed to be a function of time. For simplicity, let satisfy the boundary and initial conditions: where and .
Let the admissible control set be Associated with the wave equation (2.1) is the modified energy at terminal time : where , , and are weighing constants satisfying the condition and . The last term on the righthand side of (2.5) is a penalty term on control energy.
The optimal control problem is stated as follows: determine the optimal control function such thatsubject to (2.1) and (2.2).
3. Control Problem in Modal Space
We pose the problem at hand as a control problem for an finite system of ordinary differential equations by using modal space expansion. Let where are normalized eigenfunctions associated with eigenvalues . This implies that satisfies the eignvalueproblem It can be shown that the set forms an orthonormal set, and hence satisfies with initial conditions In view of the expansion (3.1), the performance index becomes The optimal control problem (2.6) is now modified as follows: determine such thatsubject to (3.3) and (3.4).
4. Properties of the CAS Wavelets
4.1. CAS Wavelets
Wavelets have been used by many scientists and engineers to solve several problems in areas such as signal and image processing, control problems, and stochastic problems. Wavelets are mathematical functions that are constructed using dilation and translation of a single function called the mother wavelet denoted by and must satisfy certain requirements. If the dilation parameter is and translation parameter is , then we have the following family of wavelets: Restricting and to discrete values, such as , , , and and are positive integers, gives where form a basis for . If and , then it is clear that the set forms an orthonormal basis for .
The CAS wavelets employed in this paper are defined as where The set of CAS wavelets forms an orthonormal basis for . This implies that any function defined over can be expanded as where and and are vectors given by
4.2. Operational Matrices of Integration
The integration of the function in (4.5) is given by where is an matrix, called the operational matrix, and is given by [12] in which is a zero matrix and and are matrices given by The integration of the product of two CAS function vectors is given by The product operational matrix of the CAS wavelet is given by where the matrix is given in (4.7) and is an given by [11] where , are given by
5. CASWaveletsBased Approach
To redefine the wavelet functions over the interval , we let . Then (3.3), (3.4), and (3.5), respectively, become
Using the expansion in (4.5) gives where , , , , and are vectors defined as in (4.7). Furthermore, Substituting (5.5) in (5.1) yields and hence
Using (4.12) leads to
Multiplying (5.8) by integrating, and using (4.11) give or provided that exists. Substituting equations (5.4) and (5.5) into equation (3.5) convert the performance index into a function of and hence to optimize , we solve
6. Numerical Example
Consider the wave equation with initial conditions
For the sake of illustration, the following parameters were assumed: The performance index was computed for the optimal control and compared with the performance index for the controllers , , , and . The results are summarized in Table 1.

It is observed that the proposed control is effective in significantly reducing the performance index of the problem.
7. Conclusion
A control for a wave equation where the control is a time dependent coefficient is considered. A modal space technique simplifies the optimal control of a distributed parameter system into the optimal control of a bilinear timeinvariant lumpedparameter system. A Galerkin CAS waveletbased method was developed to solve this bilinear optimal control problem. The main aspect of the proposed approach resides in converting the optimization problem into a mathematical programming problem where the necessary conditions of optimality are derived as a system of algebraic equations . A test example, which includes a variable coefficient and onedimensional hyperbolic equation, demonstrates the capability of the proposed GalerkinWavelet approach for solving optimal control problems governed by bilinear systems. Moreover, the numerical simulations show that the optimal control procedure led to a substantial damping in the bilinear system energy.
This method may be extended to treat a more general setting where the coefficients are and dependent. That is, the wave speed function is a controllable function of the form [13].
References
 V. Berdichersky, V. Jikov, and G. Papanicolcou, Eds., Homogenization, vol. 50 of Series on Advances in Mathematics for Applied Sciences, World Scientific Publishing, River Edge, NJ, USA, 1999. View at: Zentralblatt MATH  MathSciNet
 K. A. Lurie, “Control in the coefficients of linear hyperbolic equations via spaciotemporal components,” in Homogenization, vol. 50 of Series on Advances in Mathematics for Applied Sciences, pp. 285–315, World Scientific, River Edge, NJ, USA, 1999. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 K. A. Lurie, “Some new advances in the theory of dynamic materials,” Journal of Elasticity, vol. 72, no. 1–3, pp. 229–239, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Hurlebaus and L. Gaul, “Smart structure dynamics,” Mechanical Systems and Signal Processing, vol. 20, no. 2, pp. 255–281, 2006. View at: Publisher Site  Google Scholar
 A. Chambolle and F. Santosa, “Control of the wave equation by timedependent coefficient,” ESAIM: Control, Optimisation and Calculus of Variations, vol. 8, pp. 375–392, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 H. R. Joshi, “Optimal control of the convective velocity coefficient in a parabolic problem,” Nonlinear Analysis: Theory, Methods & Applications, vol. 63, no. 5–7, pp. 1383–1390, 2005. View at: Publisher Site  Google Scholar
 S. Lenhart, M. Liang, and V. Protopopescu, “Identification problem for a wave equation via optimal control,” in Control of Distributed Parameter and Stochastic Systems, vol. 141 of IFIP Conference Proceedings, pp. 79–84, Kluwer Academic Publishers, Boston, Mass, USA, 1999. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 F. Maestre, A. Münch, and P. Pedregal, “A spatiotemporal design problem for a damped wave equation,” SIAM Journal on Applied Mathematics, vol. 68, no. 1, pp. 109–132, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Tatari and M. Dehghan, “Identifying a control function in parabolic partial differential equations from over specified boundary data,” Computers & Mathematics with Applications, vol. 53, no. 12, pp. 1933–1942, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Razzaghi and S. Yousefi, “Sinecosine wavelets operational matrix of integration and its applications in the calculus of variations,” International Journal of Systems Science, vol. 33, no. 10, pp. 805–810, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Yousefi and A. Banifatemi, “Numerical solution of Fredholm integral equations by using CAS wavelets,” Applied Mathematics and Computation, vol. 183, no. 1, pp. 458–463, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Danfu and S. Xufeng, “Numerical solution of integrodifferential equations by using CAS wavelet operational matrix of integration,” Applied Mathematics and Computation, vol. 194, no. 2, pp. 460–466, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. Liang, “Bilinear optimal control for a wave equation,” Mathematical Models & Methods in Applied Sciences, vol. 9, no. 1, pp. 45–68, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2009 Taher Abualrub 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.