Research article

A particular matrix with exponential form, its inversion and some norms

  • Received: 24 December 2021 Revised: 14 February 2022 Accepted: 17 February 2022 Published: 25 February 2022
  • MSC : 15A15, 15A60

  • In this paper, we study a particular n×n matrix A=[akij]ni,j=1 and its Hadamard inverse A(1), whose entire elements are exponential form ak=e(kn)=e2πikn, where kij=min(i,j)+1. We study determinants, leading principal minor and inversions of A, A(1). Then the defined values of Euclidean norms, lp norms and spectral norms of these matrices are presented, rather than upper and lower bounds, which are different from other articles.

    Citation: Baijuan Shi. A particular matrix with exponential form, its inversion and some norms[J]. AIMS Mathematics, 2022, 7(5): 8224-8234. doi: 10.3934/math.2022458

    Related Papers:

    [1] Emrah Polatlı . On some properties of a generalized min matrix. AIMS Mathematics, 2023, 8(11): 26199-26212. doi: 10.3934/math.20231336
    [2] Hasan Gökbaş . Some properties of the generalized max Frank matrices. AIMS Mathematics, 2024, 9(10): 26826-26835. doi: 10.3934/math.20241305
    [3] Wanlin Jiang, Kezheng Zuo . Revisiting of the BT-inverse of matrices. AIMS Mathematics, 2021, 6(3): 2607-2622. doi: 10.3934/math.2021158
    [4] Qin Zhong, Ling Li . Notes on the generalized Perron complements involving inverse $ {{N}_{0}} $-matrices. AIMS Mathematics, 2024, 9(8): 22130-22145. doi: 10.3934/math.20241076
    [5] Jin Zhong, Yilin Zhang . Dual group inverses of dual matrices and their applications in solving systems of linear dual equations. AIMS Mathematics, 2022, 7(5): 7606-7624. doi: 10.3934/math.2022427
    [6] Deshu Sun . Note on error bounds for linear complementarity problems involving $ B^S $-matrices. AIMS Mathematics, 2022, 7(2): 1896-1906. doi: 10.3934/math.2022109
    [7] Maja Nedović, Dunja Arsić . New scaling criteria for $ H $-matrices and applications. AIMS Mathematics, 2025, 10(3): 5071-5094. doi: 10.3934/math.2025232
    [8] Huahai Qiu, Li Wan, Zhigang Zhou, Qunjiao Zhang, Qinghua Zhou . Global exponential periodicity of nonlinear neural networks with multiple time-varying delays. AIMS Mathematics, 2023, 8(5): 12472-12485. doi: 10.3934/math.2023626
    [9] Oe Ryung Kang, Jung Hoon Kim . The $ l_\infty $-induced norm of multivariable discrete-time linear systems: Upper and lower bounds with convergence rate analysis. AIMS Mathematics, 2023, 8(12): 29140-29157. doi: 10.3934/math.20231492
    [10] Zhigang Zhou, Li Wan, Qunjiao Zhang, Hongbo Fu, Huizhen Li, Qinghua Zhou . Exponential stability of periodic solution for stochastic neural networks involving multiple time-varying delays. AIMS Mathematics, 2024, 9(6): 14932-14948. doi: 10.3934/math.2024723
  • In this paper, we study a particular n×n matrix A=[akij]ni,j=1 and its Hadamard inverse A(1), whose entire elements are exponential form ak=e(kn)=e2πikn, where kij=min(i,j)+1. We study determinants, leading principal minor and inversions of A, A(1). Then the defined values of Euclidean norms, lp norms and spectral norms of these matrices are presented, rather than upper and lower bounds, which are different from other articles.



    Matrix theory is widely used in a variety of areas including computing science, applied mathematics, neural network nonlinear system and others. Recently, studying the determinants and some norms of particular matrix has been a hot and important topic in matrix theory [1]. From past to present, different types of matrices have been defined and studied properties of various properties, such as determinants, inverse, norms have been studied by some mathematicians. For example, Akbulak [2] studied Hadamard exponential Hankel matrix of the form eHn=[ei+j]i,j and presented lp norm, two upper bounds for spectral norm and eigenvalues of this matrix. By the identities of Gamma function, Bozkurt [3,4,5] determined lp norm of Cauchy-Toeplitz matrices and Cauchy-Hankel matrices respectively. Civciv [6] established a bound for spectral and lp norms of the Khatri-Rao products of Cauchy-Hankel matrices. Solak [7] studied the matrix of the form B=[bij], where bij=a+min(i,j)+1 and obtained its Euclidean norm and inversion of B. A few years ago, these matrices attracted the attention of C. Moler, who had experimentally discovered that most of their singular values are clustered near π. The authors in [8] studied some bounds for the spectral norm of particular matrix of the form A=[amin(i,j)]n1i,j=0, where a is a real positive number.

    Inspired by above articles, we study a particular matrix A=[akij]n1i,j=0 and its Hadamard inverse A(1), where kij=min(i,j)+1, ak=e(kn), and e(x)=e2πix. Particularly, e(n)=1, and |e(x)|=1. Hadamard inverse of A is

    A(1)=[(1e(kn))]n×n=[e(kn)]n×n.

    Two of these matrices are called Min and Hadamard inverse matrices as following, respectively, for more information, we can see reference [3],

    A=(e(1n)e(1n)e(1n)e(1n)e(1n)e(1n)e(2n)e(2n)e(2n)e(2n)e(1n)e(2n)e(3n)e(3n)e(3n)e(1n)e(2n)e(3n)e(n1n)e(n1n)e(1n)e(2n)e(3n)e(n1n)e(nn))n×n, (1.1)
    A(1)=(e(1n)e(1n)e(1n)e(1n)e(1n)e(1n)e(2n)e(2n)e(2n)e(2n)e(1n)e(2n)e(3n)e(3n)e(3n)e(1n)e(2n)e(3n)e((n1)n)e((n1)n)e(1n)e(2n)e(3n)e((n1)n)e(nn))n×n. (1.2)

    Looking at above matrices we see that these matrices are symmetric, so they are normal matrix, and the positions of the each entries make a "Γ pattern".

    For any m×n matrix A, the lp norm of A is defined by

    Ap=(mi=1nj=1|aij|p)1p.

    For p=2 this norm is called Euclidean norm showed by

    AE=(mi=1nj=1|aij|2)12.

    The spectral norm of matrix A is defined by

    A2=max1inλi(AHA),

    where λi(AHA) are the eigenvalues of matrices AHA and AH is the conjugate transpose of A.

    The spectral radius of A is defined by

    ρ(A)=maxi|λi(A)|,i=1,2,,n.

    For normal matrix, we know that the spectral norm of A is equal to spectral radius. Different from other articles, the matrix elements are complex numbers e(kn) whose modulus is 1. Therefore, based on the special identities, we can give the defined value of spectral norm rather than upper and lower bound estimates.

    First, by some methods, in the beginning of this section we show that determinants of matrix A and its Hadamard inverse A(1). Then we get some relationship of leading principle minor. Subsequently we present that these two matrices are invertible but not positive definite, and then we present that the inversion of A and A(1) are tridiagonal matrices. After that we get the defined value of Euclidean norm, lp norm and spectral norm. All definitions and statements of this paper are available in reference [1,9].

    Theorem 1. Let A be a matrix as in (1.1), n2, then

    detA=e(1n+n2)(e(1n)1)n10.

    Proof. By using elementary row operations on (1.1), we have

    detA=det(e(1n)e(1n)e(1n)e(1n)e(1n)0e(2n)e(1n)e(2n)e(1n)e(2n)e(1n)e(2n)e(1n)00e(3n)e(2n)e(3n)e(2n)e(3n)e(2n)000e(n1n)e(n2n)e(n1n)e(n2n)0000e(nn)e(n1n))n×n=e(1n+n2)(e(1n)1)n10.

    Therefore, the matrix A is invertible. However, because the matrix elements are complex numbers, we can't confirm that all leading principal minor or all eigenvalues of A are positive, and then matrix A isn't positive definite. Now we give the inversion of matrix A.

    Lemma 1. Let A is an n×n nonsigular matrix and b is an n×1 matrix, c is a real number.

    If

    M=(AbbTc),

    then the inversion of M is

    N=(A1+1lA1bbTA11lA1b1lbTA11l),

    where l=cbTA1b.

    Proof. By the definition of M and N we have MN=En+1. Thus M1=N.

    Theorem 2. Let A be as in (1.1), then A is invertible and the inversion of A is a symmetric tridiagonal matrix of the form

    A1=(c0c1000c1b2c2000c2b3c3000cn2bn1cn1000cn1cn1)n×n,

    where

    bk=e(1n)+1e(kn)(e(1n)1),ck=1e(kn)(e(1n)1).

    Proof. Using elementary row operations to solve the inversion of A, namely

    (AI)(IA1).

    By mathematical induction on n, for n=2, we have

    A=(e(1n)e(1n)e(1n)e(2n))2×2
    (e(1n)e(1n)10e(1n)e(2n)01)(e(1n)e(1n)100e(2n)e(1n)11)(111e(1n)0011e(2n)e(1n)1e(2n)e(1n))(101e(1n)+1e(2n)e(1n)1e(2n)e(1n)011e(2n)e(1n)1e(2n)e(1n)),

    so we have, for 2×2 matrix A,

    A1=(1e(1n)+1e(2n)e(1n)1e(2n)e(1n)1e(2n)e(1n)1e(2n)e(1n)).

    In this way, for 3×3 matrix A, we can get

    A1=(1e(1n)+1e(2n)e(1n)1e(2n)e(1n)01e(2n)e(1n)1e(2n)e(1n)+1e(3n)e(2n)1e(3n)e(2n)01e(3n)e(2n)1e(3n)e(2n)),

    Therefore, for n×n matrix A, we can get

    A1=(1a1+1a2a11a2a1001a2a11a2a1+1a3a21a3a2001a3a21a3a2+1a4a301an1an21an1an2+1anan11anan1001anan11anan1)n×n,

    where ak=e(kn).

    In addition,

    1e(k+1n)e(kn)=1e(kn)(e(1n)1),1e(kn)e(k1n)+1e(k+1n)e(kn)=e(1n)+1e(kn)(e(1n)1).

    In fact, by mathematical induction on n.

    The result is true for n=2, now assume that the result is true for n, that is A=[akij]n×n,A1=[akij]1n×n.

    Thus by taking b=(e(1n),e(2n),,e(nn))T,bT=(e(1n),e(2n),,e(nn)), and c=e(n+1n) along with Lemma 1 the proof is completed for n+1. Therefore the result is true for each n. For convenience, we denote

    bk=e(1n)+1e(kn)(e(1n)1),ck=1e(kn)(e(1n)1).

    Hence, we get

    A1=(c0c1000c1b2c2000c2b3c3000cn2bn1cn1000cn1cn1)n×n.

    The proof of Theorem 2 provides a method to solve the inversion of "Γ" matrix which is tridiagonal matrix.

    The next part, Euclidean norm, lp norm and spectral norms of matrix A, are considered, we can get concise results.

    Theorem 3. Let A be a matrix as in (1.1), then the lp norm of A is

    Ap=n2p.

    Proof. By definition of the Euclidean norm, by |e(kn)|=1, we have

    App=nk=1(2n2k+1)|e(kn)|p=nk=1(2n2k+1)=n2.

    Theorem 4. Let A be a matrix as in (1.1), then we get the spectral norm of A,

    A2=max(1,2sin(πn)).

    Proof. Since matrix A is Hermite matrix which are also normal matrix. So the spectral norm of A is equal to its spectral radius. By Theorem 1, we have the eigenvalues of A are respectively:

    e(1n),e(2n)e(1n),e(3n)e(2n),,e(nn)e(n1n).

    For the modulus of these eigenvalues, we have |e(1n)|=1, and

    |e(k+1n)e(kn)|=|e(kn)(e(1n)1)|=|e(1n)1|=|cos(2πn)1+isin(2πn)|=2sin(πn).

    Hence, spectral radius of A is ρ(A)=max(1,2sin(πn)).

    Namely,

    A2=ρ(A)=max(1,2sin(πn)).

    Theorem 5. Let A be a matrix as in (1.1) and n denotes the leading principal minor of A, then we have

    (1)1=e(1n),n=e(1n+n2)(e(1n)1)n1=detA;(2)k+1k=e(k+1n)e(kn);(3)nn2=e(1n)2n1;(4)12n=e(n23+n2+16)(e(1n)1)n(n1)2.

    Proof. By Theorem 1, we can get (1)

    1=e(1n),2=e(1n)e(1n)(e(1n)1),3=e(1n)e(1+2n)(e(1n)1)2,4=e(1n)e(1+2+3n)(e(1n)1)3,n2=e(1n)e(1+2+3++(n3)n)(e(1n)1)n1=e(1n+n2)(e(1n)1)n3,
    n1=e(1n)e(1+2+3++(n2)n)(e(1n)1)n1=e(1n+n2)(e(1n)1)n2,n=e(1n)e(1+2+3++(n1)n)(e(1n)1)n1=e(1n+n2)(e(1n)1)n1=detA.

    Hence, we have

    (2)k+1k=e(kn)(e(1n)1)=e(k+1n)e(kn).(3)nn2=e(5n)(e(1n)1)2n4,2n1=e(4n)(e(1n)1)2n4.

    So we get

    nn2=e(1n)2n1.

    By e(n)=1, and n(n1)2 is always integer number, that is to say e(n(n1)2)=1, then we have

    (4)12n=e(1+(1+2)+(1+2+3)++(1+2++(n1))n)(e(1n)1)n(n1)2=e(n23+n2+16)(e(1n)1)n(n1)2.

    Next we're going to discuss the determinant, inversion, Euclidean norm, lp norm and spectral norm of A(1).

    Theorem 6. Let A(1) be a matrix as in (1.2), n2, then

    detA(1)=e(n2)(1e(1n))n10.

    Proof. By using elementary row operations on (1.2), by e(12)=1, we have

    detA(1)=det(e(1n)e(1n)e(1n)e(1n)e(1n)0e(2n)e(1n)e(2n)e(1n)e(2n)e(1n)e(2n)e(1n)00e(3n)e(2n)e(3n)e(2n)e(3n)e(2n)000e(1nn)e(2nn)e(1nn)e(2nn)0000e(nn)e(1nn))=e(n2)(1e(1n))n10.

    Therefore, the matrix A(1) is invertible. Now we give the inversion of matrix A(1).

    Theorem 7. Let A(1) be a matrix as in (1.2), n2, then the inversion of A(1) is a tridiagonal matrix as following:

    [A(1)]1=(h0h1000h1f2h2000h2f3h3000hn2fn1hn1000hn1hn1)n×n,

    where

    fk=e(kn)(1+e(1n))1e(1n),hk=e(kn)1e(1n).

    Proof. In similar way, by using elementary row operations to solve the inversion of A(1), namely (A(1)I)(I[A(1)]1).

    For n=2, we have

    A(1)=(e(1n)e(1n)e(1n)e(2n))2×2
    (e(1n)e(1n)10e(1n)e(2n)01)(e(1n)e(1n)100e(2n)e(1n)11)(11e(1n)0011e(2n)e(1n)1e(2n)e(1n))(10e(1n)+1e(2n)e(1n)1e(2n)e(1n)011e(2n)e(1n)1e(2n)e(1n)),

    so for 2×2 matrix A(1),

    [A(1)]1=(e(1n)+1e(2n)e(1n)1e(2n)e(1n)1e(2n)e(1n)1e(2n)e(1n)).

    In this way, for 3×3 matrix A(1), we can get

    [A(1)]1=(e(1n)+1e(2n)e(1n)1e(2n)e(1n)01e(2n)e(1n)1e(2n)e(1n)+1e(3n)e(2n)1e(3n)e(2n)01e(3n)e(2n)1e(3n)e(2n)).

    Therefore, for n×n matrix A(1), we can get

    [A(1)]1=(1a1+1a2a11a2a1001a2a11a2a1+1a3a21a3a2001a3a21a3a2+1a4a301an1an21an1an2+1anan11anan1001anan11anan1)n×n,

    where ak=e(kn).

    Actually,

    1e(kn)e(k1n)=1e(kn)(1e(1n))=e(kn)1e(1n),1e(kn)e(k1n)+1e(k+1n)e(kn)=e(kn)(e(1n)+1)(1e(1n)).

    By mathematical induction on n. The result is true for n=2. Now assume that the result is true for n, Thus by taking b=(e(1n),e(2n),,e(nn))T,bT=(e(1n),e(2n),,e(nn)), and c=e(n+1n) along with Lemma 1 the proof is completed for n+1. Therefore the result is true for each n. For convenience, we denote

    fk=e(kn)(1+e(1n))1e(1n),hk=e(kn)1e(1n).

    Hence, we have

    [A(1)]1=(h0h1000h1f2h2000h2f3h3000hn2fn1hn1000hn1hn1).

    Theorem 8. Let A(1) be a matrix as in (1.2), then its lp norm is:

    A(1)p=n2p.

    Proof. By the definition of lp norm, we have

    A(1)pp=nk=1(2n2k+1)|e(kn)|p=nk=1(2n2k+1)=n2.

    When n=2, we can get the Euclidean normA(1)E=n.

    Theorem 9. Let A(1) be a matrix as in (1.2), n2, then the spectral norm is

    A(1)2=max(1,2sin(πn)).

    Proof. By Theorem 6, we have the eigenvalues of A(1) are respectively:

    e(1n),e(2n)e(1n),,e(nn)e(n1n).

    Since matrix A(1) is normal matrix whose spectral norm is equal to its spectral radius. So we find the modulus of above eigenvalues: |e(1n)|=1,

    |e(kn)e(k1n)|=|e(kn)(1e(1n))|=|1e(1n)|=|1cos(2πn)isin(2πn)|=2sin(πn).

    Hence,

    A(1)2=ρ(A)=max(1,2sin(πn))={1,n6;2sin(πn),n<6.

    So far, we've proved all the theorems. We can easily find that matrix A and A(1) have same Euclidean norm, lp norm and spectral norm.

    In this paper, we attempt to compute the determinant, inverse, lp-norm, and some other properties including those of its Hadamard inverse of particular matrices involving exponential forms and trigonometric functions. The computation complexity of this paper is lower than the previous work. Based on the special properties of exponential form, we get the defined values of lp norms and spectral norms of particular matrix whose entries are complex numbers e(kn). These results will expand the application range of matrix norm and enrich the system of matrix theory.

    This work is support by National Natural Science Foundation of China (NSF.11771351).

    The author declares no conflict of interest.



    [1] F. Zhang, Matrix theory: basic results and techniques, New York, NY: Springer, 2011. http://dx.doi.org/10.1007/978-1-4614-1099-7
    [2] A. Ipek, M. Akbulak, Hadamard exponential Hankel matrix, its eigenvalues and some norms, Math. Sci. Lett., 1 (2012), 81–87. http://dx.doi.org/10.12785/msl/010110 doi: 10.12785/msl/010110
    [3] D. Bozkurt, On the lp norms of almost Cauchy-Toeplitz matrices, Turk. J. Math., 20 (1996), 545–552.
    [4] D. Bozkurt, On the lp norms of Cauchy-Toeplitz matrices, Linear and Multilinear Algebra, 44 (1998), 341–346. http://dx.doi.org/10.1080/03081089808818569 doi: 10.1080/03081089808818569
    [5] D. Bozkurt, A note on the spectral norms of the matrices connected integer numbers sequence, 2011, arXiv: 1105.1724.
    [6] H. Civciv, R. Turmen, On the bounds for the spectral and lp norms of the Khatri-Rao products of Cauchy-Hankel matrices, J. Ineq. Pure and Appl. Math., 7 (2006), 195.
    [7] S. Solak, M. Bahsi, A particular matrix and its some properties, Sci. Res. Essays, 8 (2013), 1–5. http://dx.doi.org/10.5897/SRE11.410 doi: 10.5897/SRE11.410
    [8] S. H. Jafari-Petroudi, B. Pirouz, A particular matrix, its inversion and some norms, Appl. Comput. Math., 4 (2015), 47–52. http://dx.doi.org/10.11648/j.acm.20150402.13 doi: 10.11648/j.acm.20150402.13
    [9] S. H. Jafari-Petroudi, M. Pirouz, On the bounds for the spectral norm of particular matrices with Fibonacci and Lucas numbers, Int. J. Adv. Appl. Math. and Mech., 3 (2016), 82–90.
  • This article has been cited by:

    1. Baijuan Shi, Can Kızılateş, On linear algebra of r-Hankel and r-Toeplitz matrices with geometric sequence, 2024, 70, 1598-5865, 4563, 10.1007/s12190-024-02151-y
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1724) PDF downloads(67) Cited by(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog