ASYMPTOTIC EXPANSION OF SOLUTION TO SINGULARLY PERTURBED OPTIMAL CONTROL PROBLEM WITH CONVEX INTEGRAL QUALITY FUNCTIONAL WITH TERMINAL PART DEPENDING ON SLOW AND FAST VARIABLES

We consider an optimal control problem with a convex integral quality functional for a linear system with fast and slow variables in the class of piecewise continuous controls with smooth constraints on the control {(x) over dot(epsilon) = A(11)x(epsilon) + A(12)y(epsilon) + B(1)u, t is an element of [0, T], parallel to u parallel to <= 1, epsilon(y) over dot(epsilon) = A(22)y(epsilon) + B(2)u, x(epsilon)(0) = x(0), y(epsilon)(0) = y(0), del phi(2)(0) = 0, J(u) := phi(1) (x(epsilon)(T)) + phi(2)(y(epsilon)(T)) + integral(T)(0)parallel to u(t)parallel to(2) dt -> min, where x is an element of R-n, y is an element of R-m, u is an element of R-r; A(ij) and B-i, i, j = 1,2, are constant matrices of corresponding dimension, and the functions phi(1)(.), phi(2)(.) are continuously differentiable in R-n,R-m, strictly convex, and cofinite in the sense of the convex analysis. In the general case, for such problem, the Pontryagin maximum principle is a necessary and sufficient optimality condition and there exist unique vectors l(epsilon) and rho(epsilon) determining an optimal control by the formula u(epsilon)(T-t) := C-1,C-epsilon*(t)l(epsilon) + C-2,C-epsilon*(t)rho(epsilon)/S(parallel to C-1,C-epsilon*(t)l(epsilon) + C-2,C-epsilon*(t)rho(epsilon)parallel to), where C-1,C-epsilon* (t) := B-1*e(A)(11*)(t) + epsilon B--1(2)*W*(epsilon)(t), C-2(,epsilon)*(t) := epsilon B--(1)2*e(A)(22*t/epsilon), W-epsilon(t) := e(A11t )integral(t)(0) e-(A11 tau) A(12)e(A22 tau/epsilon) d tau, S(xi) := {2, 0 <= xi <= 2, xi, xi > 2. The main difference of our problem from the previous papers is that the terminal part of quality functional depends on the slow and fast variables and the controlled system is a more general form. We prove that in the case of a finite number of control change points, a power asymptotic expansion can be constructed for the initial vector of dual state lambda(epsilon) = (l(epsilon)*rho(epsilon)*)*, which determines the type of the optimal control.