Processing math: 89%
Research article

Factors Associated With Sleep Quality in Patients With Multiple Sclerosis

  • Objective: A limited number of studies with inconsistent results have assessed the factors associated with sleep quality in patients with multiple sclerosis (MS). This study aimed to evaluate sleep quality and to investigate the associations between sleep quality and demographic and health-related characteristics, anemia, fatigue and physical activity in Turkish patients with MS. Methods: A cross-sectional study was conducted in a sample of 102 patients with MS who were followed in a neurology outpatient clinic of a tertiary hospital in Turkey between March 2015 and November 2015. Data were collected by an information form, the Pittsburgh Sleep Quality Index, the Visual Analogue Scale for Fatigue and the short version of the International Physical Activity Questionnaire. Anemia was evaluated by measuring hemoglobin levels. Data analysis were performed using descriptive statistics, Mann-Whitney U-test, Kruskal-Wallis H test, Spearman’s correlation coefficients and logistic regression analysis with backward stepwise elimination. Results: The mean global Pittsburgh Sleep Quality Index score was 5.98 ± 3.94, and 52.0% of the participants reported poor sleep quality. Twenty patients (19.6%) had anemia. Patients with higher fatigue were more likely to have poor sleep quality (adjusted Odds Ratio, 1.283; 95% confidence interval, 1.031–1.596, p = 0.026). Conclusions: Poor sleep quality was common in patients with MS and higher levels of fatigue predicted poorer sleep quality. A better understanding of risk factors related to sleep quality may facilitate effective interventions that improve health outcomes.

    Citation: Belgüzar Kara, Elif Gökçe Tenekeci, Şeref Demirkaya. Factors Associated With Sleep Quality in Patients With Multiple Sclerosis[J]. AIMS Medical Science, 2016, 3(2): 203-212. doi: 10.3934/medsci.2016.2.203

    Related Papers:

    [1] Jun Moon . The Pontryagin type maximum principle for Caputo fractional optimal control problems with terminal and running state constraints. AIMS Mathematics, 2025, 10(1): 884-920. doi: 10.3934/math.2025042
    [2] Yuna Oh, Jun Moon . The infinite-dimensional Pontryagin maximum principle for optimal control problems of fractional evolution equations with endpoint state constraints. AIMS Mathematics, 2024, 9(3): 6109-6144. doi: 10.3934/math.2024299
    [3] Fátima Cruz, Ricardo Almeida, Natália Martins . A Pontryagin maximum principle for optimal control problems involving generalized distributional-order derivatives. AIMS Mathematics, 2025, 10(5): 11939-11956. doi: 10.3934/math.2025539
    [4] Chuanli Wang, Biyun Chen . An hp-version spectral collocation method for fractional Volterra integro-differential equations with weakly singular kernels. AIMS Mathematics, 2023, 8(8): 19816-19841. doi: 10.3934/math.20231010
    [5] Irmand Mikiela, Valentina Lanza, Nathalie Verdière, Damienne Provitolo . Optimal strategies to control human behaviors during a catastrophic event. AIMS Mathematics, 2022, 7(10): 18450-18466. doi: 10.3934/math.20221015
    [6] Dumitru Baleanu, Babak Shiri . Generalized fractional differential equations for past dynamic. AIMS Mathematics, 2022, 7(8): 14394-14418. doi: 10.3934/math.2022793
    [7] Ahmed M. Rajab, Saeed Pishbin, Javad Shokri . Analyzing the structure of solutions for weakly singular integro-differential equations with partial derivatives. AIMS Mathematics, 2024, 9(9): 23182-23196. doi: 10.3934/math.20241127
    [8] Regina S. Burachik, Bethany I. Caldwell, C. Yalçın Kaya . Douglas–Rachford algorithm for control- and state-constrained optimal control problems. AIMS Mathematics, 2024, 9(6): 13874-13893. doi: 10.3934/math.2024675
    [9] Xiangyun Shi, Xiwen Gao, Xueyong Zhou, Yongfeng Li . Analysis of an SQEIAR epidemic model with media coverage and asymptomatic infection. AIMS Mathematics, 2021, 6(11): 12298-12320. doi: 10.3934/math.2021712
    [10] Ahmad Alalyani, M. A. Abdou, M. Basseem . The orthogonal polynomials method using Gegenbauer polynomials to solve mixed integral equations with a Carleman kernel. AIMS Mathematics, 2024, 9(7): 19240-19260. doi: 10.3934/math.2024937
  • Objective: A limited number of studies with inconsistent results have assessed the factors associated with sleep quality in patients with multiple sclerosis (MS). This study aimed to evaluate sleep quality and to investigate the associations between sleep quality and demographic and health-related characteristics, anemia, fatigue and physical activity in Turkish patients with MS. Methods: A cross-sectional study was conducted in a sample of 102 patients with MS who were followed in a neurology outpatient clinic of a tertiary hospital in Turkey between March 2015 and November 2015. Data were collected by an information form, the Pittsburgh Sleep Quality Index, the Visual Analogue Scale for Fatigue and the short version of the International Physical Activity Questionnaire. Anemia was evaluated by measuring hemoglobin levels. Data analysis were performed using descriptive statistics, Mann-Whitney U-test, Kruskal-Wallis H test, Spearman’s correlation coefficients and logistic regression analysis with backward stepwise elimination. Results: The mean global Pittsburgh Sleep Quality Index score was 5.98 ± 3.94, and 52.0% of the participants reported poor sleep quality. Twenty patients (19.6%) had anemia. Patients with higher fatigue were more likely to have poor sleep quality (adjusted Odds Ratio, 1.283; 95% confidence interval, 1.031–1.596, p = 0.026). Conclusions: Poor sleep quality was common in patients with MS and higher levels of fatigue predicted poorer sleep quality. A better understanding of risk factors related to sleep quality may facilitate effective interventions that improve health outcomes.


    In this paper, we consider the minimization problem of

    J(x0,u())=T0l(s,x(s),u(s))ds+h(x0,x(T)), (1.1)

    subject to the Rn-valued integral-type state equation with α(0,1),

    x(t)=x0+t0f(t,s,x(s),u(s))(ts)1αds,a.e.t[0,T], (1.2)

    and the terminal state constraint

    (x0,x(T))FR2n. (1.3)

    In (1.1)–(1.3), x()Rn is the state with the initial state x0Rn, and u()U is the control with U being the control space (see Section 2). Hence, we consider the following terminal state-constrained optimal control problem:

    (P)infu()Up[0,T]J(x0,u()),subject to (1.2) and (1.3),

    where Up[0,T] is the set of admissible controls of (1.2). The problem aforementioned is referred to as (P) throughout this paper, and the precise problem statement of (P) is provided in Section 2. Note that the optimal control problems with and without state constraints capture various practical aspects of systems in science, biology, engineering, and economics [2,4,5,6,29,41,43].

    The state equation in (1.2) is known as a class of Volterra integral equations. The main feature of Volterra integral equations is the effect of memories, which does not appear in ordinary (state) differential equations. In fact, Volterra integral equations of various kinds have been playing an important role in modeling and analyzing of practical physical, biological, engineering, and other phenomena that are governed by memory effects [7,17,19,32,34]. We note that one major distinction between (1.2) and other classical Volterra integral equations is that (1.2) has a kernel f(t,s,x,u)(ts)1α, which becomes singular at s=t. In fact, α(0,1) in the singular kernel f(t,s,x,u)(ts)1α of (1.2) determines the amount of the singularity, in which the large singular behavior occurs with small α(0,1). Hence, the singular kernel in (1.2) can be applied to various areas of science, engineering, finance and economics, in which certain extraordinary phenomena have to be described. For example, the model of fluid dynamics may have some singular phenomena, leading to the deformation of the corresponding dynamic equation [17,33,36].

    The study of optimal control problems for various kinds of Volterra integral equations via the maximum principle have been studied extensively in the literature; see [3,12,13,18,20,22,23,30,35,39,40] and the references therein. Specifically, the first study on optimal control for Volterra integral equations (using the maximum principle) can be traced back to [40]. Several different formulations (with and without state constraints and/or delay) of optimal control for Volterra integral equations and their generalizations are reported in [3,7,12,18,20,30,35,39]. Some recent progress in different directions including stochastic frameworks can be found in [13,22,23,26,27,42]. We note that the above-mentioned existing works (with the exception of the stochastic case in [26,27,42]) considered the situation with nonsingular kernels only in Volterra integral equations, which corresponds to the case of f(t,s,x,u)=(ts)1αg(t,s,x,u) in (1.2). Hence, the problem settings in the earlier works can be viewed as a special case of (P).

    Recently, (P) without the terminal state constraint in (1.3) was studied in [34]. Due to the presence of the singular kernel, the technical analysis including the maximum principle (without the terminal state constraint) in [34] should be different from that of the existing works mentioned above. In particular, the proof for the well-posedness and estimates of Volterra integral equations in [34, Theorem 3.1] requires a new type of Gronwall's inequality. Furthermore, the maximum principle without the terminal state constraint in [34, Theorem 4.3] needs a different duality and variational analysis via spike perturbation. More recently, the linear-quadratic optimal control problem (without the state constraint) for linear Volterra integral equations with singular kernels was studied in [28].

    We note that the state equation in (1.2) is strongly related to classical state differential equations and (Caputo) fractional order differential equations. In particular, let DCα[x()] be the Caputo fractional derivative operator of order α(0,1) [32, Chapter 2.4]. Subsequently, by [32, Theorem 3.24 and Corollary 3.23], (1.2) becomes (Γ() denotes the gamma function)

    DCα[x()](t)=f(t,x(t),u(t))x(t)=x0+1Γ(α)t0f(s,x(s),u(s))(ts)1αds. (1.4)

    In addition, when f(t,s,x,u)=(ts)1αg(s,x,u), (1.2) is specialized to

    dx(t)dt=g(t,x(t),u(t))x(t)=x0+t0g(s,x(s),u(s))ds. (1.5)

    As (1.4) and (1.5) are special cases of our state equation in (1.2), the state equation in (1.2) is able to describe various types of differential equations including fractional and ordinary differential equations. We also mention that there are several different results on optimal control for fractional differential equations; see [1,8,14,25,31] and the references therein.

    The aim of this paper is to obtain the Pontryagin maximum principle for (P). As noted above, since [34] did not consider the state-constrained control problem, (P) can be viewed as a generalization of [34] to the terminal state-constrained control problem. The main theoretical significances and technical difficulties are the proof of the maximum principle in Theorem 3.1, which is presented in Section 5. In particular, due to (i) the inherent complex nature of the Volterra integral equation with singular kernels, (ii) the presence of the terminal state constraint, and (iii) the generalized standing assumptions of f, our maximum principle and its detailed proof must be different from those provided in the existing literature (e.g., [12,18,22,23,30,34,35,39]). Note also that as our state equation in (1.2) includes the case of Caputo fractional differential equation in (1.4), the maximum principle in Theorem 3.1 also covers [8, Theorem 3.12]. However, as the variational analysis of this paper and that of [8] are different, the proof techniques including the characterization of the adjoint equation must be different. Indeed, the maximum principle of this paper is new in the optimal control problem context and its proof requires a different technique, compared with the existing literature.

    The detailed statements of the main results of this paper are provided as follows:

    ● Regarding (ii), as mentioned above, the proof in Section 5 must be different from that of the unconstrained control problem in [34, Theorem 4.3]. Specifically, in contrast to [34, Theorem 4.3], due to the terminal state constraint in (1.3), we have to employ the Ekeland variational principle and the spike variation technique, together with the intrinsic properties of the distance function and the generalized Gronwall's inequality, to establish the duality analysis for Volterra-type forward variational and backward adjoint equations with singular kernels (see Lemma 5.2, Lemma 5.3 and Section 5.6). Note that this analysis leads to the desired necessary condition for optimality. Such a generalized proof of the maximum principle including the duality and variational analysis was not needed in [34, Theorem 4.3], as it did not consider the terminal state constraint in the corresponding optimal control problem. Note also that without the terminal state constraint, our maximum principle in Theorem 3.1 is reduced to the unconstrained case in [34, Theorem 4.3] (see Remark 3.1).

    ● As for (i), since our Volterra integral equation in (1.2) covers both singular and nonsingular kernels, the proof for the maximum principle in Section 5 should be different from that of the classical state-constrained maximum principle with the nonsingular kernel only (when f(t,s,x,u)=(ts)1αg(t,s,x,u) in (1.2)) studied in the existing literature (e.g., [12, Theorem 1], [22, Theorem 3.1] and [18,23,30,35,39]). In particular, due to the presence of the singular kernel in (1.2), we need a different proof for duality and variational analysis (see Lemma 3, Lemma 5.3 and Section 5.6). In addition, unlike the existing literature (e.g., [12, Theorem 1] and [22, Theorem 3.1]), our adjoint equation is obtained by the backward Volterra integral equation with the singular kernel, as a consequence of the new duality analysis with the singular kernel in Section 5.6 (see Theorem 3.1 and Remark 3.2).

    ● Concerning (iii), we mention that unlike the existing works for classical optimal control of Volterra integral equations with nonsingular kernels (e.g., [12,(4) and (5)], [22, page 3437 and (A3)], and [18,23,39]), our paper assumes neither the differentiability of (singular or nonsingular) kernels in (t,s,u) (time and control variables) nor the essentially boundedness of a class of admissible controls (see Remark 3.2). Hence, the detailed proof of the maximum principle provided in this paper must be different from those given in the existing literature.

    As mentioned above, the main motivation of studying (P) is to provide the general maximum principle for terminal state-constrained optimal control problems of various classes of differential and integral type equations. Specifically, as mentioned above, (1.2) covers Volterra integral equations with singular and nonsingular kernels (see Remark 3.2), ordinary differential equations (see (1.5)), and Caputo-type fractional differential equations (see (1.4)). Hence, the terminal state-constrained optimal control problems for such kind of (differential and integral) state equations can be solved via the maximum principle of this paper (see Theorem 3.1). As stated above, the maximum principle for (P) has not been presented in the existing literature.

    We mention that Volterra integral equations and terminal state-constrained optimal control problems can be applied to various practical applications in science, engineering, economics, and mathematical finance; see [2,4,5,6,7,17,19,29,32,34,41,43] and the references therein. In particular, Volterra integral equations can be used to study the so-called tautochrone problem in mechanical applications, the heat transfer problem in diffusion models, the shock wave problem, the renewal equation in the theory of industrial engineering, the investment problem in economics, and the problem of American option pricing [7,17,19,32,34]. The interacting biological population model can be described by Volterra integral equations [7]. Furthermore, as mentioned above, in a practical point of views, the singular kernel in (1.2) can be applied to those Volterra integral equations to study their singular and/or peculiar behavior. Hence, we can use the modeling framework of (1.2) for any differential and Volterra integral equations to capture the memory and/or singular effects, which can be observed in real world (e.g., the deformation of fluid dynamics due to singular phenomena [17,33,36], and the dependency of the current stock price on past investment strategies over a period of time [7,17]). We also note that the (terminal) state-constrained optimal control problems can be formulated for these examples (e.g., the Volterra-type renewal integral equation in industrial applications, the Volterra-type heat transfer integral equation in diffusion models, and the Volterra-type investment integral equation in economics). In fact, there are huge potential applications of optimal control problems, and their various different practical examples in science, engineering, biology, economics, and mathematical finance, including the optimal control for partial differential equations (PDEs), can be found in [2,4,5,6,29,41,43] and the references therein.

    In view of the preceding discussion, we believe that the maximum principle of this paper broadens both the theoretical generality and the applicability of the domain of optimal control problems.

    The rest of this paper is organized as follows. The problem statement of (P) is given in Section 2. The statement of the maximum principle for (P) is provided in Section 3. We study an example of (P) in Section 4. The proof of the maximum principle is given in Section 5. Finally, we conclude this paper in Section 6.

    Let Rn be an n-dimensional Euclidean space, where x,y:=xy is the inner product and |x|:=x,x1/2 is the norm for x,yRn. For ARm×n, A denotes the transpose of A. Let [0,T] be a time interval such that T<. Let 1S() be the indicator function of any set S. In this paper, C0 denotes a generic constant, whose value is different from line to line. For any differentiable function f:RnRl, let fx:RnRl×n be the partial derivative of f with respect to xRn. fx=[f1,xfl,x] with fj,xR1×n, and when l=1, fxR1×n. For f:Rn×RlRl, fx:Rn×RlRl×n for xRn, and fy:Rn×RlRl×l for yRl. For 1p<, define

    Lp([0,T];Rn): The space of functions ψ:[0,T]Rn such that ψ is measurable and satisfies ψ()Lp([0,T];Rn):=(T0|ψ(t)|pdt)1/p<;

    C([0,T];Rn): The space of functions ψ:[0,T]Rn such that ψ is continuous and satisfies ψ():=supt[0,T]|ψ(t)|<.

    Next, we state the precise problem statement of (P) provided in Section 1. In (1.2), α(0,1) is the parameter of the singularity, x()Rn is the state with the initial state x0Rn, and u()U is the control with U being the control space. In (1.2), f(t,s,x,u)(ts)1α denotes the singular kernel (with the singularity appearing at s=t), with f:[0,T]×[0,T]×Rn×URn being a generator. Note that α(0,1) determines the level of singularity of (1.2); see Figure 1. In (1.2), f is dependent on two time parameters, t and s. While t is the outer time variable to determine the current time, s is the inner time variable describing the path or memory of (1.2) from 0 to t. In (1.1), l:[0,T]×Rn×UR is the running cost, while h:Rn×RnR is the terminal cost.

    Figure 1.  State trajectories when x0=1 and f(t,s,x,u)=0.4sin(2πx)(ts)1αx. Note that the state trajectory exhibits further singular behavior for small α(0,1).

    Assumption 1. (i) (U,ρ) is a separable metric space. f and fx are continuous on [0,T]×[0,T] and Lipschitz continuous in (x,u)Rn×U with the Lipschitz constant K0. |f(t,s,0,u)|K for (t,s,u)[0,T]×[0,T]×U.

    (ii) l and lx are continuous on [0,T] and Lipschitz continuous in (x,u)Rn×U with the Lipschitz constant K0. |l(t,0,u)|K for (t,u)[0,T]×U. h, hx0, and hx are Lipschitz continuous in both variables with the Lipschitz constant K0.

    (iii) F is a nonempty closed convex subset of R2n.

    For p>1α and u0U, we define the space of admissible controls for (1.2) as

    Up[0,T]={u:[0,T]U|u is measurable&ρ(u(),u0)Lp([0,T];R+)}.

    Note that p>1α is needed for the well-posedness of (1.2), which is stated in the following lemma [34, Theorem 3.1 and Proposition 3.3]:

    Lemma 2.1. Let (i) of Assumption 1 hold and p>1α. For any (x0,u())Rn×Up[0,T], (1.2) admits a unique solution in C([0,T];Rn), and there is a constant C0 such that x()Lp([0,T];Rn)C(1+|x0|+ρ(u(),u0)Lp([0,T];R+)).

    Under Assumption 1, the main objective of this paper is to solve (P) via the Pontryagin maximum principle. Note that Assumption 1 is crucial for the well-posedness of the state equation in (1.2) by Lemma 2.1 as well as the maximum principle of (P). Assumptions similar to Assumption 1 have been used in various optimal control problems and their maximum principles; see [9,12,13,15,16,20,22,23,27,33,34,37,39,40,43] and the references therein. Note also that we do not need the differentiability of f in (t,s,u), which was assumed in the existing literature (e.g., [12,(4) and (5)], [22, page 3437 and (A3)], and [18,23,39]); see Remark 3.2.

    We provide the statement of the maximum principles for (P). The proof is given in Section 5.

    Theorem 3.1. Let Assumption 1 hold. Suppose that (¯u(),¯x())Up[0,T]×C([0,T];Rn) is the optimal pair for (P), where ¯x()C([0,T];Rn) is the corresponding optimal state trajectory of (1.2). Then there exists a pair (λ,ξ), where λR and ξ=(ξ1,ξ2)R2n, such that the following conditions hold:

    Nontriviality condition: (λ,ξ)0, where λ0 and ξNF(¯x0,¯x(T)) with NF(x) being the normal cone to the convex set F defined in (5.1);

    Adjoint equation: p()Lp([0,T];Rn) is the unique solution to the following backward Volterra integral equation with the singular kernel:

    p(t)=λlx(t,¯x(t),¯u(t))+Ttfx(r,t,¯x(t),¯u(t))(rt)1αp(r)dr1[0,T)(t)fx(T,t,¯x(t),¯u(t))(Tt)1α(λhx(¯x0,¯x(T))+ξ2),a.e.t[0,T];

    Transversality conditions:

    (i)0ξ1,¯x0y1+ξ2,¯x(T)y2,y=(y1,y2)F,(ii)T0p(t)dt=ξ1+ξ2+λhx0(¯x0,¯x(T))+λhx(¯x0,¯x(T));

    Hamiltonian-like maximum condition:

    λl(t,¯x(t),¯u(t))+Ttp(r)f(r,t,¯x(t),¯u(t))(rt)1αdr1[0,T)(t)(λhx(¯x0,¯x(T))+ξ2)f(T,t,¯x(t),¯u(t))(Tt)1α=maxuU{λl(t,¯x(t),u)+Ttp(r)f(r,t,¯x(t),u)(rt)1αdr1[0,T)(t)(λhx(¯x0,¯x(T))+ξ2)f(T,t,¯x(t),u)(Tt)1α},a.e.t[0,T].

    Remark 3.1. Without the terminal state constraint in (1.3), Theorem 3.1 holds with λ=1 and ξ=0. This is equivalent to [34, Theorem 4.3].

    Remark 3.2. By taking f(t,s,x,u)=(ts)1αg(t,s,x,u) in Theorem 3.1, we obtain the maximum principle for the optimal control problems of classical Volterra integral equations with nonsingular kernels. Note that Theorem 3.1 is more general than the existing maximum principles for Volterra integral equations with nonsingular kernels (e.g., [12, Theorem 1], [22, Theorem 3.1] and [18,23,30,35,39]), where Theorem 3.1 assumes neither the differentiability of the corresponding kernel with respect to time and control variables nor the essential boundedness of a class of admissible controls (e.g., [12,(4) and (5)], [22, page 3437 and (A3)]). In addition, unlike the existing literature, the adjoint equation in Theorem 3.1 must be obtained in integral form, as a consequence of the new duality analysis with singular kernels in Section 5.6.

    The motivation of this section is to demonstrate the applicability of Theorem 3.1 to the state-constrained linear-quadratic problem for Volterra integral equations with singular kernels. We consider the minimization of the objective functional of J(x0,u())=T0[x(s)+12u(s)2]ds subject to the R-valued Volterra integral equation with the singular kernel:

    x(t)=x0+t0t2u(s)(ts)1αds,a.e. t[0,T]. (4.1)

    For Case I, the terminal state constraint is given by (x0,x(T))F={x0=10,x(T)=10}R2. Here, we may assume that U is an appropriate sufficiently large compact set satisfying Assumption 1. Note that this example is closely related to the (terminal) state-constrained linear-quadratic control problem, which can be applied to study the singular aspects of various applications in science, engineering, and economics.

    Let T=3. By Theorem 3.1, the adjoint equation is p(t)=λ. From the first transversality condition in (i), we have ξR2. Then by the second transversality condition in (ii), 30p(t)dt=3λ=ξ1+ξ2. Hence, we can choose λ=1 and ξ1=ξ2=1.5. By the Hamiltonian-like maximum condition and the first-order optimality condition, the (candidate) optimal solution of Case I is as follows:

    ¯u(t)=3tt2(rt)1αdr+t21.51[0,3)(t)(3t)1α,a.e. t[0,3]. (4.2)

    We apply (4.2) to get the optimal state trajectory of (4.1). For Case I, the optimal state trajectory of (4.1) controlled by (4.2) when α=0.1 and α=0.5 is shown in Figure 2. Similar1y, when T=10, the optimal state trajectory is depicted in Figure 3.

    Figure 2.  Case I: The optimal state trajectory of (4.1) when T=3 (left: α=0.1, right: α=0.5).
    Figure 3.  Case I: The optimal state trajectory of (4.1) when T=10 (left: α=0.1, right: α=0.5).

    We now consider Case II, where the (terminal) state constraint is generalized by

    (x0,x(T))F={x0=10,x(T)[10,10]}R2.

    Similar to Case I, the adjoint equation is given by p(t)=λ. The first transversality condition in (i) implies that ξ1R and ξ2=0. From the second transversality condition in (ii), it follows that 30p(t)dt=3λ=ξ1. Hence, we can choose λ=1, ξ1=3, and ξ2=0. Then analogous to (4.2), the (candidate) optimal solution of Case II is given by ¯u(t)=3tt2(rt)1αdr, a.e. t[0,3]. Figure 4 depicts the optimal state trajectory (with α=0.1 and α=0.5) when T=3, and Figure 5 shows the case when T=10. Note that the numerical results of this section are obtained by the finite-difference method to approximate the solutions of (4.1) and (4.2).

    Figure 4.  Case II: The optimal state trajectory of (4.1) when T=3 (left: α=0.1, right: α=0.5).
    Figure 5.  Case II: The optimal state trajectory of (4.1) when T=10 (left: α=0.1, right: α=0.5).

    Remark 4.1. The purpose of the examples in this section is to show an analytic method for applying Theorem 3.1. In the future research problem, a shooting-like method has to be developed to solve (P) numerically, which can be applied to several different complex situations of (P). This requires extending the classical approach (e.g., [11,16]) to the case of Volterra integral equations with singular kernels.

    This section provides the proof of Theorem 3.1.

    Recall that F is a nonempty closed convex subsets of R2n. Let dF:R2nR+ be the standard Euclidean distance function to F defined by dF(x):=infyF|xy| for xR2n. Note that dF(x)=0 when xF. Then it follows from the projection theorem [38, Theorem 2.10] that there exists a unique PF(x)F with PF(x):R2nFR2n, the projection of xR2n onto F, such that dF(x)=infyF|xy|=|xPF(x)|. By [38, Lemma 2.11], PF(x)F is the corresponding projection if and only if xPF(x),yPF(x)0 for all yF, which leads to the characterization of PF(x). In view of [38, Definition 2.37], we have xPF(x)NF(PF(x)) for xR2n, where NF(x) is the normal cone to the convex set F at a point xF defined by

    NF(x):={yR2n|y,yx0,yF}. (5.1)

    Based on the distance function dF, (1.3) can be written by dF(¯x0,¯x(T))=0(¯x0,¯x(T))F.

    Lemma 5.1. [24, page 167] and [21, Proposition 2.5.4] The function dF(x)2 is Fréchet differentiable on R2n with the Fréchet differentiation of dF(x)2 at xR2n given by DdF(x)2(h)=2xPF(x),h for hR2n.

    Recall that the pair (¯x(),¯u())C([0,T];Rn)×Up[0,T] is the optimal pair of (P). Note that the pair (¯x0,¯x(T)) holds the terminal state constraint in (1.3). The optimal cost of (P) under (¯x(),¯u()) can be written by J(¯x0,¯u()).

    Recall the distance function dF in Section 5.1. For ϵ>0, we define the penalized objective functional as follows:

    Jϵ(x0,u())=(([J(x0,u())J(¯x0,¯u())+ϵ]+)2+dF(x0,x(T))2)12. (5.2)

    We can easily observe that Jϵ(¯x0,¯u())=ϵ>0, i.e., (¯x0,¯u())Rn×Up[0,T] is the ϵ-optimal solution of (5.2). We define the Ekeland metric ˆd:(Rn×Up[0,T])×(Rn×Up[0,T])R+:

    ˆd((x0,u()),(˜x0,˜u())):=|x0˜x0|+¯d(u(),˜u()), (5.3)

    where

    ¯d(u(),˜u()):=|{t[0,T]|u(t)˜u(t)}|,u(),˜u()Up[0,T]. (5.4)

    Note that (Rn×Up[0,T],ˆd) is a complete metric space, and Jϵ(x0,u) in (5.2) is continuous on (Rn×Up[0,T],ˆd) [33, Proposition 3.10 and Corollary 4.2,Chapter 4].

    In view of (5.2)–(5.4), we have

    {Jϵ(x0,u())>0,(x0,u())Rn×Up[0,T],Jϵ(¯x0,¯u())=ϵinf(x0,u())Rn×Up[0,T]Jϵ(x0,u())+ϵ. (5.5)

    By the Ekeland variational principle [33, Corollary 2.2,Chapter 4], there exists a pair (xϵ0,uϵ)Rn×Up[0,T] such that

    ˆd((xϵ0,uϵ()),(¯x0,¯u()))ϵ, (5.6)

    and for any (x0,u())Rn×Up[0,T],

    {Jϵ(xϵ0,uϵ())Jϵ(¯x0,¯u())=ϵ,Jϵ(xϵ0,uϵ())Jϵ(x0,u())+ϵˆd((xϵ0,uϵ()),(x0,u())). (5.7)

    We write (xϵ(),uϵ())C([0,T];Rn)×Up[0,T], where xϵ() is the state trajectory of (1.2) under (xϵ0,uϵ())Rn×Up[0,T].

    The next step is to derive the necessary condition for (xϵ0,uϵ())Rn×Up[0,T]. To this end, we employ the spike variation technique, as U does not have any algebraic structure. Hence, standard (convex) variations cannot be used.

    For δ(0,1), we define Eδ:={E[0,T]||E|=δT}, where |E| denotes the Lebesgue measure of E. For EδEδ, we introduce the spike variation associated with uϵ:

    uϵ,δ(s):={uϵ(s),s[0,T]Eδ,u(s),sEδ,

    where u()Up[0,T]. Clearly, uϵ,δ()Up[0,T]. By definition of ¯d in (5.4),

    ¯d(uϵ,δ(),uϵ())|Eδ|=δT. (5.8)

    We also consider the variation of the initial state given by xϵ0+δa, where aRn. Let xϵ,δ()C([0,T];Rn) be the state trajectory of (1.2) under (xϵ0+δa,uϵ,δ())Rn×Up[0,T]. By (5.7) and (5.8), it follows that

    ϵ(|a|+T)1δ(Jϵ(xϵ0+δa,uϵ,δ())Jϵ(xϵ0,uϵ())). (5.9)

    We now study variational analysis of (5.9).

    Lemma 5.2. The following result holds:

    supt[0,T]|xϵ,δ(t)xϵ(t)δZϵ(t)|=o(δ),

    where Zϵ is the solution to the first variational equation related to the optimal pair (xϵ0,uϵ())Rn×Up[0,T] given by

    Zϵ(t)=a+t0[fx(t,s,xϵ(s),uϵ(s))(ts)1αZϵ(s)ds+ˆf(t,s)(ts)1α]ds,a.e.t[0,T],

    with

    ˆf(t,s):=f(t,s,xϵ(s),u(s))f(t,s,xϵ(s),uϵ(s)).

    Proof. For δ(0,1), let Zϵ,δ(t):=xϵ,δ(t)xϵ(t)δ, a.e. t[0,T], where based on the Taylor expansion, Zϵ,δ can be written as (note that fϵ,δx(t,s):=10fx(t,s,xϵ(s)+r(xϵ,δ(s)xϵ(s)),uϵ,δ(s))dr)

    Zϵ,δ(t)=a+t0(fϵ,δx(t,s)(ts)1αZϵ,δ(s)+1Eδ(s)δˆf(t,s)(ts)1α)ds,a.e.t[0,T].

    Then we complete the proof by showing that limδ0supt[0,T]|Zϵ,δ(t)Zϵ(t)|=0.

    By Assumption 1, we have |ˆf(t,s)|4K+4K|xϵ(s)|+4Kρ(uϵ(s),u(s))=:˜ψ(s) and |fϵ,δx(t,s)|K. Based on Assumption 1, we can show that

    |xϵ,δ(t)xϵ(t)|b(t)+t0C(ts)1α|xϵ,δ(s)xϵ(s)|ds, (5.10)

    where b(t)=|δa|+t01Eδ(s)˜ψ(s)(ts)1αds. By [34, Lemma 2.3], it follows that

    supt[0,T]|b(t)||δa|+|t01Eδ(s)˜ψ(s)(ts)1αds|C(|δa|+|Eδ|).

    We apply the Gronwall's inequality in [34, Lemma 2.4] to get

    supt[0,T]|xϵ,δ(t)xϵ(t)|C(|δa|+|Eδ|)0,as δ0 due to (5.8). (5.11)

    From Lemma 2.1, Zϵ is a well-defined linear Volterra integral equation. Then with ˆb(t):=|a|+t0˜ψ(s)(ts)1αds and using a similar approach, we can show that

    supt[0,T]|Zϵ(t)|C(|a|+˜ψ()Lp([0,T];R)). (5.12)

    We obtain

    Zϵ,δ(t)Zϵ(t)=t0fϵ,δx(t,s)(ts)1α[Zϵ,δ(s)Zϵ(s)]ds+t0(1Eδ(s)δ1)ˆf(t,s)(ts)1αds+t0fϵ,δx(t,s)fx(t,s,xϵ(s),uϵ(s))(ts)1αZϵ(s)ds. (5.13)

    Notice that by Assumption 1, together with (5.11) and (5.12), it follows that as δ0,

    sup(t,s)Δ|(fϵ,δx(t,s)fx(t,s,xϵ(s),uϵ(s)))Zϵ(s)|C(|δa|+|Eδ|)0. (5.14)

    For convenience, let

    b(1,1)(t):=t0|fϵ,δx(t,s)fx(t,s,xϵ(s),uϵ(s))(ts)1αZϵ(s)|ds

    and

    b(1,2)(t):=|t0(1Eδ(s)δ1)ˆf(t,s)(ts)1αds|.

    By (5.14), limδ0supt[0,T]b(1,1)(t)=0. In addition, by invoking [34, Lemma 5.3.2], it follows that supt[0,T]b(1,2)(t)δ.. Then from (5.13), we have

    |Zϵ,δ(t)Zϵ(t)|b(1,1)(t)+δ+t0Zϵ,δ(s)Zϵ(s)(ts)1αds.

    By applying [34, Lemma 2.4], it follows that

    supt[0,T]|Zϵ,δ(t)Zϵ(t)|C(supt[0,T]b(1,1)(t)+δ).

    Hence, we have

    limδ0supt[0,T]|Zϵ,δ(t)Zϵ(t)|=0.

    This completes the proof.

    Based on the Taylor expansion,

    1δ(J(xϵ0+δa,uϵ,δ())J(xϵ0,uϵ()))=T0lϵ,δx(s)Zϵ,δ(s)ds+T01Eδδˆl(s)ds+hϵ,δx0(T)a+hϵ,δx(T)Zϵ,δ(T),

    where

    ˆl(s):=l(s,xϵ(s),u(s))l(s,xϵ(s),uϵ(s)),
    lϵ,δx(s):=10lx(s,xϵ(s)+r(xϵ,δ(s)xϵ(s)),uϵ,δ(s))dr,
    hϵ,δx0(T):=10hx0(xϵ0+rδa,xϵ(T)+r(xϵ,δ(T)xϵ(T)))dr,

    and

    hϵ,δx(T):=10hx(xϵ0+rδa,xϵ(T)+r(xϵ,δ(T)xϵ(T)))dr.

    Let us define

    ˆZϵ(T)=T0lx(s,xϵ(s),uϵ(s))Zϵ(s)ds+T0ˆl(s)ds+hx0(xϵ0,xϵ(T))a+hx(xϵ0,xϵ(T))Zϵ(T).

    By Lemma 2.1, ˆZϵ is a well-defined Volterra integral equation. It then follows that

    1δ(J(xϵ0+δa,uϵ,δ())J(xϵ0,uϵ()))ˆZϵ(T)=T0lϵ,δx(s)[Zϵ,δ(s)Zϵ(s)]ds+T0[lϵ,δx(s)lx(s,xϵ(s),uϵ(s))]Zϵ(s)ds+T0(1Eδδ1)ˆl(s)ds+[hϵ,δx0(T)hx0(xϵ0,xϵ(T))]a+hϵ,δx(T)[Zϵ,δ(T)Zϵ(T)]+[hϵ,δx(T)hx(xϵ0,xϵ(T))]Zϵ(T).

    Notice that limδ0supt[0,T]|Zϵ,δ(t)Zϵ(t)|=0 by Lemma 5.2. By [33, Corollary 3.9,Chapter 4], it follows that supt[0,T]|t0(1δ1Eδ(s)1)ˆl(s)ds|δ. Hence, similar to Lemma 5.2, we have

    limδ0|1δ(J(xϵ0+δa,uϵ,δ)J(xϵ0,uϵ))ˆZϵ(T)|=0. (5.15)

    From (5.9), we have

    ϵ(|a|+T)1Jϵ(xϵ0+δa,uϵ,δ())+Jϵ(xϵ0,uϵ())×1δ(([J(xϵ0+δa,uϵ,δ())J(¯x0,¯u())+ϵ]+)2([J(xϵ0,uϵ())J(¯x0,¯u())+ϵ]+)2+dF(xϵ0+δa,xϵ,δ(T))2dF(xϵ0,xϵ(T))2). (5.16)

    Let

    λϵ:=[J(xϵ0,uϵ())J(¯x0,¯u())+ϵ]+Jϵ(xϵ0,uϵ())0, (5.17)

    For ξϵR2n (ξϵ1Rn and ξϵ2Rn), let

    ξϵ:=[ξϵ1ξϵ2]:={[xϵ0xϵ(T)]PF(xϵ0,xϵ(T))Jϵ(xϵ0,uϵ())NF(PF(xϵ0,xϵ(T))),(xϵ0,xϵ(T))F,0NF(PF(xϵ0,xϵ(T))),(xϵ0,xϵ(T))F. (5.18)

    Then from (5.15), (5.17), and (5.18), as δ0, (5.16) becomes

    ϵ(|a|+T)λϵˆZϵ(T)+ξϵ1,a+ξϵ2,Zϵ(T), (5.19)

    where by (5.2), (5.17), and (5.18) (see also the discussion in Section 5.1),

    |λϵ|2+|ξϵ|2=1. (5.20)

    Lemma 5.3. For any (a,u())Rn×Up[0,T], the following results hold:

    (i) limϵ0{|xϵ0¯x0|+¯d(uϵ(),¯u())}=0;

    (ii) supt[0,T]|Zϵ(t)Z(t)|=o(1),|ˆZϵ(T)ˆZ(T)|=o(1), where

    Z(t)=a+t0[fx(t,s,¯x(s),¯u(s))(ts)1αZ(s)ds+˜f(t,s)(ts)1α]ds,ˆZ(T)=T0[lx(s,¯x(s),¯u(s))Z(s)+˜l(s)]ds+¯hx0a+¯hxZ(T),

    with ˜f(t,s):=f(t,s,¯x(s),u(s))f(t,s,¯x(s),¯u(s)), ˜l(s):=l(s,¯x(s),u(s))l(s,¯x(s),¯u(s)), ¯hx0:=hx0(¯x0,¯x(T)), and ¯hx:=hx(¯x0,¯x(T)).

    Proof. Part (i) is due to (5.6). The proof of (ii) is similar to Lemma 5.2.

    We consider the limit of ϵ0. Instead of taking the limit with respect to ϵ0, let {ϵk} be the sequence of ϵ such that ϵk0 and ϵk0 as k. We replace ϵ with ϵk. Then by (5.20), the sequences ({λϵk},{ξϵk}) are bounded for k0. From the standard compactness argument, we may extract a subsequence of {ϵk}, still denoted by {ϵk}, such that

    ({λϵk},{ξϵk})(λ0,ξ0)=:(λ,ξ),as k. (5.21)

    By (5.17) and the property of the limiting normal cones [41, page 43],

    λ0,ξNF(PF(¯x0,¯x(T)))=NF(¯x0,¯x(T)). (5.22)

    From (5.21) and (5.20), together with Lemma 5.3, as k, it follows that

    λϵkˆZϵk(T)λˆZ(T)+|ˆZϵk(T)ˆZ(T)|+|λϵkλ|ˆZ(T)λˆZ(T),ξϵk1,a=ξ1,a+ξϵk1,aξ1,aξ1,a,ξϵk2,Zϵk(T)ξ2,Z(T)+|Zϵk(T)Z(T)|+|ξϵk2ξ2||Z(T)|ξ2,Z(T).

    Therefore, as k, (5.19) becomes for any (a,u)Rn×Up[0,T],

    0λˆZ(T)+ξ1,a+ξ2,Z(T). (5.23)

    When λ>0, the nontriviality condition in Theorem 3.1 holds. When λ=0, from (5.20)–(5.22), we must have ξNF(¯x0,¯x(T)) and |ξϵk||ξ|=1. Hence, we have (λ,ξ)0, i.e., they cannot be zero simultaneously. This proves the nontriviality condition in Theorem 3.1.

    Then by using the variational equations in Lemma 5.3, (5.23) becomes

    0ξ1+λ¯hx0,a+ξ2+λ¯hx,a+T0λ˜l(s)ds+T0(λ¯hx+ξ2)1[0,T)(s)fx(T,s,¯x(s),¯u(s))(Ts)1αZ(s)ds+T0[λlx(s,¯x(s),¯u(s))Z(s)+(λ¯hx+ξ2)1[0,T)(s)˜f(T,s)(Ts)1α]ds. (5.24)

    With the adjoint equation p(\cdot) \in L^p([0, T];\mathbb{R}^n) in Theorem 3.1, (5.24) becomes

    \begin{align} 0 & \leq \langle \xi_1 + \lambda \overline{h}_{x_0}^\top, a \rangle + \langle \xi_2 + \lambda \overline{h}_x^\top, a \rangle + \int_0^T \lambda \widetilde{l}(s) \text{d} s \\ &\; \; \; + \int_0^T \Bigl [ - p(s) + \int_s^T \frac{f_x(r,s,\overline{x}(s),\overline{u}(s))^\top}{(r-s)^{1-\alpha}} p(r) \text{d} r \Bigr ]^\top Z(s) \text{d} s \\ &\; \; \; + \int_0^T \Bigl ( \lambda \overline{h}_x + \xi_2^\top \Bigr ) \mathbb{1}_{[0,T)}(s) \frac{\widetilde{f}(T,s)}{(T-s)^{1-\alpha}} \text{d} s. \end{align} (5.25)

    In (5.25), the standard Fubini's formula and Lemma 5.3 lead to

    \begin{align*} & \int_0^T \Bigl [ - p(s) + \int_s^T \frac{f_x(r,s,\overline{x}(s),\overline{u}(s))^\top}{(r-s)^{1-\alpha}} p(r) \text{d} r \Bigr ]^\top Z(s) \text{d} s \\ & = \int_0^T - p(s)^\top \Bigl [ Z(s) - \int_0^s \frac{f_x(s,r,\overline{x}(r),\overline{u}(r))}{(s-r)^{1-\alpha}} Z(r) \text{d} r \Bigr ] \text{d} s \\ & = \int_0^T - p(s)^\top \Bigl [a + \int_0^s \frac{\widetilde{f}(s,r)}{(s-r)^{1-\alpha}} \text{d} r \Bigr ] \text{d} s. \end{align*}

    Moreover, by definition of N_F in (5.1), it follows that for any y = (y_1, y_2) \in F ,

    \begin{align*} \langle \xi_1, a \rangle + \langle \xi_2, a \rangle & \leq \langle \xi_1, \overline{x}_0 - y_1 + a \rangle + \langle \xi_2, \overline{x}(T) - y_2 + a \rangle. \end{align*}

    Hence, (5.25) becomes for any (a, u) \in \mathbb{R}^n \times \mathcal{U}^p[0, T] and y = (y_1, y_2) \in F ,

    \begin{align} 0 & \leq \langle \xi_1, \overline{x}_0 - y_1 + a \rangle + \langle \xi_2, \overline{x}(T) - y_2 + a \rangle + \lambda \overline{h}_{x_0} a + \lambda \overline{h}_x a \\ &\; \; \; - \langle \int_0^T p(s) \text{d} s , a \rangle + \int_0^T - p(s)^\top \int_0^s \frac{\widetilde{f}(s,r)}{(s-r)^{1-\alpha}} \text{d} r \text{d} s \\ &\; \; \; + \int_0^T \Bigl ( \lambda \overline{h}_x + \xi_2^\top \Bigr ) \mathbb{1}_{[0,T)}(s) \frac{\widetilde{f}(T,s)}{(T-s)^{1-\alpha}} \text{d} s + \int_0^T \lambda \widetilde{l}(s) \text{d} s. \end{align} (5.26)

    We use (5.26) to prove the transversality conditions and the Hamiltonian-like maximum condition.

    In (5.26), when u = \overline{u} , for any y = (y_1, y_2) \in F , we have

    \begin{align*} 0 & \leq \langle \xi_1, \overline{x}_0 - y_1 + a \rangle + \langle \xi_2, \overline{x}(T) - y_2 + a \rangle + \lambda \overline{h}_{x_0} a + \lambda \overline{h}_x a - \langle \int_0^T p(s) \text{d} s , a \rangle. \end{align*}

    When (y_1, y_2) = (\overline{x}_0, \overline{x}(T)) , the above inequality holds for any a, -a \in \mathbb{R}^n , which implies

    \begin{align} \int_0^T p(s) \text{d} s = \xi_1 + \xi_2 + \lambda \overline{h}_{x_0}^\top +\lambda \overline{h}_x^\top. \end{align} (5.27)

    Under (5.27), the above inequality becomes 0 \leq \langle \xi_1, \overline{x}_0 - y_1 \rangle + \langle \xi_2, \overline{x}(T) - y_2 \rangle for any y = (y_1, y_2) \in F . This, together with (5.27), proves the transversality conditions in Theorem 3.1. In addition, as p (\cdot) \in L^p([0, T];\mathbb{R}^n) , the nontriviality condition implies the nontriviality of the adjoint equation in Theorem 3.1.

    We prove the Hamiltonian-like maximum condition in Theorem 3.1. Define

    \begin{align*} \Lambda(s,\overline{x}(s),u) & : = - \lambda l(s,\overline{x}(s),u) + \int_s^T p(r)^\top \frac{f(r,s,\overline{x}(s),u)}{(r-s)^{1-\alpha}} \text{d} r \\ &\; \; \; - \mathbb{1}_{[0,T)}(s) \Bigl ( \lambda \overline{h}_x + \xi_2^\top \Bigr ) \frac{f(T,s,\overline{x}(s),u)}{(T-s)^{1-\alpha}}. \end{align*}

    When (y_1, y_2) = (\overline{x}_0, \overline{x}(T)) and a = 0 in (5.26), by the Fubini's formula, (5.26) can be written as

    \begin{align} \int_0^T \Lambda(s,\overline{x}(s),u(s)) \text{d} s \leq \int_0^T \Lambda(s,\overline{x}(s),\overline{u}(s)) \text{d} s. \end{align} (5.28)

    As U is separable, there exists a countable dense set U_0 = \{u_i, \; i \geq 1\} \subset U . Moreover, there exists a measurable set S_i \subset [0, T] such that |S_i| = T and any t \in S_i is the Lebesgue point of \Lambda(t, \overline{x}(t), u(t)) , i.e., it holds that \lim_{\tau \downarrow 0} \frac{1}{2\tau} \int_{t-\tau}^{t+\tau} \Lambda(s, \overline{x}(s), u(s)) \text{d} s = \Lambda(t, \overline{x}(t), u(t)) [10, Theorem 5.6.2]. We fix u_i \in U_0 . For any t \in S_i , define

    \begin{align*} u(s) : = \begin{cases} \overline{u}(s), & s \in [0,T] \setminus (t-\tau,t+\tau),\\ u_i, & s \in (t-\tau,t+\tau). \end{cases} \end{align*}

    Then (5.28) becomes

    \begin{align*} 0 & \leq \lim\limits_{\tau \downarrow 0} \frac{1}{2\tau} \int_{t - \tau}^{t + \tau} [ \Lambda(s,\overline{x}(s),\overline{u}(s)) - \Lambda(s,\overline{x}(s),u_i) ] \text{d} s = \Lambda (t,\overline{x}(t),\overline{u}(t)) - \Lambda (t,\overline{x}(t),u_i). \end{align*}

    This is implies that

    \begin{align} \Lambda(t,\overline{x}(t),u_i) \leq \Lambda(t,\overline{x}(t),\overline{u}(t)),\; \forall u_i \in U_0,\; t \in \cap_{i \geq 1} S_i. \end{align} (5.29)

    Since \cap_{i \geq 1} S_i = [0, T] by the fact that U_0 is countable, \Lambda is continuous in u \in U , and U is separable, (5.29) implies the Hamiltonian-like maximum condition in Theorem 3.1. This completes the proof.

    In this paper, we have studied the Pontryagin maximum principle for the optimal control problem of Volterra integral equations with singular kernels. The main technical difference in the proof compared with the existing literature is the variational and duality analysis. In fact, the variational analysis in our paper needs to handle the singular kernel, and the duality analysis requires to characterize the integral-type adjoint equation with the singular kernel to get the desired Hamiltonian-like maximization condition.

    There are several important potential future research problems. One {direction} is to study (P) with additional generalized state constraints including the running state constraint. Another {direction} is to study the numerical aspects of solving (P) such as the shooting method and the Lagrangian collocation approach. Finally, one can extend (P) to the infinite-dimensional problem, which can be applied to studying optimal control of integral-type partial differential equations and integral-type systems with delays and singular kernels.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This research was supported in part by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) and the Ministry of Trade, Industry & Energy (MOTIE) of the Republic of Korea (RS-2023-00235742), in part by the National Research Foundation of Korea (NRF) Grant funded by the Ministry of Science and ICT, South Korea (NRF-2021R1A2C2094350) and in part by Institute of Information & communications Technology Planning & Evaluation (IITP) grant funded by the Korea government (MSIT) (No. 2020-0-01373, Artificial Intelligence Graduate School Program (Hanyang University)).

    The author declares no conflict of interest in this paper.

    [1] Kara B, Çelik A (2015) The relationship between risk for eating disorder and health-related quality of life in patients with multiple sclerosis. Gulhane Med J 57: 36-40.
    [2] Bøe Lunde HM, Aae TF, Indrevåg W, et al. (2012) Poor sleep in patients with multiple sclerosis. PLoS One 7: e49996.
    [3] Čarnická Z, Kollár B, Šiarnik P, et al. (2015) Sleep disorders in patients with multiple sclerosis. J Clin Sleep Med 11: 553-557.
    [4] Merlino G, Fratticci L, Lenchig C, et al. (2009) Prevalence of 'poor sleep' among patients with multiple sclerosis: an independent predictor of mental and physical status. Sleep Med 10: 26-34. doi: 10.1016/j.sleep.2007.11.004
    [5] Nagaraj K, Taly AB, Gupta A, et al. (2013) Depression and sleep disturbances in patients with multiple sclerosis and correlation with associated fatigue. J Neurosci Rural Pract 4: 387-391.
    [6] Vitkova M, Gdovinova Z, Rosenberger J, et al. (2014) Factors associated with poor sleep quality in patients with multiple sclerosis differ by disease duration. Disabil Health J 7: 466-471. doi: 10.1016/j.dhjo.2014.05.004
    [7] Sarraf P, Azizi S, Moghaddasi AN, et al. (2014) Relationship between sleep quality and quality of life in patients with multiple sclerosis. Int J Prev Med 5: 1582-1586.
    [8] Cameron MH, Peterson V, Boudreau EA, et al. (2014) Fatigue is associated with poor sleep in people with multiple sclerosis and cognitive impairment. Mult Scler Int 2014: 872732.
    [9] Najafi MR, Toghianifar N, Etemadifar M, et al. (2013) Circadian rhythm sleep disorders in patients with multiple sclerosis and its association with fatigue: a case-control study. J Res Med Sci 18: S71-73.
    [10] Motl RW, Weikert M, Suh Y., et al. (2010) Symptom cluster and physical activity in relapsing-remitting multiple sclerosis. Res Nurs Health 33: 398-412.
    [11] Pilutti LA, Dlugonski D, Sandroff BM, et al. (2014) Randomized controlled trial of a behavioral intervention targeting symptoms and physical activity in multiple sclerosis. Mult Scler 20: 594-601.
    [12] Kara B, Kılıç Ö. (2015) Predictors of poor sleep quality and excessive daytime sleepiness in Turkish adults with type 2 diabetes. J Clin Nurs 24: 1436-1439. doi: 10.1111/jocn.12710
    [13] Kara B, Tenekeci EG. (2015) Sleep quality and associated factors in older Turkish adults with hypertension: a pilot study. J Transcult Nurs (in press) doi: 10.1177/1043659615623330
    [14] Murat S, Ali U, Serdal K, et al. (2015) Assessment of subjective sleep quality in iron deficiency anaemia. Afr Health Sci 15: 621-627.
    [15] Pai MF, Hsu SP, Yang SY, et al. (2007) Sleep disturbance in chronic hemodialysis patients: the impact of depression and anemia. Ren Fail 29: 673-677. doi: 10.1080/08860220701459642
    [16] World Health Organization. (2011) Haemoglobin concentrations for the diagnosis of anemia and assessment of severity: vitamin and mineral nutrition information system. (Cited 27 August 2015). Available from: www.who.int/vmnis/indicators/haemoglobin.pdf
    [17] Ağargün MY, Kara H, Anlar Ö. (1996) Reliability and validity of the Pittsburgh Sleep Quality Index. Turk Psikiyatri Derg 7: 2-4 (in Turkish).
    [18] Yurtsever S, Beduk T. (2003) Evaluation of fatigue on hemodialysis patients. HEMAR-G 5: 3-12 (in Turkish).
    [19] Saglam M, Arikan H, Savci S, et al. (2010) International Physical Activity Questionnaire: reliability and validity of the Turkish version. Percept Motor Skill 111: 278-284. doi: 10.2466/06.08.PMS.111.4.278-284
    [20] Kara B, İşcan B. (2016) Predictors of health behaviors in Turkish female nursing students. Asian Nurs Res 10: 75-81.
  • This article has been cited by:

    1. Faïçal Ndaïrou, Delfim F. M. Torres, Pontryagin Maximum Principle for Incommensurate Fractional-Orders Optimal Control Problems, 2023, 11, 2227-7390, 4218, 10.3390/math11194218
    2. Yuna Oh, Jun Moon, The infinite-dimensional Pontryagin maximum principle for optimal control problems of fractional evolution equations with endpoint state constraints, 2024, 9, 2473-6988, 6109, 10.3934/math.2024299
    3. Andres Cardenas, Sergio Pulido, Rafael Serrano, Existence of optimal controls for stochastic Volterra equations, 2025, 31, 1292-8119, 30, 10.1051/cocv/2025020
  • Reader Comments
  • © 2016 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5003) PDF downloads(1236) Cited by(4)

Figures and Tables

Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog