This paper considers the problem of controlling the solution of an initial boundary-value problem for a wave equation with time-dependent 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 time-invariant lumped-parameter systems. A wavelet-based 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 time-varying lumped parameter system. The parameterization approach is used to approximate the state-variable and each component of the control variable using finite-term 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 [59] 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 right-hand 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 eignvalue-problem 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. CAS-Wavelets-Based 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 time-invariant lumped-parameter system. A Galerkin CAS wavelet-based 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 one-dimensional hyperbolic equation, demonstrates the capability of the proposed Galerkin-Wavelet 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].