Loading [MathJax]/jax/element/mml/optable/Latin1Supplement.js
Research article

Ergodic stationary distribution of stochastic virus mutation model with time delay

  • The virus mutation can increase the complexity of the infectious disease. In this paper, the dynamical characteristics of the virus mutation model are discussed. First, we built a stochastic virus mutation model with time delay. Second, the existence and uniqueness of global positive solutions for the proposed model is proved. Third, based on the analysis of the ergodic stationary distribution for the model, we discuss the influence mechanism between the different factors. Finally, the numerical simulation verifies the theoretical results.

    Citation: Juan Ma, Shaojuan Ma, Xinyu Bai, Jinhua Ran. Ergodic stationary distribution of stochastic virus mutation model with time delay[J]. AIMS Mathematics, 2023, 8(9): 21371-21392. doi: 10.3934/math.20231089

    Related Papers:

    [1] Xiaoyong Chen, Yating Li, Liang Liu, Yaqiang Wang . Infinity norm upper bounds for the inverse of SDD1 matrices. AIMS Mathematics, 2022, 7(5): 8847-8860. doi: 10.3934/math.2022493
    [2] Lanlan Liu, Yuxue Zhu, Feng Wang, Yuanjie Geng . Infinity norm bounds for the inverse of SDD+1 matrices with applications. AIMS Mathematics, 2024, 9(8): 21294-21320. doi: 10.3934/math.20241034
    [3] Dizhen Ao, Yan Liu, Feng Wang, Lanlan Liu . Schur complement-based infinity norm bounds for the inverse of S-Sparse Ostrowski Brauer matrices. AIMS Mathematics, 2023, 8(11): 25815-25844. doi: 10.3934/math.20231317
    [4] Yingxia Zhao, Lanlan Liu, Feng Wang . Error bounds for linear complementarity problems of SDD1 matrices and SDD1-B matrices. AIMS Mathematics, 2022, 7(7): 11862-11878. doi: 10.3934/math.2022662
    [5] Xinnian Song, Lei Gao . CKV-type B-matrices and error bounds for linear complementarity problems. AIMS Mathematics, 2021, 6(10): 10846-10860. doi: 10.3934/math.2021630
    [6] Fatih Yılmaz, Aybüke Ertaş, Samet Arpacı . Some results on circulant matrices involving Fibonacci polynomials. AIMS Mathematics, 2025, 10(4): 9256-9273. doi: 10.3934/math.2025425
    [7] Man Chen, Huaifeng Chen . On ideal matrices whose entries are the generalized kHoradam numbers. AIMS Mathematics, 2025, 10(2): 1981-1997. doi: 10.3934/math.2025093
    [8] Yuanjie Geng, Deshu Sun . Error bounds for linear complementarity problems of strong SDD1 matrices and strong SDD1-B matrices. AIMS Mathematics, 2023, 8(11): 27052-27064. doi: 10.3934/math.20231384
    [9] Baijuan Shi . A particular matrix with exponential form, its inversion and some norms. AIMS Mathematics, 2022, 7(5): 8224-8234. doi: 10.3934/math.2022458
    [10] Lanlan Liu, Pan Han, Feng Wang . New error bound for linear complementarity problem of S-SDDS-B matrices. AIMS Mathematics, 2022, 7(2): 3239-3249. doi: 10.3934/math.2022179
  • The virus mutation can increase the complexity of the infectious disease. In this paper, the dynamical characteristics of the virus mutation model are discussed. First, we built a stochastic virus mutation model with time delay. Second, the existence and uniqueness of global positive solutions for the proposed model is proved. Third, based on the analysis of the ergodic stationary distribution for the model, we discuss the influence mechanism between the different factors. Finally, the numerical simulation verifies the theoretical results.



    Let n be an integer number, N={1,2,,n}, and let Cn×n be the set of all complex matrices of order n. A matrix M=(mij)Cn×n(n2) is called a strictly diagonally dominant (SDD) matrix [1] if

    |mii|>ri(M)=nj=1,ji|mij|,      iN.

    It was shown that an SDD matrix is an H-matrix [1], where matrix M=(mij)Cn×n(n2) is called an H-matrix [1, 2, 3] if and only if there exists a positive diagonal matrix X such that MX is an SDD matrix [1, 2, 4]. H-matrices are widely applied in many fields, such as computational mathematics, economics, mathematical physics and dynamical system theory, see [1, 4, 5, 6]. Meanwhile, upper bounds for the infinity norm of the inverse matrices of H-matrices can be used in the convergence analysis of matrix splitting and matrix multi-splitting iterative methods for solving the large sparse of linear equations [7], as well as linear complementarity problems. Moreover, upper bounds of the infinity norm of the inverse for different classes of matrices have been widely studied, such as CKV-type matrices [8], S-SDDS matrices [9], DZ and DZ-type matrices [10, 11], Nekrasov matrices [12, 13, 14, 15], S-Nekrasov matrices [16], Q-Nekrasov matrices [17], GSDD1 matrices [18] and so on.

    In 2011, Peňa [19] proposed a new subclass of H-matrices called SDD1 matrices, whose definition is listed below. A matrix M=(mij)Cn×n(n2) is called an SDD1 matrix if

    |mii|>pi(M),iN1(M),

    where

    pi(M)=jN1(M){i}|mij|+jN2(M){i}rj(M)|mjj||mij|,N1(M)={i||mii|ri(M)},    N2(M)={i||mii|>ri(M)}.

    In 2023, Dai et al. [18] gave a new subclass of H-matrices named generalized SDD1 (GSDD1) matrices, which extends the class of SDD1 matrices. Here, a matrix M=(mij)Cn×n(n2) is said a GSDD1 matrix if

    ri(M)pN2(M)i(M)>0,iN2(M),

    and

    (ri(M)pN2(M)i(M))(|ajj|pN1(M)j(M))>pN1(M)i(M)pN2(M)j(M),iN2(M),jN1(M),

    where

    pN2(M)i(M)=jN2(M){i}rj(M)|mjj||mij|,pN1(M)i(M)=jN1(M){i}|mij|,iN.

    Subsequently, some upper bounds for the infinite norm of the inverse matrices of SDD matrices, SDD1 matrices and GSDD1 matrices are presented, see [18, 20, 21]. For example, the following results that will be used later are listed.

    Theorem 1. (Varah bound) [21] Let matrix M=(mij)Cn×n(n2) be an SDD matrix. Then

    ||M1||1min1in(|mii|ri(M)).

    Theorem 2. [20] Let matrix M=(mij)Cn×n(n2) be an SDD matrix. Then

    ||M1||maxiNpi(M)|mii|+εminiNZi,0<ε<miniN|mii|pi(M)ri(M),

    where

    Zi=ε(|mii|ri(M))+jN{i}(rj(M)pj(M)|mjj|)|mij|.

    Theorem 3. [20] Let matrix M=(mij)Cn×n(n2) be an SDD matrix. If ri(M)>0(iN), then

    ||M1||maxiNpi(M)|mii|miniNjN{i}rj(M)pj(M)|mjj||mij|.

    Theorem 4. [18] Let M=(mij)Cn×n be a GSDD1 matrix. Then

    ||M1||max{ε,maxiN2(M)ri(M)|mii|}min{miniN2(M)ϕi,miniN1(M)ψi},

    where

    ϕi=ri(M)jN2(M){i}rj(M)|mjj||mij|jN1(M){i}|mij|ε,ψi=|mii|εjN1(M){i}|mij|ε+jN2(M){i}rj(M)|mjj||mij|,

    and

    maxiN1(M)pN2(M)i(M)|mii|pN1(M)i(M)<ε<minjN2(M)rj(M)pN2(M)j(M)pN1(M)j(M).

    On the basis of the above articles, we continue to study special structured matrices and introduce a new subclass of H-matrices called SDDk matrices, and provide some new upper bounds for the infinite norm of the inverse matrices for SDD matrices and SDDk matrices, which improve the previous results. The remainder of this paper is organized as follows: In Section 2, we propose a new subclass of H-matrices called SDDk matrices, which include SDD matrices and SDD1 matrices, and derive some properties of SDDk matrices. In Section 3, we present some upper bounds for the infinity norm of the inverse matrices for SDD matrices and SDDk matrices, and provide some comparisons with the well-known Varah bound. Moreover, some numerical examples are given to illustrate the corresponding results.

    In this section, we propose a new subclass of H-matrices called SDDk matrices, which include SDD matrices and SDD1 matrices, and derive some new properties.

    Definition 1. A matrix M=(mij)Cn×n(n2) is called an SDDk matrix if there exists kN such that

    |mii|>p(k1)i(M),iN1(M),

    where

    p(k)i(M)=jN1(M){i}|mij|+jN2(M){i}p(k1)j(M)|mjj||mij|,p(0)i(M)=jN1(M){i}|mij|+jN2(M){i}rj(M)|mjj||mij|.

    Immediately, we know that SDDk matrices contain SDD matrices and SDD1 matrices, so

    {SDD}{SDD1}{SDD2}{SDDk}.

    Lemma 1. A matrix M=(mij)Cn×n(n2) is an SDDk(k2) matrix if and only if for iN, |mii|>p(k1)i(M).

    Proof. For iN1(M), from Definition 1, it holds that |mii|>p(k1)i(M).

    For iN2(M), we have that |mii|>ri(M). When k=2, it follows that

    |mii|>ri(M)jN1(M){i}|mij|+jN2(M){i}rj(M)|mjj||mij|=p(0)i(M)jN1(M){i}|mij|+jN2(M){i}p(0)j(M)|mjj||mij|=p(1)i(M).

    Suppose that |mii|>p(k1)i(M)(kl,l>2) holds for iN2(M). When k=l+1, we have

    |mii|>ri(M)jN1(M){i}|mij|+jN2(M){i}p(l2)j(M)|mjj||mij|=p(l1)i(M)jN1(M){i}|mij|+jN2(M){i}p(l1)j(M)|mjj||mij|=p(l)i(M).

    By induction, we obtain that for iN2(M), |mii|>p(k1)i(M). Consequently, it holds that matrix M is an SDDk matrix if and only if |mii|>p(k1)i(M) for iN. The proof is completed.

    Lemma 2. If M=(mij)Cn×n(n2) is an SDDk(k2) matrix, then M is an H-matrix.

    Proof. Let X be the diagonal matrix diag{x1,x2,,xn}, where

    (0<) xj={1,jN1(M),p(k1)j(M)|mjj|+ε,jN2(M),

    and

    0<ε<miniN|mii|p(k1)i(M)jN2(M){i}|mij|.

    If jN2(M){i}|mij|=0, then the corresponding fraction is defined to be . Next we consider the following two cases.

    Case 1: For each iN1(M), it is not difficult to see that |(MX)ii|=|mii|, and

    ri(MX)=j=1,ji|mij|xj=jN1(M){i}|mij|+jN2(M){i}(p(k1)j(M)|mjj|+ε)|mij|jN1(M){i}|mij|+jN2(M){i}(p(k2)j(M)|mjj|+ε)|mij|=jN1(M){i}|mij|+jN2(M){i}p(k2)j(M)|mjj||mij|+jN2(M){i}ε|mij|=p(k1)i(M)+εjN2(M){i}|mij|<p(k1)i(M)+|mii|p(k1)i(M)=|mii|=|(MX)ii|.

    Case 2: For each iN2(M), we can obtain that

    |(MX)ii|=|mii|(pk1i(M)|mii|+ε)=p(k1)i(M)+ε|mii|,

    and

    ri(MX)=j=1,ji|mij|xj=jN1(M){i}|mij|+jN2(M){i}(p(k1)j(M)|mjj|+ε)|mij|jN1(M){i}|mij|+jN2(M){i}(p(k2)j(M)|mjj|+ε)|mij|=jN1(M){i}|mij|+jN2(M){i}p(k2)j(M)|mjj||mij|+jN2(M){i}ε|mij|=p(k1)i(M)+εjN2(A){i}|mij|<p(k1)i(M)+ε|mii|=|(MX)ii|.

    Based on Cases 1 and 2, we have that MX is an SDD matrix, and consequently, M is an H-matrix. The proof is completed.

    According to the definition of SDDk matrix and Lemma 1, we obtain some properties of SDDk matrices as follows:

    Theorem 5. If M=(mij)Cn×n(n2) is an SDDk matrix and N1(M), then for iN1(M), ji,jN2(M)|mij|>0.

    Proof. Suppose that for iN1(M), ji,jN2(M)|mij|=0. By Definition 1, we have that p(k1)i(M)=ri(M), iN1(M). Thus, it is easy to verify that |mii|>p(k1)i(M)=ri(M)|mii|, which is a contradiction. Thus for iN1(M), ji,jN2(M)|mij|>0. The proof is completed.

    Theorem 6. Let M=(mij)Cn×n(n2) be an SDDk(k2) matrix. It holds that ji,jN2(M)|mij|>0, iN2(M). Then

    |mii|>p(k2)i(M)>p(k1)i(M)>0,iN2(M),

    and

    |mii|>p(k1)i(M)>0,iN.

    Proof. By Lemma 1 and the known conditions that for iN2(M), ji,jN2(M)|mij|>0, it holds that

    |mii|>p(k2)i(M)>p(k1)i(M)>0,iN2(M),

    and

    |mii|>p(k1)i(M),iN.

    We now prove that |mii|>p(k1)i(M)>0(iN) and consider the following two cases.

    Case 1: If N1(M)=, then M is an SDD matrix. It is easy to verify that |mii|>p(k1)i(M)>0, iN2(M)=N.

    Case 2: If N1(M), by Theorem 5 and the known condition that for iN2(M), ji,jN2(M)|mij|>0, then it is easy to obtain that |mii|>p(k1)i(M)>0(iN).

    From Cases 1 and 2, we have that |mii|>p(k1)i(M)>0(iN). The proof is completed.

    Theorem 7. Let M=(mij)Cn×n(n2) be an SDDk(k2) matrix and for iN2(M), ji,jN2(M)|mij|>0. Then there exists a diagonal matrix X=diag{x1,x2,,xn} such that MX is an SDD matrix. Elements x1,x2,,xn are determined by

    xi=p(k1)i(M)|mii|,iN.

    Proof. We need to prove that matrix MX satisfies the following inequalities:

    |(MX)ii|>ri(MX),iN.

    From Theorem 6 and the known condition that for iN2(M), ji,jN2(M)|mij|>0, it is easy to verify that

    0<p(k1)i(M)|mii|<p(k2)i(M)|mii|<1,iN2(M).

    For each iN, we have that |(MX)ii|=p(k1)i(M) and

    ri(MX)=j=1,ji|mij|xj=jN1(M){i}p(k1)j(M)|mjj||mij|+jN2(M){i}p(k1)j(M)|mjj||mij|<jN1(M){i}|mij|+jN2(M){i}p(k2)j(M)|mjj||mij|=p(k1)i(M)=|(MX)ii|,

    that is,

    |(MX)ii|>ri(MX).

    Therefore, MX is an SDD matrix. The proof is completed.

    In this section, by Lemma 2 and Theorem 7, we provide some new upper bounds of the infinity norm of the inverse matrices for SDDk matrices and SDD matrices, respectively. We also present some comparisons with the Varah bound. Some numerical examples are presented to illustrate the corresponding results. Specially, when the involved matrices are SDD1 matrices as subclass of SDDk matrices, these new bounds are in line with that provided by Chen et al. [20].

    Theorem 8. Let M=(mij)Cn×n(n2) be an SDDk(k2) matrix. Then

    ||M1||max{1,maxiN2(M)p(k1)i(M)|mii|+ε}min{miniN1(M)Ui,miniN2(M)Vi},

    where

    Ui=|mii|jN1(M){i}|mij|jN2(M){i}(p(k1)j(M)|mjj|+ε)|mij|,Vi=ε(|mii|jN2(M){i}|mij|)+jN2(M){i}(p(k2)j(M)p(k1)j(M)|mjj|)|mij|,

    and

    0<ε<miniN|mii|p(k1)i(M)jN2(M){i}|mij|.

    Proof. By Lemma 2, we have that there exists a positive diagonal matrix X such that MX is an SDD matrix, where X is defined as Lemma 2. Thus,

    ||M1||=||X(X1M1)||=||X(MX)1||||X||||(MX)1||,

    and

    ||X||=max1inxi=max{1,maxiN2(M)p(k1)i(M)|mii|+ε}.

    Notice that MX is an SDD matrix. Hence, by Theorem 1, we have

    ||(MX)1||1min1in(|(MX)ii|ri(MX)).

    Thus, for any iN1(M), it holds that

    |(MX)ii|ri(MX)=|mii|jN1(M){i}|mij|jN2(M){i}(p(k1)j(M)|mjj|+ε)|mij|=Ui.

    For any iN2(M), it holds that

    |(MX)ii|ri(MX)=p(k1)i(M)+ε|mii|jN1(M){i}|mij|jN2(M){i}(p(k1)j(M)|mjj|+ε)|mij|=jN1(M){i}|mij|+jN2(M){i}p(k2)j(M)|mjj||mij|+ε|mii|jN1(M){i}|mij|jN2(M){i}(p(k1)j(M)|mjj|+ε)|mij|=ε(|mii|jN2(M){i}|mij|)+jN2(M){i}(p(k2)j(M)p(k1)j(M)|mjj|)|mij|=Vi.

    Therefore, we get

    ||M1||max{1,maxiN2(M)p(k1)i(M)|mii|+ε}min{miniN1(M)Xi,miniN2(M)Yi}.

    The proof is completed.

    From Theorem 8, it is easy to obtain the following result.

    Corollary 1. Let M=(mij)Cn×n(n2) be an SDD matrix. Then

    ||M1||maxiNp(k1)i(M)|mii|+εminiNZi,

    where k2,

    Zi=ε(|mii|ri(M))+jN{i}(p(k2)j(M)p(k1)j(M)|mjj|)|mij|,

    and

    0<ε<miniN|mii|p(k1)i(M)ri(M).

    Example 1. Consider the n×n matrix:

    M=(421.51.5424824824824824823.54).

    Take that n=20. It is easy to verify that M is an SDD matrix.

    By calculations, we have that for k=2,

    maxiNp(1)i(M)|mii|+ε1=0.5859+ε1,miniNZi=0.4414+0.5ε1,0<ε1<0.4732.

    For k=4,

    maxiNp(3)i(M)|aii|+ε2=0.3845+ε2,miniNZi=0.2959+0.5ε2,0<ε2<0.7034.

    For k=6,

    maxiNp(5)i(M)|mii|+ε3=0.2504+ε3,miniNZi=0.1733+0.5ε3,0<ε3<0.8567.

    For k=8,

    maxiNp(7)i(M)|mii|+ε4=0.1624+ε4,miniNZi=0.0990+0.5ε4,0<ε4<0.9572.

    So, when k=2,4,6,8, by Corollary 1 and Theorem 1, we can get the upper bounds for ||M1||, see Table 1. Thus,

    ||M1||0.5859+ε10.4414+0.5ε1<2,||M1||0.3845+ε20.2959+0.5ε2<2,
    Table 1.  The bounds in Corollary 1 and Theorem 1.
    k 2 4 6 8
    Cor 1 0.5859+ε10.4414+0.5ε1 0.3845+ε20.2959+0.5ε2 0.2504+ε30.1733+0.5ε3 0.1624+ε40.0990+0.5ε4
    Th 1 2 2 2 2

     | Show Table
    DownLoad: CSV

    and

    ||M1||0.2504+ε30.1733+0.5ε3<2,||M1||0.1624+ε40.0990+0.5ε4<2.

    Moreover, when k=1, by Theorem 2, we have

    ||M1||0.7188+ε50.4844+0.5ε5,0<ε5<0.3214.

    The following Theorem 9 shows that the bound in Corollary 1 is better than that in Theorem 1 of [20] in some cases.

    Theorem 9. Let matrix M=(mij)Cn×n(n2) be an SDD matrix. If there exists k2 such that

    maxiNp(k1)i(M)|mii|miniN(|mii|ri(M))miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|,

    then

    ||M1||maxiNp(k1)i(M)|mii|+εminiNZi1min1in(|mii|ri(M)),

    where Zi and ε are defined as in Corollary 1, respectively.

    Proof. From the given condition, we have that there exists k2 such that

    maxiNp(k1)i(M)|mii|miniN(|mii|ri(M))miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|,

    then

    maxiNp(k1)i(M)|mii|miniN(|mii|ri(M))+εminiN(|mii|ri(M))miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|+εminiN(|mii|ri(M)).

    Thus, we get

    (maxiNp(k1)i(M)|mii|+ε)miniN(|mii|ri(M))miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|+εminiN(|mii|ri(M))=miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|+miniN(ε(|mii|ri(M)))miniN(ε(|mii|ri(M))+jN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|)=miniNZi.

    Since M is an SDD matrix, then

    |mii|>ri(M),Zi>0,iN.

    It's easy to verify that

    maxiNp(k1)i(M)|mii|+εminiNZi1min1in(|mii|ri(M)).

    Thus, by Corollary 1, it holds that

    ||M1||maxiNp(k1)i(M)|mii|+εminiNZi1min1in(|mii|ri(M)).

    The proof is completed.

    We illustrate Theorem 9 by the following Example 2.

    Example 2. This is the previous Example 1. For k=4, we have

    maxiNp(3)i(M)|mii|miniN(|mii|ri(M))=0.1923<0.2959=miniNjN{i}p(2)j(M)p(3)j(M)|mjj||mij|.

    Thus, by Theorem 8, we obtain that for each 0<ε2<0.7034,

    ||M1||0.3845+ε20.2959+0.5ε2<2=1min1in(|mii|ri(M)).

    However, we find that the upper bounds in Theorems 8 and 9 contain the parameter ε. Next, based on Theorem 7, we will provide new upper bounds for the infinity norm of the inverse matrices of SDDk matrices, which only depend on the elements of the given matrices.

    Theorem 10. Let M=(mij)Cn×n(n2) be an SDDk(k2) matrix and for each iN2(M), ji,jN2(M)|mij|>0. Then

    ||M1||maxiNp(k1)i(M)|mii|miniN(p(k1)i(M)jN{i}p(k1)j(M)|mjj||mij|).

    Proof. By Theorems 7 and 8, we have that there exists a positive diagonal matrix X such that MX is an SDD matrix, where X is defined as in Theorem 7. Thus, it holds that

    ||M1||=||X(X1M1)||=||X(MX)1||||X||||(MX)1||,

    and

    ||X||=max1inxi=maxiNp(k1)i(M)|mii|.

    Notice that MX is an SDD matrix. Thus, by Theorem 1, we get

    ||(MX)1||1min1in(|(MX)ii|ri(MX))=1min1in(|miixi|ri(MX))=1miniN(p(k1)i(M)jN{i}p(k1)j(M)|mjj||mij|).

    Therefore, we have that

    ||M1||maxiNp(k1)i(M)|mii|miniN(p(k1)i(M)jN{i}p(k1)j(M)|mjj||mij|).

    The proof is completed.

    Since SDD matrices are a subclass of SDDk matrices, by Theorem 10, we can obtain the following result.

    Corollary 2. Let M=(mij)Cn×n(n2) be an SDD matrix. If ri(M)>0(iN), then there exists k2 such that

    ||M1||maxiNp(k1)i(M)|mii|miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|.

    Two examples are given to show the advantage of the bound in Theorem 10.

    Example 3. Consider the following matrix:

    M=(4012120104.14620233480462023042040).

    It is easy to verify that M is not an SDD matrix, an SDD1 matrix, a GSDD1 matrix, an S-SDD matrix, nor a CKV-type matrix. Therefore, we cannot use the error bounds in [1, 8, 9, 18, 20] to estimate ||M1||. But, M is an SDD2 matrix. So by the bound in Theorem 10, we have that M10.5820.

    Example 4. Consider the tri-diagonal matrix MRn×n arising from the finite difference method for free boundary problems [18], where

    M=(b+αsin(1n)c00ab+αsin(2n)c00ab+αsin(n1n)c00ab+αsin(1)).

    Take that n=4, a=1, b=0, c=3.7 and α=10. It is easy to verify that M is neither an SDD matrix nor an SDD1 matrix. However, we can get that M is a GSDD1 matrix and an SDD3 matrix. By the bound in Theorem 10, we have

    M18.2630,

    while by the bound in Theorem 4, it holds that

    M1εmin{2.1488ε,0.3105,2.474ε3.6272},ε(1.4661,2.1488).

    The following two theorems show that the bound in Corollary 2 is better than that in Theorem 1 in some cases.

    Theorem 11. Let M=(mij)Cn×n(n2) be an SDD matrix. If ri(M)>0(iN) and there exists k2 such that

    miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|miniN(|mii|ri(M)),

    then

    ||M1||maxiNp(k1)i(M)|mii|miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|<1min1in(|mii|ri(M)).

    Proof. Since M is an SDD matrix, then N1(M)= and M is an SDDk matrix. By the given condition that ri(M)>0(iN), it holds that

    |mii|>ri(M)>jN{i}rj(M)|mjj||mij|=p(0)i(M)>0,iN,p(0)i(M)=jN{i}rj(M)|mjj||mij|>jN{i}p(0)j(M)|mjj||mij|=p(1)i(M)>0,iN.

    Similarly, we can obtain that

    |mii|>ri(M)>p(0)i(M)>>p(k1)i(M)>0,iN,

    that is,

    maxiNp(k1)i(M)|mii|<1.

    Since there exists k2 such that

    miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|miniN(|mii|ri(M)),

    then we have

    maxiNp(k1)i(M)|mii|miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|<1min1in(|mii|ri(M)).

    Thus, from Corollary 2, we get

    ||M1||maxiNp(k1)i(M)|mii|miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|<1min1in(|mii|ri(M)).

    The proof is completed.

    We illustrate the Theorem 11 by following Example 5.

    Example 5. Consider the matrix M=(mij)Cn×n(n2), where

    M=(430.91622522522522521620.934).

    Take that n=20. It is easy to check that M is an SDD matrix. Let

    lk=miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|,m=miniN(|mii|ri(M)).

    By calculations, we have

    l2=0.2692>0.1=m,l3=0.2567>0.1=m,l4=0.1788>0.1=m,l5=0.1513>0.1=m,l6=0.1037>0.1=m.

    Thus, when k=2,3,4,5,6, the matrix M satisfies the conditions of Theorem 11. By Theorems 1 and 11, we can derive the upper bounds for ||M1||, see Table 2. Meanwhile, when k=1, by Theorem 3, we get that ||M1||1.6976.

    Table 2.  The bounds in Theorem 11 and Theorem 1.
    k 2 3 4 5 6
    Th 11 1.9022 1.5959 1.8332 1.7324 2.0214
    Th 1 10 10 10 10 10

     | Show Table
    DownLoad: CSV

    From Table 2, we can see that the bounds in Theorem 11 are better than that in Theorems 1 and 3 in some cases.

    Theorem 12. Let M=(mij)Cn×n(n2) be an SDD matrix. If ri(M)>0(iN) and there exists k2 such that

    maxiNp(k1)i(M)|mii|miniN(|mii|ri(M))miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|<miniN(|mii|ri(M)),

    then

    ||M1||maxiNp(k1)i(M)|mii|miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|1min1in(|mii|ri(M)).

    Proof. By Theorem 7 and the given condition that ri(M)>0(iN), it is easy to get that

    jN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|>0,iN.

    From the condition that there exists k2 such that

    maxiNp(k1)i(M)|mii|miniN(|mii|ri(M))miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|,

    we have

    maxiNp(k1)i(M)|mii|miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|1min1in(|mii|ri(M)).

    Thus, from Corollary 2, it holds that

    ||M1||maxiNp(k1)i(M)|mii|miniNjN{i}p(k2)j(M)p(k1)j(M)|mjj||mij|1min1in(|mii|ri(M)).

    The proof is completed.

    Next, we illustrate Theorem 12 by the following Example 6.

    Example 6. Consider the tri-diagonal matrix M=(mij)Cn×n(n2), where

    M=(32.51.2422.8512.8512.8512.8511.2422.53).

    Take that n=20. It is easy to verify that M is an SDD matrix.

    By calculations, we have that for k=2,

    maxiNp(1)i(M)|mii|miniN(|mii|ri(M))=0.2686<miniNjN{i}p(0)j(M)p(1)j(M)|mjj||mij|=0.3250<0.5=miniN(|mii|ri(M)).

    For k=5, we get

    maxiNp(4)i(M)|mii|miniN(|mii|ri(M))=0.1319<miniNjN{i}p(3)j(M)p(4)j(M)|mjj||mij|=0.1685<0.5=miniN(|mii|ri(M)).

    For k=10, it holds that

    maxiNp(9)i(M)|mii|miniN(|mii|ri(M))=0.0386<miniNjN{i}p(8)j(M)p(9)j(M)|mjj||mij|=0.0485<0.5=miniN(|mii|ri(M)).

    Thus, for k=2,5,10, the matrix M satisfies the conditions of Theorem 12. Thus, from Theorems 12 and 1, we get the upper bounds for ||M1||, see Table 3. Meanwhile, when k=1, by Theorem 3, we have that ||M1||1.7170.

    Table 3.  The bounds in Theorem 12 and Theorem 1.
    k 2 5 10
    Th 12 1.6530 1.5656 1.5925
    Th 1 2 2 2

     | Show Table
    DownLoad: CSV

    From Table 3, we can see that the bound in Theorem 12 is sharper than that in Theorems 1 and 3 in some cases.

    SDDk matrices as a new subclass of H-matrices are proposed, which include SDD matrices and SDD1 matrices, and some properties of SDDk matrices are obtained. Meanwhile, some new upper bounds of the infinity norm of the inverse matrices for SDD matrices and SDDk matrices are presented. Furthermore, we prove that the new bounds are better than some existing bounds in some cases. Some numerical examples are also provided to show the validity of new results. In the future, based on the proposed infinity norm bound, we will explore the computable error bounds of the linear complementarity problems for SDDk matrices.

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

    This research is supported by Guizhou Provincial Science and Technology Projects (20191161), and the Natural Science Research Project of Department of Education of Guizhou Province (QJJ2023062, QJJ2023063).

    The authors declare that they have no competing interests.



    [1] L. Stone, R. Olinky, A. Huppert, Seasonal dynamics of recurrent epidemics, Nature, 446 (2007), 533–536. https://doi.org/10.1038/nature05638 doi: 10.1038/nature05638
    [2] M. Sinan, K. J. Ansari, A. Kanwal, K. Shah, T. Abdeljawad, B. Abdalla, et al., Analysis of the mathematical model of cutaneous leishmaniasis disease, Alex. Eng. J., 72 (2023), 117–134. https://doi.org/10.1016/j.aej.2023.03.065 doi: 10.1016/j.aej.2023.03.065
    [3] A. R. Sheergojri, P. Iqbal, P. Agarwal, N. Ozdemir, Uncertainty-based Gompertz growth model for tumor population and its numerical analysis, International Journal of Optimization and Control: Theories and Applications, 12 (2022), 137–150. https://doi.org/10.11121/ijocta.2022.1208 doi: 10.11121/ijocta.2022.1208
    [4] Y. Sabbar, Asymptotic extinction and persistence of a perturbed epidemic model with different intervention measures and standard Leˊvy jumps, Bulletin of Biomathematics, 1 (2023), 58–77. https://doi.org/10.59292/bulletinbiomath.2023004 doi: 10.59292/bulletinbiomath.2023004
    [5] F. Evirgen, U. Esmehan, U. Sümeyra, N. Ozdemir, Modelling influenza a disease dynamics under Caputo-Fabrizio fractional derivative with distinct contact rates, Mathematical Modelling and Numerical Simulation with Applications, 3 (2023), 58–73. https://doi.org/10.53391/mmnsa.1274004 doi: 10.53391/mmnsa.1274004
    [6] K. Shah, T. Abdeljawad, H. Alrabaish, On coupled system of drug therapy via piecewise equations, Fractals, 30 (2022), 2240206. https://doi.org/10.1142/S0218348X2240206X doi: 10.1142/S0218348X2240206X
    [7] K. Shah, T. Abdeljawad, F. Jarad, F. Jarad, Q. Al-Mdallal, On nonlinear conformable fractional order dynamical system via differential transform method, CMES Comp. Model. Eng., 163 (2023), 1457–1472. http://doi.org/10.32604/cmes.2023.021523 doi: 10.32604/cmes.2023.021523
    [8] S. W. Ahmad, M. Sarwar, K. Shah, A. Ahmadian, S. Salahshour, Fractional order mathematical modeling of novel corona virus (COVID‐19), Math. Methods Appl. Sci., 46 (2023), 7847–7860. https://doi.org/10.1002/mma.7241 doi: 10.1002/mma.7241
    [9] L. Liu, X. Ren, X. Liu, Dynamical behaviors of an influenza epidemic model with virus mutation, J. Biol. Syst., 26 (2018), 455–472. https://doi.org/10.1142/S0218339018500201 doi: 10.1142/S0218339018500201
    [10] B. Li, A. Deng, K. Li, Y. Hu, Z. C. Li, Y. L. Shi, et al., Viral infection and transmission in a large, well-traced outbreak caused by the SARS-CoV-2 Delta variant, Nat. Commun., 13 (2022), 460. https://doi.org/10.1038/s41467-022-28089-y doi: 10.1038/s41467-022-28089-y
    [11] Y. Liu, A. Feng, S. Zhao, W. Wang, D. He, Large-scale synchronized replacement of Alpha (B.1.1.7) variant by the Delta (B.1.617.2) variant of SARS-COV-2 in the COVID-19 pandemic, Math. Biosci. Eng., 19 (2022), 3591–3596. https://doi.org/10.3934/mbe.2022165 doi: 10.3934/mbe.2022165
    [12] R. M. Chen, Track the dynamical features for mutant variants of COVID-19 in the UK, Math. Biosci. Eng., 18 (2021), 4572–4585. https://doi.org/10.3934/mbe.2021232 doi: 10.3934/mbe.2021232
    [13] Y. Yu, Y. Liu, S. Zhao, D. He, A simple model to estimate the transmissibility of the Beta, Delta, and Omicron variants of SARS-COV-2 in South Africa, Math. Biosci. Eng., 19 (2022), 10361–10373. https://doi.org/10.3934/mbe.2022485 doi: 10.3934/mbe.2022485
    [14] S. P. Otto, T. Day, J. Arino, C. Colijn, J. Dushoff, M. Li, et al., The origins and potential future of SARS-CoV-2 variants of concern in the evolving COVID-19 pandemic, Curr. Biol., 31 (2021), R918–R929. https://doi.org/10.1016/j.cub.2021.06.049 doi: 10.1016/j.cub.2021.06.049
    [15] G. Cacciapaglia, C. Cot, A. D. Hoffer, S. Hohenegger, F. Sannino, S. Vatani, Epidemiological theory of virus variants, Phys. A, 596 (2022), 127071. https://doi.org/10.1016/j.physa.2022.127071 doi: 10.1016/j.physa.2022.127071
    [16] A. P. Dobie, Susceptible-infectious-susceptible (SIS) model with virus mutation in a variable population size, Ecol. Complex., 50 (2022), 101004. https://doi.org/10.1016/j.ecocom.2022.101004 doi: 10.1016/j.ecocom.2022.101004
    [17] U. A. de León, A. G. C. Pérez, E. Avila-Vales, Modeling the SARS-CoV-2 Omicron variant dynamics in the United States with booster dose vaccination and waning immunity, Math. Biosci. Eng., 20 (2023), 10909–10953. https://doi.org/10.3934/mbe.2023484 doi: 10.3934/mbe.2023484
    [18] Y. R. Kim, Y. J. Choi, Y. Min, A model of COVID-19 pandemic with vaccines and mutant viruses, Plos One, 17 (2022), e0275851. https://doi.org/10.1371/journal.pone.0275851 doi: 10.1371/journal.pone.0275851
    [19] G. Liu, J. Chen, Z. Liang, Z. Peng, J. Li, Dynamical analysis and optimal control for a SEIR model based on virus mutation in WSNs, Mathematics, 9 (2021), 929. https://doi.org/10.3390/math9090929 doi: 10.3390/math9090929
    [20] D. Xu, X. Xu, Y. Xie, C. Yang, Optimal control of an SIVRS epidemic spreading model with virus variation based on complex networks, Commun. Nonlinear Sci., 48 (2017), 200–210. https://doi.org/10.1016/j.cnsns.2016.12.025 doi: 10.1016/j.cnsns.2016.12.025
    [21] Q. Liu, D. Jiang, N. Shi, T. Hayat, Dynamics of a stochastic delayed SIR epidemic model with vaccination and double diseases driven by L\mathrm{\acute{e}}vy jumps, Phys. A, 492 (2018), 2010–2018. https://doi.org/10.1016/j.physa.2017.11.116 doi: 10.1016/j.physa.2017.11.116
    [22] X. Zhang, M. Liu, Dynamical analysis of a stochastic delayed SIR epidemic model with vertical transmission and vaccination, Adv. Cont. Discr. Mod., 2022 (2022), 35. https://doi.org/10.1186/s13662-022-03707-7 doi: 10.1186/s13662-022-03707-7
    [23] C. Xu, X. Li, The threshold of a stochastic delayed SIRS epidemic model with temporary immunity and vaccination, Chaos Soliton. Fract., 111 (2018), 227–234. https://doi.org/10.1016/j.chaos.2021.110772 doi: 10.1016/j.chaos.2021.110772
    [24] A. E. Koufi, The power of delay on a stochastic epidemic model in a switching environment, Complexity, 2022 (2022), 5121636. https://doi.org/10.1155/2022/5121636 doi: 10.1155/2022/5121636
    [25] B. Boukanjime, M. El-Fatini, A. Laaribi, R. Taki, K. Wang, A Markovian regime-switching stochastic hybrid time-delayed epidemic model with vaccination, Automatica, 133 (2021), 109881. https://doi.org/10.1016/j.automatica.2021.109881 doi: 10.1016/j.automatica.2021.109881
    [26] H. J. Alsakaji, F. A. Rihan, S. Kundu, O. Mohamed, Dynamics of a time-delay differential model for tumour-immune interactions with random noise, Alex. Eng. J., 61 (2022), 11913–11923. https://doi.org/10.1016/j.aej.2022.05.027 doi: 10.1016/j.aej.2022.05.027
    [27] I. Ali, S. U. Khan, Analysis of stochastic delayed SIRS model with exponential birth and saturated incidence rate, Chaos Soliton. Fract., 138 (2022), 110008. https://doi.org/10.1016/j.chaos.2020.110008 doi: 10.1016/j.chaos.2020.110008
    [28] H. J. Alsakaji, F. A. Rihan, A. Hashish, Dynamics of a stochastic epidemic model with vaccination and multiple time-delays for COVID-19 in the UAE, Complexity, 2022 (2022), 4247800. https://doi.org/10.1155/2022/4247800 doi: 10.1155/2022/4247800
    [29] A. Khan, R. Ikram, A. Din, U. W. Humphries, A. Akgul, Stochastic COVID-19 SEIQ epidemic model with time-delay, Results Phys., 30 (2021), 104775. https://doi.org/10.1016/j.rinp.2021.104775 doi: 10.1016/j.rinp.2021.104775
    [30] F. A. Rihan, H. J. Alsakaji, Analysis of a stochastic HBV infection model with delayed immune response, Math. Biosci. Eng., 18 (2021), 5194–5220. https://doi.org/10.3934/mbe.2021264 doi: 10.3934/mbe.2021264
    [31] R. Ikram, A. Khan, M. Zahri, A. Saeed, M. Yavuz, P. Kumam, Extinction and stationary distribution of a stochastic COVID-19 epidemic model with time-delay, Comput. Biol. Med., 141 (2022), 105115. https://doi.org/10.1016/j.compbiomed.2021.105115 doi: 10.1016/j.compbiomed.2021.105115
    [32] Q. Liu, D. Jiang, T. Hayat, A. Alsaedi, Dynamics of a stochastic SIR epidemic model with distributed delay and degenerate diffusion, J. Franklin I., 356 (2019), 7347–7370. https://doi.org/10.1016/j.jfranklin.2019.06.030 doi: 10.1016/j.jfranklin.2019.06.030
    [33] J. Sun, M. Gao, D. Jiang, Threshold dynamics of a non-linear stochastic viral model with time delay and CTL responsiveness, Life, 11 (2021), 766. https://doi.org/10.3390/life11080766 doi: 10.3390/life11080766
    [34] F. A. Rihan, H. J. Alsakaji, Dynamics of a stochastic delay differential model for COVID-19 infection with asymptomatic infected and interacting people: case study in the UAE, Results Phys., 28 (2021), 104658. https://doi.org/10.1016/j.rinp.2021.104658 doi: 10.1016/j.rinp.2021.104658
    [35] R. Khasminskii, Stochastic stability of differential equations, Berlin: Springer, 2011. https://doi.org/10.1007/978-3-642-23280-0
    [36] E. Buckwar, Euler-Maruyama and Milstein approximations for stochastic functional differential equations with distributed memory term, Berlin: Humboldt-Universität, 2005. http://doi.org/10.18452/3583
    [37] D. J. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations, SIAM Rev., 43 (2001), 525–546. https://doi.org/10.1137/S0036144500378302 doi: 10.1137/S0036144500378302
    [38] J. Gao, T. Zhang, Analysis on an SEIR epidemic model with logistic death rate of virus mutation, Journal of Mathematical Research with Applications, 39 (2019), 259–268. https://doi.org/10.3770/j.issn:2095-2651.2019.03.005 doi: 10.3770/j.issn:2095-2651.2019.03.005
  • This article has been cited by:

    1. Qin Li, Wenwen Ran, Feng Wang, Infinity norm bounds for the inverse of generalized {SDD_2} matrices with applications, 2024, 41, 0916-7005, 1477, 10.1007/s13160-024-00658-2
    2. Qin Li, Wenwen Ran, Feng Wang, Infinity norm bounds for the inverse of Quasi-SDD_k matrices with applications, 2024, 1017-1398, 10.1007/s11075-024-01949-y
    3. Wenwen Ran, Feng Wang, Extended SDD_1^{\dag } matrices and error bounds for linear complementarity problems, 2024, 0916-7005, 10.1007/s13160-024-00685-z
    4. Yuanjie Geng, Yuxue Zhu, Fude Zhang, Feng Wang, Infinity Norm Bounds for the Inverse of \textrm{SDD}_1-Type Matrices with Applications, 2025, 2096-6385, 10.1007/s42967-024-00457-z
    5. L. Yu. Kolotilina, SSDD Matrices and Relations with Other Subclasses of the Nonsingular ℋ-Matrices, 2025, 1072-3374, 10.1007/s10958-025-07711-6
  • 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(1469) PDF downloads(72) Cited by(0)

Figures and Tables

Figures(6)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog