PDF version | Download Matlab Code | Download Scilab Code | **GITHUB**

In this work, we address the optimal control of parameter-dependent systems. We introduce the notion of averaged control in which the quantity of interest is the average of the states with respect to the parameter family . More precisely, we are interested in solving the minimization problem

where is the space of admissible controls and the average state target. The optimization problem (1) is subject to the finite dimensional linear control system

(2)

In (2), the vector valued function is the state of the system, is a -matrix governing its free dynamics and is a -component control, with , entering and acting on the system through the control operator , a parameter-dependent matrix. Moreover, the initial datum to be controlled is independent of the parameter , but the state of the system itself depends on . The effective value of the parameter being unknown, we aim at choosing a control that would perform optimally in an averaged sense, this is, rather than controlling specific realizations of the state, the average with respect to is controlled. This allows building a control independent of the parameter and making a robust compromise of all the possible realizations of the system for the various possible values of the unknown parameter

We use the classical gradient descent method based on the adjoint methodology, and obtain the corresponding adjoint system for (2),

(3)

To minimize the functional in (1), we take the steepest descent direction given by

(4)

Hence, the new control reads as

(5)

for some small enough.

We have also used the conjugate gradient method in order to reach faster the optimal control. In order to be able to apply this method the state vector has been split as

(6)

where is the solution to the controlled system with zero initial condition,

(7)

and solves the free dynamics problem,

(8)

The functional in 1 can be expressed as

(9)

We introduce the linear operator

(10)

and its dual counterpart,

(11)

where is solution to

(12)

By doing this we can write the functional gradient as

(13)

After having defined and we can apply the conjugate gradient method to solve the control problem.

**Algorithm 1** Optimal control with Conjugate Gradient Method

**Require:** , , , , , , ,

- compute
**While****do**

###### Bibliography

**[1]** E. Zuazua (2014) *Averaged Control.* Automatica, 50 (12), p. 3077-3087.