In this paper, a class of survival red blood cells model with time-varying delays and impulsive effects is considered. First, some sufficient conditions for the persistence are derived by use of the theory on impulsive differential equations. The persistence describes the persistent survival of the mature red blood cells in the mammal under delay and impulsive perturbations. Then assuming that the coefficients in the model are ω-periodic, some criteria ensuring the existence-uniqueness and global attractivity of positive ω-periodic solution of the addressed model are obtained, which are suitable for survival red blood cells model with any ω∈R+. These global attractivity criteria describe the nonexistence of any dynamic diseases in the mammal. Moreover, our proposed results in this paper extend and improve some recent works in the literature. Finally, two examples and their computer simulations are given to show the effectiveness and advantages of the results.
Citation: Tengda Wei, Xiang Xie, Xiaodi Li. Persistence and periodicity of survival red blood cells model with time-varying delays and impulses[J]. Mathematical Modelling and Control, 2021, 1(1): 12-25. doi: 10.3934/mmc.2021002
[1] | Hui Yan, Hongxing Wang, Kezheng Zuo, Yang Chen . Further characterizations of the weak group inverse of matrices and the weak group matrix. AIMS Mathematics, 2021, 6(9): 9322-9341. doi: 10.3934/math.2021542 |
[2] | Jinyong Wu, Wenjie Shi, Sanzhang Xu . Revisiting the m-weak core inverse. AIMS Mathematics, 2024, 9(8): 21672-21685. doi: 10.3934/math.20241054 |
[3] | Xiaofei Cao, Yuyue Huang, Xue Hua, Tingyu Zhao, Sanzhang Xu . Matrix inverses along the core parts of three matrix decompositions. AIMS Mathematics, 2023, 8(12): 30194-30208. doi: 10.3934/math.20231543 |
[4] | Yang Chen, Kezheng Zuo, Zhimei Fu . New characterizations of the generalized Moore-Penrose inverse of matrices. AIMS Mathematics, 2022, 7(3): 4359-4375. doi: 10.3934/math.2022242 |
[5] | Wanlin Jiang, Kezheng Zuo . Revisiting of the BT-inverse of matrices. AIMS Mathematics, 2021, 6(3): 2607-2622. doi: 10.3934/math.2021158 |
[6] | Hongjie Jiang, Xiaoji Liu, Caijing Jiang . On the general strong fuzzy solutions of general fuzzy matrix equation involving the Core-EP inverse. AIMS Mathematics, 2022, 7(2): 3221-3238. doi: 10.3934/math.2022178 |
[7] | Almudena Campos-Jiménez, Francisco Javier García-Pacheco . The core of the unit sphere of a Banach space. AIMS Mathematics, 2024, 9(2): 3440-3452. doi: 10.3934/math.2024169 |
[8] | Yuxu Chen, Hui Kou . Core compactness of ordered topological spaces. AIMS Mathematics, 2023, 8(2): 4862-4874. doi: 10.3934/math.2023242 |
[9] | Wanlin Jiang, Kezheng Zuo . Further characterizations of the m-weak group inverse of a complex matrix. AIMS Mathematics, 2022, 7(9): 17369-17392. doi: 10.3934/math.2022957 |
[10] | Liying Yang, Jinjin Li, Yiliang Li, Qifang Li . Sub-base local reduct in a family of sub-bases. AIMS Mathematics, 2022, 7(7): 13271-13277. doi: 10.3934/math.2022732 |
In this paper, a class of survival red blood cells model with time-varying delays and impulsive effects is considered. First, some sufficient conditions for the persistence are derived by use of the theory on impulsive differential equations. The persistence describes the persistent survival of the mature red blood cells in the mammal under delay and impulsive perturbations. Then assuming that the coefficients in the model are ω-periodic, some criteria ensuring the existence-uniqueness and global attractivity of positive ω-periodic solution of the addressed model are obtained, which are suitable for survival red blood cells model with any ω∈R+. These global attractivity criteria describe the nonexistence of any dynamic diseases in the mammal. Moreover, our proposed results in this paper extend and improve some recent works in the literature. Finally, two examples and their computer simulations are given to show the effectiveness and advantages of the results.
The weak core inverse was introduced in [1] where the authors presented some characterizations and properties. In [2], the authors introduced an extension of the weak core inverse. Continuing previous research about the weak core inverse, our purpose is to present new characterizations and representations of the weak core inverse. Additionally, we also give several equivalent conditions for a matrix to be a weak core matrix.
Let Cm×n be the set of all m×n complex matrices and Z+ denotes the set of all positive integers. The symbols R(A), N(A), A∗, r(A) and In will denote the range space, null space, conjugate transpose, rank of A∈Cm×n and the identity matrix of order n. Ind(A) means the index of A∈Cn×n. Let Cn×nk be the set consisting of all n×n complex matrices with index k. The symbol dim(S) represents the dimension of a subspace S⊆Cn. PL stands for the orthogonal projection onto the subspace L. PA, PA∗ respectively denote the orthogonal projection onto R(A) and R(A∗), i.e., PA=AA†, PA∗=A†A.
We will now introduce definitions of several generalized inverses that will be used throughout the paper. The Moore-Penrose inverse of A∈Cm×n, denoted by A†, is defined as the unique matrix X∈Cn×m satisfying [3]:
(1) AXA=A, (2) XAX=X, (3) (AX)∗=AX, (4) (XA)∗=XA. |
In particular, X is an outer inverse of A which is denoted as A(2) if XAX=X. For any matrix A∈Cm×n with r(A)=r, let T⊆Cn, S⊆Cm be two subspaces such that dim(T)=t≤r and dim(S)=m−t. Then A has an outer inverse X that satisfies R(X)=T and N(X)=S if and only if AT⊕S=Cm. In that case X is unique and denoted by A(2)T,S [4].
The Drazin inverse of A∈Cn×nk, denoted by AD, is the unique matrix X∈Cn×n satisfying [5]: XAX=X, AX=XA, XAk+1=Ak.
For any matrix A∈Cn×n1, a new generalized inverse, which is called core inverse [6] was introduced. Two other generalizations of the core inverse for A∈Cn×nk such as core-EP inverse [7], DMP inverse [8] were also introduced.
In 2018, Wang and Chen [9] defined the weak group inverse of A∈Cn×nk, denoted by AⓌ, as the unique matrix X∈Cn×n such that [9]: . Moreover, it was verified that
.
Recently, Ferreyra et al. introduced a new generalization of core inverse called the weak core inverse of A∈Cn×nk, denoted by AⓌ,†(or, in short, WC inverse). It is defined as the unique matrix X∈Cn×n satisfying [1]:
XAX=X, AX=CA†, XA=ADC, |
where C=AAⓌA. Moreover, it is proved that AⓌ,†=ADCA†=AⓌAA†.
The structure of this paper is as follows: In Section 2, we give some preliminaries which will be made use of later in this paper. In Section 3, we discuss some characterizations of the WC inverse based on its range space, null space and matrix equations. In Section 4, several new representations of the WC inverse are proposed. Section 5 is devoted to deriving some properties of the WC inverse by the core-EP decomposition. Moreover, in Section 6, we present several equivalent conditions for a matrix to be a weak core matrix.
For convenience, we will use the following notations: CCMn, CEPn, CPn and COPn will denote the subsets of Cn×n consisting of core matrices, EP matrices, idempotent matrices and Hermitian idempotent matrices, respectively, i.e.,
● CCMn={A∣A∈Cn×n,r(A2)=r(A)};
● CEPn={A∣A∈Cn×n,R(A)=R(A∗)};
● CPn={A∣A∈Cn×n,A2=A};
● COPn={A∣A∈Cn×n,A2=A=A∗}.
Before giving characterizations of the WC inverse, we first present the following auxiliary lemmas which will be repeatedly used throughout this paper.
Lemma 2.1. [10] Let A∈Cn×nk. Then A can be represented as
A=U[TS0N]U∗, | (2.1) |
where T∈Ct×t is nonsingular and t=r(T)=r(Ak), N is nilpotent with index k, and U∈Cn×n is unitary.
Moreover, the representation of A given by (2.1) is unique [10,Theorem 2.4]. In that case, we have that
Ak=U[Tk˜T00]U∗, | (2.2) |
where ˜T=k−1∑j=0TjSNk−1−j.
Lemma 2.2. [1,9,10,11,12] Let A∈Cn×nk be given by (2.1). Then :
A†=U[T∗△−T∗△SN†(In−t−N†N)S∗△N†−(In−t−N†N)S∗△SN†]U∗, | (2.3) |
AD=U[T−1(Tk+1)−1˜T00]U∗, | (2.4) |
![]() |
(2.5) |
AD,†=U[T−1(Tk+1)−1˜TNN†00]U∗, | (2.6) |
A†,D=U[T∗△T∗△T−k˜T(In−t−N†N)S∗△(In−t−N†N)S∗△T−k˜T]U∗, | (2.7) |
AⓌ=U[T−1T−2S00]U∗, | (2.8) |
.
AⓌ,†=U[T−1T−2SNN†00]U∗, | (2.9) |
where ˜T=k−1∑j=0TjSNk−1−j and △=[TT∗+S(In−t−N†N)S∗]−1.
˜T and △ will be often used throughout this paper.
Lemma 2.3. [1] Let A∈Cn×nk. Then
(a) AⓌ,†=A(2)R(Ak), N((Ak)∗A2A†);
(b) AAⓌ,†=PR(Ak), N((Ak)∗A2A†);
(c) AⓌ,†A=PR(Ak), N((Ak)∗A2).
Lemma 2.4. Let A∈Cn×nk and C=AAⓌA. The following conditions hold:
(a) [1]
(b) [1] r(AⓌ,†)=r(Ak);
(c) [1] CA†C=C;
(d) [9] AⓌAk+1=Ak;
(e) Ck=AkAⓌA.
Proof. Item (e) can be directly verified by (2.1), (2.2) and (2.8).
Applying existing results for the WC inverse with respect to R(X)=R(Ak) and N(X)=N((Ak)∗A2A†), some new results can be obtained for the WC inverse in the next result.
Theorem 3.1. Let A∈Cn×nk and C=AAⓌA. The following statements are equivalent:
(a) X=AⓌ,†;
(b) N(X)=N((Ak)∗A2A†) and XAA∗=AⓌAA∗;
(c) N(X)=N((Ak)∗A2A†) and XA=AⓌA;
(d) R(X)=R(Ak) and A∗AX=A∗CA†;
(e) R(X)=R(Ak) and AX=CA†;
(f) R(X)=R(Ak) and AkX=CkA†;
(g) R(X)=R(Ak), N(X)=N((Ak)∗A2A†) and XAAⓌ=AⓌ;
(h) R(X)=R(Ak), N(X)=N((Ak)∗A2A†) and XAk+1=Ak.
Proof. (a)⇒(b). By the definition of AⓌ,†, we have that XAA∗=ADCA∗=ADAAⓌAA∗=AⓌAA∗. Hence, by (a) of Lemma 2.3, we now obtain that (b) holds.
(b)⇒(c). Postmultiplying XAA∗=AⓌAA∗ by (A†)∗, we obtain that XA=AⓌA.
(c)⇒(d). From N(X)=N((Ak)∗A2A†), we have that N(AA†)⊆N((Ak)∗A2A†)=N(X), which leads to X=XAA†. Thus we get that X=XAA†=AⓌAA†=AⓌ,† by XA=AⓌA. Hence, by the definition of AⓌ,† and (a) of Lemma 2.3, we have that (d) holds.
(d)⇒(e). Evidently.
(e)⇒(f). Since C=AAⓌA and Ck=AkAⓌA, premultiplying AX=CA† by Ak−1, we have that AkX=CkA†.
(f)⇒(g). From (2.2) and R(X)=R(Ak), we can set X=U[X1X200]U∗, where X1∈Ct×t, X2∈Ct×(n−t) and t=r(Ak). Furthermore, it follows from AkX=CkA† and (2.9) that X=AⓌ,†. Therefore, by the definition of AⓌ,† and (a) of Lemma 2.3, we obtain that (g) holds.
(g)⇒(h). It follows from AⓌAk+1=Ak and XAAⓌ=AⓌ that XAk+1=XAAⓌAk+1=AⓌAk+1=Ak.
(h)⇒(a). By R(X)=R(Ak) and XAk+1=Ak, we get that XAX=X. Hence, by (a) of Lemma 2.3, we get that X=AⓌ,†.
Now we will consider other characterizations of the WC inverse by the fact that AⓌ,†AAⓌ,†=AⓌ,†.
Theorem 3.2. Let A∈Cn×nk and C=AAⓌA. The following statements are equivalent:
(a) X=AⓌ,†;
(b) XAX=X, R(X)=R(Ak) and N(X)=N((Ak)∗A2A†);
(c) XAX=X, R(X)=R(Ak) and AX=CA†;
(d) XAX=X, AX=CA† and XAk=AⓌAk;
(e) XAX=X, XA=AⓌA and AkX=CkA†;
(f) XAX=X, XA=AⓌA and N(X)=N((Ak)∗A2A†).
Proof. (a)⇒(b). The proof can be demonstrated by (a) of Lemma 2.3.
(b)⇒(c). By the definition of AⓌ,† and (b) of Lemma 2.3, we get that AX∈CPn, R(AX)=AR(X)=R(Ak+1)=R(Ak)=R(AAⓌ,†)=R(CA†) and N(AX)=N(X)=N((Ak)∗A2A†)=N(AAⓌ,†)=N(CA†). On the other hand, Lemma 2.4 (c) implies CA†∈CPn, hence AX=CA†.
(c)⇒(d). By item (c) of Lemma 2.3, we obtain that R(X)=R(Ak)=R(AⓌ,†A). So we get that AⓌ,†AX=X, which implies that XAk=AⓌ,†AXAk=AⓌ,†CA†Ak=AⓌ,†AAⓌAA†Ak=AⓌAk.
(d)⇒(e). By conditions and AAⓌ=Ak(AⓌ)k, we can infer that X=XCA†=XAAⓌAA†=XAk(AⓌ)kAA†=AⓌAk(AⓌ)kAA†=AⓌ,†. Hence, by AⓌ,†=AⓌAA† and AkAⓌA=Ck, we obtain that (e) holds.
(e)⇒(f). Since XAX=X, XA=AⓌA, we have that R(X)=R(XA)=R(AⓌA)=R(Ak). We now obtain that X=AⓌ,† by (f) of Theorem 3.1. Hence (f) holds by (a) of Lemma 2.3.
(f)⇒(a). It follows from XAX=X that N(AX)=N(X), by conditions and (a) of Lemma 2.3. We now obtain that X=XAX=AⓌAX=AⓌAA†AX=AⓌ,†PR(AX), N(AX)=AⓌ,†.
Notice the fact that XAk+1=Ak if X=AⓌ,†. Therefore, we will characterize the WC inverse in terms of AⓌ,†Ak+1=Ak.
Theorem 3.3. Let A∈Cn×nk and C=AAⓌA. The following statements are equivalent:
(a) X=AⓌ,†;
(b) XAk+1=Ak, A∗AX=A∗CA† and r(X)=r(Ak);
(c) XAk+1=Ak, AX=CA† and r(X)=r(Ak);
(d) XAk+1=Ak, AkX=CkA† and r(X)=r(Ak).
Proof. (a)⇒(b). Since AⓌ,†=AⓌAA†, we can show that XAk+1=Ak, A∗AX=A∗CA†. Then, by (b) of Lemma 2.4, we get that (b) holds.
(b)⇒(c). Obviously.
(c)⇒(d). Premultiplying AX=CA† by Ak−1, we have that AkX=CkA† from AkAⓌA=Ck.
(d)⇒(a). It follows from XAk+1=Ak and r(X)=r(Ak) that R(X)=R(Ak). Hence, we obtain that X=AⓌ,† from (f) of Theorem 3.1.
In the following example, we show that the condition r(X)=r(Ak) in Theorem 3.3 is necessary.
Example 3.4. Let
A=[100003000], X=[100 002000]. |
Then Ind(A)=2,
A†=[10000001/30], C=[100000000] and AⓌ,†=[100000000]. |
It can be directly verified that XA3=A2, A∗AX=A∗CA† and r(X)≠r(A2), but X≠AⓌ,†. The other cases follow similarly.
By Lemma 2.3, it is clear that AX=PR(Ak), N((Ak)∗A2A†) and XA=PR(Ak), N((Ak)∗A2) if X=AⓌ,†. However, the converse is invalid as shown in the next example:
Example 3.5. Let A, X be the same as in Example 3.4. Then
AX=[100000000], XA=[100000000]and AⓌ,†=[100000000]. |
It can be directly verified that AX=PR(A2), N((A2)∗A2A†) and XA=PR(A2), N((A2)∗A2), but X≠AⓌ,†.
In the next result, we will present some new equivalent conditions for the converse implication:
Theorem 3.6. Let A∈Cn×nk and X∈Cn×n. The following statements are equivalent:
(a) X=AⓌ,†;
(b) AX=PR(Ak), N((Ak)∗A2A†),XA=PR(Ak), N((Ak)∗A2) and r(X)=r(Ak);
(c) AX=PR(Ak), N((Ak)∗A2A†),XA=PR(Ak), N((Ak)∗A2) and XAX=X;
(d) AX=PR(Ak), N((Ak)∗A2A†),XA=PR(Ak), N((Ak)∗A2) and AX2=X.
Proof. (a)⇒(b). The proof can be demonstrated by (b) and (c) of Lemma 2.3 and (b) of Lemma 2.4.
(b)⇒(c). By R(XA)=R(Ak) and r(X)=r(Ak), we obtain that R(X)=R(XA)=R(Ak), hence we further derive that XAX=X.
(c)⇒(d). By conditions and (a) of Lemma 2.3, we have that X=AⓌ,†. Therefore, by (2.9), it can be directly verified that AX2=X.
(d)⇒(a). From AX2=X, we have that X=AX2=A2X3=⋯=AkXk+1, which implies R(X)⊆R(Ak). Combining with the condition R(Ak)=R(XA)⊆R(X), we get that R(X)=R(Ak). From (2.2), we now set X=U[X1X200]U∗, where X1∈Ct×t, X2∈Ct×(n−t) and t=r(Ak). On the other hand, it follows from N(AX)=N((Ak)∗A2A†) that (Ak)∗A2A†=(Ak)∗A2X, which yields X1=T−1 and X2=T−2SNN†. Therefore, by (2.9), we obtain that X=AⓌ,†.
In [1], the authors introduced the definition of AⓌ,† with an algebraic approach. In the next result, we will consider characterization of AⓌ,† with a geometrical point of view.
Theorem 3.7. Let A∈Cn×nk. Then:
(a) AⓌ,† is the unique matrix X that satisfies:
AX=PR(Ak), N((Ak)∗A2A†), R(X)⊆R(Ak). | (3.1) |
(b) AⓌ,† is the unique matrix X that satisfies:
XA=PR(Ak), N((Ak)∗A2), N(A∗)⊆N(X). | (3.2) |
Proof. (a). Since R(AD)=R(Ak), it is a consequence of [2,Corollary 3.2] by properities of Drazin and MP inverse.
(b). Since items (c) of Lemma 2.3, AⓌ,† satisfies XA=PR(Ak), N((Ak)∗A2). Additionally, we derive that N(A∗)=N(A†)⊆N(AⓌAA†)=N(X). Now it remains to prove that X is unique.
Assume that X1, X2 satisfy (3.2), then X1A=X2A, \ N(A∗)⊆N(X1) and N(A∗)⊆N(X2). Furthermore, we get that (X1−X2)A=0 and R(X∗i)⊆R(A) for i=1,2, which further imply that A∗(X∗1−X∗2)=0 and R(X∗1−X∗2)⊆R(A). Therefore we have that R(X∗1−X∗2)⊆R(A)∩N(A∗)=R(A)∩R(A)⊥={0}. Thus, X∗1=X∗2, i.e., X1=X2.
Remark 3.8. In Theorem 3.7, R(X)⊆R(Ak) in (3.1) can be replaced by R(X)=R(Ak). However, if we replace N(A∗)⊆N(X) with N(A∗)=N(X) in (3.2), item (b) of Theorem 3.7 does not hold.
Characterizations of some generalized inverses by using its block matrices have been investigated in [13,14,15,16,17]. In [18,Theorem 3.2], the authors presented a characterization for the WC inverse using its block matrices. Next we will give another proof of it by using characterization of projection operator.
Theorem 3.9. Let A∈Cn×nk and r(Ak)=t. Then there exist a unique matrix P such that
P2=P, PAk=0, (Ak)∗A2P=0, r(P)=n−t, | (3.3) |
a unique matrix Q such that
Q2=Q, QAk=0, (Ak)∗A2A†Q=0, r(Q)=n−t, | (3.4) |
and a unique matrix X such that
r([AI−QI−PX])=r(A). | (3.5) |
Furthermore, X is the WC inverse AⓌ,† of A and
P=PN((Ak)∗A2), R(Ak), Q=PN((Ak)∗A2A†), R(Ak). | (3.6) |
Proof. It is not difficult to prove that
the condition (3.3) hold⟺(I−P)2=I−P, (I−P)Ak=Ak,(Ak)∗A2(I−P)=(Ak)∗A2, r(P)=n−t⟺I−P=PR(Ak), N((Ak)∗A2)⟺P=PN((Ak)∗A2), R(Ak). |
Similarly, we can show that (3.4) have the unique solution Q=PN((Ak)∗A2A†), R(Ak).
Furthermore, comparing (3.6) and items (b) and (c) of Lemma 2.3 immediately leads to the conclusion that
r([AI−QI−PX])=r([AAAⓌ,†AⓌ,†AX])=r(A)+r(X−AⓌ,†). |
By (3.5), we obtain that X=AⓌ,†.
In [19], Drazin introduced the (b,c)-inverse in semigroup. In [20], Benítez et al. investigated the (B,C)-inverse of A∈Cm×n, as the unique matrix X∈Cn×m satisfying [20]:
CAX=C, XAB=B, R(X)=R(B), N(X)=N(C), |
where B,C∈Cn×m. In the next result, we will show that the WC inverse is a special (B,C)-inverse.
Theorem 4.1. Let A∈Cn×nk. Then
AⓌ,†=A(Ak,(Ak)∗A2A†). |
Proof. According to Lemma 2.3, we get that
R(AⓌ,†)=R(Ak), N(AⓌ,†)=N((Ak)∗A2A†)). |
Observe that AⓌ,†AAk=AⓌAk+1=Ak and (Ak)∗A2A†AAⓌ,†=(Ak)∗A2AⓌAA†=(Ak)∗A2A†. Thus, we obtain AⓌ,†=A(Ak,(Ak)∗A2A†).
In [21], the authors introduced the Bott-Duffin inverse of A∈Cn×n when APL+PL⊥ is nonsingular, i.e., A(−1)L=PL(APL+PL⊥)−1=PL(APL+I−PL)−1. In [22], the authors showed the weak group inverse by a special Bott-Duffin inverse. Next we will show that the WC inverse of A is indeed the Bott-Duffin inverse of A2 with respect to R(Ak).
Theorem 4.2. Let A∈Cn×nk be given by (2.1). Then
AⓌ,†=(A2)(−1)(R(Ak))APA=(PAkA2PAk)†APA. | (4.1) |
Proof. It follows from (2.3) and (2.2) that
PA=U[It00NN†]U∗, | (4.2) |
PAk=U[It000]U∗. | (4.3) |
We now obtain that
(A2)(−1)(R(Ak))APA=PAk(A2PAk+I−PAk)−1APA=U[It000][T200In−t]−1[TS0N][It00NN†]U∗=U[T−1T−2SNN†00]U∗=AⓌ,†. |
Similarly, by a direct calculation, we can derive that AⓌ,†=(PAkA2PAk)†APA.
Working with the fact that P=PN((Ak)∗A2), R(Ak) and Q=PN((Ak)∗A2A†), R(Ak) in Theorem 3.9, we will consider other representations of AⓌ,† in the next theorem.
Theorem 4.3. Let A∈Cn×nk and P=PN((Ak)∗A2), R(Ak), Q=PN((Ak)∗A2A†), R(Ak). Then for any a,b≠0, we have
AⓌ,†=(A+aP)−1(I−Q)=(I−P)(A+bQ)−1. | (4.4) |
Proof. From items (b) and (c) of Lemma 2.3, it is not difficult to conclude that
(A+aP)AⓌ,†=I−Q. |
Now we only need to show the invertibility of A+aP. Assume that α=U[α1α2]∈Cn such that (A+aP)α=0, i.e., Aα=−aPα, where α1∈Cp, α2∈Cn−p. Now it follows from condition (c) of Lemma 2.3 and (3.6) that
[TS0N][α1α2]=−a[0−T−1S−T−2SN0I][α1α2], |
implying α1=0 and α2=0 since a≠0, N is nilpotent and T is non singular. Thus A+aP is nonsingular.
Analogously, we can prove that A+bQ is invertible and AⓌ,†=(I−P)(A+bQ)−1.
The limit expressions for some generalized inverses of matrices have been given in [14,15,16,17,23,24]. Similarly, the WC inverse can also be characterized as limit value as shown in the next result:
Theorem 4.4. Let A∈Cn×nk. Then:
(a) AⓌ,†=limλ→0Ak(λIn+(Ak)∗Ak+2)−1(Ak)∗A2A∗(λIn+AA∗)−1;
(b) AⓌ,†=limλ→0Ak(Ak)∗A(λIn+Ak+1(Ak)∗A)−1AA∗(λIn+AA∗)−1;
(c) AⓌ,†=limλ→0(λIn+Ak(Ak)∗A2)−1Ak(Ak)∗A2A∗(λIn+AA∗)−1;
(d) AⓌ,†=limλ→0Ak(Ak)∗A2A∗(λIn+AA∗)−1(λIn+Ak+1(Ak)∗A2A∗(λIn+AA∗)−1)−1.
Proof. According to condition (a) of Lemma 2.3, it is not hard to show that
AⓌ,†=A(2)R(Ak), N((Ak)∗A2A†)=A(2)R(Ak(Ak)∗A2A†), N(Ak(Ak)∗A2A†). |
Thus, by [25,Theorem 2.1], we have the following results:
(a) Let X=Ak, Y=(Ak)∗A2A† and by A†=limλ→0A∗(λIn+AA∗)−1. We have
AⓌ,†=limλ→0Ak(λIn+(Ak)∗Ak+2)−1(Ak)∗A2A∗(λIn+AA∗)−1. |
(b) Let X=Ak(Ak)∗A, Y=AA† and by A†=limλ→0A∗(λIn+AA∗)−1. We have
AⓌ,†=limλ→0Ak(Ak)∗A(λIn+Ak+1(Ak)∗A)−1AA∗(λIn+AA∗)−1. |
(c) Let X=In, Y=Ak(Ak)∗A2A† and by A†=limλ→0A∗(λIn+AA∗)−1. We have
AⓌ,†=limλ→0(λIn+Ak(Ak)∗A2)−1Ak(Ak)∗A2A∗(λIn+AA∗)−1. |
(d) Let X=Ak(Ak)∗A2A†, Y=In and by A†=limλ→0A∗(λIn+AA∗)−1. We have
AⓌ,†=limλ→0Ak(Ak)∗A2A∗(λIn+AA∗)−1(λIn+Ak+1(Ak)∗A2A∗(λIn+AA∗)−1)−1. |
We end up this section with three examples of computing the WC inverse of a matrix using three different expressions in Theorems 4.2–4.4.
Example 4.5. Let
A=[2210342000010000]. | (4.5) |
Then Ind(A)=2 and the weak core inverse of A is
AⓌ,†=A2(A4)†A2A†=[2−1−1/20−3/211/2000000000]. |
Firstly, using the expression (4.1) to compute the WC inverse of A. Then
(A2PA2+I−PA2)−1=[11/2−300−9/25/20000100001] and (PA2A2PA2)†=[11/2−300−9/25/20000000000]. |
After simplification, it follows that (A2)(−1)(R(A2))APA=AⓌ,† and (PA2A2PA2)†APA=AⓌ,†.
Secondly, using the expression (4.4), we obtain
(A−6P)−1=[2−1−1/2−19/12−3/217/1297/7200−1/6−1/36000−1/6] and (A+15Q)−1=[2−1−1/25/2−3/21−210005−250005]. |
Therefore, it can be directly verified (A−6P)−1(I−Q)=AⓌ,†, (I−P)(A+15Q)−1=AⓌ,†.
Finally, using the limit expressions of item (a) in Theorem 4.4.
Let B=A2(λIn+(A2)∗A4)−1(A2)∗A2A∗(λIn+AA∗)−1, then
B=A2(λIn+(A2)∗A4)−1(A2)∗A2A∗(λIn+AA∗)−1=[2(30321λ2+5361λ+580)/λ12(54629λ2+6699λ−290)/λ1(1305λ−58)/λ20(110503λ2+19405λ−870)/λ14(49773λ2+6090λ+145)/λ1(2378λ+58)/λ2000000000]. |
where λ1=λ4+38734λ3+1470569λ2+197888λ+580, λ2=λ3+38697λ2+38812λ+116.
After simplification, it follows that
limλ→0B=limλ→0A2(λIn+(A2)∗A4)−1(A2)∗A2A∗(λIn+AA∗)−1=AⓌ,†. |
The other cases in Theorem 4.4 can be similarly verified.
In this section, we discuss some properties of the WC inverse and consider the connection between the WC inverse and other known classes of matrices.
Lemma 5.1. Let A∈Cn×nk be given by (2.1). Then:
(a) A∈CEPn⇔S=0 and N=0;
(b) A∈CPn⇔T=It and N=0;
(c) A∈COPn⇔T=It, S=0 and N=0.
Proof. (a) The proof can be easily verified from (2.1) and (2.3).
(b) By (2.1), we obtain that A∈CPn is equivalent with
U[T2TS+SN0N2]U∗=U[TS0N]U∗, |
which is further equivalent with T2=T, TS+SN=S and N2=N. Hence, by non singularity of T and Nk=0, we can conclude that A∈CPn if and only if T=It and N=0.
(c) Since COPn⊆CPn, it is a direct consequence of item (b) and (2.1).
Theorem 5.2. Let A∈Cn×nk be given by (2.1). The following statements hold:
(a) AⓌ,†=0⇔A is nilpotent $;
(b) AⓌ,†=A⇔A∈CEPn and A3=A;
(c) AⓌ,†=A∗⇔A∈CEPn and AA∗=PAk;
(d) AⓌ,†=PA⇔A∈CPn;
(e) AⓌ,†=PA∗⇔A∈COPn.
Proof. (a) By (2.1) and (2.9), we directly get that
AⓌ,†=0⟺r(Ak)=t=0⟺A is nilpotent. |
(b) It follows (2.1), (2.9) and (a) of lemma 5.1 that
AⓌ,†=A⟺U[T−1T−2SNN†00]U∗=U[TS0N]U∗⟺S=0 , N=0 and T3=T⟺A∈CEPn and A3=A. |
(c) By (2.1), (2.9) and (a) of lemma 5.1, we have that
AⓌ,†=A∗⟺U[T−1T−2SNN†00]U∗=U[T∗0S∗N∗]U∗⟺S=0 , N=0 and TT∗=It⟺A∈CEPn and AA∗=PAk. |
(d) From (2.9), (4.2) and (b) of lemma 5.1, we obtain that
AⓌ,†=PA⟺U[T−1T−2SNN†00]U∗=U[It00NN†]U∗⟺T=It, N=0⟺A∈CPn. |
(e) It follows from (2.1) and (2.3) that
PA∗=A†A=U[T∗△TT∗△S(In−t−N†N)(In−t−N†N)S∗△TN†N+(In−t−N†N)S∗△S(In−t−N†N)]U∗. | (5.1) |
By (2.9) and (5.1), we now get that AⓌ,†=PA∗ is equivalent with
U[T−1T−2SNN†00]U∗=U[T∗△TT∗△S(In−t−N†N)(In−t−N†N)S∗△TN†N+(In−t−N†N)S∗△S(In−t−N†N)]U∗, |
which is further equivalent with T−1=T∗△T, (In−t−N†N)S∗△T=0 and N†N+(In−t−N†N)S∗△S(In−t−N†N)=0. Hence, by nonsingularity of △T and (c) of lemma 5.1, we can conclude that AⓌ,†=PA∗ if and only if A∈COPn.
From Lemma 2.3, we know that both AAⓌ,† and AⓌ,†A are oblique projectors. The next theorem will further discuss other characteriations for AAⓌ,† and AⓌ,†A.
Theorem 5.3. Let A∈Cn×nk be given by (2.1). The following statements hold:
(a) AAⓌ,†=PA⇔A∈CCMn; (b) AAⓌ,†=PA∗⇔A∈CEPn;
(c) AⓌ,†A=PA⇔A∈CEPn; (d) AⓌ,†A=PA∗⇔A∈CEPn.
Proof. It follows from (2.1) and (2.9) that
AAⓌ,†=U[ItT−1SNN†00]U∗, | (5.2) |
AⓌ,†A=U[ItT−1S+T−2SN00]U∗. | (5.3) |
(a) By (4.2) and (5.2), the result can be directly verified.
(b) By (5.1) and (5.2), we can show that AAⓌ,†=PA∗ if and only if (In−t−N†N)S∗△T=0 and N†N+(In−t−N†N)S∗△S(In−t−N†N)=0, which is further equivalent with S=0 and N=0, i.e., A∈CEPn.
(c) It follows from (4.2) and (5.3) that AⓌ,†A=PA is equivalent with A∈CEPn.
(d) From (5.1) and (5.3), it is similar to the proof of (b).
Recall from [6] that the core inverse is necessarily EP. The next Theorem shows that this is not the case with the WC inverse.
Theorem 5.4. Let A∈Cn×nk be given by (2.1) and t∈Z+. The following statements are equivalent:
(a) AⓌ,†∈CEPn; (b) SN=0;
(c) ; (d)
;
(e) AⓌ,†At=AtAⓌ.
Proof. (a)⇔(b). Since AⓌ,†∈CEPn is equivalent with R(AⓌ,†)=R((AⓌ,†)∗). Using (2.9), we have that AⓌ,†∈CEPn if and only if SN=0.
(c)⇔(b). By (2.5) and (5.3), it can be directly verified that if and only if SN=0.
(d)⇔(b). By (2.5) and (2.9), it follows that
![]() |
(e)⇔(b). From (2.8) and (2.9), it follows that
AⓌ,†At=AtAⓌ⟺U[Tt−1Tt−2S+T−2TtN00]U∗=U[Tt−1Tt−2S00]U∗⟺T−2TtN=0⟺SN=0. |
where Tt=t−1∑j=0TjSNt−1−j.
In [26], the authors introduced that a matrix A to be a weak group matrix if A∈CWGn, which is equivalent with SN=0. Therefore, we have that following remark:
Remark 5.5. It is worth noting that conditions (a), (c)–(e) in Theorem 5.4 are equivalent with A∈CWGn.
The next theorems provide some equivalent conditions for AⓌ,†∈CPn and AⓌ,†∈COPn.
Theorem 5.6. Let A∈Cn×nk be given by (2.1). The following statements are equivalent:
(a) AⓌ,†∈CPn; (b) T=It;
(c) AAⓌ,†=AⓌ,†; (d) AⓌ,†Ak=Ak;
(e) Ak(AⓌ,†)k=AⓌ,†; (f) A(AⓌ,†)k=(AⓌ,†)k.
Proof. (a)⇔(b). From (2.9), it is not hard to prove that AⓌ,†∈CPn is equivalent with T=It.
(c)⇔(b). From (2.9) and (5.3), it follows that AAⓌ,†=AⓌ,† if and only if T=It.
(d)⇔(b). By (2.2) and (2.9), it is easy to verify that AⓌ,†Ak=Ak if and only if T=It.
The proofs (e)⇔(b) and (f)⇔(b) are similar to the proof (d)⇔(b).
Theorem 5.7. Let A∈Cn×nk be given by (2.1). The following statements are equivalent:
(a) AⓌ,†∈COPn; (b) T=It and SN=0;
(c) AAⓌ,†=(AⓌ,†)∗; (d) AⓌ,†A=AⓌ;
(e) (AⓌ,†)kAk=AⓌ; (f) (AⓌ,†)kA=(AⓌ)k.
Proof. (a)⇔(b). From (2.9) and Theorem 5.6, we can show that AⓌ,†∈COPn is equivalent with T=It and SN=0.
(c)⇔(b). By (2.9) and (5.3), it follows from AAⓌ,†=(AⓌ,†)∗ that
[ItT−1SNN†00]=[(T−1)∗0(T−2SNN†)∗0]. |
Hence, we get that AAⓌ,†=(AⓌ,†)∗ is equivalent with T=It and SN=0.
The proofs of (d)⇔(b), (e)⇔(b) and (f)⇔(b) are similar to the proof of (c)⇔(b).
Corollary 5.8. Let A∈Cn×nk. Then A∈COPn if and only if AⓌ,†∈CPn∩CEPn.
Proof. It is a direct consequence from Theorem 5.4 and Theorem 5.6.
Working with Theorem 5.6 and Theorem 5.7, we have the following corollary.
Corollary 5.9. Let A∈Cn×nk and for any l∈N, l≥k. The following statements statements hold:
(a)A∈CPn⇔AⓌ,†∈CPn and Al=A;
(b)A∈COPn⇔AⓌ,†∈CPn and Al=A∗.
Proof. (a) The result can be easily derived by lemma 5.1 and Theorem 5.6,
(b) From Lemma 5.1 and Theorem 5.7, we can show that (b) holds.
Ferreyra et al.[1] introduced the weak core matrix. The set of all n×n weak core matrices is denoted by CWCn, that is:
CWCn={A∣A∈Cn×n,AⓌ,†=AD,†}. |
In this section, we discuss some equivalent conditions satisfied by a matrix A such that A∈CWCn using the core-EP decomposition. For convenience, we introduce a necessary lemma.
Lemma 6.1. [1] Let A∈Cn×nk be given by (2.1). Then the following statements are equivalent:
(a) A∈CWCn;
(b) SN2=0;
(c) AⓌ=AD.
Theorem 6.2. Let A∈Cn×nk be given by (2.1) and t∈Z+. The following statements are equivalent:
(a) A∈CWCn;
(b) AⓌA=ADA;
(c) AtAⓌA=AtADA;
(d) AtAⓌ,†=AtAD,†;
(e) AkAⓌ,†=AkA†;
(f) AkAⓌA=Ak.
Proof. (a)⇒(b). It is a direct consequence from condition (c) of Lemma 6.1.
(b)⇒(c). Evident.
(c)⇒(a). By condition, it follows from (2.4) and (2.8) that
[TtTt−1S+Tt−2SN00]=[TtTt−1S+Tt−k−1˜TN00], |
which implies Tt−k−1(Tk−2SN2+⋯+TSNk−1)=0. We now obtain that SN2=0 since T is invertible. By Lemma 6.1, we obtain that A∈CWCn.
(a)⇔(d). It follows form (2.6), (2.9) and Lemma 6.1 that
AtAⓌ,†=AtAD,†⟺U[Tt−1Tt−2SNN†00]U∗=U[Tt−1Tt−k−1˜TNN†00]U∗⟺Tt−2SNN†=Tt−k−1˜TNN†⟺SN2=0⟺A∈CWCn. |
(a)⇒(e). From the definition of the weak core matrix, we have that AkAⓌ,†=AkAD,†=AkA†.
(e)⇒(f). Evident.
(f)⇒(a). If AkAⓌA=Ak, by (2.2) and (2.8), we can conclude that SN2. Hence item (a) holds.
Corollary 6.3. Let A∈Cn×nk be given by (2.1). Then A∈CWCn if and only if Ak=Ck, where C=AAⓌA.
Proof. Since AkAⓌA=Ck, the result is a direct consequence of item (f) of Theorem 6.2.
Theorem 6.4. Let A∈Cn×nk be given by (2.1) and for some t∈Z+. The following statements are equivalent:
(a) A∈CWCn;
(b) A(AⓌ)tA=(AⓌ)tA2;
(c) (AⓌ)tA=(AⓌ)t+1A2;
(d) A(AⓌ)tA commutes with (AⓌ)tA2;
(e) (AⓌ)tA commutes with (AⓌ)t+1A2.
Proof. By (2.1) and (2.8), we get that
A(AⓌ)tA=U[T−t+2T−t+1S+T−tSN00]U∗, | (6.1) |
(AⓌ)tA2=U[T−t+2T−t+1S+T−tSN+T−t−1SN200]U∗. | (6.2) |
(a)⇔(b). By (6.1), (6.2) and Lemma 6.1, we get that A(AⓌ)tA=(AⓌ)tA2 if and only if A∈CWCn.
(a)⇔(c). Similar to the part (a)⇔(b).
(a)⇔(d). It follows from (6.1), (6.2) and Lemma 6.1 that
A(AⓌ)tA(AⓌ)tA2−(AⓌ)tA2A(AⓌ)tA=U[0T−2t+1SN200]U∗, |
which implies that A(AⓌ)tA commutes with (AⓌ)tA2 if and only if A∈CWCn.
(a)⇔(e). It is analogous to that of the part (a)⇔(d).
Corollary 6.5. Let A∈Cn×nk and t∈Z+. The following statements are equivalent:
(a) A∈CWCn;
Proof. Since , Corollary 6.5 can be directly verified.
In this paper, new characterizations and properties of the WC inverse are derived by using range, null space, matrix equations, respectively. Several expressions of the WC inverse are also given. Finally, we show various characterizations of the weak core matrix.
According to the current research background, more characterizations and applications for the WC inverse are worthy of further discussion which as follows:
1) Characterizing the WC inverse by maximal classes of matrices, full rank decomposition, integral expressions and so on;
2) New iterative algorithms and splitting methods for computing the WC inverse;
3) Using the WC inverse to solve appropriately constrained systems of linear equations;
4) Investigating the WC inverse of tensors.
This work was supported by the Natural Science Foundation of China under Grants 11961076. The authors are thankful to two anonymous referees for their careful reading, detailed corrections and pertinent suggestions on the first version of the paper, which enhanced the presentation of the results distinctly.
All authors read and approved the final manuscript. The authors declare no conflict of interest.
[1] | M. Wazewska-Czyzewska, A. Lasota, Mathematical problems of the dynamics of red blood cells system, Annals of the Polish Mathematical Society, Seines III, Applied Mathematics, 17 (1988), 23–40. |
[2] |
L. Duan, L. Huang, Y. Chen, Global exponential stability of periodic solutions to a delay Lasota-Wazewska model with discontinuous harvesting, P. Am. Math. Soc., 144 (2015), 561–573. doi: 10.1090/proc12714
![]() |
[3] |
J. Shao, Pseudo almost periodic solutions for a Lasota-Wazewska model with an oscillating death rate, Appl. Math. Lett., 43 (2015), 90–95. doi: 10.1016/j.aml.2014.12.006
![]() |
[4] |
Z. Yao, New results on existence and exponential stability of the unique positive almost periodic solution for Hematopoiesis model, Appl. Math. Model., 39 (2015), 7113–7123. doi: 10.1016/j.apm.2015.03.003
![]() |
[5] |
Q. Su, S. Ruan, Existence of periodic solutions in abstract semilinear equations and applications to biological models, J. Differ. Equations, 269 (2020), 11020–11061. doi: 10.1016/j.jde.2020.07.014
![]() |
[6] |
Z. Huang, S. Gong, L. Wang, Positive almost periodic solution for a class of Lasota-Wazewska model with multiple time-varying delays, Comput. Math. Appl., 61 (2011), 755–760. doi: 10.1016/j.camwa.2010.12.019
![]() |
[7] | S. Abbas, S. Dhama, M. Pinto, D. Sepúlveda, Pseudo compact almost automorphic solutions for a family of delayed population model of Nicholson type, J. Math. Anal. Appl., 495 (2020). |
[8] |
H. El-Morshedy, A. Ruiz-Herrera, Criteria of global attraction in systems of delay differential equations with mixed monotonicity, J. Differ. Equations, 268 (2020), 5945–5968. doi: 10.1016/j.jde.2019.11.016
![]() |
[9] |
S. Saker, Qualitative analysis of discrete nonlinear delay survival red blood cells model, Nonlinear Anal-Real, 9 (2008), 471–489. doi: 10.1016/j.nonrwa.2006.11.013
![]() |
[10] |
D. Fan, J. Wei, Bifurcation analysis of discrete survival red blood cells model, Commun. Nonlinear Sci., 14 (2009), 3358–3368. doi: 10.1016/j.cnsns.2009.01.015
![]() |
[11] |
S. Glasgow, Z. Perkins, N. Tai, K. Brohi, C. Vasilakis, Development of a discrete event simulation model for evaluating strategies of red blood cell provision following mass casualty events, Eur. J. Oper. Res., 270 (2018), 362–374. doi: 10.1016/j.ejor.2018.03.008
![]() |
[12] | A. Nicholson, The balance of animal population, J. Anim. Ecol., 2 (1993), 132–178. |
[13] |
K. Gopalsamy, S. Trofimchuk, Almost periodic solutions of Lasota-Wazewska type delay differential equations, J. Math. Anal. Appl., 237 (1999), 106–127. doi: 10.1006/jmaa.1999.6466
![]() |
[14] |
J. Li, Z. Wang, Existence and global attractivity of positive periodic solutions of a survival model of red blood cells, Comput. Math. Appl., 50 (2005), 41–47. doi: 10.1016/j.camwa.2005.03.003
![]() |
[15] | D. Jiang, J. Wei, Existence of positive periodic solutions for nonautonomous delay differential equations, Chinese Annals of Mathematics, Series A, 20 (1999), 715–720. |
[16] | S. Saker, S. Agarwal, Oscillation and global attractivity of a periodic survival red blood cells model, Dynamics of Continuous, Discrete and Impulsive Systems Series A: Mathematical Analysis, 12 (2005), 429–440. |
[17] |
G. Liu, A. Zhao, J. Yan, Existence and global attractivity of unique positive periodic solution for a Lasota-Wazewska model, Nonlinear Anal., 64 (2006), 1737–1746. doi: 10.1016/j.na.2005.07.022
![]() |
[18] | R. Games, J. Mawhin, Coincidence degree and nonlinear differential equations, Berlin: Springer, 1997. |
[19] |
C. Wang, R. Agarwal, Almost periodic solution for a new type of neutral impulsive stochastic Lasota-Wazewska timescale model, Appl. Math. Lett., 70 (2017), 58–65. doi: 10.1016/j.aml.2017.03.009
![]() |
[20] |
G. Stamov, I. Stamova, J. Cao, Uncertain impulsive functional differential systems of fractional order and almost periodicity, J. Franklin I., 355 (2018), 5310–5323. doi: 10.1016/j.jfranklin.2018.05.021
![]() |
[21] |
G. Stamov, On the existence of almost periodic solutions for the impulsive Lasota-Wazewska model, Appl. Math. Lett., 22 (2009), 516–520. doi: 10.1016/j.aml.2008.07.002
![]() |
[22] | Z. Yao, Existence and exponential stability of the unique positive almost periodic solution for impulsive Nicholson's blowflies model with linear harvesting term, J. Math. Anal. Appl., 39 (2015), 7124–7133. |
[23] |
J. Yan, Existence and global attractivity of positive periodic solution for an impulsive Lasota-Wazewska model, J. Math. Anal. Appl., 279 (2003), 111–120. doi: 10.1016/S0022-247X(02)00613-3
![]() |
[24] |
X. Liu, Y. Takeuchi, Periodicity and global dynamics of an impulsive delay Lasota-Wazewska model, J. Math. Anal. Appl., 327 (2007), 326–341. doi: 10.1016/j.jmaa.2006.04.026
![]() |
[25] |
X. Yang, X. Li, Q. Xi, P. Duan, Review of stability and stabilization for impulsive delayed systems, Math. Biosci. Eng., 15 (2018), 1495–1515. doi: 10.3934/mbe.2018069
![]() |
[26] | X. Li, X. Yang, T. Huang, Persistence of delayed cooperative models: Impulsive control method, Appl. Math. Comput., 342 (2019), 130–146. |
[27] |
W. Chen, Z. Ruan, W. Zheng, Stability and L2-gain analysis for impulsive delay systems: An impulse-time-dependent discretized Lyapunov functional method, Automatica, 86 (2017), 129–137. doi: 10.1016/j.automatica.2017.08.023
![]() |
[28] |
X. Yang, J. Lam, D. Ho, Z. Feng, Fixed-time synchronization of complex networks with impulsive effects via nonchattering control, IEEE T. Automat. Contr., 62 (2017), 5511–5521. doi: 10.1109/TAC.2017.2691303
![]() |
[29] |
X. Liu, K. Zhang, Synchronization of linear dynamical networks on time scales: Pinning control via delayed impulses, Automatica, 72 (2016), 147–152. doi: 10.1016/j.automatica.2016.06.001
![]() |
[30] | V. Lakshmikantham, D. Bainov, P. Simeonov, Theory of Impulsive Differential Equations, Singapore: World Scientific, 1989. |
[31] |
Z. Yang, D. Xu, Existence and exponential stability of periodic solution for impulsive delay differential equations and applications, Nonlinear Anal., 64 (2006), 130–145. doi: 10.1016/j.na.2005.06.014
![]() |
[32] | X. Li, Existence and global exponential stability of periodic solution for impulsive Cohen-Grossberg-type BAM neural networks with continuously distributed delays, Appl. Math. Comput., 215 (2009), 292–307. |
1. | Wende Li, Jianlong Chen, Yukun Zhou, CHARACTERIZATIONS AND PROPERTIES OF WEAK CORE INVERSES IN RINGS WITH INVOLUTION, 2024, 54, 0035-7596, 10.1216/rmj.2024.54.793 | |
2. | Jiaxuan Yao, Hongwei Jin, Xiaoji Liu, The weak group-star matrix, 2023, 37, 0354-5180, 7919, 10.2298/FIL2323919Y | |
3. | Dijana Mosić, Janko Marovt, Weighted MP weak group inverse, 2024, 32, 1844-0835, 221, 10.2478/auom-2024-0012 |