Skip to main content

BRIEF RESEARCH REPORT article

Front. Phys., 10 July 2020
Sec. Statistical and Computational Physics
Volume 8 - 2020 | https://doi.org/10.3389/fphy.2020.00239

A Relation Between Moore-Penrose Inverses of Hermitian Matrices and Its Application in Electrical Networks

Yujun Yang1* Dayong Wang2 Douglas J. Klein3
  • 1School of Mathematics and Information Sciences, Yantai University, Yantai, China
  • 2Business School, Hohai University, Nanjing, China
  • 3Department of Marine Sciences, Texas A&M University at Galveston, Galveston, TX, United States

A novel relation between the Moore-Penrose inverses of two nullity-1 n × n Hermitian matrices which share a common null eigenvector is established, and its application in electrical networks is illustrated by applying the result to Laplacian matrices of graphs.

1. Introduction

The Hermitian matrices are an important class of matrices arising in many contexts. A complex squared matrix is called a Hermitian matrix if it is equal to its conjugate transpose, in other words, for all i and j, its (i, j)-th element (i.e., the element in the i-th row and j-th column) is equal to the complex conjugate of its (j, i)-th element. It is widely known that all the eigenvalues of a Hermitian matrix are real. In addition, it is easily seen that Hermitian matrices contain real symmetric matrices as special cases.

Let M be an n × m matrix. An m × n matrix X is called the Moore-Penrose (generalized) inverse of M, if X satisfies the following equations:

MXM=M,XMX=X,(MX)H=MX,(XM)H=XM,

where XH represents the conjugate transpose of the matrix M. It is well-known [1] that for any matrix M, the Moore-Penrose inverse of M does exist and is unique. For this reason, the unique Moore-Penrose inverse of M is denoted by M+.

We proceed to introduce a special class of Hermitian matrices – the Laplacian matrices of graphs, which play a fundamental role in graph theory and electrical network theory. Let G = (V, E) be a connected weighted graph of order n. For each edge e of G, we assign a positive real number we to e, and we call we the weight of e. Then the adjacency matrix of G, denoted by A, is a n × n matrix such that the (i, j)-th element of A is equal to the weight of the edge ij if i and j are connected by an edge and 0 otherwise. Suppose that D is the n × n diagonal matrix such that the i-th diagonal element is equal to the sum of the weights of the edges incident to i. Then the Laplacian matrix L of G is defined as L = DA. It is easily seen that the Laplacian matrix is real and symmetric. Thus, the Laplacian matrix is a Hermitian matrix. According to the definition of the Laplacian matrix, we readily seen that the Laplacian matrix is singular and not invertible.

It is natural to consider a weighted graph G as a (resistive) electrical network N by viewing each edge e as a resistor such that the conductance of the resistor is we, where we is the weight on e. In this guise, the resistance distance [2] between any two vertices i and j of G, denoted by Ω(i, j), is defined as the net effective resistance between corresponding nodes i and j in N. It should be mentioned that resistance distance, as an important component of circuit theory, has been studied for a long time, dating back to the classical work of Kirchhoff in 1847. It is amazing that the resistance distance turns out to have many purely mathematical interpretations, although it comes from physics and engineering, among which a fundamental one is the classical result which is given via the Moore-Penrose inverse of the Laplacian matrix [2]:

Ω(i,j)=Lii+-2Lij++Ljj+,    (1.1)

where Lij+ denote the (i, j)-th element of L+. Since the identification of resistance distance as a novel distance function on graphs, the resistance distance has been extensively studied in the literature of mathematics, physics, and chemistry. For more information on resistance distances, we refer the readers to recent papers [313] and references therein.

In this paper, a relation between the Moore-Penrose inverses of two nullity-1 n × n Hermitian matrices which share a common null eigenvector is established. Then its application in electrical networks is illustrated by applying the result to Laplacian matrices of graphs.

2. A Relation Between Moore-Penrose Inverses of Two Hermitian Matrices

All the matrices considered in this section are square matrices of order n. For an invertible matrix M, we use M−1 to denote the inverse of M. Let I and O denote the identity matrix and zero matrix, respectively. This section is devoted to establish a relation between Moore-Penrose inverses of two Hermitian matrices of nullity-1 which share a common null eigenvector. To this end, we first give some properties on nullity-1 Hermitian matrices, which will be used in the later.

Lemma 2.1. Let M be a nullity-1 Hermitian n × n matrix. Suppose that 0 = λ1, λ2, …, λn are eigenvalues of M with corresponding orthonormal eigenvectors u1, u2, …, un. Then

M+=(M+u1u1H)-1-u1u1H.    (2.1)
MM+=M+M=I-u1u1H.    (2.2)
u1u1HM+=O.    (2.3)

Proof: Let U = (u1, u2, …, un) and Λ = diag{0, λ2, …, λn}. Then

M=UΛUH.

As u1u1H=Udiag{1,0,,0}UH, it follows that

M+u1u1H=UΛUH+Udiag{1,0,,0}UH       =Udiag{1,λ2,,λn}UH.

Thus M+u1u1H is invertible with

(M+u1u1H)-1=Udiag {1,1λ2,,1λn}UH.

Consequently,

(M+u1u1H)-1-u1u1H=Udiag {0,1λ2,,1λn}UH.

Thus it is easily verified by the definition of the Moore-Penrose inverse that

M+=(M+u1u1H)-1-u1u1H.

To prove Equation (2.2), note first that

MM+=UΛUHUΛ0HH=UΛΛ0HH andM+M=UΛ0UHUΛUH=UΛ0ΛUH,

where Λ0=Udiag{1,1λ2,,1λn}UH. Then, note that

ΛΛ0=Λ0Λ=diag{0,1,,1}.

Thus we have

MM+=M+M=U(diag{0,1,,1})UH    =U(I-diag{1,0,,0})UH    =UUH-Udiag{1,0,,0}UH=I-u1u1H.

For Equation (2.3), by the above arguments we have

u1u1HM+=(Udiag{1,0,,0}UH)(Udiag{0,1λ2,,1λn}UH)         =Udiag{1,0,,0}diag{0,1λ2,,1λn}UH=O,

as required.

According to the properties given in Lemma 2.1, a relation between Moore-Penrose inverses of two Hermitian matrices of nullity-1 which share a common null eigenvector could be established, as given in the following result.

Theorem 2.2. Let M and M′ be two nullity-1 Hermitian n × n matrices which share a common null eigenvector. Then

(M)+=M+[I+(M-M)M+]-1.    (2.4)

Proof. For the sake of simplicity, set Δ: = M′ − M and ∇: = (M′)+M+. Then

M(M)+=(M+Δ)(M++)=MM++M+ΔM++Δ.    (2.5)

Let u1 be the common null eigenvector shared by M and M′. Then by Lemma 2.1, we know that

M(M)+=MM+=I-u1u1H.

Thus, Equation (2.5) gives

M+ΔM++Δ=O,

that is,

M=-ΔM+.

Left-multiply both sides of the above equation by (M′)+, we get

(M)+M=-(M)+ΔM+.

Bearing in mind that (M)+M=I-u1u1H and that (M′)+ = M+ + ∇, we arrive at

(I-u1u1H)=-(M++)ΔM+,

that is,

-u1u1H=-(M++)ΔM+.    (2.6)

Since it is shown in Lemma 2.1 that

u1u1H(M)+=u1u1HM+=O,

we have

u1u1H=u1u1H[(M)+-M+]=O.

Hence Equation (2.6) becomes

=-M+ΔM+-ΔM+,

or equivalently,

(I+ΔM+)=-M+ΔM+.

So if I + ΔM+ is invertible, then by right-multiplying the above equation by (I + ΔM+)−1, we could obtain

=-M+ΔM+(I+ΔM+)-1,

which yields

(M)+=M++=M+-M+ΔM+(I+ΔM+)-1      =M+[I-ΔM+(I+ΔM+)-1]      =M+[I-(I+ΔM+)(I+ΔM+)-1+(I+ΔM+)-1]      =M+[I-I+(I+ΔM+)-1]      =M+(I+ΔM+)-1

It remains to verify that I + ΔM+ is invertible. As

M+=(M+u1u1H)-1-u1u1H,

it follows that

I+ΔM+=I+Δ[(M+u1u1H)-1-u1u1H]      =I+Δ(M+u1u1H)-1-Δu1u1H     =I+Δ(M+u1u1H)-1-(M-M)u1u1H     =I+Δ(M+u1u1H)-1-Mu1u1H+Mu1u1H.

Noticing that u1 is an 0-eigenvalue eigenvector of M and M′, it gives that

I+ΔM+=(M+u1u1H)(M+u1u1H)-1+Δ(M+u1u1H)-1     =(M+u1u1H+Δ)(M+u1u1H)-1     =(M+u1u1H)(M+u1u1H)-1.

As M+u1u1H is non-singular, by the same reason we know that M+u1u1H is non-singular, so that I + ΔM+ is invertible. The proof is complete.

Obviously, the Laplacian matrix is a Hermitian matrix. In addition, all the Laplacian matrices of connected graphs of the same order are nullity-1 and share the same eigenvector. Hence, Theorem 2.2 can be directly applied to Laplacian matrices. Let G and G′ be weighted connected graphs of order n. As a straightforward consequence of Theorem 2.2, we have

Corollary 2.3. Let G and G′ be connected weighted graphs of order n with Laplacian matrices L and L′, respectively. Then

(L)+=L+[I+(L-L)L+]-1.    (2.7)

3. An Application to Electrical Networks

The Laplacian matrix, also known as the Kirchhoff matrix, or admittance matrix, has wide applications in electrical networks. As introduced in the first section, the resistance distance could be computed in terms of the Moore-Penrose inverse of the Laplacian matrix. Actually, the computation of resistance distances is a classical problem in circuit theory and electrical network theory. Besides, this problem is relevant to a number of problems ranging from Lattice Green's functions, harmonic functions to random walks on graphs. For this reason, many researchers devote themselves to the computation of the resistance distance. With the development of more than 170 years, various formulae and techniques have been established, such as the traditional techniques like series and parallel circuits, Kirchhoff's laws and star-triangle transformation, as well as newly developed techniques like (algebraic, probabilistic, and combinatorial) formulae, local and global sum rules, recursion relations. In [14], a novel recursion formula for computing resistance distance is obtained. It turns out that resistance distances in some networks could be computed very easily by the recursion formula. In addition, the recursion formula extends the famous Rayleigh's monotonicity law by giving quantitative characterization to the law.

In this section, we use Corollary 2.3 to give a new proof to the recursion formula on resistance distances proposed in [14].

Theorem 3.1. [14] Let G and G′ be two weighted graphs which are the same except for the weights on an edge e = ij are we and we. For any two vertices p and q, denote the resistance distance between them in G and G′ by Ω(p, q) and Ω′(p, q), respectively. Then

Ω(p,q)=Ω(p,q)-δ·[Ω(p,i)+Ω(q,j)-Ω(p,j)-Ω(q,i)]24[1+δ·Ω(i,j)],    (3.1)

where δwe-we.

Proof. Denote the Laplacian matrices of G and G′ respectively by L and L′, and let e be the (column) vector of order n whose components are 0 except the i-th and j-th components are respectively 1 and −1. Then

L=L+δ·eeH.

By Corollary 2.3, we have

(L)+=L+[I+(L-L)L+]-1=L+(I+δ·eeHL+)-1.

To compute (L′)+, we first compute (I + δ · eeHL+)−1. Note that the elements of I + δ · eeHL+ are given by

[I+δ·eeHL+]kl={1,if k=li,j,δ·(Lil+-Ljl+),if k=i and li,1+δ·(Lii+-Lji+),if k=l=i,-δ·(Lil+-Ljl+),if k=j and lj,1-δ·(Lij+-Ljj+),if k=l=j,0otherwise.

Simple algebraic calculation leads to

det(I+δ·eeHL+)=1+δ·(Lii++Ljj+-2Lij+).

Then by the adjoint method, we could obtain the inverse of I + δ · eeHL+, whose elements are given by

[(I+δ·eeHL+)-1]kl=    {1,if k=li,j,-δ·(Lil+-Ljl+)1+δ·(Lii++Ljj+-2Lij+),if k=i and li,1-δ·(Lii+-Lji+)1+δ·(Lii++Ljj+-2Lij+),if k=l=i,δ·(Lil+-Ljl+)1+δ·(Lii++Ljj+-2Lij+),if k=j and lj,1+δ·(Lij+-Ljj+)1+δ·(Lii++Ljj+-2Lij+),if k=l=j0,otherwise.

Then, by algebraic calculation, we could obtain the product of L+ and (I + δ · eetL+)−1. Thus, (L′)+ is obtained, whose elements are given below. For 1 ≤ k, ln,

(L)kl+=Lkl+-δ·(Lki+-Lkj+)(Lil+-Ljl+)1+δ·(Lii++Ljj+-2Lij+).

Now we are ready to prove Equation (3.1) according to the formula given in Equation (1.1). By Equation (1.1), we have

Ω(p,q)=(L)pp++(L)qq+-2(L)pq+=Lpp++Lqq+-2Lpq+     -δ·(Lpi+-Lpj+)2+δ·(Lqi+-Lqj+)21+δ·(Lii++Ljj+-2Lij+)     -2δ·[(Lpi+-Lpj+)(Lqi+-Lqj+)]1+δ·(Lii++Ljj+-2Lij+)    =Lpp++Lqq+-2Lpq+-δ·[(Lpi+-Lpj+)2+(Lqi+-Lqj+)2]1+δ·(Lii++Ljj+-2Lij+)     -2δ·[(Lpi+-Lpj+)(Lqi+-Lqj+)]1+δ·(Lii++Ljj+-2Lij+)    =Lpp++Lqq+-2Lpq+-δ·(Lpi+-Lpj+-Lqi++Lqj+)21+δ·(Lii++Ljj+-2Lij+)    =Lpp++Lqq+-2Lpq+     -δ·[(Lpi+-Lpj+-Lqi++Lqj+)+12(Lpp+-Lpp++Lii+-Lii++Lqq+-Lqq++Ljj+-Ljj+)]1+δ·(Lii++Ljj+-2Lij+)2    =Lpp++Lqq+-2Lpq+     -δ·[(-12Lpp++Lpi+-12Lii+)+(12Lpp+-Lpj++12Ljj+)+(12Lqq+-Lqi++12Lii+)+(-12Lqq++Lqj+-12Ljj+)]1+δ·(Lii++Ljj+-Lij+-Lji+)2    =Ω(p,q)-δ·[-12Ω(p,i)+12Ω(p,j)+12Ω(q,i)-12Ω(q,j)]2[1+δ·Ω(i,j)]    =Ω(p,q)-δ·[Ω(p,i)+Ω(q,j)-Ω(p,j)-Ω(q,i)]24[1+δ·Ω(i,j)].

The proof is completed.

4. Conclusion

The Moore-Penrose inverse of the Hermitian matrix has various applications. In this paper, a relation between generalized inverses of two nullity-1 n × n Hermitian matrices which share a common null eigenvector is established, and a simple application in electrical networks is illustrated. Further applications of the relation needs to be revealed in the future.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

This research was funded by National Natural Science Foundation of China through grant number 116711347, and Natural Science Foundation of Shandong Province through grant number ZR2019YQ02.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors would like to thank the reviewers for their careful reading of the manuscript and valuable suggestions.

References

1. Ben-Israel A, Greville TN. Generalized Inverses: Theory and Application. New York, NY: Springer-Verlag (2003).

Google Scholar

2. Klein DJ, Randić M. Resistance distance. J Math Chem. (1993) 12:81–95. doi: 10.1007/BF01164627

CrossRef Full Text | Google Scholar

3. Jiang Z, Yan W. Some two-point resistances of the Sierpinski gasket network. J Stat Phys. (2018) 172:824–31. doi: 10.1007/s10955-018-2067-0

CrossRef Full Text | Google Scholar

4. Liu Q, Liu J, Wang S. Resistance distance and Kirchhoff index of the corona-vertex and the corona-edge of subdivision graph. IEEE Access. (2018) 6:55673–9. doi: 10.1109/ACCESS.2018.2871840

CrossRef Full Text | Google Scholar

5. Yang Y, Rosenfeld VR. Some spectral invariants of the neighborhood corona of graphs. Discrete Appl Math. (2018) 247:300–8. doi: 10.1016/j.dam.2018.03.052

CrossRef Full Text | Google Scholar

6. Barrett W, Evans EJ, Francis AE. Resistance distance in straight linear 2-trees. Discrete Appl Math. (2019) 258:13–34. doi: 10.1016/j.dam.2018.10.043

CrossRef Full Text | Google Scholar

7. Cao J, Liu J, Wang S. Resistance distances in corona and neighborhood corona networks based on Laplacian generalized inverse approach. J Alg Appl. (2019) 18:1950053. doi: 10.1142/S0219498819500531

CrossRef Full Text | Google Scholar

8. Carmona A, Encinas A. M., Mitjana M. Resistance distances in extended or contracted networks. Linear Alg Appl. (2019) 576:5–34. doi: 10.1016/j.laa.2018.01.026

CrossRef Full Text

9. Fan J, Zhu J, Tian L, Wang Q. Resistance distance in potting networks. Phys A. (2019) 540:123053. doi: 10.1016/j.physa.2019.123053

CrossRef Full Text | Google Scholar

10. Hong Y, Zhu Z, Luo A. Extremal graphs with diameter 2 for two indices on resistance-distance. Discrete Math. (2019) 342:487–97. doi: 10.1016/j.disc.2018.10.023

CrossRef Full Text | Google Scholar

11. Ye L, Yan W. Resistance between two vertices of almost complete bipartite graphs. Discrete Appl Math. (2019) 257:299–306. doi: 10.1016/j.dam.2018.08.030

CrossRef Full Text | Google Scholar

12. Thulasiraman K, Yadav M, Naik K. Network science meets circuit theory: resistance distance, Kirchhoff index, and Foster's theorems with generalizations and unification. IEEE Trans Circuits Syst. (2019) 66:1090–103. doi: 10.1109/TCSI.2018.2880601

CrossRef Full Text | Google Scholar

13. Yang Y, Klein DJ. Two-point resistances and random walks on stellated regular graphs. J Phys A Math Theor. (2019) 52:775201. doi: 10.1088/1751-8121/aaf8e7

CrossRef Full Text | Google Scholar

14. Yang Y, Klein DJ. A recursion formula for resistance distances and its applications. Discrete Appl Math. (2013) 161:2702–15. doi: 10.1016/j.dam.2012.07.015

CrossRef Full Text | Google Scholar

Keywords: resistance distance, electrical network, Hermitian matrix, Laplacian matrix, Moore-Penrose inverse

Citation: Yang Y, Wang D and Klein DJ (2020) A Relation Between Moore-Penrose Inverses of Hermitian Matrices and Its Application in Electrical Networks. Front. Phys. 8:239. doi: 10.3389/fphy.2020.00239

Received: 17 April 2020; Accepted: 02 June 2020;
Published: 10 July 2020.

Edited by:

Jia-Bao Liu, Anhui Jianzhu University, China

Reviewed by:

Weigen Yan, Jimei University, China
Zhongxun Zhu, South-Central University for Nationalities, China
Guihai Yu, Guizhou University of Finance and Economics, China

Copyright © 2020 Yang, Wang and Klein. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yujun Yang, yangyj@yahoo.com

Download