Processing math: 100%
Research article Special Issues

Transcranial Doppler analysis based on computer and artificial intelligence for acute cerebrovascular disease


  • Cerebrovascular disease refers to damage to brain tissue caused by impaired intracranial blood circulation. It usually presents clinically as an acute nonfatal event and is characterized by high morbidity, disability, and mortality. Transcranial Doppler (TCD) ultrasonography is a non-invasive method for the diagnosis of cerebrovascular disease that uses the Doppler effect to detect the hemodynamic and physiological parameters of the major intracranial basilar arteries. It can provide important hemodynamic information that cannot be measured by other diagnostic imaging techniques for cerebrovascular disease. And the result parameters of TCD ultrasonography such as blood flow velocity and beat index can reflect the type of cerebrovascular disease and serve as a basis to assist physicians in the treatment of cerebrovascular diseases. Artificial intelligence (AI) is a branch of computer science which is used in a wide range of applications in agriculture, communications, medicine, finance, and other fields. In recent years, there are much research devoted to the application of AI to TCD. The review and summary of related technologies is an important work to promote the development of this field, which can provide an intuitive technical summary for future researchers. In this paper, we first review the development, principles, and applications of TCD ultrasonography and other related knowledge, and briefly introduce the development of AI in the field of medicine and emergency medicine. Finally, we summarize in detail the applications and advantages of AI technology in TCD ultrasonography including the establishment of an examination system combining brain computer interface (BCI) and TCD ultrasonography, the classification and noise cancellation of TCD ultrasonography signals using AI algorithms, and the use of intelligent robots to assist physicians in TCD ultrasonography and discuss the prospects for the development of AI in TCD ultrasonography.

    Citation: Lingli Gan, Xiaoling Yin, Jiating Huang, Bin Jia. Transcranial Doppler analysis based on computer and artificial intelligence for acute cerebrovascular disease[J]. Mathematical Biosciences and Engineering, 2023, 20(2): 1695-1715. doi: 10.3934/mbe.2023077

    Related Papers:

    [1] Yong Ding, Jian-Hong Liu . The signature lncRNAs associated with the lung adenocarcinoma patients prognosis. Mathematical Biosciences and Engineering, 2020, 17(2): 1593-1603. doi: 10.3934/mbe.2020083
    [2] Rui Huang, Xiwen Liao, Qiaochuan Li . Integrative genomic analysis of a novel small nucleolar RNAs prognostic signature in patients with acute myelocytic leukemia. Mathematical Biosciences and Engineering, 2022, 19(3): 2424-2452. doi: 10.3934/mbe.2022112
    [3] Siqi Hu, Fang Wang, Junjun Yang, Xingxiang Xu . Elevated ADAR expression is significantly linked to shorter overall survival and immune infiltration in patients with lung adenocarcinoma. Mathematical Biosciences and Engineering, 2023, 20(10): 18063-18082. doi: 10.3934/mbe.2023802
    [4] Wei Niu, Lianping Jiang . A seven-gene prognostic model related to immune checkpoint PD-1 revealing overall survival in patients with lung adenocarcinoma. Mathematical Biosciences and Engineering, 2021, 18(5): 6136-6154. doi: 10.3934/mbe.2021307
    [5] Jinqi He, Wenjing Zhang, Faxiang Li, Yan Yu . Development of metastasis-associated seven gene signature for predicting lung adenocarcinoma prognosis using single-cell RNA sequencing data. Mathematical Biosciences and Engineering, 2021, 18(5): 5959-5977. doi: 10.3934/mbe.2021298
    [6] Yuan Yang, Lingshan Zhou, Xi Gou, Guozhi Wu, Ya Zheng, Min Liu, Zhaofeng Chen, Yuping Wang, Rui Ji, Qinghong Guo, Yongning Zhou . Comprehensive analysis to identify DNA damage response-related lncRNA pairs as a prognostic and therapeutic biomarker in gastric cancer. Mathematical Biosciences and Engineering, 2022, 19(1): 595-611. doi: 10.3934/mbe.2022026
    [7] Gaozhong Sun, Tongwei Zhao . Lung adenocarcinoma pathology stages related gene identification. Mathematical Biosciences and Engineering, 2020, 17(1): 737-746. doi: 10.3934/mbe.2020038
    [8] Yi Liu, Long Cheng, Xiangyang Song, Chao Li, Jiantao Zhang, Lei Wang . A TP53-associated immune prognostic signature for the prediction of the overall survival and therapeutic responses in pancreatic cancer. Mathematical Biosciences and Engineering, 2022, 19(1): 191-208. doi: 10.3934/mbe.2022010
    [9] Ji-Ming Wu, Wang-Ren Qiu, Zi Liu, Zhao-Chun Xu, Shou-Hua Zhang . Integrative approach for classifying male tumors based on DNA methylation 450K data. Mathematical Biosciences and Engineering, 2023, 20(11): 19133-19151. doi: 10.3934/mbe.2023845
    [10] Pingping Song, Jing Chen, Xu Zhang, Xiaofeng Yin . Construction of competitive endogenous RNA network related to circular RNA and prognostic nomogram model in lung adenocarcinoma. Mathematical Biosciences and Engineering, 2021, 18(6): 9806-9821. doi: 10.3934/mbe.2021481
  • Cerebrovascular disease refers to damage to brain tissue caused by impaired intracranial blood circulation. It usually presents clinically as an acute nonfatal event and is characterized by high morbidity, disability, and mortality. Transcranial Doppler (TCD) ultrasonography is a non-invasive method for the diagnosis of cerebrovascular disease that uses the Doppler effect to detect the hemodynamic and physiological parameters of the major intracranial basilar arteries. It can provide important hemodynamic information that cannot be measured by other diagnostic imaging techniques for cerebrovascular disease. And the result parameters of TCD ultrasonography such as blood flow velocity and beat index can reflect the type of cerebrovascular disease and serve as a basis to assist physicians in the treatment of cerebrovascular diseases. Artificial intelligence (AI) is a branch of computer science which is used in a wide range of applications in agriculture, communications, medicine, finance, and other fields. In recent years, there are much research devoted to the application of AI to TCD. The review and summary of related technologies is an important work to promote the development of this field, which can provide an intuitive technical summary for future researchers. In this paper, we first review the development, principles, and applications of TCD ultrasonography and other related knowledge, and briefly introduce the development of AI in the field of medicine and emergency medicine. Finally, we summarize in detail the applications and advantages of AI technology in TCD ultrasonography including the establishment of an examination system combining brain computer interface (BCI) and TCD ultrasonography, the classification and noise cancellation of TCD ultrasonography signals using AI algorithms, and the use of intelligent robots to assist physicians in TCD ultrasonography and discuss the prospects for the development of AI in TCD ultrasonography.



    The susceptible-infective-removed-susceptible (SIRS) model is one of the most popular epidemic models, in which total host population is divide into three classes called susceptible(S), infective(I) and removed(R). Since the Kermack and McKendrick pioneering work [1], a large number of studies have been carried out on the model of continuous infectious diseases attempting to gain a better understanding of diseases transmission, especially for the control policies and dynamics [2,3,4,5]. Recent years, many researchers have studied the epidemic with SIRS models such as COVID-19(2019), hand-foot-and-mouth diseases(HFMD) and other infectious diseases recently. Although scholars have flocked to study many properties of the SIRS models, there are still many aspects that deserve further study, such as the threshold dynamics of SIRS model with vaccination and optimal control using dynamic programming method.

    In studying the dynamics of the disease, one of the most important concepts is the basic reproduction number R0. It is epidemiologically defined as the expected number of secondary cases produced by a typical infectious individual during its entire infectious period in a completely susceptible host population. To investigate the global behavior of the prevalence of infectious diseases, the stability analysis of equilibrium for epidemic models have been carried out (see [6,7] and the references therein). For ordinary differential SIRS models, Mena-Lorca and Hethcote considered several kinds of SIRS epidemic models and a threshold parameters were also found in [6] to determine whether the disease dies out or approaches to an endemic equilibrium, we also note that the global stability of SIRS epidemic model have already been investigated and the threshold theorems are well obtained. In the subsequent studies, it was generally believed that the introduction of age-structure into epidemic models is necessary and reasonable. Since the fact that the age structure of a population affects the dynamics of disease transmission was recognized, various age-structured epidemic models have been investigated by many authors. Some theoretics on the transmission dynamics of age-structured epidemic models have been developed in [8,9,10,11]. He et al. [12] studied the optimal birth control of age-dependent competitive species. Yang et al. [13] proposed epidemic model with age-since-infection and diffusion, the next generation operator R is also given as the basic regeneration numbers. In paper [14], Yang et al. also studied the threshold dynamics by the method of Lyapunov functionals for an age-of-infection epidemiological model with nonlinear incidence rate. Chekroun et al. [15] showed the stability of the equilibrium by using Lyapunov functions. Lan et al. [16] just studied the impact of hospital resources and environmental perturbation to the dynamics of SIRS model. But the effect of vaccination on threshold dynamics was not considered in these articles.

    In fact, we notice that vaccination has a significant impact on the thresholds dynamics of many infectious diseases, such as COVID-19, HPV, etc. Obviously, vaccination of susceptible individuals can greatly reduce the probability that susceptible individuals become infected. Thus, decision makers hope to achieve the optimal proportion of vaccines with the lowest cost. Optimal control is a promising strategy to control disease outbreaks. The main purpose of using the control in some diseases is to search for the most effective one that reduces the infection to a minimum level while minimizing the cost of applying control measures. In papers [7,9,17], the strategies such as health promotion campaigns or lockdown policies can contribute to reducing the disease transmission. In paper [18], Mu and Zhang investigated the near-optimal control of SIRS multi-strain epidemic model with age-structure using Pontryagin maximum principle. Bolzoni et al. [19] investigated optimal control of epidemic size and duration with limited resources. Zhou et al. [20] researched optimal isolation strategies of emerging infectious diseases with limited resources. In papers [21,22,23], some measures correspond to screening and quarantining of infected were implemented to achieve the purpose of controlling the disease.

    However, we note that the most of these articles by using maximum principle to show optimal controls and only present necessary conditions. Defectively, the method of maximum principle only addresses the optimal control problem with initial value case of t0=0. In fact, vaccination can be administered at any time t=t0 in the control of the disease. Surpassingly, the dynamic programming method can give the necessary and sufficient conditions for the existence of optimal control via the HJB equation at any initial times t=t0 and initial states. In other words, the dynamic programming method is an extension of the maximum principle, namely, the dynamic programming method is equivalent to maximum principle method when t0=0.

    Overall, the threshold dynamics results about age-structured SIRS epidemic models with the proportion of the vaccine injections are comparatively scarce, as well as few articles show that dynamic programming approach to solve the optimal control problem. Thus, in this paper, taking into account for the effect of the proportion of the vaccine injections, we incorporate the vaccine injection rate in solving the threshold of age-structured SIRS epidemic models. And vaccination optimal control is also obtained through dynamic programming approach.

    The layout of this paper is as follows: in section 2, an age-structure SIRS epidemic model is reproduced and new age-structure SIRS epidemic model with vaccination is given; in section 3, we give the basic regeneration number and discuss the threshold dynamics; in section 4, we formulate and solve the corresponding optimal control problem. The existence of optimal control is proved by HJB equation and the expression of optimal control is present. in section 5, some numerical simulations are performed to demonstrate the theoretical results. Brief conclusions are given in section 6.

    We reproduce the age-structured SIRS epidemic model based on paper [24] as follows

    {dS(a,t)=[S(a,t)a+λ(a)μ(a)S(a,t)β(a)S(a,t)I(a,t)+ρ(a)η(a)I(a,t)]dt,dI(a,t)=[I(a,t)a+β(a)S(a,t)I(a,t)δ(a)I(a,t))]dt,dR(a,t)=[R(a,t)a+(1ρ(a))η(a)I(a,t)(μ(a)+γ(a))R(a,t)]dt,S(a,0)=S0(a),I(a,0)=I0(a),R(a,0)=R0(a), (2.1)

    where S(a,0)=S0(a),I(a,0)=I0(a),R(a,0)=R0(a) are initial conditions, with the following conditions:

    {S(0,t)=A0β(a)S(a,t)da,I(0,t)=k1A0β(a)I(a,t)da,R(0,t)=k2A0β(a)R(a,t)da, (2.2)

    where k1,k2[0,1] present the weight of the infected and recovered, respectively. All the parameters are related to age a, positive and bounded, These parameters and their meaning are listed in the Table 1, where η(a)=σ(a)+α(a), and δ(a)=μ(a)+ϕ(a)+σ(a)+α(a).

    Table 1.  Parameter values used in numerical simulations.
    Notation Biological meanings
    S(a,t) Densities of susceptible individuals of age a at time t
    I(a,t) Densities of infective individuals of age a at time t
    R(a,t) Densities of recovery individuals of age a at time t
    λ(a) Recruitment rate of age a
    μ(a) Mortality rate of age a
    β(a) Transmission coefficient of age a
    γ(a) Immunity loss rate of age a
    δ(a) Infected removal rate of age a
    η(a) Infected recovery rate of age a
    ρ(a) Switch between SISSIR(S) of age a
    α(a) Treatment rate of age a
    σ(a) Nature cure rate of age a
    ϕ(a) Infected-induced mortality rate of age a

     | Show Table
    DownLoad: CSV

    As we all know that the age density distribution function of the total population reaches a stable state, i.e

    S(a,t)+I(a,t)+R(a,t)=P(a). (2.3)

    We focus on threshold dynamics on an age-structured SIRS epidemic model with vaccination. The spread of many epidemics are age-related, such as Human Papilloma Virus (HPV) spread across the population after ten years olds usually. And the World Health Organization believes that the most appropriate age for HPV vaccination was 11-12 years old. Therefore, it makes sense to study the vaccination is only given for a certain age to age-structured SIRS epidemic model.

    The introduced vaccination rate u(A0) satisfies 0u(A0)1 and is only given for a certain age and susceptible individuals of age A0 is inoculated at some point. The disease does not spread among people younger than age A0, that is, there are no patients before the vaccination age A0. Therefore, we can get that the system satisfied by the susceptible, infective, and recovering individuals based on model (2.1) as following

    {S(a,t)=P(a),I(a,t)=0,R(a,t)=0,0<a<A0,dS=(Sa+λ(a,t)μ(a)Sϱ(a,t)S+ρ(a)η(a)I)dt,A0<a<A,dI=(Ia+ϱ(a,t)Sδ(a)I)dt,A0<a<A,dR=(Ra+(1ρ(a))η(a)Iμ(a)Rdt,A0<a<A,S(A0,t)=(1u(A0))S(A0,t),I(A0,t)=I(A0,t)=0,R(A0,t)=u(A0)S(A0,t),S(a,0)=S0(a),I(a,0)=I0(a),R(a.0)=R0(a),ϱ(a,t)=k(a)AA0I(ˆa,t)dˆa, (2.4)

    we call infectivity

    ϱ(a,t)=β(a)I(a,t)=k(a)AA0I(ˆa,t)dˆa<,

    where S(A0,t) and I(A0,t) are the densities of susceptible and recovering individuals in the left neighborhood of age A0, k(a) is a positive continuous function. The vaccination rate u(A0) is a constant and k(a) is a nonnegative continuous function defined on 0aA, where A is the maximum age.

    In this paper, unless otherwise specified. The H represents the Hilbert space, with inner product

    φ1,ϕ1=A0φ1(a)ϕ1(a)da,φ1,ϕ1H, (2.5)

    and we consider the state equations on space L1(0,A;H3), when φ=(φ1,φ2,φ3) and ϕ=(ϕ1,ϕ2,ϕ3), φi,ϕiH, the definition of inner product on L1(0,A;H3) as

    φ,ϕ=A0(φ1(a)ϕ1(a)+φ2(a)ϕ2(a)+φ3(a)ϕ3(a))da=A0φ(a)ϕ(a)da,φ,ϕH3. (2.6)

    The existence of time-independent equilibrium solutions for the model (2.4) is discussed. Obviously, the model (2.4) has disease-free equilibrium

    S(a,t)=P(a),I(a,t)=0,R(a,t)=0,0<a<A0,S(a,t)=(1u(A0))P(a),I(a,t)=0,R(a,t)=u(A0)P(a),A0<a<A. (3.1)

    In order to obtain the endemic disease equilibrium solution, S(a,t)=S(a),I(a,t)=I(a) and R(a,t)=R(a) are substituted into the model (2.4) obtained

    R(a)=0,0<a<A0,R(a)=u(A0)P(a),A0<a<A. (3.2)

    And we know when R0(a)=u(A0)P(a), the R(a,t)=u(A0)P(a), and when the population reaches a dynamic equilibrium P(a), Combining the equations (2.3) and (3.2), and substitute into (2.4) for the equation that satisfies I(a,t). From (2.3), it shows that

    I(a)a+I(a)t=δ(a)I(a)+k(a)(P(a)I(a)R(a))D(t), (3.3)

    where D(t)=AA0I(ˆa,t)dˆa. Let ζ(a,D(t))=δ(a)+k(a)D(t), in order to solve the (3.3), when A0<a<A, the expression of I(a,t) is

    I(a,t)=I0(at)exp(t0ζ(at+θ,D(θ))dθ)+[1u(A0)]t0exp(tτζ(at+θ,D(θ))dθ)×k(at+θ)P(at+θ)D(τ)dτ,a>t+A0,I(a,t)=[1u(A0)]aA0exp(aτζ(θ,D(ta+θ))dθ)×k(τ)P(τ)D(ta+τ)dτ,at+A0. (3.4)

    we can substitute (3.4) into the D(t) yields when A0<a<A

    D(t)=[1u(A0)]AA0aA0k(τ)P(τ)D(tτ)×exp(a+τaζ(θ,D(taτ+θ))dθ)dadτ+aA0I0(a)exp(t0ζ(a+θ,D(θ))dθ)da. (3.5)

    Since

    limtaA0I0(a)exp(t0ζ(a+θ,D(θ))dθ)daaA0I0(a)exp(t0μ(a+θ,N)dθ)da=0,

    obviously, we have D(t)=D if I(a,t)=I as t is equilibrium of model, then

    D=[1u(A0)]AA0aA0k(τ)P(τ)exp(a+τaζ(θ,D)dθ)Ddadτ+aA0I0(a)exp(t0ζ(a+θ,D)dθ)da, (3.6)

    then

    D=D[1u(A0)]AA0aA0k(τ)P(τ)exp(a+τaδ(θ)dθ)×exp(Da+τak(θ)dθ)dτda. (3.7)

    Then it expresses I=0 or

    1=[1u(A0)]AA0aA0k(τ)P(τ)exp(a+τaδ(θ)dθ)×exp(Da+τak(θ)dθ)dτda. (3.8)

    It is quite obvious the right-hand side of the equation (3.8) is the minus function about D, thus, we can define

    R0=(1u(A0))A0aA0k(τ)p(τ)exp(aτδ(θ)dθ)dτda. (3.9)

    When R0<1, (3.7) has only a unique zero solution and no positive root, when R0>1, the (3.7) formula has zero root and the only positive root I>0 satisfies (3.8) formula, thus there are only positive root E=(S,I,R) of (2.4). For ease of analysis, we converse i(a,t)=I(a,t)(1u(A0))P, equation (2.4) about I(a,t) can be transformed into

    {i(a,t)t+i(a,t)a=[1i(a,t)][1u(A0)]k(a)AA0P(ˆa)i(ˆa,t)dˆaδ(a)i(a,t),a>A0,i(A0,t)=0,t0,i(a,0)=i0(a)=I0(a)(1u(A0))P,aA0. (3.10)

    From the transformation, we can see that the existence and stability of the positive equilibrium of model (2.4) is equivalent to positive equilibrium of (3.10) model. R0 can be used as a threshold for the existence or absence of endemic diseases. It is also a threshold for the disappearance of diseases. Let =11R0 and we call vaccination elimination point.

    Similar to the [25] method, the following conclusion can be obtained. Before discussing the global stability of equilibriums, we first propose the lemma 3.1, and the similar proof process have been presented in [26].

    Lemma 3.1. Let i1(a,t) and i2(a,t) are the solutions of the (3.10) equations, satisfying i1(a,0)=i10(a)i2(a,0)=i20(a), where i10(a),i20(a)K, then there are following properties

    (1) ij(a,t)K,j=1,2;

    (2) i1(a,t)i2(a,t),t0;

    (3) When 0<ξ<1, the solution i(a,t) of (3.10) that the initial value satisfies the i(a,0)=i10(a)ξ, then it must satisfy the i1(a,t)ξi(a,t),

    where the K is denoted

    K={f(a,t)|f(a,t)C,f(0,t)=0,0f(a,t)1}.

    Theorem 3.2. When R0<1, the model (3.10) has only a disease-free equilibrium, which is globally asymptotically stable, while when R0>1, the model (3.10) also has a unique endemic equilibrium E, which is globally asymptotically stable and the disease-free equilibrium is unstable.

    Proof. We only consider the asymptotic behavior of solutions where the initial values of the model (3.10) satisfy 0i(a,0)1 in (A0,A). The theorem of the global stability is given for solutions.

    It is not difficult to see if the model (3.10) has a positive equilibrium equivalent to the (3.5) has a positive root. Since the right of the (3.5) is reduced about D(t) and tends to zero when the D(t) tends to positive infinity. There is no positive real root to the (3.5) when R0<1, and there is a unique positive real root when R0>1, i.e. the model (2.4) has only a unique positive equilibrium when R0>1 and, instead, no positive equilibrium when R0<1.

    When R0<1, in order to study the global stability of the disease-free equilibrium, expression of the solution of (3.10) about I(a,t) is obtained by using the characteristic line method

    i(a,t)=i0(at)exp(t0[δ(at+τ)+[1u(A0)]k(at+τ)AA0P(ˆa)i(ˆa,τ)dˆa]dτ)+t0[1u(A0)]k(at+τ)AA0P(ˆa)i(ˆa,τ)dˆa×exp(tτ[δ(at+θ)+[1u(A0)]k(at+θ)AA0P(ˆa)i(ˆa,θ)dˆa]dθ)dτ,at+A0, (3.11)
    i(a,t)=aA0[1u(A0)]k(τ)AA0P(ˆa)i(ˆa,ta+τ)dˆa×exp(aτ[δ(θ)+[1u(A0)]k(θ)AA0P(ˆa)i(ˆa,ta+θ)]dθ)dτ,a<t+A0. (3.12)

    When t>A, we obtain

    i(a,t)=aA0[1u(A0)]k(τ)χ(ta+τ)×exp(aτ[δ(θ)+[1u(A0)]k(θ)χ(ta+θ)]dθ)dτ,a<t+A0, (3.13)

    where χ(t)=AA0P(ˆa)i(ˆa,t)dˆa, when t<A, the χ(t) satisfy is

    χ(t)=[1u(A0)]A0P(ˆa)ˆaA0k(τ)χ(tˆa+τ)×exp(ˆaτ[δ(θ)+[1u(A0)]k(θ)χ(tˆa+θ)]dθ)dτdˆa,a<t+A0. (3.14)

    For natural numbers n>1, let χn=maxnAt(n+1)Aχ(t),Ln=max{χn1,χn}. According to the (3.14) when R0<1, we get

    χ(t)[1u(A0)]A0P(ˆa)ˆaA0χ(tˆa+τ)exp(ˆaτδ(θ)dθ)dτdˆaR0maxtAτtχ(τ). (3.15)

    So we obtain χ(t)R0max{χn1,χn} when nAt(n+1)A. Because of R0<1, we also obtain χn<Ln=χn1, i.e. χn<χn1, χn is monotonous and reduced and χn<χn1R0, from this delivery limn+χn=0. Considering above conditions and (3.13) can also derive limt+i(a,t)=0. Finally, it is concluded from (3.13) that the disease-free equilibrium of the model (3.10) is stable R0<1.

    When R0>1, substituting i(a,t)=i(a) into the (3.10), the only endemic equilibrium of (3.10) is given by

    i(a)=[1u(A0)]aA0k(a)AA0P(ˆa)i(ˆa)dˆa×exp(aτδ(θ)+[1u(A0)]k(τ)AA0P(ˆa)i(ˆa)dˆadθ)dτ. (3.16)

    When A0a<A, the i(a,t) satisfy i(a,0)>0, and χ(t)>0 can be derived from (3.11) and (3.14) when 0tA, and we also know that i(a,A)>0(A0<a<A) by (3.13). Repeat this process then the χ(t)>0 is obtained when 0t3A, let χ0=minAt2Aχ(t), χ=maxAt2Aχ(t), the following is got as

    i(a,2A)[1u(A0)]aA0χ0k(τ)×exp(aτ(δ(θ)+[1u(A0)]χk(θ))dθ)dτ. (3.17)

    Comparing the forms of (3.16) and (3.17), we can find a positive number 0<ξ<1 such that i(a,2A)ξi(a), namely, the following formula holds for the solution i(a,t) of (2.4)

    ξi(a)i(a,2A)1.

    The Lemma 3.1 shows that the solutions iξ(a,t) and i1(a,t) with initial values of ξi(a,0) and 1 respectively when t=0, must have

    iξ(a,t)i(a,t+2A)i1(a,t).

    Based on the fact that i(a) is the equilibrium of (3.10), then the ξi(a)iξ(a,t), i.e. iξ(a,t) is a monotone increasing function about t. Since i1(a,t) is a monotone subtraction function about t, and the corresponding χξ(t) and χ1(t) are monotone increasing and monotone decreasing, respectively, and there are limtχξ(t)=χ0ξlimtχ1(t)=χ01. From (3.14), the equation has only unique normal number solutions when R0>1, i.e. χ0ξ=χ01=ξ, the χ(t) corresponding to the i(a,t) also obtain limtχ(t)=ξ. Reuse (3.13) when t>A0, can be obtained limti(a,t)=i(a), i.e. the local disease equilibrium of (3.10) is globally attractive. The stability of the local disease equilibrium can be obtained by linearization. Therefore, the local disease equilibrium is globally asymptotically stable.

    Remark 3.3. 3.2 shows that the epidemic disease will die out if R0<1, while the disease will be prevalent if R0>1. This implies that R0 is the threshold of model (2.4).

    In section 3, we consider the effect of a age determined vaccine immunization of the system on threshold dynamics, which is one of the cases in which vaccines are added to all susceptible persons in the system. In this section, we focus on the optimal control of (4.1) on the basis of age-structured SIRS by implementing proportion of vaccination u(a,t) as control strategy at all ages of system, and vaccination rate u(a,t) satisfies 0u(a,t)1. We aim to minimize the infected population, while keeping the cost of applying such control strategy at the minimum level.

    Let ((a,t0),x(a,t0))[0,T)×(H×H×H), and consider the following control system over on (a,t)(0,A)×[t0,T]:

    {S(a,t)t+S(a,t)a=λ(a)μ(a)S(a,t)β(a)S(a,t)I(a,t)+ρ(a)η(a)I(a,t)+γ(a)R(a,t)u(a,t)S(a,t),I(a,t)t+I(a,t)a=β(a)S(a,t)I(a,t)δ(a)I(a,t)),R(a,t)t+R(a,t)a=(1ρ(a))η(a)I(a,t)(μ(a)+γ(a))R(a,t)+u(a,t)S(a,t),S(a,t0)=St0(a),I(a,t0)=It0(a),R(a,t0)=Rt0(a), (4.1)

    where S(a,t0)=St0(a),I(a,t0)=It0(a),R(a,t0)=Rt0(a) is varying initial time value and denote the initial state x(a,t0) as xt0(a). Here the control {u(,):(0,A)×[t0,T]U}, and u(,) is measurable. We obtain the schematic diagram of model (4.1) with introducing vaccination control (see Figure 1).

    Figure 1.  Schematic diagram of the model (4.1).

    The object is to design the optimal controller u(a,t),(a,t)(0,A)×[t0,T] for the system (4.1) which minimizes the following performance index during the finite time interval [t0,T] [27,28,29,30]. Our goal is to minimize the total number of the infected and susceptible individuals by using minimal control efforts. The cost function as follows

    J(t0,xt0(a),u(a,t))=Tt0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt+A0h(I(a,T))da, (4.2)

    where I(a,t):(a,t)[0,A]×[t0,T] and we denote x=(S,I,R),xt0(H×H×H), the terminal value is V(x,T)=A0h(I(a,T))da. Let L(t,x,u):=τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t), where τ1, τ2 and τ3 are the positive weighting factors, representing the cost per unit time of the components I(a,t), u(a,t) and u2(a,t), respectively. In particular, Tt0A0τ1I(a,t)dadt is the cost that infected individual creates for the society due to standard medical care, not including the vaccination treatment u(a,t). Tt0A0(τ2u(a,t)+12τ3u2(a,t))dadt is the cost of treating infected individuals. h:H×H×HH are measurable and h(I(a,T)) denotes the total number disease in population depend on the age at the time T. The value function as

    {V(t0,xt0(a))=infu()UJ(t0,xt0(a),u(a,t)),t0[0,T),V(T,xt0(a))=A0h(I(a,T))da. (4.3)

    The function V(,) will play an important role in obtaining optimal control, thus, we would like to study V(,) in great detail. Let (x,u) is optimal pair, where x=(S,I,R). Noting that the initial time (t=0) and the initial state (x(a,0)=x0) are fixed in the model [18]. The h(I(a,T)) satisfies

    A0h(I(a,T))h(˜I(a,T))daCI(a,T)˜I(a,T)C, (4.4)

    where C is a positive constant, represents norm on space H and I(a,T), ˜I(a,T)H. The following hypothesis are put forward and need to be satisfied.

    D1 (U,d) is a separable metric space and T>0.

    D2 The control set U is bounded convex

    D3 Different controls correspond to same terminal value.

    Then we study dynamic programming method, to give optimal strategy of model (4.1). The basic idea of this approach applied to optimal controls is to consider a family of optimal control problems with different initial times (t=t00) and initial states (x(a,t0)=xt0(a)) [see [31,32,33]], to establish relationships among these problems through the HJB equation. Firstly, we present the Bellman's principle of optimality as follows.

    Theorem 4.1. Let (D1)-(D3) hold. Then for any (t0,xt0(a))[0,T)×H×H×H,

    V(t0,xt0(a))=infu()U[t0,T]{ˆt0t0A0L(t,x(t;t0,xt0(a),u(,)),u(a,t))dt+V(ˆt0,x(ˆt0;t0,xt0(a),u(,)))},0t0ˆt0<T. (4.5)

    where L(t,x(t;t0,xt0(a),u(,)),u(a,t))=τ1It0(a)+τ2u(a,t)+12τ3u2(a,t).

    Proof. Let us denote the right-hand side of (4.5) by ¯V(t0,xt0(a)). By (4.3), we have

    V(t0,xt0(a))J(t0,xt0(a),u(a,t))=ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt+J(ˆt0,x(a,ˆt0),u(,)),u(,)U.

    Thus, taking the infimum over u(,)U we get

    V(t0,xt0(a))¯V(t0,xt0(a)). (4.6)

    Conversely, for any ε>0, there exists a uε(,)U such that

    V(t0,xt0(a))+εJ(t0,xt0(a),uε(,))ˆt0t0A0(τ1Iε(a,t)+τ2uε(a,t)+12τ3u2ε(a,t))dadt+V(ˆt0,xε(a,ˆt0))¯V(t0,xt0(a)). (4.7)

    Combing (4.6) and (4.7), we obtain (4.5).

    The equation (4.5) is dynamic programming equation, which gives a relationship among globally and locally optimal via value function.

    We would like to further study the (4.5), trying to obtain an equation for V(t0,xt0(a)) simpler form. The following results give a partial differential equation that a continuously differentiable value function ought to satisfy. We know the V(t0,xt0(a)) is a continuously differentiable value function via 4.1 and basic theorem of calculus.

    Define the Lyapunov equation (LE)

    Vx,(f(a,t,x)+g(a,t,x)u(a,t))+A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da=0.

    That is

    0=VS,(Sa+λ(a)μ(a)Sβ(a)SI+ρ(a)η(a)I+γ(a)Ru(a,t)S)+VR,(Ra+(1ρ(a))η(a)I(μ(a)+γ(a))R+u(a,t)S)+VI,(Ia+β(a)SIδ(a)I)+A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da=A0VS(Sa+λ(a)μ(a)Sβ(a)SI+ρ(a)η(a)I+γ(a)Ru(a,t)S)da+A0VR(Ra+(1ρ(a))η(a)I(μ(a)+γ(a))R+u(a,t)S)da+A0VI(Ia+β(a)SIδ(a)I)da+A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da, (4.8)

    where Vx=(VS,VI,VR), Vt=Vt. The equation LE is the HJB equation when u(a,t) is optimal control of (4.1).

    {0=A0VS(Sa+λ(a)μ(a)Sβ(a)SI+ρ(a)η(a)I+γ(a)Ru(a,t)S)da+A0VR(Ra+(1ρ(a))η(a)I(μ(a)+γ(a))R+u(a,t)S)da+A0VI(Ia+β(a)SIδ(a)I)da+A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da,Vt=T=h(x),xH×H×H. (4.9)

    In this subpart, we write the HJB equation of model (4.1) as well as its viscosity solution is discussed as following.

    Theorem 4.2. Let (D1)-(D3) and V(t,x)C1([0,T]×H×H×H). Then V(t,x) is a solution of the following terminal value problem of HJB equation:

    {Vt+infuUH(t,x,u,Vx)=0,t[0,T),V|t=T=h(x),(t,x,u)[0,T]×H×H×H×U. (4.10)

    which can be written as (4.9).

    Proof. Fix (t0,xt0(a))[0,T)×H×H×H and uU. Let x(,) be the state trajectory corresponding to the control u(,)Uω[s,T] with u(a,t)u. On the one hand, we can write according to (4.5)

    0{V(ˆt0,x(ˆt0,a))V(t0,xt0(a))}ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt. (4.11)

    Rearranging the equation (4.11) and dividing the both side of the equation by ˆt0t0 yield

    0{V(ˆt0,x(a,ˆt0))V(t0,xt0(a))}ˆt0t01ˆt0t0ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt={V(ˆt0,x(a,ˆt0))V(ˆt0,xt0(a))+V(ˆt0,x(a,ˆt0))V(t0,x(a,ˆt0))}ˆt0t01ˆt0t0ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt. (4.12)

    By (4.5) with ˆt0t0, we obtain

    0Vt(t0,xt0(a))Vx(t0,xt0(a)),dxdtA0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da=A0{VS(t0,xt0(a))(S(a,t)a+λ(a)μ(a)S(a,t)β(a)S(a,t)I(a,t)+ρ(a)η(a)I(a,t)+γ(a)R(a,t)u(a,t)S(a,t))VI(t0,xt0(a))(I(a,t)a+β(a)S(a,t)I(a,t)δ(a)I(a,t))VR(t0,xt0(a))(R(a,t)a+(1ρ(a))η(a)I(a,t)(μ(a)+γ(a))R(a,t)+u(a,t)S(a,t))}daA0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da. (4.13)

    This results in

    0Vt(t,x(a,t))+infuUA0{VS(t,x(a,t))(S(a,t)aμ(a)S(a,t)β(a)S(a,t)I(a,t)+ρ(a)η(a)I(a,t)+λ(a)+γ(a)R(a,t))+VI(t,x(a,t))(I(a,t)a+β(a)S(a,t)I(a,t)δI(a,t))+VR(t,x(a,t))(R(a,t)a+(1ρ(a))η(a)I(a,t)(μ(a)+γ(a))R(a,t))}da+A0(τ1I(a,t)τS(VS(t,x(a,t))VR(t,x(a,t)))+12τ3(S(VS(t,x(a,t))VR(t,x(a,t))))2)da. (4.14)

    On the other hand, for any ε>0, 0t0^t0<T with ˆt0t0>0 small enough, there exists a u(,)uε,ˆt0(,)U such that

    V(t0,xt0(a))+ε(ˆt0t0)V(ˆt0,x(ˆt0;t0,xt0(a),u(a,t)))+ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt. (4.15)

    Thus, it follows that as (noting VC1([0,T]×H×H×H))

    ε{V(ˆt0,x(ˆt0))V(t0,x(t0))}ˆt0t01ˆt0t0ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt1ˆt0t0ˆt0t0[Vt(t0,xt0(a))Vx(t0,xt0(a)),dxdtA0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da]dt1ˆt0t0ˆt0t0[Vt(t0,xt0(a))Vx(t0,xt0(a)),dxdtinfuUA0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))da]dt1ˆt0t0ˆt0t0Vt(t0,xt0(a))+infuUH(t,x(a,t),u(a,t),Vx(t,x(a,t)))dt. (4.16)

    Taking limit on (4.16) with respect to ˆt0t0, we have

    0Vt(t,x(a,t))+infuUA0{VS(t,x(a,t))(S(a,t)aμ(a)S(a,t)β(a)S(a,t)I(a,t)+ρ(a)η(a)I(a,t)+λ(a)+γ(a)R(a,t))+VI(t,x(a,t))(I(a,t)a+β(a)S(a,t)I(a,t)δI(a,t))+VR(t,x(a,t))(R(a,t)a+(1ρ(a))η(a)I(a,t)(μ(a)+γ(a))R(a,t))}da+A0(τ1I(a,t)τS(VS(t,x(a,t))VR(t,x(a,t)))+12τ3(S(VS(t,x(a,t))VR(t,x(a,t))))2)da. (4.17)

    Combining (4.14) and (4.17), we obtain our results.

    Definition 4.3. [34] A function vC([0,T]×H×H×H) is called a viscositysolution of (4.9) if

    v(T,x(a,T))h(x(a,T)),x(a,T)H×H×H, (4.18)

    and for any ψC1([0,T]×R), whenever vψ attains a local maximum at (t,x)[0,T)×H×H×H, we have

    ψt(t,x)+supuUH(t,x,u,ψx(t,x))0. (4.19)

    The function vC([0,T]×H×H×H) is called a viscositysolution of (4.9) if in (4.18)(4.19) the inequalities are changed to and localmaximum is changed to localminimum. In the case that v is both a viscosity subsolution and supersolution of (4.9), it is called a viscosity solution of (4.9).

    Theorem 4.4. Let (D1)-(D3) hold. Then the value function of V(t,x) satisfies

    |V(t0,xt0(a))V(ˆt0,xˆt0(a))|K{|t0ˆt0|},(t0,xt0(a)),(ˆt0,xˆt0(a))[0,T]×H×H×H, (4.20)

    for some K>0. Moreover, V(t,x(a,t)) is the only solution of (4.9) in the class C([0,T]×H×H×H).

    Proof. Through 4.2 and t0ˆt0 we obtain

    |V(t0,xt0(a))V(ˆt0,xˆt0(a))|infuU{ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt}ˆt0t0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt|t0ˆt0|maxt0tˆt0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))daK|t0ˆt0|. (4.21)

    Thus, V(t,x(a,t)) satisfies (4.21). We can know that V(t,x(a,t)) is unique viscosity of the HJB equation via the same method as [[34], Theorem 2.5].

    In this subsection, we study the existence of the optimal control pair in finite time for the system (4.1) and construct the Hamiltonian of the optimal control problem.

    Theorem 4.5. There exists an optimal control uU and a corresponding optimal state S,I,R such that

    V(t0,xt0(a))=infu()UJ(t0,xt0(a),u(a,t)),(t0,xt0(a))[0,T)×H×H×H,

    subject to the control system (4.1).

    Proof. To show the existence of the optimal control for the problem under consideration, we use the result in [35,36]. We notice that the state and control variables are non-negative, and the control set U, by definition, is closed and bounded. The optimal system is bounded, which ensures the compactness needed for the existence of an optimal control. Further, the integrand Tt0A0(τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t))dadt is convex on the control set U due to the biquadratic and quadratic nature of control variable u(a,t), which represents vaccination proportion. In addition, there exists a constant ν>1 and positive numbers κ1, κ2 such that

    τ1I(a,t)+τ2u(a,t)+12τ3u2(a,t)κ1|u(a,t)|νκ2.

    We complete the existence of the optimal control of state variables.

    Theorem 4.6. Let u(a,t) be optimal control variable, S(a,t), I(a,t) and R(a,t) be corresponding optimal state variable of model (4.1). The corresponding optimal control is given as follows:

    u(a,t)=min{max[0,1τ3S(a,t)(VR(x,t)VS(x,t))+τ],1}

    where τ2τ3=τ.

    Proof. We define the Hamiltonian function H(t,x,u(a,t),Vx),t[0,T],u(a,t)Uad,x,Vx(x,t)H×H×H and by constructing the Hamiltonian function and HJB equation to characterize this optimal control u(a,t) as follows:

    H=(Sa+λ(a)μ(a)Sβ(a)SI+ρ(a)η(a)I+γ(a)Ru(a,t)S)VS(x,t)+(Ia+β(a)SIδ(a)I)VI(x,t)τ1Iτ2u(a,t)12τ3u2(a,t)+(Ra+(1ρ(a))η(a)I(μ(a)+γ(a))R+u(a,t)S)VR(x,t). (4.22)

    Let u(a,t) be given optimal control and S(a,t), I(a,t) and R(a,t) be corresponding optimal state variable of model (4.1). By the Hamiltonian function (4.22), optimal control u(a,t)=u(a,t) is obtained

    u(a,t)=1τ3S(a,t)(VR(x,t)VS(x,t))+τ. (4.23)

    From the properties of the control set with the findings, then there is exists optimal control

    u(a,t)={0,1τ3S(VRVS)+τ<0,1τ3S(VR(x,t)VS(x,t))+τ,01τ3S(VRVS)+τ<1,1,11τ3S(VRVS)+τ. (4.24)

    Remark 4.7. If we could obtain the value function V(t,x(a,t)) by solving the HJB equation, then we could, at least formally, construct an optimal pair for each of optimal control problem for different initial values. The principle involves the following steps to solve optimal control problem

    Step.1. Solve the HJB equation (4.10) to find the V(T,x(a,t));

    Step.2. Find u(a,t) through Hamiltonian function (4.22);

    Step.3. Combine the u(a,t) to solve model (4.1) to get the optimal pair (x(a,t),u(a,t)).

    This section aims to illustrate the effectiveness of our theoretical results obtained in previous sections.

    In the numerical simulation, the parameters are given in Table 2. We assume that (a,t)[0,80]×[0,500],Δt=0.1,Δa=1, take p=1.0,ϱ=0.2,δ=0.1827,k=0.2, then R02.69>1. The parameters that we partially need to assume are reasonable according to the study of Yang et al. [14] has shown that the range of the basic reproduction number of the pandemic influenza in 1.04.0. We see from 3.2 the E is globally asymptotically stable. Actually, as showing in the following Figure 2 (a), the density of the infected individual I(a,t) tends to be a positive constant over time.

    Table 2.  Parameter values used in numerical simulations.
    Parameter Value Source of data
    λ 0.20548 [16], Lan et al.
    β 1.6805×104 [16], Lan et al.
    μ 3.4246×105 [18], Mu & Zhang
    η 0.02 [18], Mu & Zhang
    ρ 0.05 [18], Mu & Zhang
    γ 0.027 [16], Lan et al.
    δ 0.2227 Assumed
    ϱ 0.2 Assumed
    u 0.35 Assumed
    A 80 years [16], Lan et al.
    T 500 days Assumed

     | Show Table
    DownLoad: CSV
    Figure 2.  The evolution of I(a,t) of (2.4) for R02.69>1 (a) and R00.78<1 (b).

    Take ϱ=0.1,δ=0.1827,k=0.1, then R00.78<1, thus we can know the E0 is globally asymptotically stable, as showing in the following Figure 2 (b), the density of the infected individual I(a,t) tends to zero over time. These are the same conclusions as 3.2.

    We choose δ=0.1827, 0.2027 and 0.2227 to study the effect on the infected population and all the other parameters are fixed as in Table 2. We calculate the corresponding basic reproduction number R00.9427, 0.8979 and 0.7802 when R0<1, respectively. In Figure 3 (a), we can see that the extinction time of decreases with the increase of infected removal rate δ. In order to simulate the cases of R0>1, we calculate the corresponding basic reproduction number R01.3468, 2.2441 and 2.6930, respectively. The simulation result is presented in Figure 3 (b). Synthetically, it is obvious to see from Figure 3 that the peak of I(a,t) decreases significantly as infected removal rate δ of age a goes up.

    Figure 3.  (a): the path of I(a,t) under different δ when R0<1 at age a=30; (b): the path of I(a,t) under different δ when R0>1 at age a=30.

    In this section, we show the result of optimal control using a numerical example. The parameters are given in Table 2. (a,t)[0,80]×[0,500],Δt=0.1,Δa=1 by performing some numerical simulations. We take τ1=11.7, τ2=3.4 and τ3=12. The corresponding paths of control and infective populations I(a,t) are plotted in Figure 4, Figure 5 (a), Figure 6 and Figure 7 (a) show the section view at age a=20 and a=30 of the image variation tendency of I(a,t) with and without control and the tendency of corresponding optimal control u(a,t). We observe that the comprehensive use of control u(a,t) works better than no control applied during the course of disease. The effect of applying control reduces the number of infective populations I(a,t) dramatically. In short, control variable u(a,t) play an important role in the control of the disease. At the same time, we can also observe a reduced epidemic time of population infection.

    Figure 4.  When R0>1, (a): the paths of I(a,t) with and without control u(a,t) at age a=20; (b): The paths of I(a,t) with and without control u(a,t) at age a=30.
    Figure 5.  (a): when R0>1, the paths of control at age a=20 and a=30; (b): the path of control at t=500.
    Figure 6.  When R0<1, (a): the paths of I(a,t) with and without control u(a,t) at age a=20; (b): The paths of I(a,t) with and without control u(a,t) at age a=30.
    Figure 7.  (a): when R0<1, the paths of control at age a=20 and a=30; (b): the paths of control at t=500.

    We note that a significant decline in the number of infected persons I(a,t) during vaccine control of susceptible at age a=20 and a=30 from these following figures.

    When R0>1, I(a,t) trends is regard to the optimal controls, where the upper-level image is the evolution of I(a,t) without control and the lower-level image is the evolution of I(a,t) with control on in Figure 4. In Figure 5 (a) shows the optimal control pathes at ages a=20 and a=30, repectively. Eventually the disease tends to extinction with control. Compared with the path of non-vaccinated I(a,t), the epidemic time of I(a,t) with vaccination is significantly shortened, as well as the density of infected population verge to zero eventually.

    When R0<1, I(a,t) trends is regard to the optimal controls, where the upper-level image is the evolution of I(a,t) without control and the lower-level image is the evolution of I(a,t) with control in Figure 6. In Figure 7 (a) shows the optimal control pathes at ages a=20 and a=30, repectively. Eventually the disease tends to extinction with control. Compared with Compared with the path of non-vaccinated I(a,t), the epidemic time of I(a,t) with vaccination is significantly shortened. Combined with Figure 4 and the corresponding conclusions, it is clear that vaccination have a certain degree of hindrance to the transmission of epidemics. We also observe the Figure 5 (b) and Figure 7 (b) that vaccine control varies for different ages, mainly adolescents and middle age, and this control strategy is consistent with the actual control of many epidemics such as pneumonia, influenza and so on. The pneumococcal vaccine, COVID-19 vaccine and influenza vaccine are not suitable for the elderly. Overall, population age has a significant effect on the implementation of vaccine control strategies for many common diseases.

    In this paper, we recurred an age-structured SIRS epidemic model (2.1) with vaccination control, and obtained the condition for disease extinction and persistence. The research results showed that when R0<0, the disease extincted; when R0>0, the disease persisted. The conditions of disease extinction and permanence were given in 3.2. And we verified the resluts by numerical simulation. However, since the bifurcation phenomenon will occur when R0=0, this part will also be an important point for our further study of the age-structured SIRS epidemic system dynamics. Besides, we also studied the control problem of SIRS epidemic model with age structure. By proving the existence of viscosity solution, we obtain the existence of optimal control. Then, we obtained the implicit expression for optimal control and HJB equation. Finally, by using Milstein method, the optimal proportion of vaccination is obtained. In addition, we plan that in the next step, traveling wave solutions of an age-structured SIRS epidemic model will become one of our main research contents. Wu et al. [37] studied the existence and non-existence of traveling wave solutions for a diffusive age-structured SIR epidemic model. Their work provides some insights on how to deal with the high dimensional age-structured epidemic models.

    The authors thank the editor and referees for their careful reading and valuable comments. The research was supported in part by the Ningxia Natural Science Foundation (No. 2020AAC03067).

    The authors have no conflict of interest.



    [1] M. Portegies, P. Koudstaal, M. Ikram, Cerebrovascular disease, Handb. Clin. Neurol., 138 (2016), 239–261. https://doi.org/10.1016/B978-0-12-802973-2.00014-8 doi: 10.1016/B978-0-12-802973-2.00014-8
    [2] X. Tong, Q. Yang, M. D. Ritchey, M. G. George, S. L. Jackson, C. Gillespie, et al., The burden of cerebrovascular disease in the United States, Prev. Chronic. Dis., 16 (2019), 180411. https://doi.org/10.5888/pcd16.180411 doi: 10.5888/pcd16.180411
    [3] S. Yoshimura, N. Sakai, H. Yamagami, K. Uchida, M. Beppu, K. Toyoda, et al., Endovascular therapy for acute stroke with a large ischemic region, N. Engl. J. Med., 386 (2022), 1303–1313. https://doi.org/10.1056/NEJMoa2118191 doi: 10.1056/NEJMoa2118191
    [4] Q. Z. Tuo, S. T. Zhang, P. Lei, Mechanisms of neuronal cell death in ischemic stroke and their therapeutic implications, Med. Res. Rev., 42 (2022), 259–305. https://doi.org/10.1002/med.21817 doi: 10.1002/med.21817
    [5] K. N. Kadyrovich, S. K. Erkinovich, K. M. Ilhomovna, Microscopic examination of postcapillary cerebral venues in hemorrhagic stroke, Am. J. Med. Sci. Pharm. Res., 3 (2021), 69–73. https://doi.org/10.37547/TAJMSPR/Volume03Issue08-11 doi: 10.37547/TAJMSPR/Volume03Issue08-11
    [6] M. L. Osgood, Aneurysmal subarachnoid hemorrhage: Review of the pathophysiology and management strategies, Curr. Neurol. Neurosci. Rep., 21 (2021), 1–11. https://doi.org/10.1007/s11910-021-01136-9 doi: 10.1007/s11910-021-01136-9
    [7] S. N. Neifert, E. K. Chapman, M. L. Martini, W. H. Shuman, A. J. Schupper, E. K. Oermann, et al., Aneurysmal subarachnoid hemorrhage: The last decade, Transl. Stroke Res., 12 (2021), 428–446. https://doi.org/10.1007/s12975-020-00867-0 doi: 10.1007/s12975-020-00867-0
    [8] W. S. Dodd, D. Laurent, A. S. Dumont, D. M. Hasan, P. M. Jabbour, R. M. Starke, et al., Pathophysiology of delayed cerebral ischemia after subarachnoid hemorrhage: A review, J. Am. Heart Assoc., 10 (2021), e021845. https://doi.org/10.1161/jaha.121.021845 doi: 10.1161/jaha.121.021845
    [9] I. Daghals, M. Sargurupremraj, R. Danning, P. Gormley, R. Malik, P. Amouyel, et al., Migraine, stroke, and cervical arterial dissection: Shared genetics for a triad of brain disorders with vascular involvement, Neurol. Genet., 8 (2022), e653. https://doi.org/10.1212/nxg.0000000000000653 doi: 10.1212/nxg.0000000000000653
    [10] M. Marciniec, K. Sapko, M. Kulczyński, S. Popek-Marciniec, A. Szczepańska-Szerej, K. Rejdak, Non-traumatic cervical artery dissection and ischemic stroke: A narrative review of recent research, Clin. Neurol. Neurosurg., 187 (2019), 105561. https://doi.org/10.1016/j.clineuro.2019.105561 doi: 10.1016/j.clineuro.2019.105561
    [11] M. M. Ruchoux, C. A. Maurage, CADASIL: Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy, J. Neuropathol. Exp. Neurol., 56 (1997), 947–964. https://doi.org/10.1097/00005072-199709000-00001 doi: 10.1097/00005072-199709000-00001
    [12] M. Bousser, E. Tournier-Lasserve, Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy: From stroke to vessel wall physiology, J. Neurol. Neurosurg. Psychiatry, 70 (2001), 285–287. https://doi.org/10.1136/jnnp.70.3.285 doi: 10.1136/jnnp.70.3.285
    [13] E. Navarro, F. Díaz, L. Muñoz, S. Giménez-Roldán, Dominant autosomal cerebral arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL): A review, Neurologia, 17 (2002), 410–417.
    [14] L. Prado, C. Han, S. P. Oh, H. Su, Recent advances in basic research for brain arteriovenous malformation, Int. J. Mol. Sci., 20 (2019), 5324. https://doi.org/10.3390/ijms20215324 doi: 10.3390/ijms20215324
    [15] P. Pan, S. Weinsheimer, D. Cooke, E. Winkler, A. Abla, H. Kim, et al., Review of treatment and therapeutic targets in brain arteriovenous malformation, J. Cereb. Blood Flow Metab., 41 (2021), 3141–3156. https://doi.org/10.1177/0271678x211026771 doi: 10.1177/0271678x211026771
    [16] S. S. Shaligram, E. Winkler, D. Cooke, H. Su, Risk factors for hemorrhage of brain arteriovenous malformation, CNS Neurosci. Ther., 25 (2019), 1085–1095. https://doi.org/10.1111/cns.13200 doi: 10.1111/cns.13200
    [17] M. Capecchi, M. Abbattista, I. Martinelli, Cerebral venous sinus thrombosis, J. Thromb. Haemostasis, 76 (2018), 12–15. https://doi.org/10.1111/jth.14210 doi: 10.1111/jth.14210
    [18] G. Sébire, B. Tabarki, D. E. Saunders, I. Leroy, R. Liesner, C. Saint-Martin, et al., Cerebral venous sinus thrombosis in children: Risk factors, presentation, diagnosis and outcome, Brain, 128 (2005), 477–489. https://doi.org/10.1093/brain/awh412 doi: 10.1093/brain/awh412
    [19] J. M. Ferro, P. Canhão, Cerebral venous sinus thrombosis: Update on diagnosis and management, Curr. Cardiol. Rep., 16 (2014), 523. https://doi.org/10.1007/s11886-014-0523-2 doi: 10.1007/s11886-014-0523-2
    [20] H. Zhang, L. Zheng, L. Feng, Epidemiology, diagnosis and treatment of moyamoya disease, Exp. Ther. Med., 17 (2019), 1977–1984. https://doi.org/10.3892/etm.2019.7198 doi: 10.3892/etm.2019.7198
    [21] J. Li, M. Jin, X. Sun, J. Li, Y. Liu, Y. Xi, et al., Imaging of moyamoya disease and moyamoya syndrome: Current status, J. Comput. Assist. Tomogr., 43 (2019), 257–263. https://doi.org/10.1097/rct.0000000000000834 doi: 10.1097/rct.0000000000000834
    [22] S. Shang, D. Zhou, J. Ya, S. Li, Q. Yang, Y. Ding, et al., Progress in moyamoya disease, Neurosurg. Rev., 43 (2020), 371–382. https://doi.org/10.1007/s10143-018-0994-5 doi: 10.1007/s10143-018-0994-5
    [23] C. Esenwa, J. Gutierrez, Secondary stroke prevention: Challenges and solutions, Vasc. Health Risk Manag., 11 (2015), 437–450. https://doi.org/10.2147/vhrm.S63791 doi: 10.2147/vhrm.S63791
    [24] V. L. Feigin, M. H. Forouzanfar, R. Krishnamurthi, G. A. Mensah, M. Connor, D. A. Bennett, et al., Global and regional burden of stroke during 1990–2010: Findings from the global burden of disease study 2010, Lancet, 383 (2014), 245–254. https://doi.org/10.1016/s0140-6736(13)61953-4 doi: 10.1016/s0140-6736(13)61953-4
    [25] I. Nagy, D. Fabó, Clinical neurophysiological methods in diagnosis and treatment of cerebrovascular diseases, Ideggyogy. Sz., 71 (2018), 7–14. https://doi.org/10.18071/isz.71.0007 doi: 10.18071/isz.71.0007
    [26] J. A. Barios, S. Ezquerro, A. Bertomeu-Motos, J. M. Catalan, J. M. Sanchez-Aparicio, L. Donis-Barber, et al., Movement-related EEG oscillations of contralesional hemisphere discloses compensation mechanisms of severely affected motor chronic stroke patients, Int. J. Neural. Syst., 31 (2021), 2150053. https://doi.org/10.1142/s0129065721500532 doi: 10.1142/s0129065721500532
    [27] S. Shaban, B. Huasen, A. Haridas, M. Killingsworth, J. Worthington, P. Jabbour, et al., Digital subtraction angiography in cerebrovascular disease: Current practice and perspectives on diagnosis, acute treatment and prognosis, Acta Neurol. Belg., 122 (2022), 763–780. https://doi.org/10.1007/s13760-021-01805-z doi: 10.1007/s13760-021-01805-z
    [28] C. C. Young, R. H. Bonow, G. Barros, M. Mossa-Basha, L. J. Kim, M. R. Levitt, Magnetic resonance vessel wall imaging in cerebrovascular diseases, Neurosurg. Focus, 47 (2019), E4. https://doi.org/10.3171/2019.9.Focus19599 doi: 10.3171/2019.9.Focus19599
    [29] S. F. Xiang, J. T. Li, S. J. Yang, F. F. Ding, W. W. Wang, S. Huo, et al., Whole-brain volume CT angiography can effectively detect early ischemic cerebrovascular diseases, Curr. Med. Imaging, 18 (2022), 731–738. https://doi.org/10.2174/1573405617666211206091831 doi: 10.2174/1573405617666211206091831
    [30] K. Taneja, H. Lu, B. G. Welch, B. P. Thomas, M. Pinho, D. Lin, et al., Evaluation of cerebrovascular reserve in patients with cerebrovascular diseases using resting-state MRI: A feasibility study, Magn. Reson. Imaging, 59 (2019), 46–52. https://doi.org/10.1016/j.mri.2019.03.003 doi: 10.1016/j.mri.2019.03.003
    [31] L. Carnevale, G. Lembo, Innovative MRI techniques in neuroimaging approaches for cerebrovascular diseases and vascular cognitive impairment, Int. J. Mol. Sci., 20 (2019), 2656. https://doi.org/10.3390/ijms20112656 doi: 10.3390/ijms20112656
    [32] J. S. Minhas, R. B. Panerai, G. Ghaly, P. Divall, T. G. Robinson, Cerebral autoregulation in hemorrhagic stroke: A systematic review and meta-analysis of transcranial Doppler ultrasonography studies, J. Clin. Ultrasound, 47 (2019), 14–21. https://doi.org/10.1002/jcu.22645 doi: 10.1002/jcu.22645
    [33] S. Sharma, R. J. Lubrica, M. Song, R. Vandse, W. Boling, P. Pillai, The role of transcranial Doppler in cerebral vasospasm: A literature review, Subarachnoid Hemorrh., 127 (2020), 201–205. https://doi.org/10.1007/978-3-030-04615-6_32 doi: 10.1007/978-3-030-04615-6_32
    [34] N. Fatima, A. Shuaib, T. S. Chughtai, A. Ayyad, M. Saqqur, The role of transcranial Doppler in traumatic brain injury: A Systemic review and Meta-analysis, Asian J. Neurosurg., 14 (2019), 626–633. https://doi.org/10.4103/ajns.AJNS_42_19 doi: 10.4103/ajns.AJNS_42_19
    [35] V. K. Sharma, K. S. Wong, A. V. Alexandrov, Transcranial Doppler, Front. Neurol. Neurosci., 40 (2016), 124–140. https://doi.org/10.1159/000448309 doi: 10.1159/000448309
    [36] C. C. Bishop, S. Powell, D. Rutt, N. L. Browse, Transcranial Doppler measurement of middle cerebral artery blood flow velocity: A validation study, Stroke, 17 (1986), 913–915. https://doi.org/10.1161/01.str.17.5.913 doi: 10.1161/01.str.17.5.913
    [37] S. Purkayastha, F. Sorond, Transcranial Doppler ultrasound: Technique and application, Semin. Neurol., 32 (2012), 411–420. https://doi.org/10.1055/s-0032-1331812 doi: 10.1055/s-0032-1331812
    [38] C. K. Willie, F. L. Colino, D. M. Bailey, Y. C. Tzeng, G. Binsted, L. W. Jones, et al., Utility of transcranial Doppler ultrasound for the integrative assessment of cerebrovascular function, J. Neurosci. Methods., 196 (2011), 221–237. https://doi.org/10.1016/j.jneumeth.2011.01.011 doi: 10.1016/j.jneumeth.2011.01.011
    [39] G. Tsivgoulis, A. V. Alexandrov, M. A. Sloan, Advances in transcranial Doppler ultrasonography, Curr. Neurol. Neurosci. Rep., 9 (2009), 46–54. https://doi.org/10.1007/s11910-009-0008-7 doi: 10.1007/s11910-009-0008-7
    [40] F. A. Rasulo, E. De Peri, A. Lavinio, Transcranial Doppler ultrasonography in intensive care, Eur. J. Anaesthesiol. Suppl., 42 (2008), 167–173. https://doi.org/10.1017/s0265021507003341 doi: 10.1017/s0265021507003341
    [41] T. Montrief, S. Alerhand, C. Jewell, J. Scott, Incorporation of transcranial Doppler into the ED for the neurocritical care patient, Am. J. Emerg. Med., 37 (2019), 1144–1152. https://doi.org/10.1016/j.ajem.2019.03.003 doi: 10.1016/j.ajem.2019.03.003
    [42] J. H. Fetzer, What is artificial intelligence, in Artificial Intelligence: Its Scope and Limits (eds J. H. Fetzer), (1990), 3–27. https://doi.org/10.1007/978-94-009-1900-6_1
    [43] M. Haenlein, A. Kaplan, A brief history of artificial intelligence: On the past, present, and future of artificial intelligence, Calif. Manage. Rev., 61 (2019), 5–14. https://doi.org/10.1177/0008125619864925 doi: 10.1177/0008125619864925
    [44] C. Zhang, Y. Lu, Study on artificial intelligence: The state of the art and future prospects, J. Ind. Inf. Integr., 23 (2021), 100224. https://doi.org/10.1016/j.jii.2021.100224 doi: 10.1016/j.jii.2021.100224
    [45] Y. K. Chan, Y. F. Chen, T. Pham, W. Chang, M. Y. Hsieh, Artificial intelligence in medical applications, J. Healthcare Eng., 2018 (2018), 1–2. https://doi.org/10.1155/2018/4827875 doi: 10.1155/2018/4827875
    [46] J. Stewart, P. Sprivulis, G. Dwivedi, Artificial intelligence and machine learning in emergency medicine, Emerg. Med. Australas., 30 (2018), 870–874. https://doi.org/10.1111/1742-6723.13145 doi: 10.1111/1742-6723.13145
    [47] A. Kirubarajan, A. Taher, S. Khan, S. Masood, Artificial intelligence in emergency medicine: A scoping review, J. Am. Coll. Emerg. Phys. Open., 1 (2020), 1691–1702. https://doi.org/10.1002/emp2.12277 doi: 10.1002/emp2.12277
    [48] Y. Berlyand, A. S. Raja, S. C. Dorner, A. M. Prabhakar, J. D. Sonis, R. V. Gottumukkala, et al., How artificial intelligence could transform emergency department operations, Am. Jo. Emerg. Med., 36 (2018), 1515–1517. https://doi.org/10.1016/j.ajem.2018.01.017 doi: 10.1016/j.ajem.2018.01.017
    [49] A. J. Myrden, A. Kushki, E. Sejdić, A. M. Guerguerian, T. Chau, A brain-computer interface based on bilateral transcranial Doppler ultrasound, PLoS One, 6 (2011), e24170. https://doi.org/10.1371/journal.pone.0024170 doi: 10.1371/journal.pone.0024170
    [50] N. Birbaumer, Breaking the silence: Brain-computer interfaces (BCI) for communication and motor control, Psychophysiology, 43 (2006), 517–532. https://doi.org/10.1111/j.1469-8986.2006.00456.x doi: 10.1111/j.1469-8986.2006.00456.x
    [51] H. Uğuz, A. Öztürk, R. Saraçoğlu, A. Arslan, A biomedical system based on fuzzy discrete hidden Markov model for the diagnosis of the brain diseases, Expert Syst. Appl., 35 (2008), 1104–1114. https://doi.org/10.1016/j.eswa.2007.08.006 doi: 10.1016/j.eswa.2007.08.006
    [52] I. Güler, F. Hardalaç, M. Kaymaz, Comparison of FFT and adaptive ARMA methods in transcranial Doppler signals recorded from the cerebral vessels, Comput. Biol. Med., 32 (2002), 445–453. https://doi.org/10.1016/s0010-4825(02)00036-7 doi: 10.1016/s0010-4825(02)00036-7
    [53] N. Karaboga, F. Latifoglu, Elimination of noise on transcranial Doppler signal using ⅡR filters designed with artificial bee colony—ABC-algorithm, Digital Signal Process., 23 (2013), 1051–1058. https://doi.org/10.1016/j.dsp.2012.09.015 doi: 10.1016/j.dsp.2012.09.015
    [54] T. Koza, S. Koçkana, N. Karaboğa, F. Latifoğlu, SSA analysis of noise eliminated transcranial doppler signals with iIR filters, in 2015 23nd Signal Processing and Communications Applications Conference (SIU), (2015), 2021–2024. https://doi.org/10.1109/SIU.2015.7130262
    [55] Y. J. Mei, R. T. Hu, J. Lin, H. Y. Xu, L. Y. Wu, H. P. Li, et al., Diagnosis of middle cerebral artery stenosis using transcranial Doppler images based on convolutional neural network, World Neurosurg., 161 (2022), e118–e125. https://doi.org/10.1016/j.wneu.2022.01.068 doi: 10.1016/j.wneu.2022.01.068
    [56] A. Baig, C. Manion, V. Iyer, W. Khawar, B. Donnelly, A. Monteiro, et al., E-142 robotic transcranial doppler with artificial intelligence to identify cerebral emboli during transcatheter aortic valve replacement - a novel neuromonitoring tool, J. NeuroInterventional Surg., 14 (2022), A152. https://doi.org/10.1136/neurintsurg-2022-SNIS.253 doi: 10.1136/neurintsurg-2022-SNIS.253
    [57] S. Esmaeeli, C. M. Hrdlicka, A. Brenes Bastos, J. Wang, S. Gomez-Paz, K. A. Hanafy, et al., Robotically assisted transcranial Doppler with artificial intelligence for assessment of cerebral vasospasm after subarachnoid hemorrhage, J. Neurocrit. Care, 13 (2020), 32–40. https://doi.org/10.18700/jnc.200002 doi: 10.18700/jnc.200002
    [58] H. White, B. Venkatesh, Applications of transcranial Doppler in the ICU: A review, Intensive Care Med., 32 (2006), 981–994. https://doi.org/10.1007/s00134-006-0173-y doi: 10.1007/s00134-006-0173-y
    [59] Y. Pan, W. Wan, M. Xiang, Y. Guan, Transcranial Doppler ultrasonography as a diagnostic tool for cerebrovascular disorders, Front. Hum. Neurosci., 16 (2022), 841809. https://doi.org/10.3389/fnhum.2022.841809 doi: 10.3389/fnhum.2022.841809
    [60] J. C. Granry, Transcranial Doppler in anesthesia and intensive care, Ann. Fr. Anesth. Reanim., 10 (1991), 127–136. https://doi.org/10.1016/s0750-7658(05)80453-8 doi: 10.1016/s0750-7658(05)80453-8
    [61] M. P. Spencer, D. Whisler, Transorbital Doppler diagnosis of intracranial arterial stenosis, Stroke, 17 (1986), 916–921. https://doi.org/10.1161/01.str.17.5.916 doi: 10.1161/01.str.17.5.916
    [62] L. Thomassen, U. Waje-Andreassen, H. Naess, J. Aarseth, D. Russell, Doppler ultrasound and clinical findings in patients with acute ischemic stroke treated with intravenous thrombolysis, Eur. J. Neurol., 12 (2005), 462–465. https://doi.org/10.1111/j.1468-1331.2005.01008.x doi: 10.1111/j.1468-1331.2005.01008.x
    [63] J. Allendoerfer, M. Goertler, G. M. von Reutern, Prognostic relevance of ultra-early doppler sonography in acute ischaemic stroke: A prospective multicentre study, Lancet Neurol., 5 (2006), 835–840. https://doi.org/10.1016/s1474-4422(06)70551-8 doi: 10.1016/s1474-4422(06)70551-8
    [64] A. Mattioni, S. Cenciarelli, P. Eusebi, M. Brazzelli, T. Mazzoli, M. Del Sette, et al., Transcranial Doppler sonography for detecting stenosis or occlusion of intracranial arteries in people with acute ischaemic stroke, Cochrane Database Syst. Rev., 2 (2020), Cd010722. https://doi.org/10.1002/14651858.CD010722.pub2 doi: 10.1002/14651858.CD010722.pub2
    [65] N. Samagh, H. Bhagat, K. Jangra, Monitoring cerebral vasospasm: How much can we rely on transcranial Doppler, J. Anaesthesiol. Clin. Pharmacol., 35 (2019), 12–18. https://doi.org/10.4103/joacp.JOACP_192_17 doi: 10.4103/joacp.JOACP_192_17
    [66] R. Aaslid, Transcranial Doppler assessment of cerebral vasospasm, Eur. J. Ultrasound, 16 (2002), 3–10. https://doi.org/10.1016/s0929-8266(02)00045-9 doi: 10.1016/s0929-8266(02)00045-9
    [67] D. D. Li, J. Y. Chang, C. X. Zhou, J. B. Cui, Clinical diagnosis of cerebral vasospasm after subarachnoid hemorrhage by using transcranial Doppler sonography, Eur. Rev. Med. Pharmacol. Sci., 22 (2018), 2029–2035. https://doi.org/10.26355/eurrev_201804_14732 doi: 10.26355/eurrev_201804_14732
    [68] C. E. Scherle Matamoros, E. A. Samaniego, K. Sam, J. A. Roa, J. P. Nellar, D. R. Rodríguez, Prediction of symptomatic vasospasm in patients with aneurysmal subarachnoid hemorrhage using early transcranial Doppler, J. Vasc. Interv. Neurol., 11 (2020), 19–26.
    [69] R. Behrouz, The rise and fall of transcranial Doppler ultrasonography for the diagnosis of vasospasm in aneurysmal subarachnoid hemorrhage, J. Neurosurg. Anesthesiol., 31 (2019), 79–80. https://doi.org/10.1097/ana.0000000000000490 doi: 10.1097/ana.0000000000000490
    [70] B. P. D. Inusa, L. Sainati, C. MacMahon, R. Colombatti, M. Casale, S. Perrotta, et al., An educational study promoting the delivery of transcranial Doppler ultrasound screening in paediatric sickle cell disease: A european multi-centre perspective, J. Clin. Med., 9 (2019), 44. https://doi.org/10.3390/jcm9010044 doi: 10.3390/jcm9010044
    [71] J. Kanter, S. Phillips, A. M. Schlenz, M. Mueller, M. Dooley, L. Sirline, et al., Transcranial doppler screening in a current cohort of children with sickle cell anemia: Results from the displace study, J. Pediatr. Hematol. Oncol., 43 (2021), e1062–e1068. https://doi.org/10.1097/mph.0000000000002103 doi: 10.1097/mph.0000000000002103
    [72] S. Thurn, K. Kleinschmidt, I. Kovacic, C. Wendl, R. A. Linker, S. Corbacioglu, et al., Transcranial Doppler sonography and the effect of haematopoietic stem cell transplantation in sickle cell disease, Neurol. Res. Pract., 4 (2022), 12. https://doi.org/10.1186/s42466-022-00175-y doi: 10.1186/s42466-022-00175-y
    [73] A. Pedicelli, M. Bartocci, E. Lozupone, F. D'Argento, A. Alexandre, G. Garignano, et al., The role of cervical color Doppler ultrasound in the diagnosis of brain death, Neuroradiology, 61 (2019), 137–145. https://doi.org/10.1007/s00234-018-2111-4 doi: 10.1007/s00234-018-2111-4
    [74] Y. Zurynski, N. Dorsch, I. Pearson, R. Choong, Transcranial Doppler ultrasound in brain death: Experience in 140 patients, Neurol. Res., 13 (1991), 248–252. https://doi.org/10.1080/01616412.1991.11740000 doi: 10.1080/01616412.1991.11740000
    [75] D. Escudero, J. Otero, B. Quindós, L. Viña, Transcranial Doppler ultrasound in the diagnosis of brain death. Is it useful or does it delay the diagnosis, Med. Intensiva, 39 (2015), 244–250. https://doi.org/10.1016/j.medin.2014.11.005 doi: 10.1016/j.medin.2014.11.005
    [76] D. Grosset, Aneurysm re-rupture: Doppler evidence of first phase vasospasm, J. Neurol. Neurosurg. Psychiatry, 53 (1990), 361. https://doi.org/10.1136/jnnp.53.4.361-b doi: 10.1136/jnnp.53.4.361-b
    [77] J. De Léan, N. Thomas, Neurological Complications of Pregnancy Advances in Neurology, VOLUME 64, 1993, Edited by O. Devinsky, E. Feldman and B. Hainline, Published by Raven Press, 286 pages, $C117.00, Can. J. Neurol. Sci., 24 (1977), 88. https://doi.org/10.1017/S031716710002120X
    [78] R. A. Purdy, Neurology in Clinical Practice, VOLUME 1 & 2, 1995, Edited by W. G. Bradley, R. B. Daroff, G. M. Fenichel, C. D. Marsden, Published by Butterworth-Heinemann, 3264 pages, $C439.00, Can. J. Neurol. Sci., 24 (1977), 88–89. https://doi.org/10.1017/S0317167100021223
    [79] R. Aaslid, T. M. Markwalder, H. Nornes, Noninvasive transcranial Doppler ultrasound recording of flow velocity in basal cerebral arteries, J. Neurosurg., 57 (1982), 769–774. https://doi.org/10.3171/jns.1982.57.6.0769 doi: 10.3171/jns.1982.57.6.0769
    [80] A. H. Ropper, S. M. Kehne, L. Wechsler, Transcranial Doppler in brain death, Neurology, 37 (1987), 1733–1733. https://doi.org/10.1212/wnl.37.11.1733 doi: 10.1212/wnl.37.11.1733
    [81] N. K. Altinbas, E. Ustuner, H. Ozcan, S. Bilgic, T. Sancak, E. Dusunceli, Effect of carotid artery stenting on ophthalmic artery flow patterns, J. Ultrasound Med., 33 (2014), 629–638. https://doi.org/10.7863/ultra.33.4.629 doi: 10.7863/ultra.33.4.629
    [82] J. Sun, Y. Liu, J. Zhang, X. Chen, Z. Lin, S. Nie, et al., Electroacupuncture improves cerebral vasospasm and functional outcome of patients with aneurysmal subarachnoid hemorrhage, Front. Neurosci., 12 (2018), 724. https://doi.org/10.3389/fnins.2018.00724 doi: 10.3389/fnins.2018.00724
    [83] D. B. Berry, E. K. Englund, S. Chen, L. R. Frank, S. R. Ward, Medical imaging of tissue engineering and regenerative medicine constructs, Biomater. Sci., 9 (2021), 301–314. https://doi.org/10.1039/d0bm00705f doi: 10.1039/d0bm00705f
    [84] O. M. Pinillos, C. N. Rodríguez, R. Hakimi, Transcranial Doppler ultrasound pulsatility index: Utility and clinical interpretation, in Neurosonology in Critical Care: Monitoring the Neurological Impact of the Critical Pathology (eds C. N. Rodríguez and C. Baracchini), (2022), 357–376. https://doi.org/10.1007/978-3-030-81419-9_21
    [85] N. G. Rainov, J. B. Weise, W. Burkert, Transcranial Doppler sonography in adult hydrocephalic patients, Neurosurg. Rev., 23 (2000), 34–38. https://doi.org/10.1007/s101430050029 doi: 10.1007/s101430050029
    [86] L. Wang, Y. Xing, Y. Li, K. Han, J. Chen, Evaluation of flow velocity in unilateral middle cerebral artery stenosis by Transcranial Doppler, Cell Biochem. Biophys., 70 (2014), 823–830. https://doi.org/10.1007/s12013-014-9986-4 doi: 10.1007/s12013-014-9986-4
    [87] M. B. McCarville, Comparison of duplex and nonduplex transcranial Doppler ultrasonography, Ultrasound Q, 24 (2008), 167–171. https://doi.org/10.1097/RUQ.0b013e3181862b67 doi: 10.1097/RUQ.0b013e3181862b67
    [88] L. Duan, L. D. Xu, Business Intelligence for Enterprise Systems: A Survey, IEEE Trans. Ind. Inf., 8 (2012), 679–687. https://doi.org/10.1109/tii.2012.2188804 doi: 10.1109/tii.2012.2188804
    [89] Y. Lu, L. D. Xu, Internet of things (IoT) cybersecurity research: A review of current research topics, IEEE Int. Things J., 6 (2019), 2103–2115. https://doi.org/10.1109/jiot.2018.2869847 doi: 10.1109/jiot.2018.2869847
    [90] P. Hamet, J. Tremblay, Artificial intelligence in medicine, Metabolism, 69 (2017), S36–S40. https://doi.org/10.1016/j.metabol.2017.01.011 doi: 10.1016/j.metabol.2017.01.011
    [91] V. Kaul, S. Enslin, S. A. Gross, History of artificial intelligence in medicine, Gastrointest. Endosc., 92 (2020), 807–812. https://doi.org/10.1016/j.gie.2020.06.040 doi: 10.1016/j.gie.2020.06.040
    [92] E. H. Shortliffe, R. Davis, S. G. Axline, B. G. Buchanan, C. C. Green, S. N. Cohen, Computer-based consultations in clinical therapeutics: Explanation and rule acquisition capabilities of the MYCIN system, Comput. Biomed. Res., 8 (1975), 303–320. https://doi.org/10.1016/0010-4809(75)90009-9 doi: 10.1016/0010-4809(75)90009-9
    [93] B. G. Buchanan, E. A. Feigenbaum, The stanford heuristic programming project: Goals and activities, AI Mag., 1 (2017), 25. https://doi.org/10.1609/aimag.v1i1.89 doi: 10.1609/aimag.v1i1.89
    [94] S. M. Weiss, C. A. Kulikowski, S. Amarel, A. Safir, A model-based method for computer-aided medical decision-making, Artif. Intell., 11 (1978), 145–172. https://doi.org/10.1016/0004-3702(78)90015-2 doi: 10.1016/0004-3702(78)90015-2
    [95] E. Strickland, IBM Watson, heal thyself: How IBM overpromised and underdelivered on AI health care, IEEE Spectr., 56 (2019), 24–31. https://doi.org/10.1109/mspec.2019.8678513 doi: 10.1109/mspec.2019.8678513
    [96] M. Christ, F. Grossmann, D. Winter, R. Bingisser, E. Platz, Modern triage in the emergency department, Dtsch. Arztebl. Int., 107 (2010), 892–898. https://doi.org/10.3238/arztebl.2010.0892 doi: 10.3238/arztebl.2010.0892
    [97] C. I. Price, L. Shaw, S. Islam, M. Javanbakht, A. Watkins, P. McMeekin, et al., Effect of an enhanced paramedic acute stroke treatment assessment on thrombolysis delivery during emergency stroke care: A cluster randomized clinical trial, JAMA Neurol., 77 (2020), 840–848. https://doi.org/10.1001/jamaneurol.2020.0611 doi: 10.1001/jamaneurol.2020.0611
    [98] M. Kanglie, S. Bipat, I. A. H. van den Berk, T. S. R. van Engelen, M. G. W. Dijkgraaf, J. M. Prins, et al., OPTimal IMAging strategy in patients suspected of non-traumatic pulmonary disease at the emergency department: Chest X-ray or ultra-low-dose chest CT (OPTIMACT) trial-statistical analysis plan, Trials, 21 (2020), 407. https://doi.org/10.1186/s13063-020-04343-w doi: 10.1186/s13063-020-04343-w
    [99] G. Vingerhoets, N. Stroobant, Lateralization of cerebral blood flow velocity changes during cognitive tasks. A simultaneous bilateral transcranial Doppler study, Stroke, 30 (1999), 2152–2158. https://doi.org/10.1161/01.str.30.10.2152 doi: 10.1161/01.str.30.10.2152
    [100] T. Wessels, J. U. Harrer, C. Jacke, U. Janssens, C. Klö tzsch, The prognostic value of early transcranial Doppler ultrasound following cardiopulmonary resuscitation, Ultrasound Med. Biol., 32 (2006), 1845–1851. https://doi.org/10.1016/j.ultrasmedbio.2006.06.023 doi: 10.1016/j.ultrasmedbio.2006.06.023
    [101] A. Myrden, A. Kushki, E. Sejdić, T. Chau, Towards increased data transmission rate for a three-class metabolic brain-computer interface based on transcranial Doppler ultrasound, Neurosci. Lett., 528 (2012), 99–103. https://doi.org/10.1016/j.neulet.2012.09.030 doi: 10.1016/j.neulet.2012.09.030
    [102] A. Goyal, A. A. Samadani, A. M. Guerguerian, T. Chau, An online three-class Transcranial Doppler ultrasound brain computer interface, Neurosci. Res., 107 (2016), 47–56. https://doi.org/10.1016/j.neures.2015.12.013 doi: 10.1016/j.neures.2015.12.013
    [103] A. Khalaf, M. Sybeldon, E. Sejdic, M. Akcakaya, A brain-computer interface based on functional transcranial doppler ultrasound using wavelet transform and support vector machines, J. Neurosci. Methods, 293 (2018), 174–182. https://doi.org/10.1016/j.jneumeth.2017.10.003 doi: 10.1016/j.jneumeth.2017.10.003
    [104] S. Serhatlioğlu, F. Hardalaç, I. Güler, Classification of transcranial Doppler signals using artificial neural network, J. Med. Syst., 27 (2003), 205–214. https://doi.org/10.1023/a:1021821229512 doi: 10.1023/a:1021821229512
    [105] A. Ozturk, A. Arslan, F. Hardalac, Comparison of neuro-fuzzy systems for classification of transcranial Doppler signals with their chaotic invariant measures, Expert Syst. Appl., 34 (2008), 1044–1055. https://doi.org/10.1016/j.eswa.2006.12.006 doi: 10.1016/j.eswa.2006.12.006
    [106] H. Uğuz, A hybrid system based on information gain and principal component analysis for the classification of transcranial Doppler signals, Comput. Methods Prog. Biomed., 107 (2011), 598–609. https://doi.org/10.1016/j.cmpb.2011.03.013 doi: 10.1016/j.cmpb.2011.03.013
    [107] H. Uğuz, A. Arslan, A new approach based on discrete hidden Markov model using Rocchio algorithm for the diagnosis of the brain diseases, Digital Signal Process., 20 (2010), 923–934. https://doi.org/10.1016/j.dsp.2009.11.001 doi: 10.1016/j.dsp.2009.11.001
    [108] M. Seera, C. P. Lim, K. S. Tan, W. S. Liew, Classification of transcranial Doppler signals using individual and ensemble recurrent neural networks, Neurocomputing, 249 (2017), 337–344. https://doi.org/10.1016/j.neucom.2016.05.11 doi: 10.1016/j.neucom.2016.05.11
    [109] N. Karaboga, F. Latifoglu, Adaptive filtering noisy transcranial Doppler signal by using artificial bee colony algorithm, Eng. Appl. Artif. Intell., 26 (2013), 677–684. https://doi.org/10.1016/j.engappai.2012.10.011 doi: 10.1016/j.engappai.2012.10.011
    [110] N. Karaboga, F. Latifoglu, T. Koza, SSA analysis of transcranial Doppler signal using ⅡR filters designed with ABC algorithm, Curr. Opin. Biotechnol., 22 (2011), 5. https://doi.org/10.1016/j.copbio.2011.05.159 doi: 10.1016/j.copbio.2011.05.159
    [111] B. Kamişlıoğlu, T. Koza, S. Koçkanat, N. Karaboğa, Noise cancellation on mitral valve Doppler signal with ⅡR digital filter using harmony search algorithm, in 2013 IEEE INISTA, (2013), 1–5. https://doi.org/10.1109/INISTA.2013.6577620
    [112] K. Clare, A. Stein, N. Damodara, E. Feldstein, H. Alshammari, S. Ali, et al., Safety and efficacy of a novel robotic transcranial doppler system in subarachnoid hemorrhage, Sci. Rep., 12 (2022), 1–6. https://doi.org/10.1038/s41598-021-04751-1 doi: 10.1038/s41598-021-04751-1
    [113] T. Wang, Y. Chen, H. Du, Y. Liu, L. Zhang, M. Meng, Monitoring of neuroendocrine changes in acute stage of severe craniocerebral injury by transcranial Doppler ultrasound image features based on artificial intelligence algorithm, Comput. Math. Methods Med., 2021 (2021), 3584034. https://doi.org/10.1155/2021/3584034 doi: 10.1155/2021/3584034
    [114] M. L. P. Portegies, P. J. Koudstaal, M. A. Ikram, Cerebrovascular disease, in Handbook of Clinical Neurology (eds M. J. Aminoff and F. Boller), (2016), 239–261. https://doi.org/10.1016/B978-0-12-802973-2.00014-8
    [115] R. Capildeo, S. Haberman, F. C. Rose, The definition and classification of stroke: A new approach, QJM Int. J. Med., 47 (1978), 177–196. https://doi.org/10.1093/oxfordjournals.qjmed.a067535 doi: 10.1093/oxfordjournals.qjmed.a067535
  • 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(3930) PDF downloads(300) Cited by(7)

Figures and Tables

Figures(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog