Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 301971 | 7 pages | https://doi.org/10.1155/2012/301971

A Fast Algorithm of Moore-Penrose Inverse for the Loewner-Type Matrix

Academic Editor: Tadeusz Kaczorek
Received25 Jul 2011
Accepted08 Oct 2011
Published04 Jan 2012

Abstract

In this paper, we present a fast algorithm of Moore-Penrose inverse for Loewner-type matrix with full column rank by forming a special block matrix and studying its inverse. Its computation complexity is , but it is by using .

1. Introduction

Loewner matrix was first studied by Loewner in 1934 in [1]. He studied the relations of various Loewner matrices via the characteristic of the monotone matrix function and the problem of rational interpolation at that time. Since then, more further studies had been carried out by many scientists in [25], such as Belevitch, Donoghue Jr., Fiedler, Chen and Zhang. From their mentioned works, we can find various properties of Loewner matrix and its application in the rational interpolation. In 1984, Vavrín presented a fast algorithm for the inverse of Loewner matrix in [6], then the fast algorithm for the Loewner system was got accordingly. Rost and Vavrín presented a fast algorithm for the system whose coefficient is a Loewner-Vandermonde matrix from 1995 to 1996 in [7, 8]. Lu gave a fast triangular factorization algorithm for the symmetrical Loewner-type matrix in 2003 in [9]. Xu et al. and so forth gave a fast triangular factorization algorithm for the inverse of symmetrical Loewner-type matrix in 2003 in [10]. In 2009, Tong et al. gave a fast algorithm of the Moore-Penrose inverse for symmetrical Loewner-type matrix [11]. In this paper, we generalize the fast algorithm of symmetrical Loewner-type matrix to the Loewner-type matrix. The theory and computation of generalized inverse matrix arise in various applications such as optimization theory, control theory, computation mathematics, and mathematical statistics. Therefore, there are very important theoretical and practical significance when we study the fast algorithm of Moore-Penrose inverse for the Loewner-type matrix.

This paper is organized as follows: in Section 1, we present some preliminaries results. The fast algorithm of Moore-Penrose inverse for Loewner-type matrix is driven in Section 2. We give some numerical examples to illustrate the fast algorithm obtained in Section 3.

2. Preliminaries

A Loewner matrix is a matrix of the form , where () are given numbers and . A Loewner-type matrix is a matrix which satisfies where , , , and . Loewner matrix is a special case of Loewner-type matrix, and it satisfies .

Let be the rank of Loewner-type matrix . Forming an matrix we obtain

It is obvious that we may obtain the Moore-Penrose inverse of Loewner-type matrix by (2.3) if we can get .

Now let us begin to find .

The following result will be useful for getting .

Lemma 2.1 (see [12]). Let all the leading principal submatrices of matrix whose order is be invertible. Linear system be given, where . Note that , and let , be solution vectors of linear systems , differently, then where .

3. Fast Algorithm of the Moore-Penrose Inverse for Loewner-Type Matrix

By using (2.1), we know that satisfies

Let all the leading principal submatrices ?? of be invertible. If , by virtue of (3.1), we have

Let , be solution vectors of linear systems differently, then using , we obtain , . By virtue of Lemma 2.1 we have where , , and is a solution vector of . Multiplying (3.2) by on the left and on the right differently, we obtain Multiplying (3.5) by and noting that , we have Substituting (3.4) in (3.6) and noting that , , we have and hence, Now, let us look for . Choosing the th of , namely and using (3.8), we have

Note that , where is the th column of , then Letting in (3.5) and multiplying it by on the right, we have

So, we obtain the fast algorithm of Moore-Penrose inverse of Loewner-type matrix by (2.3), (3.4), and (3.8)~(3.12).

Algorithm
Step 1.
For ,
Step 2.
For , then

The algorithm requires multiplication and division operations and addition and subtraction operations, and the computation complexity is , but it is by using .

4. Numerical Examples

We get the Moore-Penrose inverse matrix of Loewner-type matrix with Fortran program in computer, and the results of the partial numerical examples are given as follows in Table 1 (the error is measured by the 2-norm of vector, and the time is measured by second).


Orders of matrixFast algorithm
Err.Tim.Err.Tim.

10000207.254 −1628.109 −1411
20000202.398 −1522.788 −1431
30000202.019 −1552.106 −1489
40000209.108 −1481.001 −14201
60000204.901 −14363.760 −14705

Example 1. Consider

From above examples and many more examples not given, we find that the stability of the fast algorithm is very good, and it needs shorter time than that of .

Acknowledgment

This work is being supported by the National Natural Science Foundation Grants of China (no. 60974082).

References

  1. K. Loewner, “Uber mototone matrixfunktionen,” Math, vol. 38, pp. 177–216, 1934. View at: Google Scholar
  2. V. Belevitch, “Interpolation matrices,” Philips Research Report, vol. 25, no. 5, pp. 337–369, 1970. View at: Google Scholar
  3. W. F. Donoghue Jr., Monotone Matrix Functions and Analytic Continuation, Springer, New York, NY, USA, 1974.
  4. M. Fiedler, “Hankel and Loewner matrices,” Linear Algebra and Its Applications, vol. 58, pp. 75–95, 1984. View at: Google Scholar
  5. Z. Vavřín, “Inverses of Loewner matrices,” Linear Algebra and Its Applications, vol. 63, pp. 227–236, 1984. View at: Publisher Site | Google Scholar
  6. G. N. Chen and H. P. Zhang, “More on Loewner matrices,” Linear Algebra and Its Applications, vol. 203-204, pp. 265–300, 1994. View at: Google Scholar
  7. K. Rost and Z. Vavřín, “Recursive solution of Loewner-Vandermonde systems of equations. II,” Linear Algebra and Its Applications, vol. 223-224, pp. 597–617, 1995. View at: Google Scholar
  8. K. Rost and Z. Vavřín, “Recursive solution of Loewner-Vandermonde systems of equations. I,” Linear Algebra and Its Applications, vol. 233, pp. 51–65, 1996. View at: Publisher Site | Google Scholar
  9. Q. Lu, “A fast algorithm for the triangular factorization of the symmetric Loewner type matrix,” Journal of Engineering Mathematics, vol. 20, no. 2, pp. 139–142, 2003 (Chinese). View at: Google Scholar
  10. M. Xu, Z. Xu, and Z.-K. Shi, “A fast triangular factorization algorithm for the inverse of symmetric loewner-type matrix,” Journal of China Three Gorges University, vol. 6, pp. 552–554, 2003 (Chinese). View at: Google Scholar
  11. Q. Tong, Q. Lu, S. Liu, and J. Chai, “A fast algorithm of Moore-Penrose inverse for the symmetric Loewner-type matrix,” in Proceedings of the International Conference on Information Engineering and Computer Science (ICIECS '09), pp. 2155–2158, December 2009. View at: Publisher Site | Google Scholar
  12. Z. Xu, K.-Y. Zhang, and Q. Lu, Fast Algorithm of Toeplitz-Like Matrix(in Chinese), Northwest Polytechnical University Press, Xi'an, China, 1999.

Copyright © 2012 QiuJuan Tong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

657 Views | 1761 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.