Loading [MathJax]/jax/output/SVG/jax.js
Research article

Epilepsy EEG classification method based on supervised locality preserving canonical correlation analysis


  • Received: 09 July 2021 Accepted: 15 November 2021 Published: 18 November 2021
  • Existing epileptic seizure automatic detection systems are often troubled by high-dimensional electroencephalogram (EEG) features. High-dimensional features will not only bring redundant information and noise, but also reduce the response speed of the system. In order to solve this problem, supervised locality preserving canonical correlation analysis (SLPCCA), which can effectively use both sample category information and nonlinear relationships between features, is introduced. And an epileptic signal classification method based on SLPCCA is proposed. Firstly, the power spectral density and the fluctuation index of the frequency slice wavelet transform are extracted as features from the EEG fragments. Next, SLPCCA obtains the optimal projection direction by maximizing the weight correlation between the paired samples in the class and their neighbors. And the projection combination of original features in the optimal direction is the fusion feature. The fusion features are then input into LS-SVM for training and testing. This method is verified on the Bonn dataset and the CHB-MIT dataset and gets good results. On various classification tasks of Bonn data set, the proposed method achieves an average classification accuracy of 99.16%. On the binary classification task of the inter-seizure and seizure epileptic EEG of the CHB-MIT dataset, the proposed method achieves an average accuracy of 97.18%. The experimental results show that the algorithm achieves excellent results compared with several state-of-the-art methods. In addition, the parameter sensitivity of SLPCCA and the relationship between the dimension of the fusion features and the classification results are discussed. Therefore, the stability and effectiveness of the method are further verified.

    Citation: Hongming Liu, Yunyuan Gao, Jianhai Zhang, Juanjuan Zhang. Epilepsy EEG classification method based on supervised locality preserving canonical correlation analysis[J]. Mathematical Biosciences and Engineering, 2022, 19(1): 624-642. doi: 10.3934/mbe.2022028

    Related Papers:

    [1] Mengjun Yu, Kun Li . A data-driven reduced-order modeling approach for parameterized time-domain Maxwell's equations. Networks and Heterogeneous Media, 2024, 19(3): 1309-1335. doi: 10.3934/nhm.2024056
    [2] Andreas Hiltebrand, Siddhartha Mishra . Entropy stability and well-balancedness of space-time DG for the shallow water equations with bottom topography. Networks and Heterogeneous Media, 2016, 11(1): 145-162. doi: 10.3934/nhm.2016.11.145
    [3] Leqiang Zou, Yanzi Zhang . Efficient numerical schemes for variable-order mobile-immobile advection-dispersion equation. Networks and Heterogeneous Media, 2025, 20(2): 387-405. doi: 10.3934/nhm.2025018
    [4] Patrick Henning, Mario Ohlberger . The heterogeneous multiscale finite element method for advection-diffusion problems with rapidly oscillating coefficients and large expected drift. Networks and Heterogeneous Media, 2010, 5(4): 711-744. doi: 10.3934/nhm.2010.5.711
    [5] JinJun Yong, Changlun Ye, Xianbing Luo . A fully discrete HDG ensemble Monte Carlo algorithm for a heat equation under uncertainty. Networks and Heterogeneous Media, 2025, 20(1): 65-88. doi: 10.3934/nhm.2025005
    [6] Liuchao Xiao, Wenbo Li, Leilei Wei, Xindong Zhang . A fully discrete local discontinuous Galerkin method for variable-order fourth-order equation with Caputo-Fabrizio derivative based on generalized numerical fluxes. Networks and Heterogeneous Media, 2023, 18(2): 532-546. doi: 10.3934/nhm.2023022
    [7] Salim Meddahi, Ricardo Ruiz-Baier . A new DG method for a pure–stress formulation of the Brinkman problem with strong symmetry. Networks and Heterogeneous Media, 2022, 17(6): 893-916. doi: 10.3934/nhm.2022031
    [8] Eric Chung, Yalchin Efendiev, Ke Shi, Shuai Ye . A multiscale model reduction method for nonlinear monotone elliptic equations in heterogeneous media. Networks and Heterogeneous Media, 2017, 12(4): 619-642. doi: 10.3934/nhm.2017025
    [9] Timothy Blass, Rafael de la Llave . Perturbation and numerical methods for computing the minimal average energy. Networks and Heterogeneous Media, 2011, 6(2): 241-255. doi: 10.3934/nhm.2011.6.241
    [10] Yue Tai, Xiuli Wang, Weishi Yin, Pinchao Meng . Weak Galerkin method for the Navier-Stokes equation with nonlinear damping term. Networks and Heterogeneous Media, 2024, 19(2): 475-499. doi: 10.3934/nhm.2024021
  • Existing epileptic seizure automatic detection systems are often troubled by high-dimensional electroencephalogram (EEG) features. High-dimensional features will not only bring redundant information and noise, but also reduce the response speed of the system. In order to solve this problem, supervised locality preserving canonical correlation analysis (SLPCCA), which can effectively use both sample category information and nonlinear relationships between features, is introduced. And an epileptic signal classification method based on SLPCCA is proposed. Firstly, the power spectral density and the fluctuation index of the frequency slice wavelet transform are extracted as features from the EEG fragments. Next, SLPCCA obtains the optimal projection direction by maximizing the weight correlation between the paired samples in the class and their neighbors. And the projection combination of original features in the optimal direction is the fusion feature. The fusion features are then input into LS-SVM for training and testing. This method is verified on the Bonn dataset and the CHB-MIT dataset and gets good results. On various classification tasks of Bonn data set, the proposed method achieves an average classification accuracy of 99.16%. On the binary classification task of the inter-seizure and seizure epileptic EEG of the CHB-MIT dataset, the proposed method achieves an average accuracy of 97.18%. The experimental results show that the algorithm achieves excellent results compared with several state-of-the-art methods. In addition, the parameter sensitivity of SLPCCA and the relationship between the dimension of the fusion features and the classification results are discussed. Therefore, the stability and effectiveness of the method are further verified.



    In this article, we prove the non-existence of solutions to the following quasilinear elliptic problem which has degenerate coercivity in their principal part by approximation,

    {div(a(x,u,u))+|u|q1u=λ,xΩ,u=0,xΩ, (1)

    where 1<p<N,q>1 and λ is a Radon measure. Ω is a bounded smooth subset of RN(N>2). a(x,t,ξ):Ω×R×RNRN is the Carathéodory function (i.e: a(x,t,ξ) is measure on Ω for every (t,ξ) in R×RN, and a(,t,ξ) is continuous on R×RN for almost every x in Ω), such that the following assumptions hold,

    a(x,t,ξ)ξc|ξ|p(1+|t|)θ(p1), (2)
    |a(x,t,ξ)|c0(|ξ|p1+b(x)), (3)
    [a(x,t,ξ)a(x,t,ξ)][ξξ]>0, (4)

    for almost every xΩ,tR,ξ,ξRN with ξξ, where 0θ<1, c and c0 are two positive constants, bLp(Ω) is a non-negative function, p is the conjugate Hölder exponent of p.

    It is well-known that[3,9], problem Δu+|u|q1u=δ0 has no distributional solution if qNN2. On the other hand, if q<NN2, then there exists a unique solution to

    {Δu+|u|q1u=δ0,xΩ,u=0,xΩ.

    In the famous work [9], Brezis proved that if {un} is sequence of solution to the nonlinear elliptic problem

    {Δun+|un|q1un=fn,xΩ,un=0,xΩ, (5)

    with q>NN2, and fnL(Ω) is a sequence functions such that, for any ϱ>0,

    limnΩBϱ(0)|fnf|=0.

    Then un converges to the unique solution u to the following equation

    {Δu+|u|q1u=f,xΩ,u=0,xΩ.

    This fact shows that Bϱ(0) is a removable singularity set of solution to equation (5) provided q>NN2. Orsina and Prignet[24] extended the result of [9] to more general operator div(a(x,u,u)), where a(x,u,u) satisfies (2)-(4) with θ=0. The main results of [24] shown that problem (1) with θ=0 has a solution for every given bounded measure λ if q<r(p1)rp. Some other related results see [12,6,10,8,14,23,26,27,21,19,16] and references therein.

    The main goal of this paper is to study the non-existence of solutions to problem (1). More precisely, consider the limit of approximating equation (9)(see Theorem 1.2 below), our main task is to understand which is the limit of solutions to (9) and what equation it satisfies. A point worth emphasizing is that, even if p=2, the convergence of solutions is not true if the right hand side are distributions weakly converging in W1,2(Ω), see [5] for some counterexamples.

    In order to state the main results of this paper, we need some definitions.

    Let K be a compact subset of Ω, r>1 is a real number. The r capacity of K respect to Ω is defined as

    capr(K,Ω)=inf{urW1,r0:uCc(Ω),uχK},

    where χK is the characteristic function of K.

    Let λ be a bounded measure on Ω, we say that λ is concentrated on a set E if λ(B)=λ(BE) for every Borel subset B of Ω. Thanks to the Hahn decomposition, λ can be decomposed as the difference of two nonnegative mutually singular measure, that is λ=λ+λ.

    If λ is concentrated on a set E, as a consequence of the fact that λ+ and λ are mutually singular, we have that λ+ and λ concentrated a set E+ and E respectively and E+E=.

    Let λ=λ+λ be a measure, fn=f+nfn approximations of λ in the following way:

    limn+Ωf+nφdx=Ωφdλ+,limn+Ωfnφdx=Ωφdλ, (6)

    for every function φ, which is continuous and bounded on Ω, where {f+n} and {fn} are sequences of nonnegative L(Ω) functions. We not assume that f+n and fn are the positive and negative part of fn. Observe that choosing φ1 in (6), we obtain

    f+nL1(Ω)C,fnL1(Ω)C. (7)

    For all k>0,sR, define

    Tk(s)=max{k,min{k,s}},Gk(s)=sTk(s).

    Firstly we stale the existence result.

    Theorem 1.1. Let Ω be a bounded smooth subset of RN(N>2), 1<p<N, gL1(Ω) and (2)-(4) hold. Then there exists a unique entropy solution uW1,p0(Ω) to problem

    {div(a(x,u,u))+|u|q1u=g,xΩ,u=0,xΩ. (8)

    if

    q<N(1θ)N(1+θ(p1)).

    Moreover,

    uMp1(Ω),|u|Mp2(Ω),

    where Mp1,Mp2 represents the Marcinkiewicz space with exponent

    p1=N(p1)(1θ)Np,p2=N(p1)(1θ)N(1+θ(p1)).

    Remark 1. The previous result gives existence and uniqueness of the entropy solution uW1,p0(Ω) to (8) for every 1<p<N and 0<θ<1. If θ=0, the same result for (8) can be proved by the same techniques of [2].

    Our main results are following:

    Theorem 1.2. Let 1<p<rN and λ=λ+λ be a bounded Radon measure which is concentrated on a set E with zero r capacity. Let fn=f+nfn be a sequence of L(Ω) functions which converge to λ in the sense of (6). gL1(Ω) and let gn is a sequence of L(Ω) functions which converge to g weakly in L1(Ω). Suppose unW1,p0(Ω) is the solution to problem:

    {div(a(x,un,un))+|un|q1un=fn+gn,xΩ,un=0,xΩ. (9)

    Then |un|p1 strong converges to |u|p1 in Lσ(Ω) as n for

    σ<pq(q+1+θ(p1))(p1),

    if

    q>r(p1)[1+θ(p1)]rp, (10)

    where u is unique solution of (8). Moreover,

    limn+Ω|un|q1unφdx=Ω|u|q1uφdx+Ωφdλ,φC(Ω). (11)

    Remark 2. The above theorem shows that there is not a solution to problem (1) can be obtained by approximation, if q is large enough and the measure λ is concentrated on a set with zero r capacity.

    Remark 3. Boccardo et.al [7] considered the non-existence result to the following problem

    {div(a(x,u)(1+u)γ)+u=μ,xΩ,u=0,xΩ, (12)

    where γ>1 and μ is a non-negative Radon measure, concentrated on a set E with zero harmonic capacity, a(x,ξ) satisfies (2)-(4) with θ=0, p=2 and b(x)=0. While in Theorem 1.2, λ is a bounded Radon measure concentrated on a set E with zero r capacity with p<rN, instead of p capacity. Therefore Theorem 1.2 is not a triviality extend the results of Theorem 4.1 of [7]. Furthermore, in Theorem 1.2, θ(p1)(0,p1) since θ(0,1). Note that, in problem (12), they required that γ>1. It is worth pointing out that different ranges of γ have an important impact on the behavior of solutions to problem (12), more details see [25,18,1,17,13,4].

    The structure of this paper is as follows: Section 2 mainly gives some lemmas which play a important role in the process of proof of the main theorem. The proof of theorem 1.1 and 1.2 are given in Section 3.

    In the following, C is a constant and its value may changes from line to line.

    In order to prove Theorem 1.1 and 1.2, the following basic lemmas and definitions are required.

    Lemma 2.1. (see Lemma 2.1 of [22]) Let K+ and K be two disjoin compact subsets of Ω with zero r capacity, λ=λ+λ be a measure which is concentrated on a set with zero r capacity with 1<rN, Then there exist two functions ψ+δ and ψδ in Cc(Ω), such that

    0ψ+δ1,0ψδ1,Ω|ψ+δ|rdxδ,Ω|ψδ|rdxδ,0Ω(1ψ+δ)dλ+δ,0Ω(1ψδ)dλδ,0Ωψδdλ+δ,0Ωψ+δ)dλδ,ψ+δ1,xK+,ψ+δ1,xK, (13)

    for every δ>0.

    Definition 2.2. Let u be an measurable function on Ω such that Tk(u)W1,p0(Ω) for every k>0. Then there exist a unique measurable function v:ΩRN such that

    Tk(u)=vχ{|u|k},a.einΩandforeveryk>0.

    Define the gradient of u as the function v and denote it by v=u.

    Definition 2.3. Let fL1(Ω), q>0 and (2)-(4) hold. A measurable function u is an entropy solution to problem (8), if Tk(u)W1,p0(Ω) for every k>0, |u|qL1(Ω) and

    Ωa(x,u,u)Tk(uφ)dx+Ω|u|q1uTk(uφ)dxΩgTk(uφ)dx,

    for every φW1,p0(Ω)L1(Ω).

    Definition 2.4. Marcinkiewicz space Ms(Ω)(s>0) is the space composed of all the measurable functions v that satisfy

    |{|υ|k}|Cks,

    for any k>0, where the constant C>0.

    If |Ω| is bounded and 0<ε<s1, then the following embedding relationship hold:

    Ls(Ω)Ms(Ω)Lsε(Ω).

    Lemma 2.5. Let uMs(Ω) with s>0. If there exist a constant ρ>0, such that for any k>0,

    Ω|Tk(u)|pdxCkρ,

    for some positive constant C. Then

    |u|Mpss+ρ(Ω).

    Proof. Let σ be a fixed positive real number, for every k>0,

    |{|u|>σ}|=|{|u|>σ,|u|k}|+|{|u|>σ,|u|>k}||{|Tk(u)|>σ}|+|{|u|>k}|. (14)

    Moreover,

    |{|Tk(u)|>σ}|1σpΩ|Tk(u)|pdxCkρσp. (15)

    Since uMs(Ω), by Definition 2.4, there exist a constant C such that

    |{|u|>k}|Cks. (16)

    Combining (14)-(16), we have

    |{|u|>σ}|Ckρσp+CksCkpss+ρ.

    Therefore, by Definition 2.4, we get |u|Mpss+ρ.

    Lemma 2.6. Let {un} be a sequence in W1,p0(Ω) and assume that there exist positive constants ρ and C with p>ρ, such that

    Ω|Tk(un)|pdxCkρ,

    for any k and n. Then there exists a subsequence, still denoted by {un}, which converges to a measurable function v almost everywhere in Ω.

    Lemma 2.7. Let u be an entropy solution to (8), then

    {k<|u|<k+h}|u|pdxCkθ(p1).

    Proof. For any given h and k>0,sR, define

    Tk,h(s)=Th(sTk(s))={sksgn(s),k|s|<k+h,h,|s|k+h,0,|s|k.

    Take Tk,h(u) as test function in (8), we have

    {k<|u|<k+h}(a(x,u,u)u)dx+Ω|u|q1uTk,h(u)dx=ΩgTk,h(u)dx. (17)

    Since uTk,h(u)0, we find

    {k<|u|<k+h}(a(x,u,u)u)dxΩgTk,h(u)dx, (18)

    and

    ΩgTk,h(u)dxh{|u|>k}|g|dxC. (19)

    According to the assumption (2) and (17)-(19), we get,

    {k<|u|<k+h}|u|pdxCkθ(p1).

    Proposition 1. Let uW1,p0(Ω) be an entropy solution to (8) and satisfy

    {|u|<k}|u|pdxCkρ (20)

    for every k>0 and p>ρ. Then uMp1(Ω), where p1=N(pρ)/(Np). More precisely, there exists C=C(N,p,θ)>0 such that

    |{|u|>k}|Ckp1.

    Proof. For every k>0, by the Sobolev embedding theorem and (20),

    Tk(u)pC(N,p,θ)Tk(u)pCkρp,

    where p=NpNp. For 0<ηk, we have

    {|u|η}={|Tk(u)η|}.

    Hence

    |{|u|>η}|Tk(u)ppηpC(kρ)ppηp.

    Setting η=k, we obtain

    |{|u|>k}|CkN(pρ)Np.

    This fact shows that uMp1(Ω) with p1=N(pρ)/(Np).

    Proposition 2. Assume that uW1,p0(Ω) is an entropy solution to (8), which satisfies (20) for every k. Then uMp2(Ω), where p2=N(pρ)/(Nρ), that is there exists C=C(N,p,θ)>0 such that

    |{|u|>h}|Chp2,

    for every h>0.

    Proof. For k,λ>0, set

    ψ(k,λ)=|{|u|p>λ,|u|>k}|.

    Using the fact that the function λψ(k,λ) is nonincreasing, we get, for k,λ>0,

    ψ(0,λ)=|{|u|p>λ}|1λλ0ψ(0,s)dsψ(k,0)+1λλ0ψ(0,s)ψ(k,s)ds. (21)

    By Proposition 1,

    ψ(k,0)Ckp1, (22)

    where p1=N(pρ)/(Np). Since ψ(0,s)ψ(k,s)=|{|u|p>s,|u|<k}|, thanks to (20), we have

    0ψ(0,s)ψ(k,s)ds={|u|<k}|u|pdxCkρ. (23)

    Combining (21)-(23), we arrive at

    ψ(0,λ)Ckρλ+Ckp1. (24)

    Let Ckρλ=Ckp1 and λ=hp, (24) implies that

    |{|u|>h}|ChN(pρ)Nρ.

    That is uMp2(Ω) with p2=N(pρ)/(Nρ).

    In this section we prove Theorem 1.1 and 1.2 combining the results of Sections 2.

    In the proofs of Theorem 1.1 and 1.2, ω(n,m,δ) will denote any quantity (depending on n,m and δ) such that

    limδ0+limm+limn+ω(n,m,δ)=0.

    If the quantity does not depend on one or more of the three parameters n,m and δ, we will omit the dependence from it in ω. For example, ω(n,δ) is any quantity such that

    limδ0+limn+ω(n,δ)=0.

    The proof of Theorem 1.1 will be divided in several steps.

    Proof. (1)Uniqueness: Let u1 and u2 be two entropy solutions to equation (8). The proof of the fact that u1=u2 will follow from the following four steps.

    Step 1. Assume that giL1(Ω),(i=1,2). Choosing Tk(u1Thu2) and Tk(u2Thu1) as test function in (8) respectively, we get

    I:=Ωa(x,u1,u1)Tk(u1Thu2)dx+Ωa(x,u2,u2)Tk(u2Thu1)dx=Ω|u1|q1u1Tk(u1Thu2)dxΩ|u2|q1u2Tk(u2Thu1)dx+Ωg1Tk(u1Thu2)dx+Ωg2Tk(u2Thu1)dx. (25)

    Step 2. Denote

    A0={xΩ:|u1u2|<k,|u1|<h,|u2|<h},A1={xΩ:|u1Thu2|<k,|u2|h},A2={xΩ:|u1Thu2|<k,|u2|<h,|u1|h}.

    For xA0,

    Tk(u1Thu2)=(u1u2)

    and

    Tk(u2Thu1)=Tk(u2u1).

    Thus, for every xA0,

    Ωa(x,u1,u1)Tk(u1Thu2)dx+Ωa(x,u2,u2)Tk(u2Thu1)dx=A0[a(x,u1,u1)a(x,u2,u2)](u1u2)dx:=I0. (26)

    For xA1, Tk(u1Thu2)=(u1h)=u1. By (2), we get

    Ωa(x,u1,u1)Tk(u1Thu2)dx=A1a(x,u1,u1)u1dx0. (27)

    For xA2, Tk(u1Thu2)=(u1u2). Thus

    Ωa(x,u1,u1)Tk(u1Thu2)dxA2a(x,u1,u1)u2dx. (28)

    Similarly, denote

    A1={xΩ:|u2Thu1|<k,|u1|h},A2={xΩ:|u2Thu1|<k,|u1|<h,|u2|h}.

    Then for xA1, Tk(u2Thu1)=(u2h)=u2. By (2), we get

    Ωa(x,u2,u2)Tk(u2Thu1)dx=A1a(x,u2,u2)u2dx0. (29)

    For xA2, Tk(u2Thu1)=(u2u1). Thus

    Ωa(x,u2,u2)Tk(u2Thu1)dxA2a(x,u2,u2)u1dx. (30)

    Summing up (26)-(30) in the form II0I1, where

    I1=A2a(x,u1,u1)u2dx+A2a(x,u2,u2)u1dx:=I11+I12.

    Now, we estimate I11. By the Hölder inequality and (3), we have

    I11a(x,u1,u1)Lp({h|u1|h+k})u2Lp({hk|u2|h})c0(u1p1Lp({h|u1|h+k})+b(x)Lp({|u1|h}))u2Lp({hk|u2|h}).

    Therefore, by Lemma 2.7 and Proposition 2, I110 as h for every k>0. I120 as h for every k>0 can be obtained in the same way.

    Hence, we find

    Ωa(x,u1,u1)Tk(u1Thu2)dx+Ωa(x,u2,u2)Tk(u2Thu1)dx=A0[a(x,u1,u1)a(x,u2,u2)](u1u2)dx+ε(h). (31)

    Step 3. Now estimate the terms on the right hand side of (25). Denote

    B0={xΩ:|u1|<h,|u2|<h},B1={xΩ:|u1|h},B2={xΩ:|u2|h}.

    For xB0, since Tk(u1Thu2)=Tk(u1u2) and Tk(u2Thu1)=Tk(u2u1), we arrive at

    Ω|u1|q1u1Tk(u1Thu2)dx+Ω|u2|q1u2Tk(u2Thu1)dx=B0(|u1|q1u1|u2|q1u2)Tk(u1u2)dx0, (32)

    and

    Ωg1Tk(u1Thu2)dx+Ωg2Tk(u2Thu1)dx=B0(g1g2)Tk(u1u2)dx0. (33)

    For xB1, since Tk(u2Thu1)=Tk(u2h). Then

    Ω|u1|q1u1Tk(u1Thu2)dx+Ω|u2|q1u2Tk(u2Thu1)dxkB1(|u1|q1u1+|u2|q1u2)dx:=J1,

    and

    Ωg1Tk(u1Thu2)dx+Ωg2Tk(u2Thu1)dxkB1(|g1|+|g2|)dx:=J2.

    For xB2, since Tk(u1Thu2)=Tk(u1h), we get

    Ω|u1|q1u1Tk(u1Thu2)dx+Ω|u2|q1u2Tk(u2Thu1)dxkB2(|u1|q1u1+|u2|q1u2)dx:=J1,

    and

    Ωg1Tk(u1Thu2)dx+Ωg2Tk(u2Thu1)dxkB2(|g1|+|g2|)dx:=J2.

    According to |B1|0,|B2|0 as h and |u|qL1(Ω) for fixed k>0, we get

    J1+J2+J1+J20ash. (34)

    Step 4. Combining (25) and (31)-(34), we have

    A0[a(x,u1,u1)a(x,u2,u2)](u1u2)dxε(h),

    where ε(h)0 as h. Since A0 converges to {xΩ:|u1u2|<k} by measure as h for fixed k>0, we conclude that

    {|u1u2|<k}[a(x,u1,u1)a(x,u2,u2)](u1u2)dx0,

    for all k>0. This fact, combine with (4), implies that u1=u2 a.e in Ω. Then we get u1=u2 a.e in Ω.

    (2) Existence:

    Step 1. Let

    F(x,u)=g(x)β(u),

    where β(u)=|u|q1u, which is continuous with respect to u. Then g(x)=F(x,0)L1(RN) and β is monotonous nondecreasing with respect to u with β(0)=0 and β(u)u0.

    Let gnC0, such that gn converges to g in L1(Ω), with gnL1(Ω)gL1(Ω) for every n1. Define βn(s)=Tn(β). In this way, |βn(s)||β(s)| for every sR and xΩ. Finally we take

    γn(s)=βn(s)+1n|s|p2s.

    Then by [20], there exists unW1,p0(Ω) such that

    {diva(x,un,un)+γn(x,un)=gn,xΩ,un=0,xΩ, (35)

    holds in the sense of distributions in Ω.

    By density arguments, we can take Th(unTk(un)) and Tk(un) as the test function in (35) respectively, we have

    {k|un|<k+h}a(x,un,un)undx+{|un|>k}γnTh(unTk(un))dx={|un|>k}gnTh(unTk(un))dx, (36)

    and

    {|un|>k}a(x,un,un)undx+ΩγnTk(un)dx=ΩgnTk(un)dx. (37)

    Combine (36) with (2) (fix the ellipticity constant c=1) and γnTh(unTk(un))0, we get,

    {k<|un|<k+h}|un|pdxhkθ(p1){|un|>k}gndxhkθ(p1)gnL1(Ω)=Ckθ(p1). (38)

    Since a(x,un,un)un0 by (2), we have

    {|un|>k}|γn(un)|dx{|un|>k}|gn|dxgnL1(Ω)C. (39)

    Combine (37) with γnTk(un)0, we have

    {|un|<k}|un|pdxCk1+θ(p1). (40)

    Step 2. Convergence. Using (38) and Proposition 1, we have |{|un|>k}| is bounded uniformly for every k>0. Thanks to (40), we see that {Tk(un)} is bounded in Lploc(Ω) for every k>0.

    Next we prove that unu locally in measure.

    For t,ϵ>0, we have

    {|unum|>t}{|un|>k}{|um|>k}{|Tk(un)Tk(um)|>t}.

    Thus

    |{|unum|>t}||{|un|>k}|+|{|um|>k}|+|{|Tk(un)Tk(um)|>t}|.

    Choosing k large enough such that |{|un|>k}|<ϵ and |{|um|>k}|<ϵ. Since {Tk(un)}n is bounded in Lp(Ω) and Tk(un)W1,p0(Ω) for every k>0. Assume that {Tk(un)} is a Cauchy sequence in Lq(ΩBR) for any q<pN/(Np) and any R>0,

    Tk(un)Tk(u)inLploc(Ω)anda.einΩ.

    Then

    |{|Tk(un)Tk(um)|>t}BR|tqΩBR|Tk(un)Tk(um)|qdxϵ,

    for all n,mn0(k,t,R). This show that {un} is a Cauchy sequence in BR. Hence that unu locally.

    Now to prove that un converges to some function v locally. We need to prove that {un} is a Cauchy sequence in any ball BR. Let t,ϵ>0 again, then

    {|unum|>t}BR{|unum|k,|un|l,|um|l,|unum|>t}{|un|>l}{|um|>l}({|unum|>k}BR).

    Choose l large enough such that |{|un|>l}|ϵ for all nN. If a is a continuous function independent of x, then by (4), there exists a μ>0, such that |ξ|<l,|ξ|<l and |ξξ|>t means

    [a(x,t,ξ)a(x,t,ξ)][ξξ]μ.

    This is a consequence of continuity and strict monotonicity of a. Set

    dn=gnγn(x,un). (41)

    Taking Tk(unum) as the test function of (35) and by (37), (41), we have

    {|unum|<k}[a(x,un,un)a(x,um,um)](unum)dx=Ω(dndm)Tk(unum)dxCk1+θ(p1).

    Then

    {|unum|k,|un|l,|um|l,|unum|>t}1μ{|unum|<k}[a(x,un,un)a(x,um,um)](unum)dx1μCk1+θ(p1)ϵ,

    if k is small enough such that k1+θ(p1)μϵ/C.

    Since l and k have been confirmed, if n0 large enough, we have |({|unum|>k}BR)|ϵ for n,mn0. Then we get |{|unum|>t}BR|4ϵ. This prove that un converges to some function v locally.

    Finally, since {Tk(un)}nLp(Ω) for every k>0, it converges weakly to {Tk(u)} in Lploc(Ω). We have uW1,p0(Ω) and u=v a.e in Ω.

    Step 3. In order to prove the existence of the solution completely, we still need to prove that sequence {a(x,u,u)}n is bounded in Lqloc(Ω) for all

    q(1,N(1θ)N(1+θ(p1))).

    Indeed, by Proposition 2, |un|p1MN(1θ)N(1+θ(p1))Lqloc(Ω). And by (3), we have |a(x,un,un)|Lp(Ω)Lqloc(Ω). According to the Nemitskii's theorem, unu implies that

    a(x,un,un)a(x,u,u).

    It follows that

    a(x,u,u)MN(1θ)N(1+θ(p1))Lqloc(Ω),

    for all q(1,N(1θ)N(1+θ(p1))).

    In this subsection, we give the proof of Theorem 1.2 following some ideas in [11,22].

    Proof. Step 1 (A priori estimates). Firstly, choosing Tk(un)(1φδ)s as test function in the weak formulation of (9), where s=ηηp+1 and η will be given in (48), we have

    Ωa(x,un,un)Tk(un)(1φδ)sdx+Ω|un|q1unTk(un)(1φδ)sdx=sΩa(x,un,un)φδTk(un)(1φδ)s1dx+ΩgnTk(un)(1φδ)sdx+Ωf+nTk(un)(1φδ)sdx+ΩfnTk(un)(1φδ)sdx. (42)

    By (2), we get

    Ωa(x,un,un)Tk(un)dμcΩ|Tk(un)|p(1+|Tk(un)|)θ(p1)dμ, (43)

    here dμ:=(1φδ)sdx.

    Since unTk(un)0,

    Ω|un|q1unTk(un)(1φδ)sdx{|un|k}|un|q1unTk(un)dμkq+1μ({|un|k}). (44)

    Using (3) and the Young inequality, we find

    Ω|a(x,un,un)φδTk(un)(1φδ)s1|dxc0kΩ(|un|p1+b(x))(|φ+δ|+|φ+δ|)(1φδ)s1dxCkΩ(|un|(p1)r+|b(x)|r)(1φδ)(s1)rdx+CkΩ(|φ+δ|r+|φ+δ|r)dxCk(Ω(|un|(p1)r+|b(x)|r)(1φδ)(s1)rdx+δ). (45)

    Combine (42)-(45), by (7) and {gn}L1(Ω), bLp(Ω), we have

    Ω|Tk(un)|p(1+|Tk(un)|)θ(p1)dμ+kq+1μ({|un|k})Ck(Ω|un|(p1)r(1φδ)(s1)rdx+δ+μ(Ω). (46)

    For a fixed σ0, thanks to (46), we get

    μ({|un|>σ})=μ({|un|>σ,|un|<k})+μ({|un|>σ,|un|k})1σpΩ|Tk(un)|pdμ+μ({|u|>k})(1+k)θ(p1)σpΩ|Tk(un)|p(1+|Tk(un)|)θ(p1)dμ+μ({|u|>k})C(Ω|un|(p1)r(1φδ)(s1)rdx+δ+μ(Ω))((1+k)1+θ(p1)σp+1kq),

    which implies

    μ|{|un|>σ}|Cσpqq+1+θ(p1)(Ω|un|(p1)r(1φδ)(s1)rdx+δ+μ|Ω|). (47)

    Let

    (p1)r<η<pqq+1+θ(p1). (48)

    Clearly, such η exists by (10). In view of (47)-(48), we have

    Ω|un|ηdμC(Ω|un|(p1)r(1φδ)(s1)rdx+δ+μ(Ω)).

    By the Holder's inequality,

    Ω|un|(p1)r(1φδ)(s1)rdxC(Ω|un|ηdμ)(p1)rηC(Ω|un|(p1)r(1φδ)(s1)rdx+δ+μ|Ω|)(p1)rη.

    By Lemma 2.1, 1φδ is zero both on a neighbourhood of K+ and K. Hence

    Ω|un|(p1)r(1φδ)(s1)rdxC(δ+μ|Ω|)C(δ). (49)

    Using (46) and (49), we conclude that

    Ω|Tk(un)|pdxCk1+θ(p1). (50)

    According to Lemma 2.5, we have |un|Ms(Ω), where s=pqq+1+θ(p1).

    By (50) and Lemma 2.6, there exists a subsequence, still denoted by un, which converges to a measurable function u almost everywhere in Ω. So Tk(un)Tk(u) in Ω for every k>0.

    Since Tk(un)W1,p0(Ω), by the weak lower semi-continuity of the norm, Tk(u)W1,p0(Ω) for every k>0. Thus u has an gradient u in the sense of Definition 2.2, as a consequence of the a priori estimates on un and (4), we have

    a(x,un,un)a(x,u,u)stronglyin(Ls(Ω))N, (51)

    for every s<pq(q+1+θ(p1))(p1).

    Step 2 (Energy estimates). Let ψδ=ψ+δ+ψδ, where ψ+δ and ψδ are as in Lemma 2.1. Then

    {un>2m}uqn(1ψδ)dx=ω(n,m,δ), (52)

    and

    {un<2m}|un|q(1ψδ)dx=ω(n,m,δ). (53)

    Choose βm(un)(1ψδ) as test function in the weak formulation of (9), where

    βm(s)={sm1,m<s2m,1,s>2m,0,sm.

    We obtain

    1m{m<un<2m}a(x,un,un)un(1ψδ)dx(A)Ωa(x,un,un)ψδβm(un)dx(B)+Ω|un|q1unβm(un)(1ψδ)dx(C)=Ωf+nβm(un)(1ψδ)dx(D)Ωfnβm(un)(1ψδ)dx(E)+Ωgnβm(un)(1ψδ)dx.(F)

    Since (A) and (E) are non-negative, we can get rid of them. And since βm(um) converges to βm(u) almost everywhere in Ω and in the weaktopology of L(Ω), βm(un) converges to zero in the weaktopology of L(Ω) as m, we have

    (B)=Ωa(x,u,u)ψδβm(u)dx+ω(n)=ω(n,m),

    and

    (C){un>2m}uqn(1ψδ)dx.

    By ψδ=ψ+δ+ψδ and (6),

    (D)Ωf+n(1ψδ)dx=Ω(1ψ+δ)dλ+Ωψδdλ+ω(n)=ω(n,δ),

    and

    (F)=ω(n,m).

    We get (52), the proof of (53) is identical.

    Step 3 (Passing to the limit). Now we show that u is an entropy solution to (8) with datum g. Choose Tk(unφ)(1ψδ) as test function in the weak formulation of (9), we get

    Ωa(x,un,un)Tk(unφ)(1ψδ)dx(A)Ωa(x,un,un)ψδTk(unφ)dx(B)+Ω|un|q1unTk(unφ)(1ψδ)dx(C)=Ωf+nTk(unφ)(1ψδ)dx(D)
    ΩfnTk(unφ)(1ψδ)dx(E)+ΩgnTk(unφ)(1ψδ)dx.(F)

    By (13),

    (A)={|unφ|<k}a(x,un,un)un(1ψδ)dx{|unφ|<k}a(x,un,un)φ(1ψδ)dx,

    while

    {|unφ|<k}a(x,un,un)φ(1ψδ)dx={|uφ|<k}a(x,u,u)φdx+ω(n,δ).

    The Fatou lemma implies

    {|uφ|<k}a(x,u,u)udxlimninf{|unφ|<k}a(x,un,un)undx.

    Using (13), (51), we have

    (B)=Ωa(x,u,u)ψδTk(uφ)dx+ω(n)=ω(n,δ).

    While

    (F)=ΩgTk(uφ)dx+ω(n,δ),

    and

    |(D)|+|(E)|=Ω(f+n+fn)Tk(unφ)(1ψδ)dxkΩ(f+n+fn)(1ψδ)dx=ω(n,δ).

    So that we only need to deal with (C). Let m>k+φL(Ω) be fixed,

    (C)={2mun2m}|un|q1unTk(unφ)(1ψδ)dx(G)+k{un>2m}uqn(1ψδ)dx+k{un<2m}|un|q(1ψδ)dx.(H)

    By (52) and (53), we get

    (H)=ω(n,m,δ),

    and

    (G)=Ω|u|q1uTk(uφ)(1ψδ)dx+ω(n,m)=Ω|u|q1uTk(uφ)dx+ω(n,m,δ).

    Summing up the result of (A)-(H), we have

    Ωa(x,u,u)Tk(uφ)dx+Ω|u|q1uTk(uφ)dxΩgTk(uφ)dx.

    Thus u is the entropy solution of (8).

    Finally we prove (10). Choose φCc(Ω) as test function in the weak formulation of (9), we get

    Ωa(x,un,un)φdx+Ω|un|q1unφdx=Ω(fn+gn)φdx.

    Thanks to the assumptions of fn, gn and by (51),

    limn+Ω|un|q1unφdx=Ωa(x,u,u)φdx+Ωgφdx+Ωφdλ. (54)

    Since the entropy solution of (8) is also a distributional solution of the same problem, for the same φ,

    Ωa(x,u,u)φdx+Ω|u|q1uφdx=Ωgφdx. (55)

    Together with (54) and (55), we find

    limn+Ω|un|q1unφdx=Ω|u|q1uφdx+Ωφdλ.

    Thus (11) holds for every φCc(Ω). Since |un|q1un is bounded in L1(Ω), (11) can be extended by density to the functions in Cc(Ω).

    The authors also would like to thank the anonymous referees for their valuable comments which has helped to improve the paper.



    [1] A. J. Durnford, W. Rodgers, F. J. Kirkham, M. A. Mullee, A. Whitney, M. Prevett, et al., Very good inter-rater reliability of engel and ilae epilepsy surgery outcome classifications in a series of 76 patients, Seizure, 20 (2011), 809–812. doi: 10.1016/j.seizure.2011.08.004. doi: 10.1016/j.seizure.2011.08.004
    [2] E. Howell, Epilepsy stigma: Moving from a global problem to a global solution, Seizure-Eur. J. Epilepsy, 19 (2010), 628–629. doi: 10.1016/j.seizure.2010.10.016. doi: 10.1016/j.seizure.2010.10.016
    [3] M. Yildiz, E. Bergİl, The investigation of channel selection effects on epileptic analysis of eeg signals, Balk. J. Electr. Comput. Eng., 3 (2015), 236–241.
    [4] A. Matin, R. A. Bhuiyan, S. R. Shafi, A. K. Kundu, M. U. Islam, A hybrid scheme using pca and ica based statistical feature for epileptic seizure recognition from eeg signal, in 2019 Joint 8th International Conference on Informatics, Electronics & Vision (ICIEV) and 2019 3rd International Conference on Imaging, Vision & Pattern Recognition (icIVPR), IEEE, (2019), 301–306. doi: 10.1109/ICIEV.2019.8858573.
    [5] Y. Liu, B. Jiang, J. Feng, J. Hu, H. Zhang, Classification of EEG signals for epileptic seizures using feature dimension reduction algorithm based on LPP, Multimedia Tools Appl., 80 (2020), 30261-–30282. doi: 10.1007/s11042-020-09135-7. doi: 10.1007/s11042-020-09135-7
    [6] J. Birjandtalab, M. B. Pouyan, M. Nourani, Nonlinear dimension reduction for eeg-based epileptic seizure detection, in 2016 IEEE-EMBS International Conference on Biomedical and Health Informatics (BHI), IEEE, (2016), 595–598. doi: 10.1109/BHI.2016.7455968.
    [7] Q. Hou, Y. Liu, J. Liu, S. Sun, Epilepsy detection using random forest classification based on locally linear embedding algorithm, in 2020 5th International Conference on Control, Robotics and Cybernetics (CRC), IEEE, (2020), 202–206. doi: 10.1109/CRC51253.2020.9253455.
    [8] K. C. Chua, V. Chandran, R. Acharya, C. Lim, Automatic identification of epilepsy by hos and power spectrum parameters using eeg signals: A comparative study, in 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, IEEE, (2008), 3824–3827. doi: 10.1109/IEMBS.2008.4650043.
    [9] Z. Yan, A. Miyamoto, Z. Jiang, Frequency slice wavelet transform for transient vibration response analysis, Mech. Syst. Signal Process., 23 (2009), 1474–1489. doi: 10.1016/j.ymssp.2009.01.008. doi: 10.1016/j.ymssp.2009.01.008
    [10] Z. Yan, T. Tao, Z. Jiang, H. Wang, Discrete frequency slice wavelet transform, Mech. Syst. Signal Process., 96 (2017), 385–392. doi: 10.1016/j.ymssp.2017.04.019. doi: 10.1016/j.ymssp.2017.04.019
    [11] H. Hotelling, Relations between two sets of variates, in Breakthroughs in statistics, Springer, (1992), 162–190. doi: 10.1007/978-1-4612-4380-9_14.
    [12] Q. S. Sun, S. G. Zeng, Y. Liu, P. A. Heng, D. S. Xia, A new method of feature fusion and its application in image recognition, Pattern Recognition, 38 (2005), 2437–2448. doi: 10.1016/j.patcog.2004.12.013. doi: 10.1016/j.patcog.2004.12.013
    [13] T. Sun, S. Chen, J. Yang, P. Shi, A supervised combined feature extraction method for recognition, in Procedings of the IEEE International Conference on Data Mining, Pisa, Italy, Citeseer, (2008), 1043–1048.
    [14] T. Melzer, M. Reiter, H. Bischof, Appearance models based on kernel canonical correlation analysis, Pattern Recognit., 36 (2003), 1961–1971. doi: 10.1016/S0031-3203(03)00058-X. doi: 10.1016/S0031-3203(03)00058-X
    [15] T. Sun, S. Chen, Locality preserving cca with applications to data visualization and pose estimation, Image Vision Comput., 25 (2007), 531–543. doi: 10.1016/j.imavis.2006.04.014. doi: 10.1016/j.imavis.2006.04.014
    [16] S. D. Hou, Q. S. Sun, D. S. Xia, Supervised locality preserving canonical correlation analysis algorithm, Pattern Recognit. Artif. Intell., 2012, 143–149.
    [17] R. Andrzejak, K. Lehnertz, F. Mormann, C. Rieke, P. David, C. Elger, Indications of nonlinear deterministic and finite-dimensional structures in time series of brain electrical activity: Dependence on recording region and brain state, Phys. Rev. E, 2001, 061907. doi: 10.1103/PhysRevE.64.061907. doi: 10.1103/PhysRevE.64.061907
    [18] A. Shoeb, J. Guttat, Application of machine learning To epileptic seizure detection, in International Conference on Machine Learning, 2010.
    [19] R. W. Homan, J. Herman, P. Purdy, Cerebral location of international 10–20 system electrode placement, Electroencephalogr. Clin. Neurophysiol., 66 (1987), 376–382. doi: 10.1016/0013-4694(87)90206-9. doi: 10.1016/0013-4694(87)90206-9
    [20] R. J. Martis, U. R. Acharya, J. H. Tan, A. Petznick, L. Tong, C. K. Chua, et al., Application of intrinsic time-scale decomposition (itd) to eeg signals for automated seizure prediction, Int. J. Neural Syst., 23 (2013), 1350023. doi: 10.1142/S0129065713500238. doi: 10.1142/S0129065713500238
    [21] N. Ahammad, T. Fathima, P. Joseph, Detection of epileptic seizure event and onset using EEG, BioMed. Res. Int., 2014. doi: 10.1155/2014/450573. doi: 10.1155/2014/450573
    [22] R. J. Martis, J. H. Tan, C. K. Chua, T. C. Loon, S. W. J. YEO, L. Tong, Epileptic eeg classification using nonlinear parameters on different frequency bands, J. Mech. Med. Biol., 15 (2015), 1550040. doi: 10.1142/S0219519415500402. doi: 10.1142/S0219519415500402
    [23] N. S. Tawfik, S. M. Youssef, M. Kholief, A hybrid automated detection of epileptic seizures in eeg records, Comput. Electr. Eng., 53 (2016), 177–190. doi:10.1016/j.compeleceng.2015.09.001. doi: 10.1016/j.compeleceng.2015.09.001
    [24] E. Kabir, Y. Zhang, Epileptic seizure detection from eeg signals using logistic model trees, Brain Inf., 3 (2016), 93–100. doi: 10.1007/s40708-015-0030-2. doi: 10.1007/s40708-015-0030-2
    [25] R. R. Sharma, R. B. Pachori, Time-frequency representation using ievdhm-ht with application to classification of epileptic eeg signals, IET Sci., Meas. Technol., 12 (2018), 72–82.
    [26] V. Gupta, R. B. Pachori, Epileptic seizure identification using entropy of fbse based eeg rhythms, Biomed. Signal Process. Control, 53 (2019), 101569. doi: 10.1016/j.bspc.2019.101569. doi: 10.1016/j.bspc.2019.101569
    [27] H. Al-Hadeethi, S. Abdulla, M. Diykh, R. C. Deo, J. H. Green, Adaptive boost ls-svm classification approach for time-series signal classification in epileptic seizure diagnosis applications, Expert Syst. Appl., 161 (2020), 113676. doi: 10.1016/j.eswa.2020.113676. doi: 10.1016/j.eswa.2020.113676
    [28] N. Rafiuddin, Y. U. Khan, O. Farooq, Feature extraction and classification of eeg for automatic seizure detection, in 2011 International Conference on Multimedia, Signal Processing and Communication Technologies, IEEE, (2011), 184–187. doi: 10.1109/MSPCT.2011.6150470.
    [29] Y. U. Khan, N. Rafiuddin, O. Farooq, Automated seizure detection in scalp eeg using multiple wavelet scales, in 2012 IEEE international conference on signal processing, computing and control, IEEE, (2012), 1–5. doi: 10.1109/ISPCC.2012.6224361.
    [30] M. Zabihi, S. Kiranyaz, A. B. Rad, A. K. Katsaggelos, M. Gabbouj, T. Ince, Analysis of high-dimensional phase space via poincaré section for patient-specific seizure detection, IEEE Trans. Neural Syst. Rehabil. Eng., 24 (2015), 386–398. doi: 10.1109/TNSRE.2015.2505238. doi: 10.1109/TNSRE.2015.2505238
    [31] P. Thodoroff, J. Pineau, A. Lim, Learning robust features using deep learning for automatic seizure detection, in Machine learning for healthcare conference, PMLR, (2016), 178–190.
    [32] M. Z. Ahmad, A. M. Kamboh, S. Saleem, A. A. Khan, Mallat's scattering transform based anomaly sensing for detection of seizures in scalp EEG, IEEE Access, 5 (2017), 16919–16929. doi: 10.1109/ACCESS.2017.2736014. doi: 10.1109/ACCESS.2017.2736014
    [33] D. Chen, S. Wan, J. Xiang, F. S. Bao, A high-performance seizure detection algorithm based on discrete wavelet transform (dwt) and EEG, PloS one, 12 (2017), e0173138. doi: 10.1371/journal.pone.0173138. doi: 10.1371/journal.pone.0173138
    [34] J. Bonello, L. Garg, G. Garg, E. E. Audu, Effective data acquisition for machine learning algorithm in eeg signal processing, in Soft Computing: Theories and Applications, Springer, (2018), 233–244. doi: 10.1007/978-981-10-5699-4_23.
    [35] K. M. Tsiouris, S. Markoula, S. Konitsiotis, D. D. Koutsouris, D. I. Fotiadis, A robust unsupervised epileptic seizure detection methodology to accelerate large EEG database evaluation, Biomed. Signal Process. Control, 40 (2018), 275–285. doi: 10.1016/j.bspc.2017.09.029. doi: 10.1016/j.bspc.2017.09.029
    [36] M. Zhou, C. Tian, R. Cao, B. Wang, Y. Niu, T. Hu, et al., Epileptic seizure detection based on EEG signals and CNN, Front. Neuroinf., 12 (2018), 95. doi: 10.3389/fninf.2018.00095. doi: 10.3389/fninf.2018.00095
    [37] M. B. Ahmadi, A. Craik, H. F. Azgomi, J. T. Francis, J. L. Contreras-Vidal, R. T. Faghih, Real-time seizure state tracking using two channels: A mixed-filter approach, in 2019 53rd Asilomar Conference on Signals, Systems, and Computers, IEEE, (2019), 2033–2039. doi: 10.1109/IEEECONF44664.2019.9048990.
    [38] J. Wu, T. Zhou, T. Li, Detecting epileptic seizures in eeg signals with complementary ensemble empirical mode decomposition and extreme gradient boosting, Entropy, 22 (2020), 140. doi: 10.3390/e22020140. doi: 10.3390/e22020140
    [39] S. Chakrabarti, A. Swetapadma, A. Ranjan, P. K. Pattnaik, Time domain implementation of pediatric epileptic seizure detection system for enhancing the performance of detection and easy monitoring of pediatric patients, Biomed. Signal Process. Control, 59 (2020), 101930. doi: 10.1016/j.bspc.2020.101930. doi: 10.1016/j.bspc.2020.101930
    [40] L. A. Moctezuma, M. Molinas, EEG channel-selection method for epileptic-seizure classification based on multi-objective optimization, Front. Neurosci., 14 (2020), 593. doi: 10.3389/fnins.2020.00593. doi: 10.3389/fnins.2020.00593
    [41] L. A. Moctezuma, M. Molinas, Classification of low-density eeg for epileptic seizures by energy and fractal features based on emd, J. Biomed. Res., 34 (2020), 180–190. doi: 10.7555/JBR.33.20190009. doi: 10.7555/JBR.33.20190009
  • This article has been cited by:

    1. A.M. Elaiw, A.J. Alsaedi, A.D. Hobiny, S. Aly, Stability of a delayed SARS-CoV-2 reactivation model with logistic growth and adaptive immune response, 2023, 616, 03784371, 128604, 10.1016/j.physa.2023.128604
    2. D. Burini, N. Chouhad, Virus models in complex frameworks: Towards modeling space patterns of SARS-CoV-2 epidemics, 2022, 32, 0218-2025, 2017, 10.1142/S0218202522500476
    3. N. Bellomo, F. Brezzi, M. A. J. Chaplain, New trends of mathematical sciences towards modeling virus pandemics in a globally connected world, 2022, 32, 0218-2025, 1923, 10.1142/S0218202522010011
    4. A. M. Elaiw, A. S. Shflot, A. D. Hobiny, Stability analysis of SARS-CoV-2/HTLV-I coinfection dynamics model, 2022, 8, 2473-6988, 6136, 10.3934/math.2023310
    5. A. M. Elaiw, N. H. AlShamrani, E. Dahy, A. A. Abdellatif, Aeshah A. Raezah, Effect of Macrophages and Latent Reservoirs on the Dynamics of HTLV-I and HIV-1 Coinfection, 2023, 11, 2227-7390, 592, 10.3390/math11030592
    6. Ahmed M. Elaiw, Abdullah J. Alsaedi, Afnan Diyab Al Agha, Aatef D. Hobiny, Global Stability of a Humoral Immunity COVID-19 Model with Logistic Growth and Delays, 2022, 10, 2227-7390, 1857, 10.3390/math10111857
    7. Giulia Bertaglia, Chuan Lu, Lorenzo Pareschi, Xueyu Zhu, Asymptotic-Preserving Neural Networks for multiscale hyperbolic models of epidemic spread, 2022, 32, 0218-2025, 1949, 10.1142/S0218202522500452
    8. Ahmed M. Elaiw, Raghad S. Alsulami, Aatef D. Hobiny, Modeling and Stability Analysis of Within-Host IAV/SARS-CoV-2 Coinfection with Antibody Immunity, 2022, 10, 2227-7390, 4382, 10.3390/math10224382
    9. Ryan Weightman, Anthony Sbarra, Benedetto Piccoli, Coupling compartmental models with Markov chains and measure evolution equations to capture virus mutability, 2022, 32, 0218-2025, 2099, 10.1142/S0218202522500506
    10. Jiying Ma, Shasha Ma, Dynamics of a stochastic hepatitis B virus transmission model with media coverage and a case study of China, 2022, 20, 1551-0018, 3070, 10.3934/mbe.2023145
    11. Marcelo Bongarti, Luke Diego Galvan, Lawford Hatcher, Michael R. Lindstrom, Christian Parkinson, Chuntian Wang, Andrea L. Bertozzi, Alternative SIAR models for infectious diseases and applications in the study of non-compliance, 2022, 32, 0218-2025, 1987, 10.1142/S0218202522500464
    12. N. Bellomo, M. Esfahanian, V. Secchini, P. Terna, What is life? Active particles tools towards behavioral dynamics in social-biology and economics, 2022, 43, 15710645, 189, 10.1016/j.plrev.2022.10.001
    13. Ahmed M. Elaiw, Abdulsalam S. Shflot, Aatef D. Hobiny, Global Stability of Delayed SARS-CoV-2 and HTLV-I Coinfection Models within a Host, 2022, 10, 2227-7390, 4756, 10.3390/math10244756
    14. Henrique A. Tórtura, José F. Fontanari, The synergy between two threats: Disinformation and COVID-19, 2022, 32, 0218-2025, 2077, 10.1142/S021820252250049X
    15. Juan Pablo Agnelli, Bruno Buffa, Damián Knopoff, Germán Torres, A Spatial Kinetic Model of Crowd Evacuation Dynamics with Infectious Disease Contagion, 2023, 85, 0092-8240, 10.1007/s11538-023-01127-6
    16. A.M. Elaiw, A.J. Alsaedi, A.D. Hobiny, Global stability of a delayed SARS-CoV-2 reactivation model with logistic growth, antibody immunity and general incidence rate, 2022, 61, 11100168, 12475, 10.1016/j.aej.2022.05.034
    17. A. D. Al Agha, A. M. Elaiw, Global dynamics of SARS-CoV-2/malaria model with antibody immune response, 2022, 19, 1551-0018, 8380, 10.3934/mbe.2022390
    18. Diletta Burini, Nadia Chouhad, Nicola Bellomo, Waiting for a Mathematical Theory of Living Systems from a Critical Review to Research Perspectives, 2023, 15, 2073-8994, 351, 10.3390/sym15020351
    19. Ahmed M. Elaiw, Afnan D. Al Agha, Global Stability of a Reaction–Diffusion Malaria/COVID-19 Coinfection Dynamics Model, 2022, 10, 2227-7390, 4390, 10.3390/math10224390
    20. A. M. Elaiw, Raghad S. Alsulami, A. D. Hobiny, Global dynamics of IAV/SARS-CoV-2 coinfection model with eclipse phase and antibody immunity, 2022, 20, 1551-0018, 3873, 10.3934/mbe.2023182
    21. Ali Algarni, Afnan D. Al Agha, Aisha Fayomi, Hakim Al Garalleh, Kinetics of a Reaction-Diffusion Mtb/SARS-CoV-2 Coinfection Model with Immunity, 2023, 11, 2227-7390, 1715, 10.3390/math11071715
    22. Matthew O. Adewole, Farah A. Abdullah, Majid K. M. Ali, 2024, 3203, 0094-243X, 030007, 10.1063/5.0225272
    23. Elsayed Dahy, Ahmed M. Elaiw, Aeshah A. Raezah, Hamdy Z. Zidan, Abd Elsattar A. Abdellatif, Global Properties of Cytokine-Enhanced HIV-1 Dynamics Model with Adaptive Immunity and Distributed Delays, 2023, 11, 2079-3197, 217, 10.3390/computation11110217
    24. Nicola Bellomo, Jie Liao, Annalisa Quaini, Lucia Russo, Constantinos Siettos, Human behavioral crowds review, critical analysis and research perspectives, 2023, 33, 0218-2025, 1611, 10.1142/S0218202523500379
    25. Diletta Burini, Damian A. Knopoff, Epidemics and society — A multiscale vision from the small world to the globally interconnected world, 2024, 34, 0218-2025, 1567, 10.1142/S0218202524500295
    26. Ahmed M. Elaiw, Aeshah A. Raezah, Matuka A. Alshaikh, Global Dynamics of Viral Infection with Two Distinct Populations of Antibodies, 2023, 11, 2227-7390, 3138, 10.3390/math11143138
    27. Nicola Bellomo, Massimo Egidi, From Herbert A. Simon’s legacy to the evolutionary artificial world with heterogeneous collective behaviors, 2024, 34, 0218-2025, 145, 10.1142/S0218202524400049
    28. Ryan Weightman, Benedetto Piccoli, 2024, Chapter 7, 978-3-031-56793-3, 157, 10.1007/978-3-031-56794-0_7
    29. Luca Serena, 2023, Methodological Aspects of Multilevel Modeling and Simulation, 979-8-3503-3784-6, 111, 10.1109/DS-RT58998.2023.00025
    30. Aeshah A. Raezah, A.M. Elaiw, M.A. Alshaikh, Global stability of secondary DENV infection models with non-specific and strain-specific CTLs, 2024, 10, 24058440, e25391, 10.1016/j.heliyon.2024.e25391
    31. Nicola Bellomo, Raluca Eftimie, Guido Forni, What is the in-host dynamics of the SARS-CoV-2 virus? A challenge within a multiscale vision of living systems, 2024, 19, 1556-1801, 655, 10.3934/nhm.2024029
    32. A. M. Elaiw, E. A. Almohaimeed, A. D. Hobiny, Stability of HHV-8 and HIV-1 co-infection model with latent reservoirs and multiple distributed delays, 2024, 9, 2473-6988, 19195, 10.3934/math.2024936
    33. Yunfeng Xiong, Chuntian Wang, Yuan Zhang, Tom Britton, Interacting particle models on the impact of spatially heterogeneous human behavioral factors on dynamics of infectious diseases, 2024, 20, 1553-7358, e1012345, 10.1371/journal.pcbi.1012345
    34. Ahmed M. Elaiw, Ghadeer S. Alsaadi, Aatef D. Hobiny, Global co-dynamics of viral infections with saturated incidence, 2024, 9, 2473-6988, 13770, 10.3934/math.2024671
    35. Matthew O. Adewole, Taye Samuel Faniran, Farah A. Abdullah, Majid K.M. Ali, COVID-19 dynamics and immune response: Linking within-host and between-host dynamics, 2023, 173, 09600779, 113722, 10.1016/j.chaos.2023.113722
    36. Giulia Bertaglia, Andrea Bondesan, Diletta Burini, Raluca Eftimie, Lorenzo Pareschi, Giuseppe Toscani, New trends on the systems approach to modeling SARS-CoV-2 pandemics in a globally connected planet, 2024, 34, 0218-2025, 1995, 10.1142/S0218202524500301
    37. Nicola Bellomo, Seung-Yeal Ha, Jie Liao, Wook Yoon, Behavioral swarms: A mathematical theory toward swarm intelligence, 2024, 34, 0218-2025, 2305, 10.1142/S0218202524500490
    38. Mohamed Zagour, 2024, Chapter 6, 978-3-031-56793-3, 127, 10.1007/978-3-031-56794-0_6
    39. Nisrine Outada, A forward look to perspectives, 2023, 47, 15710645, 133, 10.1016/j.plrev.2023.10.011
    40. Luca Serena, Moreno Marzolla, Gabriele D’Angelo, Stefano Ferretti, A review of multilevel modeling and simulation for human mobility and behavior, 2023, 127, 1569190X, 102780, 10.1016/j.simpat.2023.102780
    41. Christian Parkinson, Weinan Wang, Analysis of a Reaction-Diffusion SIR Epidemic Model with Noncompliant Behavior, 2023, 83, 0036-1399, 1969, 10.1137/23M1556691
    42. D. Burini, N. Chouhad, Cross-diffusion models in complex frameworks from microscopic to macroscopic, 2023, 33, 0218-2025, 1909, 10.1142/S0218202523500458
    43. Ahmed M. Elaiw, Raghad S. Alsulami, Aatef D. Hobiny, Global properties of SARS‐CoV‐2 and IAV coinfection model with distributed‐time delays and humoral immunity, 2024, 47, 0170-4214, 9340, 10.1002/mma.10074
    44. B. Bellomo, M. Esfahanian, V. Secchini, P. Terna, From a mathematical science of living systems to biology and economics, 2023, 47, 15710645, 264, 10.1016/j.plrev.2023.11.002
    45. Ahmed M. Elaiw, Amani S. Alsulami, Aatef D. Hobiny, Global properties of delayed models for SARS-CoV-2 infection mediated by ACE2 receptor with humoral immunity, 2024, 9, 2473-6988, 1046, 10.3934/math.2024052
    46. Bishal Chhetri, Krishna Kiran Vamsi Dasu, Stability and bifurcation analysis of a nested multi-scale model for COVID-19 viral infection, 2024, 12, 2544-7297, 10.1515/cmb-2024-0006
    47. Vinicius V. L. Albani, Jorge P. Zubelli, Stochastic transmission in epidemiological models, 2024, 88, 0303-6812, 10.1007/s00285-023-02042-z
    48. Hyeong-Ohk Bae, Seung Yeon Cho, Jane Yoo, Seok-Bae Yun, Mathematical modeling of trend cycle: Fad, fashion and classic, 2024, 01672789, 134500, 10.1016/j.physd.2024.134500
    49. Juan Pablo Agnelli, Claudio Armas, Damián A. Knopoff, Spatial Kinetic Modeling of Crowd Evacuation: Coupling Social Behavior and Infectious Disease Contagion, 2025, 17, 2073-8994, 123, 10.3390/sym17010123
    50. Gabriel Benedetti, Ryan Weightman, Benedetto Piccoli, Optimizing overlapping non-pharmaceutical interventions with a socio-demographic model, 2025, 1972-6724, 10.1007/s40574-025-00477-4
    51. Jorge P Zubelli, Jennifer Loria, Vinicius V L Albani, On the estimation of the time-dependent transmission rate in epidemiological models, 2025, 41, 0266-5611, 065001, 10.1088/1361-6420/add55b
  • Reader Comments
  • © 2022 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(3092) PDF downloads(127) Cited by(12)

Figures and Tables

Figures(2)  /  Tables(7)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog