Loading [MathJax]/jax/element/mml/optable/GreekAndCoptic.js
Research article Special Issues

A direct integral pseudospectral method for solving a class of infinite-horizon optimal control problems using Gegenbauer polynomials and certain parametric maps

  • We present a novel direct integral pseudospectral (PS) method (a direct IPS method) for solving a class of continuous-time infinite-horizon optimal control problems (IHOCs). The method transforms the IHOCs into finite-horizon optimal control problems in their integral forms by means of certain parametric mappings, which are then approximated by finite-dimensional nonlinear programming problems (NLPs) through rational collocations based on Gegenbauer polynomials and Gegenbauer-Gauss-Radau (GGR) points. The paper also analyzes the interplay between the parametric maps, barycentric rational collocations based on Gegenbauer polynomials and GGR points and the convergence properties of the collocated solutions for IHOCs. Some novel formulas for the construction of the rational interpolation weights and the GGR-based integration and differentiation matrices in barycentric-trigonometric forms are derived. A rigorous study on the error and convergence of the proposed method is presented. A stability analysis based on the Lebesgue constant for GGR-based rational interpolation is investigated. Two easy-to-implement pseudocodes of computational algorithms for computing the barycentric-trigonometric rational weights are described. Three illustrative test examples are presented to support the theoretical results. We show that the proposed collocation method leveraged with a fast and accurate NLP solver converges exponentially to near-optimal approximations for a coarse collocation mesh grid size. The paper also shows that typical direct spectral/PS and IPS methods based on classical Jacobi polynomials and certain parametric maps usually diverge as the number of collocation points grow large if the computations are carried out using floating-point arithmetic and the discretizations use a single mesh grid, regardless of whether they are of Gauss/Gauss-Radau type or equally spaced.

    Citation: Kareem T. Elgindy, Hareth M. Refat. A direct integral pseudospectral method for solving a class of infinite-horizon optimal control problems using Gegenbauer polynomials and certain parametric maps[J]. AIMS Mathematics, 2023, 8(2): 3561-3605. doi: 10.3934/math.2023181

    Related Papers:

    [1] Yousef Jawarneh, Humaira Yasmin, M. Mossa Al-Sawalha, Rasool Shah, Asfandyar Khan . Fractional comparative analysis of Camassa-Holm and Degasperis-Procesi equations. AIMS Mathematics, 2023, 8(11): 25845-25862. doi: 10.3934/math.20231318
    [2] Zheng Dou, Kexin Luo . Global weak solutions of nonlinear rotation-Camassa-Holm model. AIMS Mathematics, 2023, 8(7): 15285-15298. doi: 10.3934/math.2023781
    [3] Yunxi Guo, Ying Wang . The Cauchy problem to a gkCH equation with peakon solutions. AIMS Mathematics, 2022, 7(7): 12781-12801. doi: 10.3934/math.2022707
    [4] A. K. M. Kazi Sazzad Hossain, M. Ali Akbar . Solitary wave solutions of few nonlinear evolution equations. AIMS Mathematics, 2020, 5(2): 1199-1215. doi: 10.3934/math.2020083
    [5] Mahmoud A. E. Abdelrahman, S. Z. Hassan, R. A. Alomair, D. M. Alsaleh . Fundamental solutions for the conformable time fractional Phi-4 and space-time fractional simplified MCH equations. AIMS Mathematics, 2021, 6(6): 6555-6568. doi: 10.3934/math.2021386
    [6] Zhe Ji, Yifan Nie, Lingfei Li, Yingying Xie, Mancang Wang . Rational solutions of an extended (2+1)-dimensional Camassa-Holm- Kadomtsev-Petviashvili equation in liquid drop. AIMS Mathematics, 2023, 8(2): 3163-3184. doi: 10.3934/math.2023162
    [7] M. Ali Akbar, Norhashidah Hj. Mohd. Ali, M. Tarikul Islam . Multiple closed form solutions to some fractional order nonlinear evolution equations in physics and plasma physics. AIMS Mathematics, 2019, 4(3): 397-411. doi: 10.3934/math.2019.3.397
    [8] Umair Ali, Sanaullah Mastoi, Wan Ainun Mior Othman, Mostafa M. A Khater, Muhammad Sohail . Computation of traveling wave solution for nonlinear variable-order fractional model of modified equal width equation. AIMS Mathematics, 2021, 6(9): 10055-10069. doi: 10.3934/math.2021584
    [9] Ying Wang, Yunxi Guo . Blow-up solution and analyticity to a generalized Camassa-Holm equation. AIMS Mathematics, 2023, 8(5): 10728-10744. doi: 10.3934/math.2023544
    [10] M. TarikulIslam, M. AliAkbar, M. Abul Kalam Azad . Traveling wave solutions in closed form for some nonlinear fractional evolution equations related to conformable fractional derivative. AIMS Mathematics, 2018, 3(4): 625-646. doi: 10.3934/Math.2018.4.625
  • We present a novel direct integral pseudospectral (PS) method (a direct IPS method) for solving a class of continuous-time infinite-horizon optimal control problems (IHOCs). The method transforms the IHOCs into finite-horizon optimal control problems in their integral forms by means of certain parametric mappings, which are then approximated by finite-dimensional nonlinear programming problems (NLPs) through rational collocations based on Gegenbauer polynomials and Gegenbauer-Gauss-Radau (GGR) points. The paper also analyzes the interplay between the parametric maps, barycentric rational collocations based on Gegenbauer polynomials and GGR points and the convergence properties of the collocated solutions for IHOCs. Some novel formulas for the construction of the rational interpolation weights and the GGR-based integration and differentiation matrices in barycentric-trigonometric forms are derived. A rigorous study on the error and convergence of the proposed method is presented. A stability analysis based on the Lebesgue constant for GGR-based rational interpolation is investigated. Two easy-to-implement pseudocodes of computational algorithms for computing the barycentric-trigonometric rational weights are described. Three illustrative test examples are presented to support the theoretical results. We show that the proposed collocation method leveraged with a fast and accurate NLP solver converges exponentially to near-optimal approximations for a coarse collocation mesh grid size. The paper also shows that typical direct spectral/PS and IPS methods based on classical Jacobi polynomials and certain parametric maps usually diverge as the number of collocation points grow large if the computations are carried out using floating-point arithmetic and the discretizations use a single mesh grid, regardless of whether they are of Gauss/Gauss-Radau type or equally spaced.



    Nonlinear evolution equations have much achieved importance and recognition recently. The differential equations have proved to be valuable tools to the modeling of many physical phenomena. The study of the traveling wave solutions for nonlinear differential equations (NLDEs) play an important role in various types of fields, such as, quantum mechanics, electricity, plasma physics, chemical kinematics, optical fibers, biological model, electromagnetic field, viscoelasticity, electrochemistry, physics, signal processing, optical fibers, control theory, fluid mechanics and population model [1,2,3,4], etc.

    Now-a-days NLDEs have been creating important opportunity for the researchers. As a result, in the past few years different efficient and significant methods have been established and enhanced to obtain exact solutions to evolution equations, such as, the (G/G)-expansion method [5,6,7,8], the variational iteration method [9,10,11], the modified decomposition method [12,13], the Jacobi elliptic function method [14], the modified Kudryashov method [15], the fractional sub-equation method [16], the homotopy perturbation method [17,18], the MSE method [19,20,21,22], the tanh-function method [23], the finite difference method [24], the exp-function method [25,26,27], the F-expansion method [28,29], the operational method [30], the sine-cosine method [31], the first integral method [32], etc.

    In 2006, Wazwaz [33] studies a family of important physical equations which is known as modified α - equation and he provides the following form of that modified α - equation:

    utuxxt+(α+1)u2uxαuxuxxuuxxx=0, (1.1)

    where, α is a positive integer. By taking α=2, Wazwaz [33] reduces Eq. (1.1) to the following modified Camassa-Holm (MCH) equation:

    utuxxt+3u2ux2uxuxxuuxxx=0. (1.2)

    Simplifying the MCH Eq. (1.2), we attain the following equation which is known as simplified modified Camassa-Holm (SMCH) equation with involving parameters kandβ :

    ut+2kuxuxxt+βu2ux=0.

    The study of the SMCH equation has explained many significant and sensible ambiences. In present literature, the SMCH equation has investigated through the exp(ϕ(η))-expansion method [34], the exp-function method [35], the improved (G/G)-expansion method [36], etc.

    The aim of this article is to examine the fresh, useful and more general traveling wave solutions to the SMCH equation by aid of the MSE method. The MSE method is a recently developed effective, straightforward, powerful and rising method to investigate traveling wave solutions to the NLDEs. We also have shown the graphical representation of the attained solutions of its definite values of the involved parameters through depicting 3D and 2D figures.

    The rest of the article is patterned as follows: In Section 2, we have described the MSE method. Section 3 is used to determine the fresh exact traveling wave solutions to the equation. In Section 4, we have shown the graphical representations and physical significance and in Section 5, we have represented the comparison and discussed the results. Finally, in Section 6, we have concluded this article.

    Now, we expressed the MSE method which is used by many researchers to observe the solitary wave solutions to nonlinear evolution equations (NLEEs). Let us consider the nonlinear evolution equations (NLEE) in two independent variables x and t, in the following form:

    R(u,ut,ux,uxx,utt,)=0 (2.1)

    where u=u(x,t) is an unknown function, R is a polynomial of u(x,t) and its derivatives in which the highest order derivatives and nonlinear terms are involved and its partial derivatives is ut=ut,ux=ux,utt=2ut2,uxx=2ux2.

    To solve the NLEE (2.1) by aid of the MSE method, we represent the main steps in the following form:

    Step 1: We consider the real variables x and t by the traveling wave variable ξ, in the following form:

    ξ=xwt;u(x,t)=F(ξ) (2.2)

    where w denotes the velocity of the traveling wave.

    By using the traveling wave transformation Eq. (2.2) in Eq. (2.1) then reduces an ordinary differential equation (ODE) written in the form:

    G(F,F,F,) (2.3)

    where G is a polynomial of F(ξ) and its derivatives while the highest order derivatives and nonlinear terms are connected and the superscripts indicate the ordinary derivatives with respect to ξ, in which F(ξ)=dFdξ.

    Step 2: We assume that the Eq. (2.3) has the traveling wave solution in the following form:

    F(ξ)=Nj=0aj{s(ξ)s(ξ)}j, (2.4)

    where N be the positive integer, aj(j=0,1,2,,N) are unknown constants to be determined, in which aN0, and the unknown function s(ξ) to be evaluated, in which s(ξ)=dsdξ0.

    Step 3: The positive integer N showing in Eq. (2.4) can be determined by considering the homogeneous balance between the highest order derivatives and the nonlinear terms appearing in Eq. (2.3).

    Moreover, if the degree of F(ξ) is deg[F(ξ)]=N, therefore, the degree of any other expressions will be as follows:

    deg[dmF(ξ)dξm]=N+M and deg[Fm{dqF(ξ)dξq}p]=mN+p(N+q) (2.5)

    Step 4: Substitute Eq. (2.4) into Eq. (2.3), we calculate all the necessary derivatives F,F,F and then estimate the function s(ξ). Accordingly of this substitution, we get a polynomial of (1s(ξ))j,j= 0.1.2. In the resultant polynomial, we equate all the coefficients of (1s(ξ))j,j=0.1.2. to zero, which yields a set of algebraic and differential equations for aj(j=0,1,2,,N), s(ξ) and other needful parameters.

    Step 5: The solutions of the set of algebraic and differential equations obtained in Step 4, provide the values of s(ξ) and other needful parameters.

    In this section, we examine the several new travelling wave solutions to the SMCH equation by the implementation of the MSE method. Let us consider the equation in of the form:

    ut+2kuxuxxt+βu2ux=0, (3.1)

    where β>0 and k is arbitrary constant. It arises as a description of gravity water waves in the long-wave regime, water wave mechanics, turbulent motion, etc. Using the traveling wave transformation Eq. (2.2), the Eq. (3.1) converts into a nonlinear ODE

    wu+2ku+wu+βu2u=0, (3.2)

    Integrating Eq. (3.2) once and choose the integrating constant to zero, yields

    (2kw)u+wu+β3u3=0, (3.3)

    Balancing the linear term of the highest order derivative u and the nonlinear term of the highest order u3 occurring in Eq. (3.3), gives N=1. Therefore, the shape of the solution of Eq. (3.3) becomes:

    u=a0+a1(s(ξ)s(ξ)), (3.4)

    where a0 and a1 are constants to be determined, such that a10 and s(ξ) is an unknown function to be determined. Substituting the values of u,u3 and u using Eq. (3.4) into Eq. (3.3) and then equating the coefficients of s0,s1,s2,s3 to zero, these become

    2ka0wa0+13βa30=0 (3.5)
    2ka1swa1s+βa20a1s+wa1s=0 (3.6)
    βa0a21(s)23wa1ss=0 (3.7)
    2wa1(s)3+13βa31(s)3=0 (3.8)

    From Eq. (3.5) and Eq. (3.8), we determine

    a0=0,a0=±3(w2k)β and a1=±i6wβ, since a10.

    Using (3.6) and (3.7), we also obtain

    s(ξ)=c3+β2a21a20c1e3(w2kβa20)ξβa0a19(w2kβa20)2 (3.9)

    Case 1:

    When a0=0 then the solution is discarded.

    Case 2:

    When a0=3(w2k)β,a1=±i6wβ and s(ξ)=c3+c1β2a20a21e3(2k+wθa20)ξβa0a19(w2kβa20)2, the solution (3.4) becomes below:

    u(ξ)=3(w2k)β23c1w(w2k)e±i(4k2w)ξ2w(w2k)2(w2k)β(c32c1w(w2k)e±i(4k2w)ξ2w(w2k)(2k+w3(2k+w))2) (3.10)

    The solution in Eq. (3.10) is useful and further general travelling wave solution which does not exist in the previous literature. If we set different special values of c1 and c3 consistently further numerous solutions can be determined through this method, but for simplicity the remaining solutions have not been written down here.

    When w>2k, using the obtained results a0,a1 and s(ξ) and if we choose, c1 = b k(k+m) and c3=b2k2m2 and also by means of the computer algebra like Mathematica, converting the solution (3.10) from exponential to trigonometric function, we achieve the solution is simplified as the following soliton form:

    u(ξ)=±i3(w2k)βcot((w2k)ξ2k) (3.11)

    Also, when w>2k if we choose, c1 = b k(k+m) and c3=b2k2m2, the above solution (3.2) reduces to the following soliton:

    u(ξ)=±i3(w2k)βtan((w2k)ξ2k) (3.12)

    Solutions in Eq. (3.11) and Eq. (3.12) are the well-known periodic wave solutions.

    Furthermore, when w<2k, using the obtained results a0,a1 and s(ξ) and if we choose, c1 = b k(k+m) and c3=b2k2m2 and also by means of the computer algebra like Mathematica, converting the solution in Eq. (3.10) from exponential to trigonometric function, we achieve the solution is simplified as the following soliton form:

    u(ξ)=±3(2kw)βcoth((2kw)ξ2k) (3.13)

    Also, for w<2k if we choose, c1 = b k(k+m) and c3=b2k2m2, the above solution in Eq. (3.10) reduces to the following soliton:

    u(ξ)=±3(2kw)βtanh((2kw)ξ2k) (3.14)

    Solutions in Eq. (3.13) and Eq. (3.14) are the well-known singular periodic wave solutions.

    Case 3:

    When a0=3(w2k)β,a1=±i6wβ and s(ξ)=c3+c1β2a20a21e3(2k+wθa20)ξβa0a19(w2kβa20)2, the solution (3.4) reduces:

    u(ξ)=3(w2k)β+23c1w(w2k)e±i(2(w2k))ξ2w(w2k)2(w2k)β(c32c1w(w2k)e±i(2(w2k))ξ2w(w2k)(2k+w3(2k+w))2) (3.15)

    The solution in Eq. (3.15) is useful and further general travelling wave solution which does not exist in the previous literature. If we set different special values of c1 and c3 consistently further numerous solutions can be determined through this method, but for simplicity the remaining solutions have not been written down here.

    When w>2k using the obtained results a0,a1 and s(ξ) and if we choose, c1 = b k(k+m) and c3=b2k2m2 and also by means of the computer algebra like Mathematica, converting the solution in Eq. (3.15) from exponential to trigonometric function, we achieve the solution is simplified as the following soliton form:

    u(ξ)=±i3(w2k)βcot((w2k)ξ2k) (3.16)

    And, when w>k if we choose, c1 = b k(k+m) and c3=b2k2m2, the above solution in Eq. (3.15) reduces to the following soliton:

    u(ξ)=±i3(w2k)βtan((w2k)ξ2k) (3.17)

    Solutions in Eq. (3.16) and in Eq. (3.17) are the well-known periodic wave solutions.

    On the other hand, when w<k using the obtained results a0,a1 and s(ξ) and if we choose, c1 = b k(k+m) and c3=b2k2m2 and also by means of the computer algebra like Mathematica, converting the solution in Eq. (3.15) from exponential to trigonometric function, we achieve the solution is simplified as the following soliton form:

    u(ξ)=±3(2kw)βcoth((2kw)ξ2k) (3.18)

    Also, when w<k and if we choose, c1 = b k(k+m) and c3=b2k2m2, the above solution in Eq. (3.15) reduces to the following soliton:

    u(ξ)=±3(2kw)βtanh((2kw)ξ2k) (3.19)

    Solutions in Eq. (3.18) and in Eq. (3.19) are the well-known singular periodic wave solutions.

    Therefore, summarizing the solutions between Eqs. (3.11), (3.12), (3.16), (3.17) and also summarizing between Eqs. (3.14), (3.13), (3.18), (3.19) we attain the following subsequent solutions:

    u(ξ)=±i3(w2k)βcot((w2k)ξ2k) (3.20)
    u(ξ)=±i3(w2k)βtan((w2k)ξ2k) (3.21)
    u(ξ)=±3(2kw)βcoth((2kw)ξ2k) (3.22)
    u(ξ)=±3(2kw)βtanh((2kw)ξ2k) (3.23)

    Now making use of the wave variable of Eq. (2.2), the solutions in Eqs. (3.20), (3.21), (3.22) and (3.23) become as follows:

    u(1,2)(ξ)=±i3(w2k)βcot((w2k)ξ2k) (3.24)
    u(3,4)(ξ)=±i3(w2k)βtan((w2k)ξ2k) (3.25)
    u(5,6)(ξ)=±3(2kw)βcoth((2kw)ξ2k) (3.26)
    u(7,8)(ξ)=±3(2kw)βtanh((2kw)ξ2k) (3.27)

    The solutions in Eqs. (3.24), (3.25), (3.26) and (3.27) are fresh, important and further general to the SMCH equation can be extracted.

    Therefore, it is remarkable to observe that the obtained traveling wave solutions arise to be suitable to seek for gravity water wave in the long-wave regime, the water wave mechanics, the turbulent motion, the driving a flow of fluid and singular periodic waves, etc.

    In this section, we show the graphical representations of the obtained solutions and discuss the physical significance to the MSCH equation is arranged below:

    The shapes of the figures of the obtained solutions in Eq. (3.24) (positive and negative values) for the definite values of the parameters are given below:

    Figure 1.  3D and 2D modulus plot of solution in Eq. (3.24) which is periodic wave when w=1 and k = 1.5, {\rm{ \mathsf{ β} }} = 1 within the interval 0 \le x, t \le 10 for 3D and t = 0 for 2D.
    Figure 2.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{5}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 3.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{15}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 4.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{25}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 5.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{30}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 6.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{35}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 7.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{40}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 8.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{50}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 9.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{60}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 10.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{65}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 11.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{70}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 12.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{75}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.
    Figure 13.  3D and 2D modulus plot of solution in Eq. (3.24) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{80}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1, \alpha = 0.5 within the interval 0 \le x, t \le 4 for 3D and t = 0 for 2D.

    The shape of the figures of the obtained solution in Eq. (3.27) (positive and negative values) for the definite values of the parameters is shown below:

    Figure 14.  3D and 2D modulus plot of solution in Eq. (3.27) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{1}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1 within the interval - 10 \le x, t \le 10 for 3D and t = 0 for 2D.
    Figure 15.  3D and 2D modulus plot of solution in Eq. (3.27) when \mathit{\boldsymbol{w}}{\rm{ = }}{\bf{2}}, k = 1.5, {\rm{ \mathsf{ β} }} = 1 within the interval - 10 \le x, t \le 10 for 3D and t = 0 for 2D.
    Figure 16.  3D and 2D modulus plot of solution in Eq. (3.27) when \mathit{\boldsymbol{w}} = {\bf{3}}, \mathit{\boldsymbol{k}} \ne {\bf{1}}{\bf{.5}}, {\rm{ \mathsf{ β} }} = 1 within the interval - 10 \le x, t \le 10 for 3D and t = 0 for 2D.
    Figure 17.  3D and 2D modulus plot of solution in Eq. (3.27) which is singular periodic wave when \mathit{\boldsymbol{w}} = {\bf{5}}, \mathit{\boldsymbol{k}} = {\bf{1}}{\bf{.5}}, {\rm{ \mathsf{ β} }} = 1 within the interval - 10 \le x, t \le 10 for 3D and t = 0 for 2D.
    Figure 18.  3D and 2D modulus plot of solution in Eq. (3.27) which is singular periodic wave when \mathit{\boldsymbol{w}} = {\bf{7}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1 within the interval - 10 \le x, t \le 10 for 3D and t = 0 for 2D.
    Figure 19.  3D and 2D modulus plot of solution in Eq. (3.27)) which is singular periodic wave when \mathit{\boldsymbol{w}} = {\bf{35}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1 within the interval - 10 \le x, t \le 10 for 3D and t = 0 for 2D.
    Figure 20.  3D and 2D modulus plot of solution in Eq. (3.27) which is singular periodic wave when \mathit{\boldsymbol{w}} = {\bf{70}} and k = 1.5, {\rm{ \mathsf{ β} }} = 1 within the interval - 10 \le x, t \le 10 for 3D and t = 0 for 2D.

    In this sub-section, we have described the physical significance of the accomplished solutions through figures. We have obtained eight solutions where both the solutions in Eq. (3.24) and in Eq. (3.25) contain two solutions (positive and negative sign) which are in the trigonometry form and both the solutions in Eq. (3.27) and in Eq. (3.27) also contains two solutions (positive and negative sign) which are in the trigonometric hyperbolic form.

    From the above graphical representations, we may assert that the solutions might be useful to inspect the physical significance in order to know the internal mechanisms of the related physical phenomena. We find different shapes of figures for the solutions in Eq. (3.24) and in Eq. (3.25), within intervals 0 \leq x, t \leq 10 and for the solutions in Eq. (3.26) and in Eq. (3.27) within interval -10 \leq x, t \leq 10 for the different definite values of involved parameters.

    Plotting the figures of the solutions in Eq. (3.24) and in Eq. (3.25, ) we have observed that, for the different values of the velocity w of the traveling wave, the figures of the obtained solutions become different shapes. Plotting the figures of the solutions in Eq. (3.24), we watch that when w = 1, the shape of the wave length is large and when w = 5, the shape of the wave length is small. Consequently when w = 15, the shape of the wave length larger than w = 5 and again when w = 20, the shape of the wave length is decreased. It is important to notice that when w = 35, the shape of the wave length becomes maximum i.e., larger than other values of the velocity of the travelling wave w. We also notice that for larger values of w = 35, the shape of the wave lengths decreasing and again increasing and when velocity arises to w = 70 then the shape of the figure becomes like w = 35. Furthermore, when larger values than w = 70, the shapes of the wave lengths decreasing and again increasing and when velocity arises to w = 105 then the shape of the figures become same to w = 70.

    From the above discussion we have claimed that, the solutions in Eq. (3.24) to the MSCH equation are periodic wave solution and its period is w = 35.

    Similarly, plotting the figures of the solutions in Eq. (3.25), we have found the same figures (periodic wave) like the shape of the figures in Eq. (3.24) and the same situation arises and also observed that its period is w = 35. Therefore, for minimalism, the shapes of the figures of the solutions in Eq. (3.25) have not shown here.

    Furthermore, plotting the figures in the case of the solutions in Eq. (3.26) and in Eq. (3.27) we have observed that, when w < 3, k = 1.5, {\rm{\beta }} = 1 and w = 3, k\ne 1.5, {\rm{\beta }} = 1 and within interval -10\le x, t\le 10 the shapes of the figures are not singular periodic wave. But for other definite values of the including parameters within interval -10\le x, t\le 10, the shapes of the figures are singular periodic wave and for simplicity, only we have depicted the shapes of the figures of the solutions in Eq. (3.27) and the shapes of the figures of the solutions have omitted here.

    In this section, we notice that, the solutions investigated by the exp (-\phi (\eta)) -expansion method, Ali et al. [34], obtained only five (05) wave solutions which were ordinary solutions (see Appendix-1) for the SMCH equation. But by using the MSE method, we have obtained eight (08) wave solutions which have not been found by the exp (-\phi (\eta)) -expansion method which are useful, efficient and more general and give well known shapes like periodic wave solutions and singular periodic wave solutions. The different choices of the integral constants from solutions in Eq. (3.10) and in Eq. (3.15), it might be obtained huge amount of exact wave solutions. Moreover, the attained solutions might be useful to analyze the physical significance. For the definite values of the parameters, different types of periodic and singular periodic solutions are derived in our article.

    Therefore, comparing between the obtained solutions and the solutions obtained by the exp (-\phi (\eta)) -expansion method [34], we might claim that our obtained solutions are more general and useful.

    In this article, we have examined the new, significant and further general travelling wave solutions which are in the form of trigonometric and hyperbolic function to the SMCH equation by means of the efficient technique known as MSE method. The attained solutions are in general form and the definite values of the included parameters yield diverse known soliton solutions and the solutions might be useful to analyze many potential applications in mathematical physics and engineering. Furthermore, the MSE method has established better solutions than exp (-\phi (\eta)) -expansion method. We also have shown the graphical representation of the obtained solutions through depicting figures in order to interpret the physical phenomena of the tangible incidents. The established results have shown that the MSE method is straightforward, further general, efficient and more powerful and can be used to obtain exact traveling wave solutions for many other nonlinear differential equations.

    We would like to thanks the honorable reviewers for their valuable comments and suggestions, which are really helpful to enrich the quality of our paper. Furthermore, we are grateful to the journal authority for their proper judgements and kind considerations.

    The authors declare that they have no any competing interests.

    In this section, we have shown the list of the solutions to the simplified modified Camassa-Holm (SMCH) equation investigated by the exp (-\phi (\eta)) method (Ali et al., [34]) solutions, are arranged below:

    \begin{array}{l}{u_{1}(\eta) = a_{0}-\frac{2 \sqrt{6} \mu}{\left(\sqrt{\lambda^{2}-4 \mu} tanh \left(\frac{\sqrt{\lambda^{2}-4 \mu}}{2}\left(\eta+c_{1}\right)-\lambda\right)\right)}} \end{array} (A1.1)
    {u_{2}(\eta) = a_{0}-\frac{2 \sqrt{6} \mu}{\left(\sqrt{-\lambda^{2}+4 \mu} tan \left(\frac{\sqrt{-\lambda^{2}+4 \mu}}{2}\left(\eta+c_{1}\right)-\lambda\right)\right)}} (A1.2)
    \begin{array}{c}{u_{3}(\eta) = a_{0}-\frac{\sqrt{6} \lambda}{exp \left(\left(\eta+c_{1}\right)^{\lambda}-1\right)}} \end{array} (A1.3)
    {u_{4}(\eta) = a_{0}-\frac{\sqrt{6}\left(\eta+c_{1}\right) \lambda^{2}}{\left(2\left(\eta+c_{1}\right)^{\lambda}+2\right)}} (A1.4)
    {u_{5}(\eta) = a_{0}-\frac{\sqrt{6}}{\left(\eta+c_{1}\right)}} (A1.5)

    where \eta = x-v t.



    [1] S. A. Orszag, Accurate solution of the Orr–Sommerfeld stability equation, J. Fluid Mech., 50 (1971), 689–703. https://doi.org/10.1017/S0022112071002842 doi: 10.1017/S0022112071002842
    [2] G. S. Patterson, S. A. Orszag, Spectral calculations of isotropic turbulence: efficient removal of aliasing interactions, Phys. Fluids, 14 (1971), 2538–2541. https://doi.org/10.1063/1.1693365 doi: 10.1063/1.1693365
    [3] W. Kang, N. Bedrossian, Pseudospectral optimal control theory makes debut flight, saves NASA $ 1M in under three hours, SIAM News, 40 (2007), 1–3.
    [4] Z. Liu, S. Li, K. Zhao, Extended multi-interval Legendre-Gauss-Radau pseudospectral method for mixed-integer optimal control problem in engineering, Int. J. Syst. Sci., 52 (2021), 928–951. https://doi.org/10.1080/00207721.2020.1849862 doi: 10.1080/00207721.2020.1849862
    [5] K. T. Elgindy, A high-order embedded domain method combining a Predictor-Corrector-Fourier-Continuation-Gram method with an integral Fourier pseudospectral collocation method for solving linear partial differential equations in complex domains, J. Comput. Appl. Math., 361 (2019), 372–395. https://doi.org/10.1016/j.cam.2019.03.032 doi: 10.1016/j.cam.2019.03.032
    [6] M. Nazari, M. Nazari, M. H. N. Skandari, Pseudo-spectral method for controlling the drug dosage in cancer, IET Syst. Biol., 14 (2020), 261–270. https://doi.org/10.1049/iet-syb.2020.0054 doi: 10.1049/iet-syb.2020.0054
    [7] K. T. Elgindy, B. Karasözen, High-order integral nodal discontinuous Gegenbauer-Galerkin method for solving viscous Burgers' equation, Int. J. Comput. Math., 96 (2019), 2039–2078. https://doi.org/10.1080/00207160.2018.1554860 doi: 10.1080/00207160.2018.1554860
    [8] B. Fornberg, D. M. Sloan, A review of pseudospectral methods for solving partial differential equations, Acta Numer., 3 (1994), 203–267. https://doi.org/10.1017/S0962492900002440 doi: 10.1017/S0962492900002440
    [9] B. Fornberg, A practical guide to pseudospectral methods, Cambridge University Press, 1996. https://doi.org/10.1017/CBO9780511626357
    [10] J. S. Hesthaven, S. Gottlieb, D. Gottlieb, Spectral methods for time-dependent problems, Cambridge University Press, 2007. https://doi.org/10.1017/CBO9780511618352
    [11] C. Canuto, M. Y. Hussaini, A. Quarteroni, T. A. Zang, Spectral methods: fundamentals in single domains, Springer Berlin, Heidelberg, 2006. https://doi.org/10.1007/978-3-540-30726-6
    [12] D. Garg, Advances in global pseudospectral methods for optimal control, Ph.D. Thesis, Gainesville, FL: University of Florida, 2011.
    [13] G. T. Huntington, Advancement and analysis of a Gauss pseudospectral transcription for optimal control problems, Ph.D. Thesis, Massachusetts: MIT, 2007.
    [14] J. T. Betts, Practical methods for optimal control using nonlinear programming, Philadelphia, PA: SIAM, 2001.
    [15] C. L. Darby, hp-Pseudospectral method for solving continuous-time nonlinear optimal control problems, Ph.D. Thesis, Florida: University of Florida, 2011.
    [16] O. von Stryk, R. Bulirsch, Direct and indirect methods for trajectory optimization, Ann. Oper. Res., 37 (1992), 357–373. https://doi.org/10.1007/BF02071065 doi: 10.1007/BF02071065
    [17] C. C. Francolin, Costate estimation for optimal control problems using orthogonal collocation at Gaussian quadrature points, Ph.D. Thesis, Florida: University of Florida, 2013.
    [18] K. T. Elgindy, K. A. Smith-Miles, Fast, accurate, and small-scale direct trajectory optimization using a Gegenbauer transcription method, J. Comput. Appl. Math., 251 (2013), 93–116. https://doi.org/10.1016/j.cam.2013.03.032 doi: 10.1016/j.cam.2013.03.032
    [19] K. T. Elgindy, Gegenbauer collocation integration methods: advances in computational optimal control theory, Bull. Aust. Math. Soc., 89 (2014), 168–170. https://doi.org/10.1017/S0004972713001044 doi: 10.1017/S0004972713001044
    [20] K. T. Elgindy, K. A. Smith-Miles, Optimal Gegenbauer quadrature over arbitrary integration nodes, J. Comput. Appl. Math., 242 (2013), 82–106. https://doi.org/10.1016/j.cam.2012.10.020 doi: 10.1016/j.cam.2012.10.020
    [21] K. T. Elgindy, B. Karasözen, Distributed optimal control of viscous Burgers' equation via a high-order, linearization, integral, nodal discontinuous Gegenbauer-Galerkin method, Optim. Control Appl. Methods, 41 (2020), 253–277. https://doi.org/10.1002/oca.2541 doi: 10.1002/oca.2541
    [22] C. J. Kim, S. Sung, A comparative study of transcription techniques for nonlinear optimal control problems using a pseudo-spectral method, Int. J. Aeronaut. Space Sci., 16 (2015), 264–277. https://doi.org/10.5139/IJASS.2015.16.2.264 doi: 10.5139/IJASS.2015.16.2.264
    [23] K. T. Elgindy, S. A. Dahy, High-order numerical solution of viscous Burgers' equation using a Cole-Hopf barycentric Gegenbauer integral pseudospectral method, Math. Methods Appl. Sci., 41 (2018), 6226–6251. https://doi.org/10.1002/mma.5135 doi: 10.1002/mma.5135
    [24] K. T. Elgindy, H. M. Refat, High-order shifted Gegenbauer integral pseudo-spectral method for solving differential equations of Lane–Emden type, Appl. Numer. Math., 128, (2018), 98–124. https://doi.org/10.1016/j.apnum.2018.01.018
    [25] C. W. Clenshaw, A. R. Curtis, A method for numerical integration on an automatic computer, Numer. Math., 2 (1960), 197–205. https://doi.org/10.1007/bf01386223 doi: 10.1007/bf01386223
    [26] S. E. El-Gendi, Chebyshev solution of differential, integral and integro-differential equations, Comput. J., 12 (1969), 282–287. https://doi.org/10.1093/comjnl/12.3.282 doi: 10.1093/comjnl/12.3.282
    [27] X. Gao, T. Li, Q. Shan, Y. Xiao, L. Yuan, Y. Liu, Online optimal control for dynamic positioning of vessels via time-based adaptive dynamic programming, J. Ambient Intell. Human. Comput., 2019, 1–13. https://doi.org/10.1007/s12652-019-01522-9
    [28] D. Wang, M. Ha, M. Zhao, The intelligent critic framework for advanced optimal control, Artif. Intell. Rev., 55 (2022), 1–22. https://doi.org/10.1007/s10462-021-10118-9 doi: 10.1007/s10462-021-10118-9
    [29] B. Pang, L. Cui, Z. P. Jiang, Human motor learning is robust to control-dependent noise, Biol. Cybern., 116 (2022), 307–325. https://doi.org/10.1007/s00422-022-00922-z doi: 10.1007/s00422-022-00922-z
    [30] R. F. Baum, Existence theorems for Lagrange control problems with unbounded time domain, J. Optim. Theory Appl., 19 (1976), 89–116. https://doi.org/10.1007/BF00934054 doi: 10.1007/BF00934054
    [31] G. R. Bates, Lower closure and existence theorems for optimal control problems with infinite horizon, J. Optim. Theory Appl., 24 (1978), 639–649. https://doi.org/10.1007/BF00935304 doi: 10.1007/BF00935304
    [32] A. Haurie, Existence and global asymptotic stability of optimal trajectories for a class of infinite-horizon, nonconvex systems, J. Optim. Theory Appl., 31 (1980), 515–533. https://doi.org/10.1007/BF00934475 doi: 10.1007/BF00934475
    [33] D. A. Carlson, A. Haurie, Infinite horizon optimal control: theory and applications, Springer Berlin, Heidelberg, 1987. https://doi.org/10.1007/978-3-662-02529-1
    [34] E. J. Balder, An existence result for optimal economic growth problems, J. Math. Anal. Appl., 95 (1983), 195–213. https://doi.org/10.1016/0022-247x(83)90143-9 doi: 10.1016/0022-247x(83)90143-9
    [35] D. A. Carlson, Existence of finitely optimal solutions for infinite-horizon optimal control problems, J. Optim. Theory Appl., 51 (1986), 41–62. https://doi.org/10.1007/BF00938602 doi: 10.1007/BF00938602
    [36] L. Wang, Existence and uniqueness of solutions for a class of infinite-horizon systems derived from optimal control, Int. J. Math. Math. Sci., 2005 (2005), 837–843. https://doi.org/10.1155/IJMMS.2005.837 doi: 10.1155/IJMMS.2005.837
    [37] S. Pickenhain, Infinite horizon optimal control problems in the light of convex analysis in hilbert spaces, Set-Valued Var. Anal., 23 (2015), 169–189. https://doi.org/10.1007/s11228-014-0304-5 doi: 10.1007/s11228-014-0304-5
    [38] K. O. Besov, On Balder's existence theorem for infinite-horizon optimal control problems, Math. Notes, 103 (2018), 167–174. https://doi.org/10.1134/s0001434618010182 doi: 10.1134/s0001434618010182
    [39] A. V. Dmitruk, N. V. Kuz'kina, Existence theorem in the optimal control problem on an infinite time interval, Math. Notes, 78 (2005), 466–480. https://doi.org/10.1007/s11006-005-0147-3 doi: 10.1007/s11006-005-0147-3
    [40] S. M. Aseev, An existence result for infinite-horizon optimal control problem with unbounded set of control constraints, IFAC-PapersOnLine, 51 (2018), 281–285. https://doi.org/10.1016/j.ifacol.2018.11.396 doi: 10.1016/j.ifacol.2018.11.396
    [41] V. Basco, H. Frankowska, Hamilton–Jacobi–Bellman equations with time-measurable data and infinite horizon, Nonlinear Differ. Equ. Appl., 26 (2019), 7. https://doi.org/10.1007/s00030-019-0553-y doi: 10.1007/s00030-019-0553-y
    [42] H. Halkin, Necessary conditions for optimal control problems with infinite horizons, Econometrica, 42 (1974), 267–272. https://doi.org/10.2307/1911976 doi: 10.2307/1911976
    [43] D. Garg, W. Hager, A. V. Rao, Gauss pseudospectral method for solving infinite-horizon optimal control problems, In: AIAA guidance, navigation, and control conference, Toronto, Ontario, Canada: AIAA, 2012, 1–9. https://doi.org/10.2514/6.2010-7890
    [44] D. Garg, W. W. Hager, A. V. Rao, Pseudospectral methods for solving infinite-horizon optimal control problems, Automatica, 47 (2011), 829–837. https://doi.org/10.1016/j.automatica.2011.01.085 doi: 10.1016/j.automatica.2011.01.085
    [45] D. Garg, M. A. Patterson, C. Francolin, C. L. Darby, G. T. Huntington, W. W. Hager, et al., Direct trajectory optimization and costate estimation of finite-horizon and infinite-horizon optimal control problems using a Radau pseudospectral method, Comput. Optim. Appl., 49 (2011), 335–358. https://doi.org/10.1007/s10589-009-9291-0 doi: 10.1007/s10589-009-9291-0
    [46] X. Tang, J. Chen, Direct trajectory optimization and costate estimation of infinite-horizon optimal control problems using collocation at the flipped Legendre-Gauss-Radau points, IEEE/CAA J. Autom. Sin., 3 (2016), 174–183. https://doi.org/10.1109/JAS.2016.7451105 doi: 10.1109/JAS.2016.7451105
    [47] M. Shahini, M. A. Mehrpouya, Transformed Legendre spectral method for solving infinite horizon optimal control problems, IMA J. Math. Control Inf., 35 (2018), 341–356. https://doi.org/10.1093/imamci/dnw051 doi: 10.1093/imamci/dnw051
    [48] D. Gottlieb, C. W. Shu, On the Gibbs phenomenon IV: recovering exponential accuracy in a subinterval from a Gegenbauer partial sum of a piecewise analytic function, Math. Comput., 64 (1995), 1081–1095. https://doi.org/10.2307/2153484 doi: 10.2307/2153484
    [49] D. Gottlieb, C. W. Shu, On the Gibbs phenomenon and its resolution, SIAM Rev., 39 (1997), 644–668. https://doi.org/10.1137/S0036144596301390 doi: 10.1137/S0036144596301390
    [50] J. R. Kamm, T. O. Williams, J. S. Brock, S. Li, Application of Gegenbauer polynomial expansions to mitigate Gibbs phenomenon in Fourier–Bessel series solutions of a dynamic sphere problem, Int. J. Numer. Methods Biomed. Eng., 26 (2010), 1276–1292. https://doi.org/10.1002/cnm.1207 doi: 10.1002/cnm.1207
    [51] K. T. Elgindy, Optimal control of a parabolic distributed parameter system using a fully exponentially convergent barycentric shifted Gegenbauer integral pseudospectral method, J. Ind. Manag. Optim., 14 (2018), 473–496. https://doi.org/10.3934/jimo.2017056 doi: 10.3934/jimo.2017056
    [52] W. A. Light, A comparison between Chebyshev and ultraspherical expansions, IMA J. Appl. Math., 21 (1978), 455–460. https://doi.org/10.1093/imamat/21.4.455 doi: 10.1093/imamat/21.4.455
    [53] J. P. Boyd, Orthogonal rational functions on a semi-infinite interval, J. Comput. Phys., 70 (1987), 63–88. https://doi.org/10.1016/0021-9991(87)90002-7 doi: 10.1016/0021-9991(87)90002-7
    [54] C. Canuto, M. Y. Hussaini, A. Quarteroni, T. A. Zang, Spectral methods in fluid dynamics, Springer Berlin, Heidelberg, 1988. https://doi.org/10.1007/978-3-642-84108-8
    [55] F. Fahroo, I. M. Ross, Pseudospectral methods for infinite-horizon nonlinear optimal control problems, J. Guid. Control Dynam., 31 (2008), 927–936. https://doi.org/10.2514/1.33117 doi: 10.2514/1.33117
    [56] G. Szegö, Orthogonal polynomials, American Mathematical Society, 1939.
    [57] H. Wang, D. Huybrechs, S. Vandewalle, Explicit barycentric weights for polynomial interpolation in the roots or extrema of classical orthogonal polynomials, Math. Comput., 83 (2014), 2893–2914. https://doi.org/10.1090/S0025-5718-2014-02821-4 doi: 10.1090/S0025-5718-2014-02821-4
    [58] K. T. Elgindy, High-order adaptive Gegenbauer integral spectral element method for solving non-linear optimal control problems, Optimization, 66 (2017), 811–836. https://doi.org/10.1080/02331934.2017.1298597 doi: 10.1080/02331934.2017.1298597
    [59] J. P. Berrut, Linear rational interpolation of continuous functions over an interval, In: W. Gautschi, Mathematics of computation 1943–1993: a half-century of computational mathematics, Proceedings of Symposia in Applied Mathematics, Vancouver, British Columbia: AMS, 1994,261–264. https://doi.org/10.1090/psapm/048/1314853
    [60] J. P. Berrut, H. D. Mittelmann, Lebesgue constant minimizing linear rational interpolation of continuous functions over the interval, Comput. Math. Appl., 33 (1997), 77–86. https://doi.org/10.1016/S0898-1221(97)00034-5 doi: 10.1016/S0898-1221(97)00034-5
    [61] J. M. Carnicer, Weighted interpolation for equidistant nodes Carnicer, Numer. Algor., 55 (2010), 223–232. https://doi.org/10.1007/s11075-010-9399-4 doi: 10.1007/s11075-010-9399-4
    [62] Q. Wang, P. Moin, G. Iaccarino, A rational interpolation scheme with superpolynomial rate of convergence, SIAM J. Numer. Anal., 47 (2010), 4073–4097. https://doi.org/10.1137/080741574 doi: 10.1137/080741574
    [63] L. Bos, S. De Marchi, K. Hormann, J. Sidon, Bounding the Lebesgue constant for Berrut's rational interpolant at general nodes, J. Approx. Theory, 169 (2013), 7–22. https://doi.org/10.1016/j.jat.2013.01.004 doi: 10.1016/j.jat.2013.01.004
    [64] J. P. Berrut, Rational functions for guaranteed and experimentally well-conditioned global interpolation, Comput. Math. Appl., 15 (1988), 1–16. https://doi.org/10.1016/0898-1221(88)90067-3 doi: 10.1016/0898-1221(88)90067-3
    [65] L. Bos, S. De Marchi, K. Hormann, On the Lebesgue constant of Berrut's rational interpolant at equidistant nodes, J. Comput. Appl. Math., 236 (2011), 504–510. https://doi.org/10.1016/j.cam.2011.04.004 doi: 10.1016/j.cam.2011.04.004
    [66] K. T. Elgindy, High-order, stable, and efficient pseudospectral method using barycentric Gegenbauer quadratures, Appl. Numer. Math., 113 (2017), 1–25. https://doi.org/10.1016/j.apnum.2016.10.014 doi: 10.1016/j.apnum.2016.10.014
    [67] M. R. Hestenes, Multiplier and gradient methods, J. Optim. Theory Appl., 4 (1969), 303–320. https://doi.org/10.1007/BF00927673 doi: 10.1007/BF00927673
    [68] M. J. D. Powell, A method for nonlinear constraints in minimization problems, Optimization, 1969,283–298.
    [69] K. T. Elgindy, Optimization via Chebyshev polynomials, J. Appl. Math. Comput., 56 (2018), 317–349. https://doi.org/10.1007/s12190-016-1076-x doi: 10.1007/s12190-016-1076-x
    [70] P. E. Murray, W. Murray, M. A. Saunders, SNOPT: an SQP algorithm for large-scale constrained optimization, SIAM J. Optim., 12 (2002), 979–1006. https://doi.org/10.1137/S1052623499350013 doi: 10.1137/S1052623499350013
    [71] P. E. Murray, W. Murray, M. A. Saunders, SNOPT: an SQP algorithm for large-scale constrained optimization, SIAM Rev., 47 (2005), 99–131. https://doi.org/10.1137/S0036144504446096 doi: 10.1137/S0036144504446096
    [72] D. E. Kirk, Optimal control theory: an introduction, Englewood Cliffs, N.J.: Prentice-Hall, 1970.
    [73] K. Mamehrashi, A. Nemati, A new approach for solving infinite horizon optimal control problems using Laguerre functions and Ritz spectral method, Int. J. Comput. Math., 97 (2020), 1529–1544. https://doi.org/10.1080/00207160.2019.1628949 doi: 10.1080/00207160.2019.1628949
    [74] H. S. Nik, P. Rebelo, M. S. Zahedi, Solution of infinite horizon nonlinear optimal control problems by piecewise Adomian decomposition method, Math. Model. Anal., 18 (2013), 543–560. https://doi.org/10.3846/13926292.2013.841598 doi: 10.3846/13926292.2013.841598
  • This article has been cited by:

    1. Abdul Ghaffar, Ayyaz Ali, Sarfaraz Ahmed, Saima Akram, Moin-ud-Din Junjua, Dumitru Baleanu, Kottakkaran Sooppy Nisar, A novel analytical technique to obtain the solitary solutions for nonlinear evolution equation of fractional order, 2020, 2020, 1687-1847, 10.1186/s13662-020-02751-5
    2. Haci Mehmet Baskonus, Muzaffer Ercan, Extraction Complex Properties of the Nonlinear Modified Alpha Equation, 2021, 5, 2504-3110, 6, 10.3390/fractalfract5010006
    3. Thilagarajah Mathanaranjan, Solitary wave solutions of the Camassa–Holm-Nonlinear Schrödinger Equation, 2020, 19, 22113797, 103549, 10.1016/j.rinp.2020.103549
    4. Onur Alp Ilhan, M. Nurul Islam, M. Ali Akbar, Construction of Functional Closed Form Wave Solutions to the ZKBBM Equation and the Schrödinger Equation, 2020, 2228-6187, 10.1007/s40997-020-00358-5
    5. Abdulla - Al - Mamun, Samsun Nahar Ananna, Tianqing An, Nur Hasan Mahmud Shahen, , Periodic and solitary wave solutions to a family of new 3D fractional WBBM equations using the two-variable method, 2021, 3, 26668181, 100033, 10.1016/j.padiff.2021.100033
    6. M. Al-Amin, M. Nurul Islam, M. Ali Akbar, Adequate wide-ranging closed-form wave solutions to a nonlinear biological model, 2021, 4, 26668181, 100042, 10.1016/j.padiff.2021.100042
    7. Jiahua Fang, Muhammad Nadeem, Hanan A. Wahash, Arzu Akbulut, A Semianalytical Approach for the Solution of Nonlinear Modified Camassa–Holm Equation with Fractional Order, 2022, 2022, 2314-4785, 1, 10.1155/2022/5665766
    8. Md. Tarikul Islam, Mst. Armina Akter, J. F. Gómez-Aguilar, Md. Ali Akbar, Novel and diverse soliton constructions for nonlinear space–time fractional modified Camassa–Holm equation and Schrodinger equation, 2022, 54, 0306-8919, 10.1007/s11082-022-03602-1
    9. S M Rayhanul Islam, S M Yiasir Arafat, Hanfeng Wang, Abundant closed-form wave solutions to the simplified modified Camassa-Holm equation, 2022, 24680133, 10.1016/j.joes.2022.01.012
    10. Aniqa Zulfiqar, Jamshad Ahmad, Exact solitary wave solutions of fractional modified Camassa-Holm equation using an efficient method, 2020, 59, 11100168, 3565, 10.1016/j.aej.2020.06.002
    11. Ismail Onder, Melih Cinar, Aydin Secer, Mustafa Bayram, Analytical solutions of simplified modified Camassa-Holm equation with conformable and M-truncated derivatives: A comparative study, 2022, 24680133, 10.1016/j.joes.2022.06.012
    12. Mounirah Areshi, Aly R. Seadawy, Asghar Ali, Abdulrahman F. AlJohani, Weam Alharbi, Amal F. Alharbi, Construction of Solitary Wave Solutions to the (3 + 1)-Dimensional Nonlinear Extended and Modified Quantum Zakharov–Kuznetsov Equations Arising in Quantum Plasma Physics, 2023, 15, 2073-8994, 248, 10.3390/sym15010248
    13. M. Al-Amin, M. Nurul Islam, M. Ali Akbar, The closed-form soliton solutions of the time-fraction Phi-four and (2+1)-dimensional Calogero–Bogoyavlenskii–Schiff model using the recent approach, 2022, 5, 26668181, 100374, 10.1016/j.padiff.2022.100374
    14. M. Ayesha Khatun, Mohammad Asif Arefin, M. Hafiz Uddin, Mustafa Inc, Luigi Rodino, Abundant Explicit Solutions to Fractional Order Nonlinear Evolution Equations, 2021, 2021, 1563-5147, 1, 10.1155/2021/5529443
    15. M. Nurul Islam, Onur Alp İlhan, M. Ali Akbar, Fatma Berna Benli, Danyal Soybaş, Wave propagation behavior in nonlinear media and resonant nonlinear interactions, 2022, 108, 10075704, 106242, 10.1016/j.cnsns.2021.106242
    16. Noha M. Rasheed, Mohammed O. Al-Amr, Emad A. Az-Zo’bi, Mohammad A. Tashtoush, Lanre Akinyemi, Stable Optical Solitons for the Higher-Order Non-Kerr NLSE via the Modified Simple Equation Method, 2021, 9, 2227-7390, 1986, 10.3390/math9161986
    17. M. Ayesha Khatun, Mohammad Asif Arefin, M. Ali Akbar, M. Hafiz Uddin, Numerous explicit soliton solutions to the fractional simplified Camassa-Holm equation through two reliable techniques, 2023, 20904479, 102214, 10.1016/j.asej.2023.102214
    18. M. Al-Amin, M. Nurul Islam, Onur Alp İlhan, M. Ali Akbar, Danyal Soybaş, Firdous A. Shah, Solitary Wave Solutions to the Modified Zakharov–Kuznetsov and the (2 + 1)-Dimensional Calogero–Bogoyavlenskii–Schiff Models in Mathematical Physics, 2022, 2022, 2314-4785, 1, 10.1155/2022/5224289
    19. Minzhi Wei, Liping He, Existence of periodic wave for a perturbed MEW equation, 2023, 8, 2473-6988, 11557, 10.3934/math.2023585
    20. Thitthita Iatkliang, Supaporn Kaewta, Nguyen Minh Tuan, Sekson Sirisubtawee, Novel Exact Traveling Wave Solutions for Nonlinear Wave Equations with Beta-Derivatives via the sine-Gordon Expansion Method, 2023, 22, 2224-2880, 432, 10.37394/23206.2023.22.50
    21. Melike Kaplan, Rubayyi T. Alqahtani, Nadiyah Hussain Alharthi, Wave Propagation and Stability Analysis for Ostrovsky and Symmetric Regularized Long-Wave Equations, 2023, 11, 2227-7390, 4030, 10.3390/math11194030
    22. Sujoy Devnath, Shahansha Khan, M. Ali Akbar, Exploring solitary wave solutions to the simplified modified camassa-holm equation through a couple sophisticated analytical approaches, 2024, 59, 22113797, 107580, 10.1016/j.rinp.2024.107580
    23. YANZHI MA, ZENGGUI WANG, BIFURCATION AND EXACT SOLUTIONS OF SPACE-TIME FRACTIONAL SIMPLIFIED MODIFIED CAMASSA–HOLM EQUATION, 2023, 31, 0218-348X, 10.1142/S0218348X23500858
    24. Md. Abde Mannaf, Rajandra Chadra Bhowmik, Mst. Tania Khatun, Md. Ekramul Islam, Udoy S. Basak, M. Ali Akbar, Optical solitons of SMCH model in mathematical physics: impact of wind and friction on wave, 2024, 56, 0306-8919, 10.1007/s11082-023-05641-8
    25. Mrutyunjaya Sahoo, S. Chakraverty, 2024, 9780443154041, 227, 10.1016/B978-0-44-315404-1.00019-9
    26. Ghazala Akram, Maasoomah Sadaf, Saima Arshed, Muhammad Abdaal Bin Iqbal, Simulations of exact explicit solutions of simplified modified form of Camassa–Holm equation, 2024, 56, 1572-817X, 10.1007/s11082-024-06940-4
    27. Md. Tarikul Islam, Shahariar Ryehan, Farah Aini Abdullah, J.F. Gómez-Aguilar, The effect of Brownian motion and noise strength on solutions of stochastic Bogoyavlenskii model alongside conformable fractional derivative, 2023, 287, 00304026, 171140, 10.1016/j.ijleo.2023.171140
    28. Aly R. Seadway, Asghar Ali, Ahmet Bekir, Adem C. Cevikel, Analysis of the(3+1)-Dimensional Fractional Kadomtsev–Petviashvili–Boussinesq Equation: Solitary, Bright, Singular, and Dark Solitons, 2024, 8, 2504-3110, 515, 10.3390/fractalfract8090515
    29. Faisal Yasin, Muhammad Arshad, Ghulam Farid, Mohammad Ali Hoseinzadeh, Hadi Rezazadeh, W-shape and abundant of other solitary wave solutions of the positive Gardner Kadomtsov–Petviashivilli dynamical model with applications, 2024, 56, 1572-817X, 10.1007/s11082-024-06922-6
    30. Aly R. Seadawy, Asghar Ali, Ahmet Bekir, Novel solitary waves solutions of the extended cubic(3+1)-dimensional Schr\ddot{o}dinger equation via applications of three mathematical methods, 2024, 56, 1572-817X, 10.1007/s11082-024-06528-y
    31. Jamilu Sabi’u, Mayssam Tarighi Shaayesteh, Ali Taheri, Hadi Rezazadeh, Mustafa Inc, Ali Akgül, New exact solitary wave solutions of the generalized (3 + 1)-dimensional nonlinear wave equation in liquid with gas bubbles via extended auxiliary equation method, 2023, 55, 0306-8919, 10.1007/s11082-023-04870-1
    32. M. Ashikur Rahman, M. Al-Amin, Mst. Kamrunnaher, M. Abul Kawser, Rajaul Haque, M. Ali Akbar, M. Nurul Islam, Mathematical analysis of some new adequate broad-ranging soliton solutions of nonlinear models through the recent technique, 2024, 9, 26668181, 100634, 10.1016/j.padiff.2024.100634
    33. Hassan Almusawa, Musawa Yahya Almusawa, Adil Jhangeer, Zamir Hussain, A Comprehensive Study of Dynamical Behavior and Nonlinear Structures of the Modified α Equation, 2024, 12, 2227-7390, 3809, 10.3390/math12233809
    34. Md. Asaduzzaman, Farhana Jesmin, Construction of the Closed Form Wave Solutions for TFSMCH and (1 1) Dimensional TFDMBBM Equations via the EMSE Technique+, 2025, 9, 2504-3110, 72, 10.3390/fractalfract9020072
    35. Yimin Liu, Xiaoshan Zhao, Application of homotopy analysis method to solve a class of time-fractional order mCH equations, 2025, 2964, 1742-6588, 012043, 10.1088/1742-6596/2964/1/012043
    36. Kang-Jia Wang, Guo-Dong Wang, Feng Shi, Xiao-Lian Liu, Hong-Wei Zhu, Variational principle, Hamiltonian, bifurcation analysis, chaotic behaviors and the diverse solitary wave solutions of the simplified modified Camassa–Holm equation, 2025, 22, 0219-8878, 10.1142/S0219887825500136
    37. H. A. Ashi, Noufe H. Aljahdaly, Investigation of the damped geophysical KdV equation using the explicit exponential time differencing method, 2025, 81, 1573-0484, 10.1007/s11227-025-07289-5
  • Reader Comments
  • © 2023 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(2656) PDF downloads(117) Cited by(5)

Figures and Tables

Figures(20)  /  Tables(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog