\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2014 (2014), No. 26, pp. 1--16.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2014 Texas State University - San Marcos.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2014/26\hfil Relaxation in control systems]
{Relaxation in control systems of fractional semilinear evolution
equations}

\author[X. Liu, X. Fu \hfil EJDE-2014/26\hfilneg]
{Xiaoyou Liu, Xi Fu}  % in alphabetical order

\address{Xiaoyou Liu \newline
School of Mathematics and Physics, University of South China,
Hengyang 421001, Hunan Province, China}
\email{liuxiaoyou2002@hotmail.com Phone +86 13548962352}

\address{Xi Fu \newline
Department of Mathematics, Shaoxing University,
Shaoxing 312000, Zhejiang Province,  China}
\email{fuxi1984@hotmail.com}


\thanks{Submitted October 15, 2013. Published January 14, 2014.}
\subjclass[2000]{34G20, 93C20, 34A08}
\keywords{Fractional semilinear evolution equation;
relaxation property; \hfill\break\indent feedback control; nonconvex constraint;
mild solution}

\begin{abstract}
 We consider a control system described by fractional semilinear evolution
 equations with a mixed multivalued control constraint whose values are
 nonconvex  closed sets. Along with the original system, we consider
 the system in which the constraint on the control is the closed convex
 hull of the original  constraint. We obtain existence results for the
 control systems and  study relations between the solution sets
 of the two systems. An example is given to illustrate the abstract results.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{definition}[theorem]{Definition}
\allowdisplaybreaks

\section{Introduction}

Let $J=[0,b]$ and $0<\alpha<1$. In this paper, we consider a control
system described by fractional semilinear evolution equations of the form
\begin{equation}\label{mainequation}
\begin{gathered}
^CD_t^{\alpha}x(t)=Ax(t)+h(t,x(t))+g(t)u(t),\quad t\in J,\\
x(0)=x_0,
\end{gathered}
\end{equation}
with the mixed nonconvex constraint on the control
\begin{equation}\label{ucontrol}
u(t)\in U(t,x(t))\quad\text{a.e. on }J,
\end{equation}
where $^CD_t^{\alpha}$ is the Caputo fractional derivative of order $\alpha$,
$A$ is the infinitesimal generator of a strongly continuous semigroup
$\{T(t),t\geq0\}$ in a separable reflexive Banach space $X$,
$g:J\to \mathcal{L}(Y,X)$ ($\mathcal{L}(Y,X)$ is the space of continuous
linear operators from $Y$ into $X$), $h:J\times X\to X$ is a nonlinear
function and $U:J\times X\to 2^Y\backslash\{\emptyset\}$ is a multivalued
map with closed values (not necessarily convex). The space $Y$ is a separable,
 reflexive Banach space modeling the control space.

Along with the constraint \eqref{ucontrol} on the control, we also consider
the constraint
\begin{equation}\label{convexucontrol}
u(t)\in\operatorname{\overline{co}}U(t,x(t))\quad\text{a.e. on }J
\end{equation}
on the control. Here $\operatorname{\overline{co}}$ stands for the closed
convex hull of a set.

The solutions to the control systems considered in this paper are in
the mild sense and the precise definition will be given in 
Definition \ref{solutionofintegralform} below.

We denote by $\mathcal{R}_U$, $\mathcal{T}r_U$
($\mathcal{R}_{\overline{\mathrm{co}}U}$,
$\mathcal{T}r_{\overline{\mathrm{co}}U}$) the sets of all solutions,
all admissible trajectories of the control system
\eqref{mainequation}, \eqref{ucontrol}
(with the control system \eqref{mainequation}, \eqref{convexucontrol}, respectively).

The main results obtained in this paper are that:
 $\mathcal{T}r_{\overline{\mathrm{co}}U}$ is a compact set in $C(J,X)$
 and the relaxation property
\begin{equation}\label{importantresults}
\mathcal{T}r_{\overline{\mathrm{co}}U}=\overline{\mathcal{T}r_U}
\end{equation}
holds, where the bar stands for the closure in $C(J,X)$.

Recently, fractional calculus and differential equations have been proved
to be valuable tools in the modeling of many phenomena in various fields
of science and engineering. We can find its numerous applications
in viscoelasticity, electrochemistry, control, porous media,
electromagnetic, etc., see \cite{gaul,glockle,hilfer} for example.
There has been a great deal of interest in the existence of solutions
of fractional differential equations. One can see the monographs of Kilbas
et al \cite{kilbas}, Miller et al \cite{miller}, the survey of
Agarwal et al \cite{agarwal1,agarwal2}, Liu et al \cite{liuzhenhai1,liuzhenhai2}
and the references therein.

Abstract fractional semilinear differential equations represent a class
of fractional partial differential equations. For the study of
their existence results, we can refer to
Zhou and Jiao \cite{yongzhou1,yongzhou2}, Wang and Zhou \cite{wangjinrong2}
and the references therein. For control systems governed by fractional
semilinear differential equations, many literatures were devoted to give
sufficient conditions for their (approximate) controllability and optimal
control theory. For instance, Kumar and Sukavanam \cite{surendra},
Sakthivel et al
 \cite{sakganantamc2013,sakganrenantcnsns2013,sakthivela1,sakthivela},
Ganesh et al \cite{gansakmahantjam2013} (approximate controllability).
Wang and Zhou \cite{wangjinrong1} (optimal control theory).

Relaxation property, such as \eqref{importantresults}, if true,
has important ramifications in control theory, since it implies that
every trajectory of the convexified (full) system can be approximated
in $C(J,X)$ norm, with arbitrary degree of accuracy, by trajectories
of the original system. There are many papers dealing with the verification
of the relaxation property for various classes of control systems,
for instance, Tolstonogov \cite{tolstonogov12} of control systems
of subdifferential type, Mig\'orski \cite{migorski1,migorski2},
Tolstonogov\cite{tolstonogov2}, Tolstonogov et al \cite{tolstonogov11},
 Denkowski et al \cite{Denkowski} (c.f. Section 7.4)
of nonlinear evolution inclusions or equations.

In this paper, we study the relaxation property for control systems
described by a class of fractional semilinear evolution equations.
Please note that the control systems studied here are closed-loop
systems (feedback control systems) while the ones considered in papers
related to this work cited above were concerned with open-loop systems.

The rest of the paper is organized as follows:
In section \ref{preliminaries}, we introduce some useful preliminaries
and give the assumptions on the data of our problems. Some auxiliary results
needed in the proof of the main results are given in
section \ref{uuxiliaryresults}. Section \ref{existenceresults} deals with
the existence of solutions for the control systems.
The main results are presented in section \ref{secmainresults}.
An example and some concluding remarks are given in
sections \ref{sectionexample}.


\section{Preliminaries and assumptions}
\label{preliminaries}

Let $J=[0,b]$ be the closed interval of the real line with the Lebesgue
measure $\mu$ and the $\sigma$-algebra $\Sigma$ of $\mu$ measurable sets.
The norm of the space $X$ (or $Y$) will be denoted by $\|\cdot\|_X$
(or $\|\cdot\|_Y$). We denote by $C(J,X)$ the space of all continuous
functions from $J$ into $X$ with the supnorm given by
$\|x\|_{C}=\sup_{t\in J}\|x(t)\|_X$ for $x\in C(J,X)$.
For any Banach space $V$, the symbol $\omega$-$V$ stands for $V$ equipped
with the weak $\sigma(V,V^*)$ topology. The same notation will be used
for subsets of $V$. In all other cases, we assume that $V$ and its subsets
are equipped with the strong (normed) topology.

We first recall the following known definitions from the theory of
fractional calculus. For more details, please see \cite{kilbas,miller}.

\begin{definition}\label{fractionalintegral} \rm
The fractional integral of order $\alpha$ with the lower limit zero for
a function $f$ is defined as
\[
I^{\alpha}f(t)=\frac{1}{\Gamma(\alpha)}
\int_0^t\frac{f(s)}{(t-s)^{1-\alpha}}ds,\quad t>0,\;\alpha>0,
\]
provided the right hand side is point-wise defined on $[0,\infty)$,
where $\Gamma(\cdot)$ is the gamma function.
\end{definition}

\begin{definition}\label{riemannliouville} \rm
The Riemann-Liouville derivative of order $\alpha$ with the lower limit
zero for a function $f$ is defined as
\[
^LD^{\alpha}f(t)=\frac{1}{\Gamma(n-\alpha)}\frac{d^n}{dt^n}
\int_0^t\frac{f(s)}{(t-s)^{\alpha+1-n}}ds,\quad t>0,\; n-1<\alpha<n.
\]
\end{definition}

\begin{definition} \rm
The Caputo derivative of order $\alpha$ with the lower limit zero for
a function $f$ is defined as
\[
^CD^{\alpha}f(t)=\, ^LD^{\alpha}\Big(f(t)
-\sum_{k=0}^{n-1}\frac{t^k}{k!}f^{(k)}(0)\Big),\quad t>0,\; n-1<\alpha<n.
\]
\end{definition}

If $f$ is an abstract function with values in $X$, then integrals which appear
in Definitions \ref{fractionalintegral} and \ref{riemannliouville}
are taken in Bochner's sense.

We now proceed to some basic definitions and results from multivalued analysis.
For more details on multivalued analysis, see the books \cite{aubin,hu}.

We use the following symbols: $P_{f}(Y)$ is the set of all nonempty closed
subsets of $Y$, $P_{bf}(Y)$ is the set of all nonempty, closed and bounded
subsets of $Y$.

On $P_{bf}(Y)$, we have a metric known as the ``Hausdorff metric" and defined by
\[
h(A,B)=\max\big\{\sup_{a\in A}d(a,B),\,\sup_{b\in B}d(b,A)\big\},
\]
where $d(x,C)$ is the distance from a point $x$ to a set $C$.
We say a multivalued map is $h$-continuous if it is continuous in
the Hausdorff metric $h(\cdot,\cdot)$.

We say that a multivalued map $F:J\to P_f(Y)$ is measurable if
$F^{-1}(E)=\{t\in J:F(t)\cap E\neq\emptyset\}\in \Sigma$ for every closed
set $E\subseteq Y$. If $F:J\times X\to P_f(Y)$, then the measurability
of $F$ means that $F^{-1}(E)\in\Sigma\otimes\mathcal{B}_{X}$,
where $\Sigma\otimes\mathcal{B}_{X}$ is the $\sigma$-algebra of subsets in
 $J\times X$ generated by the sets $A\times B$, $A\in\Sigma$,
 $B\in\mathcal{B}_{X}$, and $\mathcal{B}_{X}$ is the $\sigma$-algebra of
the Borel sets in $X$.

Suppose $V$, $Z$ are two Hausdorff topological spaces and
$F: V\to 2^Z \backslash\{\emptyset\}$. We say that $F$ is lower semicontinuous
in the sense of Vietoris (l.s.c. for short) at a point $x_0\in V$,
if for any open set $W\subseteq Z$, $F(x_0)\cap W\neq\emptyset$, there
is a neighborhood $O(x_0)$ of $x_0$ such that $F(x)\cap W\neq\emptyset$
for all $x\in O(x_0)$. $F$ is said to be upper semicontinuous in the sense
of Vietoris (u.s.c. for short) at a point $x_0\in V$, if for any open set
 $W\subseteq Z$, $F(x_0)\subseteq W$, there is a neighborhood $O(x_0)$ of $x_0$
such that $F(x)\subseteq W$ for all $x\in O(x_0)$. For the properties of
l.s.c and u.s.c, please refer to the book \cite{hu}.

Besides the standard norm on $L^q(J,Y)$ (here $Y$ is a separable, reflexive
Banach space ), $1<q<\infty$, we also consider the so called weak norm
\begin{equation}\label{weaknorm}
\|u(\cdot)\|_{\omega}=\sup_{0\leq t_1\leq t_2\leq b}
\big\|\int_{t_1}^{t_2}u(s)ds\big\|_Y, \text{ for }u\in L^q(J,Y).
\end{equation}
The space $L^q(J,Y)$ furnished with this norm will be denoted by
$L_{\omega}^q(J,Y)$. The following result establishes a relation between
convergence in $\omega$-$L^q(J,Y)$ and convergence in $L_{\omega}^q(J,Y)$.

\begin{lemma}[\cite{tolstonogov2}] \label{weaknormforl2ty}
If a sequence $\{u_n\}_{n\geq1}\subseteq L^q(J,Y)$ is bounded and converges
to $u$ in $L_{\omega}^q(J,Y)$, then it converges to $u$ in $\omega$-$L^q(J,Y)$.
\end{lemma}

We assume the following assumptions on the data of our problems in the
whole paper.
\begin{itemize}
\item[(H1)]: The operator $A$ generates a strongly continuous semigroup
$T(t)$, $t\geq0$ in $X$, and there exists a constant $M_A\geq1$ such
that $\sup_{t\in[0,\infty)}\|T(t)\|\leq M_A$. For any $t>0$,
$T(t)$ is compact.

\item[(H2)] The operator $g:J\to\mathcal{L}(Y,X)$ is such that:
\begin{itemize}
\item[(1)] the map $t\to g(t)u$ is measurable for any $u\in Y$;
\item[(2)] for a.e. $t\in J$,
\begin{equation}
\|g(t)\|_{\mathcal{L}(Y,X)}\leq d, \text{ with }d>0.
\end{equation}
\end{itemize}

\item[(H3)] The function $h:J\times X\to X$ satisfies the following:
\begin{itemize}
\item[(1)] $t\to h(t,x)$ is measurable for all $x\in X$;
\item[(2)] there exists a function $l\in L^{\infty}(J,\mathbb{R}^+)$
such that for a.e. $t\in J$ and all $x$, $y\in X$,
\begin{equation}
\|h(t,x)-h(t,y)\|_X\leq l(t)\|x-y\|_X;
\end{equation}
\item[(3)] there exists a constant $0<\beta<\alpha$ such that for a.e.
$t\in J$, and all $x\in X$, $\|h(t,x)\|_X\leq a_1(t)+c_1\|x\|_X$,
where $a_1\in L^{1/\beta}(J,\mathbb{R}^+)$ and $c_1>0$.
\end{itemize}

\item[(H4)] The multivalued map $U:J\times X\to P_{f}(Y)$ is such that:
\begin{itemize}
\item[(1)] for all $x\in X$, $t\to U(t,x)$ is measurable;
\item[(2)] $h(U(t,x),U(t,y))\leq k_1(t)\|x-y\|_X$ a.e. on $J$,
 with $k_1\in L^{\infty}(J,\mathbb{R}^+)$;
\item[(3)] for a.e. $t\in J$, and all $x\in X$,
$\|U(t,x)\|_Y=\sup\{\|v\|_Y:v\in U(t,x)\}\leq a_2(t)+c_2\|x\|_X$,
 where $a_2\in L^{1/\beta}(J,\mathbb{R}^+)$ and $c_2>0$.
\end{itemize}
\end{itemize}

\begin{definition}[\cite{yongzhou1,yongzhou2}]\label{solutionofintegralform} \rm
A pair of functions $(x,u)$ is a solution (mild solution) of the control
system \eqref{mainequation}, \eqref{ucontrol}, if $x(0)=x_0$,
$x\in C(J,X)$ and there exists $u\in L^1(J,Y)$ such that $u(t)\in U(t,x(t))$ a.e.
 $t\in J$ and
\begin{equation}
x(t)=P_{\alpha}(t)x_0+\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)
\big(g(s)u(s)+h(s,x(s))\big)ds.
\end{equation}
\end{definition}

A similar definition can be introduced for the system \eqref{mainequation},
\eqref{convexucontrol}. Here
\begin{gather*}
P_{\alpha}(t)=\int_0^{\infty}\xi_{\alpha}(\theta)T(t^{\alpha}\theta)d\theta,\quad
Q_{\alpha}(t)=\alpha\int_0^{\infty}\theta\xi_{\alpha}(\theta)T(t^{\alpha}\theta)
 d\theta,
\\
\xi_{\alpha}(\theta)=\frac{1}{\alpha}\theta^{-1-\frac{1}{\alpha}}
 \varpi_{\alpha}\big(\theta^{-\frac{1}{\alpha}}\big)\geq0,
\\
\varpi_{\alpha}(\theta)=\frac{1}{\pi}\sum_{n=1}^{\infty}(-1)^{n-1}
 \theta^{-n\alpha-1} \frac{\Gamma(n\alpha+1)}{n!}\sin(n\pi\alpha),\quad
 \theta\in(0,\infty),
\end{gather*}
and $\xi_{\alpha}$ is a probability density function defined on $(0,\infty)$;
that is,
\[
\xi_{\alpha}(\theta)\geq0,\quad \theta\in(0,\infty),\quad
\int_0^{\infty}\xi_{\alpha}(\theta)d\theta=1.
\]
It is not difficult to verify that
\begin{equation}\label{asimplecalculation}
\int_0^{\infty}\theta\xi_{\alpha}(\theta)d\theta=\frac{1}{\Gamma(1+\alpha)}.
\end{equation}

\begin{lemma}[\cite{yongzhou1,yongzhou2}]\label{propertyforpq}
Let {\rm (H1)} hold. Then the operators $P_{\alpha}$ and $Q_{\alpha}$
have the following properties:
\begin{itemize}
\item[(1)] For any fixed $t\geq0$, $P_{\alpha}(t)$ and $Q_{\alpha}(t)$
are linear and bounded operators, i.e., for any $x\in X$,
\[
\|P_{\alpha}(t)x\|_X\leq M_A\|x\|_X,\quad \|Q_{\alpha}(t)x\|_X
\leq\frac{\alpha M_A}{\Gamma(1+\alpha)}\|x\|_X;
\]
\item[(2)] $\{P_{\alpha}(t),t\geq0\}$ and $\{Q_{\alpha}(t),t\geq0\}$
 are strongly continuous;
\item[(3)] For every $t>0$, $P_{\alpha}(t)$ and $Q_{\alpha}(t)$
are compact operators.
\end{itemize}
\end{lemma}

The proof of the above lemma can be found in \cite{yongzhou1}.


\section{Auxiliary results}\label{uuxiliaryresults}

In this section, we shall give some auxiliary results needed in the proof
of the main results. We begin with the a prior estimation of the trajectory
of the control systems.

\begin{lemma}\label{boundednessforsolution}
For any admissible trajectory $x$ of the control system
\eqref{mainequation}, \eqref{convexucontrol};
 i.e., $x\in\mathcal{T}r_{\overline{\mathrm{co}}U}$, there is a constant $L$
such that
\begin{equation}\label{boundedfortraj}
\|x\|_C\leq L.
\end{equation}
\end{lemma}

\begin{proof}
Let any $x\in\mathcal{T}r_{\overline{\mathrm{co}}U}$.
From Definition \ref{solutionofintegralform}, we have that there exists
a $u(t)\in\operatorname{\overline{co}}U(t,x(t))$ a.e. $t\in J$ and
\[
x(t)=P_{\alpha}(t)x_0+\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)\big(g(s)u(s)
+h(s,x(s))\big)ds.
\]
Then by Lemma \ref{propertyforpq}, we obtain
\begin{equation}\label{lemboundedforx1}
\begin{aligned}
\|x(t)\|_X
&\leq M_A\|x_0\|_X+\frac{\alpha M_A}{\Gamma(1+\alpha)}
\int_0^t(t-s)^{\alpha-1}\|h(s,x(s))\|_Xds \\
&\quad +\frac{\alpha M_A}{\Gamma(1+\alpha)}
\int_0^t(t-s)^{\alpha-1}\|g(s)u(s)\|_Xds.
\end{aligned}
\end{equation}
From (H3)(2), (H3)(3) and the H\"older inequality, we have
\begin{equation}\label{lemboundedforx2}
\begin{aligned}
&\int_0^t(t-s)^{\alpha-1}\|h(s,x(s))\|_Xds \\
&\leq \int_0^t(t-s)^{\alpha-1}\|h(s,x(s))-h(s,0)\|_Xds
 +\int_0^t(t-s)^{\alpha-1}\|h(s,0)\|_Xds \\
&\leq \int_0^t(t-s)^{\alpha-1}l(s)\|x(s)\|_Xds
 +\int_0^t(t-s)^{\alpha-1}a_1(s)ds \\
&\leq \Big[\frac{1-\beta}{\alpha-\beta}b^{\frac{\alpha-\beta}{1-\beta}}
 \Big]^{1-\beta}\|a_1\|_{L^{1/\beta}(J)}
+\|l\|_{L^{\infty}(J)}\int_0^t(t-s)^{\alpha-1}\|x(s)\|_Xds.
\end{aligned}
\end{equation}
Similarly, by (H2)(2) and (H4)(3), we obtain
\begin{equation}\label{lemboundedforx3}
\begin{aligned}
&\int_0^t(t-s)^{\alpha-1}\|g(s)u(s)\|_Xds \\
&\leq d\int_0^t(t-s)^{\alpha-1}\big(a_2(s)+c_2\|x(s)\|_X\big)ds \\
&\leq d\Big[\frac{1-\beta}{\alpha-\beta}b^{\frac{\alpha-\beta}{1-\beta}}
 \Big]^{1-\beta}\|a_2\|_{L^{1/\beta}(J)}
 +dc_2\int_0^t(t-s)^{\alpha-1}\|x(s)\|_Xds.
\end{aligned}
\end{equation}
Combining \eqref{lemboundedforx2}, \eqref{lemboundedforx3} with
\eqref{lemboundedforx1}, we obtain
\begin{align*}
\|x(t)\|_X
&\leq M_A\|x_0\|_X+\frac{\alpha M_A}{\Gamma(1+\alpha)}
\big(dc_2+\|l\|_{L^{\infty}(J)}\big)\int_0^t(t-s)^{\alpha-1}\|x(s)\|_Xds \\
&\quad +\frac{\alpha M_A}{\Gamma(1+\alpha)}
\Big[\frac{1-\beta}{\alpha-\beta}b^{\frac{\alpha-\beta}{1-\beta}}\Big]^{1-\beta}
\Big(\|a_1\|_{L^{1/\beta}(J)}+d\|a_2\|_{L^{1/\beta}(J)}\Big).
\end{align*}
From the above inequality, using the well-known singular-version Gronwall
inequality (see \cite[Theorem 3.1]{dixon}), we can deduce that there
exists a constant $L>0$ such that $\|x\|_C\leq L$.
\end{proof}

Let $\operatorname{pr}_L:X\to X$ be the L-radial retraction; i.e.,
\[
\operatorname{pr}_L(x)=\begin{cases}
x, & \|x\|_X\leq L,\\
\frac{Lx}{\|x\|_X}, & \|x\|_X>L.
\end{cases}
\]
This map is Lipschitz continuous. We define $U_1(t,x)=U(t,\operatorname{pr}_Lx)$.
Evidently, $U_1$ satisfies (H4)(1) and (H4)(2). Moreover, by the properties
of $\operatorname{pr}_L$, we have, for a.e. $t\in J$, all $x\in X$ and
all $u\in U_1(t,x)$ such that
\[
\|u\|_Y\leq a_2(t)+c_2L\,\,\textrm {and}\,\,\|u\|_Y\leq a_2(t)+c_2\|x\|_X.
\]
Hence, Lemma \ref{boundednessforsolution} is still valid with $U(t,x)$
substituted by $U_1(t,x)$. Consequently, henceforth we assume without any
loss of generality that, for a.e. $t\in J$ and all $x\in X$,
\begin{equation}\label{boundedforu}
\sup\{\|v\|_Y:v\in U(t,x)\}\leq\varphi(t)=a_2(t)+c_2L,\quad
\text{with }\varphi\in L^{1/\beta}(J,\mathbb{R}^+).
\end{equation}
Let $\varphi$ be defined by \eqref{boundedforu}, we put
\begin{gather}\label{boundedcontrolspace}
Y_{\varphi}=\{u\in L^{1/\beta}(J,Y):\|u(t)\|_Y\leq\varphi(t)\text{ a.e. }t\in J\},\\
\label{trajectoryrhdspace}
X_{\varphi}=\{f\in L^{1/\beta}(J,X):\|f(t)\|_X\leq d\varphi(t)+a_1(t)+c_1L\text{ a.e. }t\in J\}.
\end{gather}

In accordance with (H2) and (H3), for any $x\in C(J,X)$ and
$u\in L^{1/\beta}(J,Y)$, the function $t\to g(t)u(t)+h(t,x(t))$ is an element
of the space $L^{1/\beta}(J,X)$. Hence, we can consider an operator
$\mathcal{A}:C(J,X)\times L^{1/\beta}(J,Y)\to L^{1/\beta}(J,X)$ defined by
\begin{equation}\label{definofoperatorA}
\mathcal{A}(x,u)(t)=g(t)u(t)+h(t,x(t)).
\end{equation}

\begin{lemma}\label{propertyofoperatorA}
The map $(x,u)\to\mathcal{A}(x,u)$ is sequentially continuous from
$C(J,X)\times\omega$-$L^{1/\beta}(J,Y)$ into $\omega$-$L^{1/\beta}(J,X)$.
\end{lemma}

\begin{proof}
Suppose that $x_n\to x$ in $C(J,X)$ and $u_n\to u$ in $\omega$-$L^{1/\beta}(J,Y)$.
Let any $h\in L^{1/(1-\beta)}(J,X^*)$ be fixed.
Now we may assume that $\|x_n\|_C\leq M$ for some constant $M>0$ and $n\geq1$.
Then from (H2) and (H3), we can have the following facts
\begin{gather}\label{lemproperofoperA1}
h(t,x_n(t))\to h(t,x(t))\quad \text{in $X$ a.e. }t\in J, \\
\label{lemproperofoperA2}
\|h(t,x_n(t))\|_X\leq a_1(t)+c_1M, \\
\label{lemproperofoperA3}
\int_J\langle g^*(t)h(t),u_n(t)\rangle dt\to\int_J\langle g^*(t)h(t),u(t)
\rangle dt,
\end{gather}
where $g^*(t)$ is the operator adjoint to $g(t)$.
From \eqref{lemproperofoperA1} and \eqref{lemproperofoperA2},
using Lebesgue's dominated convergence theorem, we obtain
\begin{equation}\label{lemproperofoperA4}
h(t,x_n(t))\to h(t,x(t))\quad \text{in }L^{1/\beta}(J,X).
\end{equation}
Since $\langle h(t),g(t)u(t)\rangle=\langle g^*(t)h(t),u(t)\rangle$
and $h\in L^{1/(1-\beta)}(J,X^*)$ is arbitrary, by \eqref{lemproperofoperA3},
we deduce that
\[
g(t)u_n(t)\to g(t)u(t)\quad \text{in }\omega\text{-}L^{1/\beta}(J,X).
\]
This together and \eqref{lemproperofoperA4} imply
\[
\mathcal{A}(x_n,u_n)\to\mathcal{A}(x,u)\quad \text{in }\omega\text{-}L^{1/\beta}(J,X).
\]
The lemma is proved.
\end{proof}

Now we consider the  auxiliary problem:
\begin{equation}\label{middleequation}
\begin{gathered}
^CD_t^{\alpha}x(t)=Ax(t)+f(t),\,\,t\in J=[0,b],\\
x(0)=x_0.
\end{gathered}
\end{equation}
It is clear that, for every $f\in L^{1/\beta}(J,X)$, equation
\eqref{middleequation} has a unique mild solution $S(f)\in C(J,X)$ which
is given by
\[
S(f)(t)=P_{\alpha}(t)x_0+\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)f(s)ds.
\]

The following lemma concerns with the property of the solution map
$S$ which is crucial in our investigation.

\begin{lemma}\label{continuousofsoluoperator}
The solution map $S:X_{\varphi}\to C(J,X)$ is continuous from
$\omega$-$X_{\varphi}$ into $C(J,X)$.
\end{lemma}

\begin{proof}
Consider the operator $H:L^{1/\beta}(J,X)\to C(J,X)$ defined by
\[
H(f)(t)=\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)f(s)ds.
\]
We know $H$ is linear. From simple calculation, one has
\begin{equation}\label{propertyofoperatorh}
\|H(f)\|_C\leq\frac{\alpha M_A}{\Gamma(1+\alpha)}
\Big[\frac{1-\beta}{\alpha-\beta}b^{\frac{\alpha-\beta}{1-\beta}}\Big]^{1-\beta}\|f\|_{L^{1/\beta}(J,X)};
\end{equation}
i.e., the operator $H$ is continuous from $L^{1/\beta}(J,X)$ to $C(J,X)$,
hence $H$ is also continuous from $\omega$-$L^{1/\beta}(J,X)$ to
$\omega$-$C(J,X)$.

Let  $C\in P_{b}(L^{1/\beta}(J,X))$ and suppose that for any
$f\in C$, $\|f\|_{L^{1/\beta}(J,X)}\leq K$ ($K>0$ is a constant).
Next we will show that $H$ is completely continuous.

(a) From \eqref{propertyofoperatorh}, we know that $\|H(f)(t)\|_X$ is
uniformly bounded for any $t\in J$ and $f\in C$.

(b) $H$ is equicontinuous on $C$. Let $0\leq t_1<t_2\leq b$.
For any $f\in C$, we obtain
\begin{align*}
&\|H(f)(t_2)-H(f)(t_1)\|_X \\
&=\big\|\int_0^{t_2}(t_2-s)^{\alpha-1}Q_{\alpha}(t_2-s)f(s)ds-
\int_0^{t_1}(t_1-s)^{\alpha-1}Q_{\alpha}(t_1-s)f(s)ds\big\|_X \\
&\leq \big\|\int_{t_1}^{t_2}(t_2-s)^{\alpha-1}Q_{\alpha}(t_2-s)f(s)ds
\big\|_X \\
&\quad +\big\|\int_0^{t_1}\Big((t_2-s)^{\alpha-1}-(t_1-s)^{\alpha-1}\Big)
 Q_{\alpha}(t_2-s)f(s)ds\big\|_X  \\
&\quad +\big\|\int_0^{t_1}(t_1-s)^{\alpha-1}\Big(Q_{\alpha}(t_2-s)
 -Q_{\alpha}(t_1-s)\Big)f(s)ds\big\|_X\\
&=: I_1+I_2+I_3.
\end{align*}
By using analogous arguments as in Lemma \ref{boundednessforsolution}, 
we  have
\begin{gather*}
I_1\leq \frac{\alpha M_A}{\Gamma(1+\alpha)}
\Big[\frac{1-\beta}{\alpha-\beta}\Big]^{1-\beta}K(t_2-t_1)^{\alpha-\beta},
 \\
\begin{aligned}
I_2&\leq \frac{\alpha M_A}{\Gamma(1+\alpha)}\Big(\int_0^{t_1}
\big((t_1-s)^{\alpha-1}-(t_2-s)^{\alpha-1}\big)^{1/(1-\beta)}ds\Big)^{1-\beta}K
 \\
&\leq \frac{\alpha M_A}{\Gamma(1+\alpha)}\Big(\int_0^{t_1}
 \big((t_1-s)^{\frac{\alpha-1}{1-\beta}}-(t_2-s)^{\frac{\alpha-1}{1-\beta}}\big)ds\Big)^{1-\beta}K
 \\
&= \frac{\alpha M_A}{\Gamma(1+\alpha)}
 \Big[\frac{1-\beta}{\alpha-\beta}\Big]^{1-\beta}
 \Big(t_1^{\frac{\alpha-\beta}{1-\beta}}
 -t_2^{\frac{\alpha-\beta}{1-\beta}}+(t_2-t_1)^{\frac{\alpha-\beta}{1-\beta}}
 \Big)^{1-\beta}K \\
&\leq \frac{2\alpha M_A}{\Gamma(1+\alpha)}
 \Big[\frac{1-\beta}{\alpha-\beta}\Big]^{1-\beta}
 \big(t_2-t_1\big)^{\alpha-\beta}K.
\end{aligned}
\end{gather*}
For $t_1=0$, $0<t_2\leq b$, it is easy to see that $I_3=0$. For $t_1>0$ and 
$\epsilon>0$ be small enough, we have
\begin{align*}
I_3 &\leq \big\|\int_0^{t_1-\epsilon}(t_1-s)^{\alpha-1}\Big(Q_{\alpha}(t_2-s)
 -Q_{\alpha}(t_1-s)\Big)f(s)ds\big\|_X \\
&\quad +\big\|\int_{t_1-\epsilon}^{t_1}(t_1-s)^{\alpha-1}\Big(Q_{\alpha}(t_2-s)
 -Q_{\alpha}(t_1-s)\Big)f(s)ds\big\|_X \\
&\leq \sup_{s\in[0,t_1-\epsilon]}\|Q_{\alpha}(t_2-s)-Q_{\alpha}(t_1-s)\|\Big[\frac{1-\beta}{\alpha-\beta}\Big]^{1-\beta}
\big(t_1^{\frac{\alpha-\beta}{1-\beta}}-\epsilon^{\frac{\alpha-\beta}{1-\beta}}\big)^{1-\beta}K \\
&\quad +\frac{2\alpha M_A}{\Gamma(1+\alpha)}
\Big[\frac{1-\beta}{\alpha-\beta}\Big]^{1-\beta}\epsilon^{\alpha-\beta}K.
\end{align*}
Combining the estimations for $I_1$, $I_2$, $I_3$, and letting  
$t_2\to t_1$ and $\epsilon\to 0$ in $I_3$, we know that $H$ is equicontinuous. 
For more details, please see \cite{yongzhou2}.

(c) The set $\Pi(t)=\{H(f)(t):f\in C\}$ is relatively compact in $X$. 
Clearly, $\Pi(0)=\{0\}$ is compact, and hence, it is only necessary to 
consider $t>0$. For each $h\in(0,t)$, $t\in(0,b]$, $f\in C$, and $\delta>0$ 
being arbitrary, we define
\[
\Pi_{h,\delta}(t)=\{H_{h,\delta}(f)(t):f\in C\},
\]
where
\begin{align*}
H_{h,\delta}(f)(t)
&= \alpha\int_0^{t-h}\int_{\delta}^{\infty}\theta(t-s)^{\alpha-1}\xi_{\alpha}
(\theta) T((t-s)^{\alpha}\theta)f(s)d\theta ds \\
&= \alpha\int_0^{t-h}\int_{\delta}^{\infty}\theta(t-s)^{\alpha-1}
\xi_{\alpha}(\theta)
T(h^{\alpha}\delta)T\big((t-s)^{\alpha}\theta-h^{\alpha}\delta\big)f(s)d\theta ds
  \\
&= \alpha T(h^{\alpha}\delta)\int_0^{t-h}\int_{\delta}^{\infty}
 \theta(t-s)^{\alpha-1}\xi_{\alpha}(\theta)
T\big((t-s)^{\alpha}\theta-h^{\alpha}\delta\big)f(s)d\theta ds.
\end{align*}
From the compactness of $T(h^{\alpha}\delta)$ ($h^{\alpha}\delta>0$), we obtain 
that the set $\Pi_{h,\delta}(t)$ is relatively compact in $X$ for any 
$h\in(0,t)$ and $\delta>0$. Moreover, we have
\begin{align*}
&\|H(f)(t)-H_{h,\delta}(f)(t)\|_X \\
&= \alpha\Big\|\int_0^t\int_0^{\delta}\theta(t-s)^{\alpha-1}
 \xi_{\alpha}(\theta) T((t-s)^{\alpha}\theta)f(s)d\theta ds \\
&\quad +\int_0^t\int_{\delta}^{\infty}\theta(t-s)^{\alpha-1}\xi_{\alpha}(\theta)
 T((t-s)^{\alpha}\theta)f(s)d\theta ds \\
&\quad -\int_0^{t-h}\int_{\delta}^{\infty}\theta(t-s)^{\alpha-1}
 \xi_{\alpha}(\theta) T((t-s)^{\alpha}\theta)f(s)d\theta ds\Big\|_X 
\\
&\leq \alpha\Big\|\int_0^t\int_0^{\delta}\theta(t-s)^{\alpha-1}\xi_{\alpha}
 (\theta) T((t-s)^{\alpha}\theta)f(s)d\theta ds\Big\|_X \\
&\quad +\alpha\Big\|\int_{t-h}^t\int_{\delta}^{\infty}\theta(t-s)^{\alpha-1}
 \xi_{\alpha}(\theta) T((t-s)^{\alpha}\theta)f(s)d\theta ds\Big\|_X \\
&\leq M_A\alpha\Big(\int_0^t(t-s)^{\frac{\alpha-1}{1-\beta}}ds\Big)^{1-\beta}
 \|f\|_{L^{1/\beta}(J,X)} \int_0^{\delta}\theta\xi_{\alpha}(\theta)d\theta \\
&\quad +M_A\alpha\Big(\int_{t-h}^t(t-s)^{\frac{\alpha-1}{1-\beta}}ds
 \Big)^{1-\beta}\|f\|_{L^{1/\beta}(J,X)}
\int_{\delta}^{\infty}\theta\xi_{\alpha}(\theta)d\theta \\
&\leq M_AK\alpha\Big[\frac{1-\beta}{\alpha-\beta}\Big]^{1-\beta}
\Big(b^{\alpha-\beta}\int_0^{\delta}\theta\xi_{\alpha}(\theta)d\theta
 +\frac{1}{\Gamma(1+\alpha)}h^{\alpha-\beta}\Big).
\end{align*}
By \eqref{asimplecalculation}, the last term of the preceding inequality 
tends to zero as $h\to 0$ and $\delta\to 0$. Therefore, there are relatively 
compact sets arbitrarily close to the set $\Pi(t)$, $t>0$. Hence the set 
$\Pi(t)$, $t>0$ is also relatively compact in $X$.

Since $X_{\varphi}$ is a convex compact metrizable subset of
 $\omega$-$L^{1/\beta}(J,X)$, it suffices to prove the sequential continuity
 of the map $S$. Now let $\{f_n\}_{n\geq1}\subseteq X_{\varphi}$ such that
\begin{equation}\label{lemsolution1}
f_n\to f\quad \text{in $\omega\text{-}L^{1/\beta}(J,X)$, $f\in X_{\varphi}$}.
\end{equation}
By the property of the operator $H$, we have $H(f_n)\to H(f)$ in 
$\omega$-$C(J,X)$. Since $\{f_n\}_{n\geq1}$ is bounded, there is a 
subsequence $\{f_{n_k}\}_{k\geq1}$ of the sequence $\{f_n\}_{n\geq1}$ such 
that $H(f_{n_k})\to z$ in $C(J,X)$ for some $z\in C(J,X)$. From the facts that
\[
H(f_n)\to H(f)\quad\text{in }\omega\text{-}C(J,X),\quad \text{and}\quad
H(f_{n_k})\to z\text{ in }C(J,X),
\]
we obtain that $z=H(f)$ and $H(f_n)\to H(f)$ in $C(J,X)$.

By the definitions of the operators $S$ and $H$, we have that 
$S(f)(t)=P_{\alpha}(t)x_0+H(f)(t)$. Then due to the arguments above, 
we have $S(f_n)\to S(f)$ in $C(J,X)$. This completes the proof of the lemma.
\end{proof}


\section{Existence results for control systems}
\label{existenceresults}

In this section, we shall prove the existence of solutions for the control systems 
\eqref{mainequation}, \eqref{ucontrol} and 
\eqref{mainequation}, \eqref{convexucontrol}.

Let $\Lambda=S(X_\varphi)$. From Lemma \ref{continuousofsoluoperator}, we 
have $\Lambda$ is a compact subset of $C(J,X)$. 
It follows from \eqref{boundedforu} and \eqref{trajectoryrhdspace} 
that $\mathcal{T}r_U\subseteq\mathcal{T}
r_{\overline{\mathrm{co}}U}\subseteq\Lambda$. 
Let $\overline{U}:C(J,X)\to 2^{L^{1/\beta}(J,Y)}$ be defined by
\begin{equation}\label{uforsolutions}
\overline{U}(x)=\{h:J\to Y\text{ measurable}:h(t)\in U(t,x(t))
\text{ a.e.}\},\,x\in C(J,X).
\end{equation}

\begin{theorem}\label{existenceforcontrolsys}
The set $\mathcal{R}_U$ is nonempty and the set 
$\mathcal{R}_{\overline{\mathrm{co}}U}$ is a compact subset of the space 
$C(J,X)\times\omega$-$L^{1/\beta}(J,Y)$.
\end{theorem}

\begin{proof}
By the hypotheses (H4)(1) and (H4)(2), we have that for any measurable 
function $x:J\to X$, the map $t\to U(t,x(t))$ is measurable and has closed 
values \cite[Proposition 2.7.9]{hu}. 
Therefore it has measurable selectors \cite{Himmelberg}. 
So the operator $\overline{U}$ is well defined and its values are closed 
decomposable subsets of $L^{1/\beta}(J,Y)$. We claim that $x\to\overline{U}(x)$
is l.s.c. Let $x_*\in C(J,X)$, $h_*\in\overline{U}(x_*)$ and let 
$\{x_n\}_{n\geq1}\subseteq C(J,X)$ be a sequence converging to $x_*$. 
It follows from \cite[Lemma 3.2]{zhu1} that there is a sequence 
$h_n\in\overline{U}(x_n)$ such that
\begin{equation}\label{thmexistforconsys1}
\|h_*(t)-h_n(t)\|_Y\leq d_Y(h_*(t),U(t,x_n(t)))+\frac{1}{n},\quad
\text{a.e. }t\in J.
\end{equation}
Since the map $y\to U(t,y)$ is $h$-continuous a.e. $t\in J$ ((H4)(2)), 
then for a.e. $t\in J$, the map $y\to U(t,y)$ is l.s.c. 
\cite[Proposition 1.2.66]{hu}. Hence by Proposition 1.2.26 in \cite{hu}, 
the function $y\to d_Y(h_*(t),U(t,y))$ is u.s.c. for a.e. $t\in J$. 
It follows from \eqref{thmexistforconsys1} that, for a.e. $t\in J$,
\begin{align*}
\lim_{n\to\infty}\|h_*(t)-h_n(t)\|_Y
&\leq \limsup_{n\to\infty}d_Y(h_*(t),U(t,x_n(t)))\\
&\leq  d_Y(h_*(t),U(t,x_*(t)))=0.
\end{align*}
This together with \eqref{boundedforu} implies that $h_n\to h_*$ in 
$L^{1/\beta}(J,Y)$. Therefore the map $x\to\overline{U}(x)$ is l.s.c.
By \cite[Proposition 2.2]{tolstonogov4} (also see \cite[Theorem 2.8.7]{hu}),
 there exists a continuous function $m:\Lambda\to L^{1/\beta}(J,Y)$ such that
\begin{equation}\label{thmexistforconsys2}
m(x)\in\overline{U}(x),\quad \text{for all }x\in \Lambda.
\end{equation}
Consider the map $\mathcal{P}:L^{1/\beta}(J,X)\to L^{1/\beta}(J,Y)$
defined by $\mathcal{P}(f)=m(S(f))$. 
Thanks to Lemma \ref{continuousofsoluoperator} and the continuity of $m$, 
the map $\mathcal{P}$ is continuous from $\omega$-$X_{\varphi}$ into
 $L^{1/\beta}(J,Y)$. Then by Lemma \ref{propertyofoperatorA}, we deduce
that the map $f\to\mathcal{A}(S(f),\mathcal{P}(f))$ is continuous from 
$\omega$-$X_{\varphi}$ into $\omega$-$L^{1/\beta}(J,X)$. It follows
from \eqref{boundedforu}, \eqref{trajectoryrhdspace} and 
\eqref{definofoperatorA} that $\mathcal{A}(S(f),\mathcal{P}(f))\in X_{\varphi}$ 
for every $f\in X_{\varphi}$. Therefore, the map 
$f\to\mathcal{A}(S(f),\mathcal{P}(f))$ is continuous from $\omega$-$X_{\varphi}$ 
into $\omega$-$X_{\varphi}$. Since $\omega$-$X_{\varphi}$ is a convex metrizable 
compact set in $\omega$-$L^{1/\beta}(J,X)$, Schauder's fixed point theorem implies
that this map has a fixed point $f_*\in X_{\varphi}$; i.e., 
$f_*=\mathcal{A}(S(f_*),\mathcal{P}(f_*))$. Let $u_*=\mathcal{P}(f_*)$ and 
$x_*=S(f_*)$, then we have $u_*=m(x_*)$ and $f_*=\mathcal{A}(x_*,u_*)$. 
That is to say we have
\begin{gather*}
x_*(t)=S(f_*)(t)=P_{\alpha}(t)x_0+\int_0^t(t-s)^{\alpha-1}
Q_{\alpha}(t-s)\big(g(s)u_*(s)+h(s,x_*(s))\big)ds, \\
u_*(t)\in U(t,x_*(t))\,\,\text{a.e.}\,\,t\in J.
\end{gather*}
These imply that $(x_*(\cdot),u_*(\cdot))$ is a solution of the control 
system \eqref{mainequation}, \eqref{ucontrol}. Hence $\mathcal{R}_U$ 
is nonempty.

It is easy to see that 
$\mathcal{R}_{\overline{\mathrm{co}}U}\subseteq \Lambda\times Y_{\varphi}$.
 Since $\Lambda$ is compact in $C(J,X)$ and $Y_{\varphi}$ is metrizable convex 
compact in $\omega$-$L^{1/\beta}(J,Y)$, we have that
$\mathcal{R}_{\overline{\mathrm{co}}U}$ is relatively compact in 
$C(J,X)\times\omega$-$L^{1/\beta}(J,Y)$. Hence to complete the proof of this
theorem, it is sufficient to prove that $\mathcal{R}_{\overline{\mathrm{co}}U}$ 
is sequentially closed in $C(J,X)\times\omega$-$L^{1/\beta}(J,Y)$.

Let $\{(x_n(\cdot),u_n(\cdot))\}_{n\geq1}\subseteq
\mathcal{R}_{\overline{\mathrm{co}}U}$ be a sequence converging to
 $(x(\cdot),u(\cdot))$ in the space $C(J,X)\times\omega$-$L^{1/\beta}(J,Y)$. Denote
\begin{gather*}
f_n(t)=g(t)u_n(t)+h(t,x_n(t)), \\
f(t)=g(t)u(t)+h(t,x(t)).
\end{gather*}
According to Lemma \ref{propertyofoperatorA}, $f_n\to f$ in 
$\omega$-$L^{1/\beta}(J,X)$. Since $f_n\in X_{\varphi}$ and $x_n=S(f_n)$,
$n\geq1$, Lemma \ref{continuousofsoluoperator} implies that
\[
x=S(f).
\]
Hence, to prove that 
$(x(\cdot),u(\cdot))\in\mathcal{R}_{\overline{\mathrm{co}}U}$, 
we only need to verify that $u(t)\in \operatorname{\overline{co}}U(t,x(t))$ 
a.e. $t\in J$.

Since $u_n\to u$ in $\omega$-$L^{1/\beta}(J,Y)$, by Mazur's theorem, we have
\begin{equation}\label{thmexistforconsys3}
u(t)\in\cap_{n=1}^{\infty}\operatorname{\overline{co}}\big(\cup_{k=n}^{\infty}
u_k(t)\big),\quad \text{for a.e. }t\in J.
\end{equation}
By (H4)(2) and the fact that 
$h(\operatorname{\overline{co}}A,\operatorname{\overline{co}}B)\leq h(A,B)$ 
for sets $A,B$, the map $x\to\operatorname{\overline{co}}U(t,x)$ is 
$h$-continuous. Then from Proposition 1.2.86 in \cite{hu}, the map 
$x\to\operatorname{\overline{co}}U(t,x)$ has property Q. Therefore we have
\begin{equation}\label{thmexistforconsys4}
\cap_{n=1}^{\infty}\operatorname{\overline{co}}
\Big(\cup_{k=n}^{\infty}\operatorname{\overline{co}}U(t,x_k(t))\Big)
\subseteq\operatorname{\overline{co}}U(t,x(t)),\quad
\text{for a.e. }t\in J.
\end{equation}
By \eqref{thmexistforconsys3} and \eqref{thmexistforconsys4}, we obtain 
that $u(t)\in\operatorname{\overline{co}}U(t,x(t))$ a.e. $t\in J$. 
This means that $\mathcal{R}_{\overline{\mathrm{co}}U}$ is compact in 
$C(J,X)\times\omega$-$L^{1/\beta}(J,Y)$. The proof is complete.
\end{proof}


\section{Main results}\label{secmainresults}

Now we are in a position to obtain our main results.

\begin{theorem}\label{maintheorem1}
For any $(x_*(\cdot),u_*(\cdot))\in\mathcal{R}_{\overline{\mathrm{co}}U}$, 
we have that there exists a sequence $(x_n(\cdot),u_n(\cdot))\in\mathcal{R}_U$, 
$n\geq1$, such that
\begin{gather}\label{jielun1}
x_n\to x_*\quad\text{in }C(J,X), \\
\label{jielun2}
u_n\to u_*\quad\text{in } L^{1/\beta}_{\omega}(J,Y)\quad
\text{and}\quad \omega\text{-}L^{1/\beta}(J,Y).
\end{gather}
Moreover, we have
\begin{equation}\label{jielun3}
\overline{\mathcal{T}r_{U}}=\mathcal{T}r_{\overline{\mathrm{co}}U},
\end{equation}
where the bar stands for the closure in the space $C(J,X)$.
\end{theorem}

\begin{proof}
Let any $(x_*(\cdot),u_*(\cdot))\in\mathcal{R}_{\overline{\mathrm{co}}U}$, 
then we have $u_*(t)\in\operatorname{\overline{co}}U(t,x_*(t))$ a.e. $t\in J$. 
It follows from (H4)(1), (H4)(2) and \eqref{boundedforu} that the map 
$t\to U(t,x_*(t))$ is measurable and integrally bounded. Hence by using 
\cite[Theorem 2.2]{tolstonogov5}, we have that, for any $n\geq1$, there 
exists a measurable selection $v_n(t)$ of the multivalued map $t\to U(t,x_*(t))$ 
such that
\begin{equation}\label{mainthm1}
\sup_{0\leq t_1\leq t_2\leq b}\big\|\int_{t_1}^{t_2}(u_*(s)-v_n(s))ds
\big\|_Y\leq\frac{1}{n}.
\end{equation}
For each fixed $n\geq1$, by (H4)(2), we have that, for any $x\in X$ and a.e. 
$t\in J$, there exists a $v\in U(t,x)$ such that
\begin{equation}\label{mainthm3}
\|v_n(t)-v\|_Y<k_1(t)\|x_*(t)-x\|_X+\frac{1}{n}.
\end{equation}
Let a map $H_n:J\times X\to 2^Y$ be defined by
\begin{equation}\label{mainthm4}
H_n(t,x)=\{v\in Y:v\text{ satisfies inequality \eqref{mainthm3}}\}.
\end{equation}
It follows from \eqref{mainthm3} that $H_n(t,x)$ is well defined for a.e.
 on $J$ and all $x\in X$, and its values are open sets.
 Using \cite[Corollary 2.1]{tolstonogov3} 
(since we can assume without loss of generality that $U(t,x)$ is 
$\Sigma\otimes\mathcal{B}_X$ measurable, see \cite[Proposition 2.7.9]{hu}), 
we obtain that, for any $\epsilon>0$, there is a compact set 
$J_{\epsilon}\subseteq J$ with $\mu(J\backslash J_{\epsilon})\leq\epsilon$, 
such that the restriction of $U(t,x)$ to $J_{\epsilon}\times X$ is l.s.c and 
the restrictions of $v_n(t)$ and $k_1(t)$ to $J_{\epsilon}$ are continuous. 
So \eqref{mainthm3} and \eqref{mainthm4} imply that the graph of the 
restriction of $H_n(t,x)$ to $J_{\epsilon}\times X$ is an open set in 
$J_{\epsilon}\times X\times Y$. Let a map $H:J\times X\to 2^Y$ be defined by
\begin{equation}\label{mainthm5}
H(t,x)=H_n(t,x)\cap U(t,x).
\end{equation}
It is obvious that, for a.e. $t\in J$ and all $x\in X$, $H(t,x)\neq\emptyset$. 
Due to the arguments above and Proposition 1.2.47 in \cite{hu}, we know that 
the restriction of $H(t,x)$ to $J_{\epsilon}\times X$ is l.s.c. and 
so does $\overline{H}(t,x)=\overline{H(t,x)}$, here the bar stands for the 
closure of a set in $Y$.

Now we consider the system \eqref{mainequation} with the  constraint 
on the control
\begin{equation}\label{ucapfcontrol}
u(t)\in\overline{H}(t,x(t))\quad \text{a.e. on } J.
\end{equation}
Since $\overline{H}(t,x)\subseteq U(t,x)$, the a priori estimate 
Lemma \ref{boundednessforsolution} also holds in this situation. 
Repeating the proof of Theorem \ref{existenceforcontrolsys}, 
we obtain that there is a solution $(x_n(\cdot),u_n(\cdot))$ 
of the control system \eqref{mainequation}, \eqref{ucapfcontrol}. 
The definition of $\overline{H}$ implies that 
$(x_n(\cdot),u_n(\cdot))\in\mathcal{R}_U$ and
\begin{equation}\label{mainthm7}
\|v_n(t)-u_n(t)\|_Y\leq k_1(t)\|x_*(t)-x_n(t)\|_X+\frac{1}{n}.
\end{equation}
Since $(x_n(\cdot),u_n(\cdot))\in\mathcal{R}_U$, $n\geq1$, and 
$(x_*(\cdot),u_*(\cdot))\in\mathcal{R}_{\overline{\mathrm{co}}U}$, we have
\begin{equation}\label{mainthm8}
x_*(t)=P_{\alpha}(t)x_0+\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)
\big(g(s)u_*(s)+h(s,x_*(s))\big)ds
\end{equation}
and
\begin{equation}\label{mainthm9}
x_n(t)=P_{\alpha}(t)x_0+\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)
\big(g(s)u_n(s)+h(s,x_n(s))\big)ds.
\end{equation}
Theorem \ref{existenceforcontrolsys} and 
$\{(x_n(\cdot),u_n(\cdot))\}_{n\geq1}\subseteq\mathcal{R}_U\subseteq
\mathcal{R}_{\overline{\mathrm{co}}U}$ imply that we can assume, 
possibly up to a subsequence, that the sequence 
$(x_n(\cdot),u_n(\cdot))\to(\overline{x}(\cdot),\overline{u}
(\cdot))\in\mathcal{R}_{\overline{\mathrm{co}}U}$ in 
$C(J,X)\times\omega$-$L^{1/\beta}(J,Y)$. Subtracting \eqref{mainthm9}
from \eqref{mainthm8}, and using (H3)(2), (H2)(2) and \eqref{mainthm7}, we have
\begin{align}
&\|x_*(t)-x_n(t)\|_X \nonumber \\
&=\Big\|\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)\big(g(s)u_*(s)-g(s)u_n(s)\big)ds
  \nonumber \\
&\quad +\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)\big(h(s,x_*(s))-h(s,x_n(s))
 \big)ds\Big\|_X  \nonumber\\
&\leq \big\|\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)g(s)\big(u_*(s)-v_n(s)\big)
  ds\big\|_X  \nonumber\\
&\quad +\big\|\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)g(s)
  \big(v_n(s)-u_n(s)\big)ds\big\|_X  \nonumber\\
&\quad +\big\|\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)\big(h(s,x_*(s))
 -h(s,x_n(s))\big)ds\big\|_X  \nonumber\\
&\leq \big\|\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)g(s)\big(u_*(s)
 -v_n(s)\big)ds\big\|_X  \nonumber\\
&\quad +\frac{\alpha M_Ad}{\Gamma(1+\alpha)}\int_0^t(t-s)^{\alpha-1}
 \Big(\frac{1}{n}+k_1(s)\|x_*(s)-x_n(s)\|_X\Big)ds  \nonumber\\
&\quad +\frac{\alpha M_A\|l\|_{L^{\infty}}}{\Gamma(1+\alpha)}
 \int_0^t(t-s)^{\alpha-1}\|x_*(s)-x_n(s)\|_X ds  \nonumber\\
&\leq \big\|\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)g(s)
 \big(u_*(s)-v_n(s)\big)ds\big\|_X
+\frac{\alpha M_Adb^{\alpha}}{n\alpha\Gamma(1+\alpha)}  \nonumber\\
&\quad +\frac{\alpha M_A(d\|k_1\|_{L^{\infty}}
 +\|l\|_{L^{\infty}})}{\Gamma(1+\alpha)}\int_0^t(t-s)^{\alpha-1}
 \|x_*(s)-x_n(s)\|_Xds. \label{mainthm10}
\end{align}
Due to \eqref{mainthm1}, one has $v_n\to u_*$ in $\omega$-$L^{1/\beta}(J,Y)$
 by Lemma \ref{weaknormforl2ty}. Then it is easy to show that 
$g(t)v_n(t)\to g(t)u_*(t)$ in $\omega$-$L^{1/\beta}(J,X)$. By the property
of the operator $H$ defined in the proof of Lemma \ref{continuousofsoluoperator},
 we have that, for any $t\in J$,
\[
\big\|\int_0^t(t-s)^{\alpha-1}Q_{\alpha}(t-s)g(s)\big(u_*(s)
 -v_n(s)\big)ds\big\|_X\to0,\quad \text{as }n\to\infty.
\]
Since $\|x_*(t)\|_X\leq L$, $\|x_n(t)\|_X\leq L$ for any $n$, 
$t\in J$ and $x_n\to\overline{x}$ in $C(J,X)$, letting $n\to\infty$ 
in \eqref{mainthm10}, we obtain
\[
\|x_*(t)-\overline{x}(t)\|_X\leq\frac{\alpha M_A(d\|k_1\|_{L^{\infty}}
+\|l\|_{L^{\infty}})}{\Gamma(1+\alpha)}\int_0^t(t-s)^{\alpha-1}\|x_*(s)
-\overline{x}(s)\|_Xds.
\]
Then by \cite[Theorem 3.1]{dixon}, we obtain $x_*=\overline{x}$; 
i.e., we have $x_n\to x_*$ in $C(J,X)$. Hence from $\eqref{mainthm7}$, 
we have $(v_n-u_n)\to 0$ in $L^{1/\beta}(J,Y)$.
Therefore, $u_n=u_n-v_n+v_n\to u_*$ in $\omega$-$L^{1/\beta}(J,Y)$ and
$L^{1/\beta}_{\omega}(J,Y)$, i.e., \eqref{jielun1} and \eqref{jielun2} hold.

Since it is clear that $\mathcal{T}r_U\subseteq
\mathcal{T}r_{\overline{\mathrm{co}}U}$ and 
$\mathcal{T}r_{\overline{\mathrm{co}}U}$ is compact in $C(J,X)$ by 
Theorem \ref{existenceforcontrolsys}, then from the proof of the first 
part of this theorem, we have
\[
\overline{\mathcal{T}r_U}=\mathcal{T}r_{\overline{\mathrm{co}}U},
\]
where the bar stands for the closure in $C(J,X)$. This completes the proof.
\end{proof}


\section{An example}
\label{sectionexample}

In this section, we present an example of control systems governed by 
fractional partial differential equations. In particular, to illustrate 
the abstract results of this paper, we provide the following example which 
do not aim at generality but indicate how our theorems can be applied to 
concrete problems. Since the hypotheses on the operator $g$ and the function 
$h$ are very common, we mainly pay attention to the operator $A$ and the 
multivalued map $U$ here.

Let $J=[0,1]$ and $\Omega=[0,\pi]$. Put $X=Y=L^2(\Omega)$. We consider the 
fractional control system
\begin{equation}\label{sectionexample1}
\begin{gathered}
^CD_t^{\alpha}x(t,z)=\partial^2_zx(t,z)+\bar{h}(t,z,x(t,z))
 +\bar{b}(t)\bar{u}(t,z),\quad t\in J,\; z\in\Omega,\\
x(t,0)=x(t,\pi)=0,\\
x(0,z)=x_0(z),\\
\bar{u}(t,z)\in \bar{U}(t,z,x(t,z)),\quad \text{a.e. in } J\times\Omega,
\end{gathered}
\end{equation}
where $^CD_t^{\alpha}$ is the Caputo fractional derivative of order 
$0<\alpha<1$, $\bar{h}$, $\bar{b}$ are suitable functions, $\bar{U}$ 
is a multivalued function which will be given below.

Define the operator $A$ by $A\omega=\omega''$ with $D(A)$ consisting of 
all $\omega\in X$ with $\omega$, $\omega'$ are absolutely continuous, 
$\omega''\in X$ and $\omega(0)=\omega(\pi)=0$. Then
\[
A\omega=-\sum_{n=1}^{\infty}n^2\langle\omega,e_n\rangle e_n,\quad\omega\in D(A),
\]
where $e_n(z)=(2/\pi)^{\frac{1}{2}}\sin(nz)$, $z\in\Omega$, $n=1,2,3,\cdots$, 
is the orthogonal set of eigenfunctions of $A$ and $\langle\cdot,\cdot\rangle$ 
denotes the $L^2$ inner product. It is clear that $A$ is the infinitesimal 
generator of a strongly continuous semigroup $\{T(t),\,t\geq0\}$ in $X$ and
 $T(t)$, $t>0$ is also compact, which is given by
\[
T(t)\omega=\sum_{n=1}^{\infty}e^{-n^2t}\langle\omega,e_n\rangle e_n,\quad 
\omega\in X.
\]
Hence the assumption (H1) is satisfied.
\begin{itemize}
\item[(H5)]  $\bar{U}:J\times\Omega\times\mathbb{R}\to \mathbb{R}$ 
is a multivalued function with closed values satisfying the following conditions:
\begin{itemize}
\item[(1)] the map $(t,z)\to\bar{U}(t,z,x)$ is measurable;
\item[(2)] $h(\bar{U}(t,z,x_1),\bar{U}(t,z,x_2))\leq\bar{k}_1(t)|x_1-x_2|$ 
  a.e. in $t\in J\times\Omega$ with $\bar{k}_1$ in $L^{\infty}_+(J)$;
\item[(3)] $|\bar{U}(t,z,x)|\leq\bar{a}_2(t,z)+\bar{c}_2(t,z)|x|$ a.e. in
  $J\times\Omega$ with $\bar{a}_2\in L^{1/\beta}(J,L^2_{+}(\Omega))$,
  $0<\beta<\alpha$ and $\bar{c}_2\in L^{\infty}_{+}(J\times\Omega)$.
\end{itemize}
\end{itemize}
Put $x(t)=x(t,\cdot)$; that is $x(t)(z)=x(t,z)$, $t\in J$, $z\in\Omega$.
 Define a multivalued map $U:J\times X\to 2^Y$ by
\[ %\label{exampleu}
U(t,x)=\{u\in Y:u(z)\in\bar{U}(t,z,x(z))\text{ a.e. in }\Omega\},\quad x\in X.
\]
Suppose assumption (H5) holds, then it is easy to verify that 
(H4)(1) and (H4)(2) are satisfied. Moveover, we have
\[
\sup\{\|u\|_Y:u\in U(t,x)\}\leq|\bar{a}_2(t)|_2
 +\|\bar{c}_2\|_{L^{\infty}}\|x\|_X,
\]
where $|\bar{a}_2(\cdot)|_2\in L^{1/\beta}(J,\mathbb{R}^+)$, 
$|\bar{a}_2(t)|_2=(\int_{\Omega}a^2(t,z)dz)^{\frac{1}{2}}$. 
This means that (H4)(3) holds.

Let $h(t,x)(z)=\bar{h}(t,z,x(t)(z))$ and $g(t)=\bar{b}(t)$. 
With $A$ and $U$ defined above, the fractional control system 
\eqref{sectionexample1} can be rewritten to our abstract form 
\eqref{mainequation}, \eqref{ucontrol}. Hence the abstract results obtained 
in the previous sections can be applied to the control system 
\eqref{sectionexample1}.


\subsection*{Conclusions} %\label{sectionconclusion}
Existence results and relaxation property of a class of fractional feedback 
control systems in Banach spaces have been investigated. With some 
auxiliary results provided in section \ref{uuxiliaryresults}, we 
obtained the existence results of the control systems by Schauder's 
fixed point theorem. To get the relaxation property, we used some tools 
from multivalued analysis.

Our future work will be devoted to study the following fractional control 
problems: systems with Riemann-Liouville fractional derivative, time 
optimal control, optimal control of Lagrange type and relaxed control systems 
by using other convexification techniques.

\subsection*{Acknowledgments}
The authors would like to thank the anonymous reviewers and the editor 
for their valuable comments and suggestions. 
This work was supported in part by the NNSFC (grants 10671211, 11126170)
and  Hunan Provincial Natural Science Foundation (grant 11JJ4006).


\begin{thebibliography}{00}


\bibitem{agarwal1} Ravi P. Agarwal, M. Belmekki, M. Benchohra;
\emph{A survey on semilinear differential equations and inclusions involving 
Riemann-Liouville fractional derivative}, Adv. Difference Equ. (2009) 
Article ID 981728, 47pp.

\bibitem{agarwal2} Ravi P. Agarwal, M. Benchohra, S. Hamani; 
\emph{A survey on existence results for boundary value problems of nonlinear 
fractional differential equations and inclusions}, 
Acta Appl. Math. 109 (2010) 973-1033.

\bibitem{aubin} J. P. Aubin, A. Cellina;
\emph{Differential Inclusions: Set-Valued Maps and Viability Theory}, 
Springer-Verlag, Berlin, New York, Tokyo, 1984.

\bibitem{Denkowski} Z. Denkowski, S. Mig\'orski, N. S. Papageorgiou;
\emph{An Introduction to Nonlinear Analysis: Applications},
 Kluwer Academic/Plenum Publishers, Boston, Dordrecht, London, New York, 2003.

\bibitem{dixon} J. Dixon, S. McKee;
\emph{Weakly Singular Discrete Gronwall Inequalities}, 
ZAMM. Z. angew. Math. Mech. 68(11) (1986) 535-544.

\bibitem{gansakmahantjam2013} R. Ganesh, R. Sakthivel, N. I. Mahmudov, 
S. M. Anthoni;
\emph{Approximate controllability of fractional integrodifferential 
evolution equations}, J. Appl. Math. Volume 2013 (2013) Article ID 291816, 
7 pages.

\bibitem{gaul} L. Gaul, P. Klein, S. Kempfle;
\emph{Damping description involving fractional operators}, Mech. Syst. 
Signal Process. 5 (1991) 81-88.

\bibitem{glockle} W. G. Glockle, T. F. Nonnenmacher;
\emph{A fractional calculus approach of self-similar protein dynamics}, 
Biophys. J. 68 (1995) 46-53.

\bibitem{Himmelberg} C. J. Himmelberg;
\emph{Measurable relations}, Fund. Math. 87 (1975) 53-72.

\bibitem{hilfer} R. Hilfer;
\emph{Applications of Fractional Calculus in Physics}, 
World Scientific, Singapore, 2000.

\bibitem{hu} S. C. Hu, N. S. Papageorgiou;
\emph{Handbook of multivalued Analysis: Volume I: Theory}, 
Kluwer Academic Publishers, Dordrecht Boston, London, 1997.

\bibitem{kilbas} A. A. Kilbas, H. M. Srivastava, J. J. Trujillo;
\emph{Theory and Applications of Fractional Differential Equations}, 
in: North-Holland Mathematics Studies, vol. 204, Elsevier Science B.V, 
Amsterdam, 2006.

\bibitem{surendra} S. Kumar, N. Sukavanam;
\emph{Approximate controllability of fractional order semilinear systems 
with bounded delay}, J. Differ. Equat. 252(11) (2012) 6163-6174.

\bibitem{liuzhenhai1} Zhenhai Liu, Jihua Sun;
\emph{Nonlinear boundary value problems of fractional functional 
integro-differential equations}, Comput. Math. Appl. 64(10) (2012) 3228-3234.

\bibitem{liuzhenhai2} Zhenhai Liu, Jihua Sun;
\emph{Nonlinear boundary value problems of fractional differential systems}, 
Comput. Math. Appl. 64(4) (2012) 463-475.

\bibitem{migorski1} S. Mig\'orski;
\emph{Existence and relaxation results for nonlinear evolution inclusions
 revisited}, J. Appl. Math. Stoch. Anal. 8(2) (1995) 143-149.

\bibitem{migorski2} S. Mig\'orski;
\emph{Existence and relaxation results for nonlinear second order evolution 
inclusions}, Discuss. Math. Differential Inclusions 15 (1995) 129-148.

\bibitem{miller} K. S. Miller, B. Ross;
\emph{An Introduction to the Fractional Calculus and Differential Equations}, 
John Wiley, New York, 1993.

\bibitem{sakganantamc2013} R. Sakthivel, R. Ganesh, S. M. Anthoni;
\emph{Approximate controllability of fractional nonlinear differential 
inclusions}, Appl. Math. Comput. 225 (2013) 708-717.

\bibitem{sakganrenantcnsns2013} R. Sakthivel, R. Ganesh, Yong Ren, S. M. Anthoni;
\emph{Approximate controllability of nonlinear fractional dynamical systems}, 
Commun. Nonlinear Sci. Numer. Simul. 18 (2013) 3498-3508.

\bibitem{sakthivela1} R. Sakthivel, N. I. Mahmudovb, Juan J. Nieto;
\emph{Controllability for a class of fractional-order neutral evolution control 
systems}, Appl. Math. Comput. 218 (2012) 10334-10340.

\bibitem{sakthivela} R. Sakthivel, Yong Ren, N. I. Mahmudov;
\emph{On the approximate controllability of semilinear fractional 
differential systems}, Comput. Math. Appl. 62 (2011) 1451-1459.

\bibitem{tolstonogov12} A. A. Tolstonogov;
\emph{Relaxation in nonconvex optimal control problems with subdifferential 
operators}, Journal of Mathematical Sciences 140(6) (2007) 850-872.

\bibitem{tolstonogov2} A. A. Tolstonogov;
\emph{Relaxation in non-convex control problems described by first-order 
evolution equations}, Mat. Sb. 190:11 (1999) 135-160;
 English transl., Sb. Math. 190 (1999) 1689-1714.

\bibitem{tolstonogov3} A. A. Tolstonogov;
\emph{Scorza-Dragoni's theorem for multi-valued mappings with variable 
domain of definition}, Mat. Zametki 48:5 (1990) 109-120; English transl., 
Math. Notes 48 (1990) 1151-1158.

\bibitem{tolstonogov11} A. A. Tolstonogov, D. A. Tolstonogov;
\emph{On the ``bang-bang" principle for nonlinear evolution inclusions},
 Nonlinear differ. equ. appl. 6 (1999) 101-118.

\bibitem{tolstonogov4} A. A. Tolstonogov, D. A. Tolstonogov;
\emph{Lp-continuous extreme selectors of multifunctions with decomposable
 values: Existence theorems}, Set-Valued Anal. 4 (1996) 173-203.

\bibitem{tolstonogov5} A. A. Tolstonogov, D. A. Tolstonogov;
\emph{Lp-continuous extreme selectors of multifunctions with decomposable 
values: Relaxation theorems}, Set-Valued Anal. 4 (1996) 237-269.

\bibitem{wangjinrong2} JinRong Wang, Yong Zhou;
\emph{Existence and controllability results for fractional semilinear 
differential inclusions}, Nonlinear Anal. 12(6) (2011) 3642-3653.

\bibitem{wangjinrong1} JinRong Wang, Yong Zhou;
\emph{A class of fractional evolution equations and optimal controls}, 
Nonlinear Anal. 12 (2011) 262-272.

\bibitem{yongzhou1} Yong Zhou, Feng Jiao;
\emph{Existence of mild solutions for fractional neutral evolution equations}, 
Comput. Math. Appl. 59 (2010) 1063-1077.

\bibitem{yongzhou2} Yong Zhou, Feng Jiao;
\emph{Nonlocal Cauchy problem for fractional evolution equations}, 
Nonlinear Anal. 11 (2010) 4465-4475.

\bibitem{zhu1} J. Zhu;
\emph{On the solution set of differential inclusions in Banach space}, 
J. Differ. Equat. 93(2) (1991) 213-237.

\end{thebibliography}

\end{document}
