Mathematical Problems in Engineering

Volume 2015, Article ID 924757, 7 pages

http://dx.doi.org/10.1155/2015/924757

## A New Derivation and Recursive Algorithm Based on Wronskian Matrix for Vandermonde Inverse Matrix

^{1}School of Information Science and Engineering, Hunan International Economics University, Changsha, Hunan 410205, China^{2}College of Science, National University of Defense Technology, Changsha, Hunan 410073, China

Received 26 June 2015; Accepted 27 September 2015

Academic Editor: Xinkai Chen

Copyright © 2015 Qun Zhou 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.

#### Abstract

For an analytical expression of Vandermonde inverse matrix, a new derivation process based on Wronskian matrix and Lagrange interpolation polynomial basis is presented. Recursive formula and implementation cases for the direct formula of Vandermonde inverse matrix are given based on deriving the unified formula of Wronskian inverse matrix. For the calculation of symbol-type Vandermonde inverse matrix, the direct formula and recursive method are verified to be more efficient than Mathematica which is good at symbolic computation by comparing the computing time in Mathematica. The process and steps of recursive algorithm are relatively simple. The derivation process and idea both have very important values in theory and practice of Vandermonde and generalized Vandermonde inverse matrix.

#### 1. Introduction

Vandermonde matrix and generalized Vandermonde matrix and their inverses are always widely concerned in many research fields, such as numerical analysis, data and signal processing, and control theory [1–8]. There are many effective ways to calculate the Vandermonde inverse matrix. Fox example, Tou [1] and Reis [2] obtained the matrix formula from the coefficients of polynomial terms. Neagoe [3] deduced an analytic formula of complex-type Vandermonde inverse matrix based on symmetric polynomials. Eisinberg and Fedele [4] presented a general explicit formula for the elements of inverse matrix and two different algorithms were deduced. In this paper, based on results in [1, 2], a new derivation process of an analytical expression of Vandermonde inverse matrix is presented based on Wronskian matrix and Lagrange interpolation polynomial basis, and recursive formula and implementation cases for the direct formula of Vandermonde inverse matrix are presented based on deducing the unified formula of Wronskian inverse matrix.

The rest of this paper is arranged as follows. Firstly, we summarize the main idea and analytical formulas presented in [3] and give a new derivation process for the formula in Section 2. In Section 3, a recursive formula of Vandermonde inverse matrix is deduced. In Section 4, numerical simulations are presented with some cases in Mathematica and conclusion and prospect are given by comparing and analyzing the numerical results in the last section.

#### 2. A New Derivation of Vandermonde Inverse Matrix

If are different numbers, that is to say, they are not equal to each other, then we define the -order matrixwhere as a Vandermonde matrix. For the inverse of this matrix, Neagoe [3] obtained the calculation formula (2) according to relations between the determinant of the matrix (1) and the determinant of the matrix which was obtained by adding a row whose elements were and a column whose elements were to (1) and based on the calculation formula of inverse matrix by using the matrix determinant:where are column subscripts and are row subscripts and

Here we will use a new method entirely different from [3] to deduce formula (2).

For a group of linearly independent functions , where has the ()th derivative, Wronskian matrix is defined asIf there is another group of functions which satisfythen the inverse matrix of (4) can be obtained [9] aswhere and .

If we set , then we can verify that the Lagrange interpolation polynomial basiswhere satisfy (5); that is to say,For Wronskian matrix ,From (6), the inverse matrix of iswhere are row subscripts and are column subscripts in the above two matrices.

Denote by the diagonal matrix composed of ; that is,It is easily proved thatDenoted byusing inverse matrix formula of matrix multiplication, we haveFrom (10) and (14), we can get

Let ; then (12) changes toIts inverse matrix isFor the values of in (2), we specify them by the following case. For example, taking , from (3), we can getThen we havewhere are row subscripts and are column subscripts, and from (17) and (19), we can getThus, according to the corresponding relations of elements between and and and , we can obtain (2); that is, where are row subscripts and are column subscripts.

#### 3. Recursive Algorithm of Vandermonde Inverse Matrix

By the above method of function matrix, we denote asAccording to the property of Vandermonde determinant, we have

For the convenience of description, we introduce the following notations.

denotes the algebraic complement of the element of whose row subscript is and column subscript is . denotes th derivative on of and

It is easily obtained thatFrom (25) and (26), we can get

Calculating th derivative on at both sides of (26), we can obtainThey can lead toIt is easy to know thatDenoted by , thenFrom (29)–(31), we can getAnd set when , , , , and

Let in (32); the recursion formula of elements of Vandermonde inverse matrix can be obtained.

For simplicity, denote . Because , all elements do not contain , so we denote -order Vandermonde matrix composed of asThus, each element of Vandermonde inverse matrix can be determined by the following recursive formula:where and set when , , , , and

When , and .

When , .

#### 4. Application Case and Simulation

We will take the 3-order Vandermonde inverse matrix as an example to illustrate the recursive steps of (34) in detail.

LetFrom (34), we have

Taking and as the initial values to calculate each element of , then we can getSince , thenFrom (34) again, we can getSince , thenSubstituting to (40), we haveAs in , so from (36), this leads toMeanwhile, we can also obtain

For , there isAs in , so , and from (35), we can obtainAs in , so and andFrom (45) and (36), ; then we can get ; that is,Meanwhile, we can also get

The other elements of , , , can be obtained by a similar process. Finally, we can getMeanwhile, we can getSo, the 3-order Vandermonde inverse matrix of (35) isand the 2-order Vandermonde inverse matrix is

We choose different order Vandermonde matrices to compare the computing time spent on calculating their inverse matrices in mathematical software Mathematica 10. We divide the calculations into two types: symbolic inverse matrix and numerical inverse matrix. For the limitation of computing time, we take only the order as and the total time of each calculation is the time spent by doing 100 repeated calculations; its unit is second. When we calculate the numerical Vandermonde inverse matrix, we take . We obtain the results as Tables 1 and 2.