Next Article in Journal
On a Matrix Formulation of the Sequence of Bi-Periodic Fibonacci Numbers
Previous Article in Journal
Combining Glimm’s Scheme and Operator Splitting for Simulating Constrained Flows in Porous Media
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Totally Positive Wronskian Matrices and Symmetric Functions

Department of Applied Mathematics, University Research Institute of Mathematics and Its Applications (IUMA), University of Zaragoza, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(9), 589; https://doi.org/10.3390/axioms13090589
Submission received: 5 July 2024 / Revised: 23 August 2024 / Accepted: 28 August 2024 / Published: 29 August 2024
(This article belongs to the Section Algebra and Number Theory)

Abstract

:
The elements of the bidiagonal decomposition (BD) of a totally positive (TP) collocation matrix can be expressed in terms of symmetric functions of the nodes. Making use of this result, and studying the relation between Wronskian and collocation matrices of a given TP basis of functions, we can express the entries of the BD of Wronskian matrices as the values of certain symmetric functions evaluated at a single node. Moreover, in the case of polynomial bases, we obtain compact formulae for the entries of the BD of their Wronskian matrices. Interesting examples illustrate the applications of the obtained formulae.

1. Introduction

The class of totally positive matrices possesses rich mathematical properties and numerous applications. It has been extensively studied (cf. [1,2,3,4,5,6]), and attracts significant interest across various fields of mathematics, including approximation theory, combinatorics, computer-aided geometric design, and economics. Let us recall that a matrix is TP if all its minors are nonnegative, and strictly totally positive (STP) if all its minors are positive.
Different characterizations of the total positivity property of a matrix exist in terms of the sign of certain collections of its minors. They imply a reduction in the number of determinants that need to be analyzed to determine whether a matrix is TP or STP. For example, a matrix is STP if and only if all its minors with consecutive rows and columns are positive (see [7,8]). This characterization may be refined, as pointed out in [4], where it is shown that only minors with consecutive rows and columns that include either the first row or the first column need to be checked. These minors are usually referred to as initial minors. The above observations significantly reduce the complexity of the tests for STP matrices. For later use in this paper, we quote this fundamental result in the next theorem.
Theorem 1 
(Theorem 4.1 of [4]). Let A be a real l × n matrix. Then, A is STP if and only if every initial minor of A is positive.
Considering the Cauchy–Binet formula for determinants, it was established that the product of the TP matrices results in another TP matrix (see Theorem 3.1 of [1]). This property opened the possibility of factorizing TP matrices into products of simpler TP matrices. This endeavor has been the focus of extensive literature [9,10,11,12,13,14]. In this regard, the most definite fact is that a nonsingular TP matrix always admits a bidiagonal decomposition that can be fully determined by its initial minors [6].
In [5,6], it is shown that a nonsingular TP matrix A R n × n can be written as
A = F n 1 F 1 D G 1 G 2 G n 1 ,
where D is a diagonal matrix D = diag ( p i , i ) 1 i n whose diagonal entries are called pivots, and F i R n × n (respectively, G i R n × n ), i = 1 , , n 1 , are the TP, lower (respectively, upper) triangular bidiagonal matrices of the following form:
F i = 1 1 m i + 1 , 1 1 m n , n i 1 , G i T = 1 1 m ˜ i + 1 , 1 1 m ˜ n , n i 1 .
The off-diagonal entries m i , j are usually called multipliers, and satisfy m i , j = p i , j / p i 1 , j , with 
p i , 1 = a i , 1 , 1 i n , p i , j = det A [ i j + 1 , , i 1 , , j ] det A [ i j + 1 , , i 1 1 , , j 1 ] , 1 < j i n ,
where A [ i 1 , , i r | j 1 , , j s ] denotes the submatrix from A formed with rows i 1 , , i r and columns j 1 , , j s . In the following, we shall denote A [ i 1 , , i r ] : = A [ i 1 , , i r | i 1 , , i r ] . Similarly, m ˜ i , j = p ˜ i , j / p ˜ i 1 , j , where p ˜ i , 1 = a 1 , i and p ˜ i , j can be obtained by means of the quotient of minors in (2) for the transposition A t of the matrix A.
The bidiagonal decomposition (BD) (1) of nonsingular TP matrices is important for theoretical reasons and practical numerical applications, ensuring great accuracy in linear algebra computations involving matrices whose BD can be computed to high relative accuracy (HRA). There is rich literature delivering algorithms for resolving algebraic problems related to TP matrices to HRA through their BD (cf. [9,10,11,12,13,14,15]).
Given a basis F = ( f 1 , f 2 , , f n ) of functions defined on an interval I R and a set of parameters X = { x 1 , x 2 , , x n } with x 1 < x 2 < < x n within I, the collocation matrix of the basis F at X is defined as:
M F ( X ) : = f j ( x i ) 1 i , j n .
The basis F is TP (STP) if, for any X = { x 1 , x 2 , , x n } with x 1 < x 2 < < x n , the matrix M F ( X ) is TP (STP).
In this paper, we will show that the collocation matrices M F ( X ) connect the realm of total positivity and the field of symmetric functions. This relation happens because the initial minors of a collocation matrix
det M F ( X ) [ i j + 1 , , i 1 , , j ] , det M F ( X ) t [ i j + 1 , , i 1 , , j ]
are antisymmetric functions of the variables x 1 , x 2 , , x n and, in turn, can be expressed as the product of a Vandermonde determinant and a symmetric function of x 1 , , x n . Consequently, each initial minor of an n × n collocation matrix can be computed by evaluating a symmetric function. An interesting problem is finding the n 2 symmetric functions that encode the BD of the collocation matrices (3).
As a precedent, it was found in [16] that, surprisingly, the elements of the BD decomposition of the Cauchy-polynomial–Vandermonde matrices could be expressed in terms of Schur functions. Later, a systematic line of research was initiated in [15], where explicit formulas in terms of Schur polynomials were provided for the initial minors of collocation matrices of arbitrary polynomial bases. As a result, the straightforward computation of the BD for any collocation matrix of a polynomial basis can be achieved by evaluating the Schur polynomials at the nodes. The formulas provided in [15] were also used to determine the maximal interval, for which the polynomial bases are STP. Moreover, the techniques developed in [15] can be extended to any non-polynomial basis of functions, provided that the corresponding symmetric functions associated with the initial minors of the collocation matrices can be identified.
Intriguingly, the connection between TP bases and symmetric functions may be extended to various mathematical objects, provided that they are related to collocation matrices through certain limits. They include Wronskian matrices, which are the main focus of this work. In this paper, the initial minors of Wronskian matrices will be expressed as the limit of finite differences of certain collocation matrices. This observation allows us to apply the findings of [15], and derive a concise formula for the initial minors of Wronskian matrices.
The paper is organized as follows. Section 2 introduces the necessary concepts to make the article as self-contained as possible. Section 3 shows that any minor with consecutive rows of a Wronskian matrix can be expressed in terms of the limits of determinants of collocation matrices at equally spaced nodes. As a consequence, conditions guaranteeing the total positivity of Wronskian matrices are derived. Moreover, their minors with consecutive rows and columns are expressed in terms of symmetric functions. In particular, for the polynomial case, these minors are written in terms of Schur polynomials in Section 4. The applicability of the achieved formula is illustrated in Section 5 for Bernstein bases and recursive polynomial bases, such as Jacobi, Laguerre, Hermite, and Bessel bases. We conclude with an appendix containing the pseudocode of an algorithm for the computation of the minors (4) in the case of polynomial bases.

2. Preliminary Results

Let us consider the monomial polynomial basis ( m 1 , , m n ) with
m i ( x ) : = x i 1 , i = 1 , , n .
The corresponding collocation matrix at X = { x 1 , x 2 , , x n } is the well-known Vandermonde matrix
V X = V x 1 , , x n : = 1 x 1 x 1 n 1 1 x n x n n 1 ,
and its determinant satisfies
det V X = i < j ( x j x i ) .
Given h > 0 , and the equally spaced sequence
x k : = x + ( k 1 ) h , k = 1 , , n ,
the Vandermonde determinant satisfies
det V X = h n ( n 1 ) 2 r = 1 n 1 r ! .
This result can be extended to initial minors of V X as
det V x k , x k + 1 , , x l = h ( l k ) ( l k + 1 ) 2 r = 1 l k r ! .
Vandermonde matrices are STP if 0 < x 1 < < x n . Thus, we can say that the monomial polynomial basis is STP on the interval ( 0 , ) .
Now, let us recall that f ( x 1 , x 2 , , x n ) is a symmetric function if
f ( x 1 , x 2 , , x n ) = f ( x σ ( 1 ) , x σ ( 2 ) , , x σ ( n ) ) ,
for any permutation σ of the indices { 1 , 2 , , n } . On the other hand, f ( x 1 , x 2 , , x n ) is antisymmetric if
f ( x σ ( 1 ) , x σ ( 2 ) , , x σ ( n ) ) = sgn ( σ ) f ( x 1 , x 2 , , x n ) ,
where sgn ( σ ) is the signature of σ , taking the value + 1 if σ is even and 1 if σ is odd.
Note that any minor of the collocation matrix M F ( X ) can be considered as an antisymmetric function of the nodes involved. Moreover, since the Vandermonde determinant is nonzero for different values of the nodes, it is always possible to express any minor of the collocation matrix as the product of a Vandermonde determinant and a symmetric function.
Given a partition λ : = ( λ 1 , λ 2 , , λ p ) of size | λ | : = λ 1 + + λ p and length l ( λ ) : = p , such that λ 1 λ 2 λ p > 0 , the Jacobi definition of the corresponding Schur polynomial in n variables is, via the Weyl’s formula,
s λ ( x 1 , x 2 , , x n ) : = det ( x i λ j + n j ) 1 i , j n det ( x i n j ) 1 i , j n ,
and, by convention, s Ø; ( x 1 , x 2 , , x n ) : = 1 for the empty partition Ø.
Schur polynomials are symmetric functions in their arguments. In addition, we now list other well-known properties that will be used in the following sections:
(i)
s λ ( x 1 , , x n ) > 0 for positive values of x i , i = 1 , , n . Additionally,
s λ ( 1 , , 1 j ) = D λ , j ,
with
D λ , j = 1 k < l j λ k λ l + l k l k 0 .
(ii)
s λ ( x 1 , , x n ) = 0 if l ( λ ) > n .
(iii)
s λ ( x 1 , , x n ) is a homogeneous function of degree | λ | , i.e.,
s λ ( α x 1 , α x 2 , , α x n ) = α | λ | s λ ( x 1 , x 2 , , x n ) .
(iv)
As running over all the partitions of size | λ | , the corresponding Schur polynomials provide a basis for the space of symmetric homogeneous polynomials of degree | λ | . When considering all partitions, Schur polynomials furnish a basis of symmetric functions.
For more details, interested readers are referred to [17].
We finish this section by defining some symmetric functions that will be used in the following:
g F , i , j ( x i j + 1 , , x i ) : = det M F ( X ) [ i j + 1 , , i | 1 , , j ] det V x i j + 1 , , x i , 1 j i n , g F , j , i ( x 1 , , x j ) : = det M F t ( X ) [ i j + 1 , , i | 1 , , j ] det V x 1 , , x j , 1 j < i n .

3. Initial Minors and Total Positivity of Wronskian Matrices

Let F = ( f 1 , f 2 , , f n ) be a system of C n 1 functions on I. The Wronskian matrix of F at x I is defined as:
W F ( x ) : = f 1 ( x ) f 2 ( x ) f n ( x ) f 1 ( x ) f 2 ( x ) f n ( x ) f 1 ( n 1 ) ( x ) f 2 ( n 1 ) ( x ) f n ( n 1 ) ( x ) ,
where f ( x ) and f ( k ) ( x ) , k = 2 , , n 1 , denote the first and k-th derivative of f at x.
For a given F = ( f 1 , , f n ) , we shall denote
F 0 : = F , F m : = ( f 1 ( m ) , f 2 ( m ) , , f n ( m ) ) , m = 1 , , n 1 .
Note that the system F m needs to not be a basis, since it is not guaranteed that the functions are linearly independent.
Let us recall that the forward finite-difference approximation of the derivative of a function f at x = x 0 is given by:
f ( x 0 ) f ( x 0 + h ) f ( x 0 ) h ,
where h > 0 is a small step size. Then, we define the forward finite-difference of f as:
Δ h f ( x ) : = f ( x + h ) f ( x ) ,
and recursively, for  n N , the higher-order difference of f:
Δ h n f ( x ) : = i = 0 n ( 1 ) n i n i f ( x + i h ) .
The relationship of the higher-order differences with the respective derivatives is straightforward, and can be expressed as follows:
f ( n ) ( x ) = lim h 0 1 h n Δ h n f ( x ) .
The next result demonstrates that the determinant of a Wronskian matrix can be expressed as a limit of collocation matrices.
Proposition 1. 
Let F = ( f 1 , f 2 , , f n ) be a basis of C n 1 functions on I R , and  W F ( x ) be the Wronskian matrix of F at x I . Then, we have
det W F ( x ) = lim h 0 + h n ( 1 n ) / 2 f 1 ( x 1 ) f 2 ( x 1 ) f n ( x 1 ) f 1 ( x 2 ) f 2 ( x 2 ) f n ( x 2 ) f 1 ( x n ) f 2 ( x n ) f n ( x n ) ,
for x k : = x + ( k 1 ) h , k = 1 , , n .
Proof. 
By applying elementary properties of determinants, we can express
f 1 ( x 1 ) f n ( x 1 ) f 1 ( x 2 ) f n ( x 2 ) f 1 ( x n ) f n ( x n ) = f 1 ( x 1 ) f n ( x 1 ) Δ h f 1 ( x 1 ) Δ h f n ( x 1 ) Δ h n 1 f 1 ( x 1 ) Δ h n 1 f n ( x 1 ) ,
and then,
f 1 ( x 1 ) f n ( x 1 ) 1 h Δ h f 1 ( x 1 ) 1 h Δ h f n ( x 1 ) 1 h n 1 Δ h n 1 f 1 ( x 1 ) 1 h n 1 Δ h n 1 f n ( x 1 ) = h n ( 1 n ) / 2 f 1 ( x 1 ) f n ( x 1 ) f 1 ( x 2 ) f n ( x 2 ) f 1 ( x n ) f n ( x n ) .
The Formula (11) follows as we take the limit h 0 in (12) and consider (10).    □
With similar reasoning, the previous result can be extended to any minor of W F ( x ) with consecutive rows.
Proposition 2. 
Let F = ( f 1 , f 2 , , f n ) be the basis of C n 1 functions on I R , and  W F ( x ) be the Wronskian matrix of F at x I . Then, for any 1 j i n and 1 α 1 < < α j n ,
det W F ( x ) [ i j + 1 , , i | α 1 , , α j ] = lim h 0 + h j ( 1 j ) / 2 f α 1 ( i j ) ( x 1 ) f α j ( i j ) ( x 1 ) f α 1 ( i j ) ( x 2 ) f α j ( i j ) ( x 2 ) f α 1 ( i j ) ( x j ) f α j ( i j ) ( x j ) ,
where x k : = x + ( k 1 ) h , k = 1 , , n .
Using the previous result, conditions to guarantee the total positivity property of a Wronskian matrix can be derived.
Theorem 2. 
Let F = ( f 1 , f 2 , , f n ) be the basis of C n 1 functions on ( a , b ) and F k ; k = 1 , , n 1 ; the systems are defined in (9). If the Wronskian matrix W F ( x ) at x ( a , b ) is nonsingular, and F k is TP on ( a , b ) , for  k = 0 , , n 1 , then W F ( x ) is TP at x ( a , b ) .
Proof. 
According to Theorem 2.3 of [18], we only have to check that all minors of W F ( x ) with consecutive rows are nonnegative. Let us consider ( x δ , x + δ ) ( a , b ) and h < δ j 1 . The condition that F i j is TP on ( a , b ) for i j = 0 , , n 1 implies that, for any pair ( i , j ) such that 1 j i n , the determinants in (13) at the nodes x k = x + ( k 1 ) h are nonnegative, and remain so when h 0 + .    □
Now, we shall express the initial minors of W F ( x ) in terms of symmetric functions.
Theorem 3. 
Let F = ( f 1 , f 2 , , f n ) be a basis of C n 1 functions on I R , and  W F ( x ) be the Wronskian matrix of F at x I . Then,
det W F ( x ) [ i j + 1 , , i | 1 , , j ] = g F i j , j , j ( x , , x ) r = 1 j 1 r ! , det W F t ( x ) [ i j + 1 , , i | 1 , , j ] = g F , j , i ( x , , x ) r = 1 j 1 r ! ,
where g F i j , j , j and g F , j , i , 1 j i n , are the symmetric functions in (8).
Proof. 
By applying (8) in Formula (13) for the sequence α k = k , and taking the limit h 0 + , we deduce that
det W F ( x ) [ i j + 1 , , i | 1 , , j ] = lim h 0 + h j ( 1 j ) 2 f 1 ( i j ) ( x 1 ) f j ( i j ) ( x 1 ) f 1 ( i j ) ( x 2 ) f j ( i j ) ( x 2 ) f 1 ( i j ) ( x j ) f j ( i j ) ( x j ) = lim h 0 + h j ( 1 j ) 2 det M F i j ( X ) [ 1 , , j ] = lim h 0 + h j ( 1 j ) 2 | V x 1 , , x j | g F i j , j , j ( x 1 , , x j ) = r = 1 j 1 r ! g F i j , j , j ( x , , x ) ,
where in the last equality, we have used (5). Similarly,
det W F t ( x ) [ i j + 1 , , i | 1 , , j ] = f i j + 1 ( x ) f i j + 1 ( j 1 ) ( x ) f i j + 2 ( x ) f i j + 2 ( j 1 ) ( x ) f i ( x ) f i ( j 1 ) ( x ) = lim h 0 + h j ( 1 j ) 2 f i j + 1 ( x 1 ) f i j + 1 ( x j ) f i j + 2 ( x 1 ) f i j + 2 ( x j ) f i ( x 1 ) f i ( x j ) = lim h 0 + h j ( 1 j ) 2 det M F t ( X ) [ i j + 1 , , i | 1 , , j ] = r = 1 j 1 r ! g F , j , i ( x , , x ) .
   □
Starting from this point, the discussion will be narrowed to the polynomial basis.

4. Initial Minors and Total Positivity of Wronskian Matrices for Polynomial Bases

Let P n ( I ) be the space of polynomials of a degree not greater than n, defined on I R , and P = ( p 1 , , p n ) is a basis of P n 1 ( I ) , such that
p i ( x ) = j = 1 n a i , j x j 1 , x I , i = 1 , , n .
Let A = ( a i , j ) 1 i , j n be the change of the basis matrix from the monomial basis M = ( m 1 , , m n ) , with  m i ( x ) = x i 1 ; i = 1 , , n , satisfying
( p 1 , p 2 , , p n ) T = A ( m 1 , m 2 , , m n ) T .
The first derivatives of P form a system P 1 = ( p 1 , p 2 , , p n ) of polynomials of a degree not greater than n 2 are given by
p i ( x ) = j = 1 n 1 ( A 1 ) i , j x j 1 , x I , i = 1 , , n ,
where
( A 1 ) i , j : = j a i , j + 1 , 1 i n , 1 j n 1 , 0 , 1 i n , j = n .
Higher derivative sets P m = ( p 1 ( m ) , p 2 ( m ) , , p n ( m ) ) can be defined in a similar fashion from the matrix A,
p i ( m ) ( x ) = j = 1 n m ( j + m 1 ) ! ( j 1 ) ! a i , j + m x j 1 = j = 1 n m ( A m ) i , j x j 1 , x I , i = 1 , , n ,
with
( A m ) i , j : = ( j + m 1 ) ! ( j 1 ) ! a i , j + m , 1 i n , 1 j n m , 0 , 1 i n , 1 j = n m + 1 , , n .
In [15], it was shown how to express the BD of the collocation matrix M P ( x ) in terms of Schur polynomials and some minors of A. Specifically,
det M P ( x ) [ i j + 1 , , i | 1 , , j ] = | V x i j + 1 , , x i | l 1 > > l j det A [ 1 , , j | l j , , l 1 ] s ( l 1 j , , l j 1 ) ( x i j + 1 , , x i ) .
Schur polynomials are naturally labeled by partitions, and their product rule and other properties are easily stated in terms of them. Thus, for operativeness, it is convenient to write the linear combination appearing in (17) in terms of partitions. We shall consider the partitions λ = ( λ 1 , , λ j ) , where λ r = l r + r j 1 , for  r = 1 , , j . Note that, since l 1 > > l j , λ is a well-defined partition. So, for the minors of A, we shall use the notation
A [ i , λ ] : = det A [ i j + 1 , , i | l j , , l 1 ] = det A [ i j + 1 , , i | λ j + 1 , , λ 1 + j ] .
Be aware that as the dummy variables satisfy l 1 > > l j and l k n , for  1 k n , the partitions they correspond to will have j parts of maximal length n j each. In other words, the sum in (17) must be over all Young diagrams that fit in a j × ( n j ) rectangular box. This fact can be expressed as
l ( λ ) j , λ 1 n j .
Taking into account this notation, we can write
l j < < l 1 det A [ 1 , , j | l j , , l 1 ] s ( l 1 j , , l j 1 ) ( x i j + 1 , , x i ) = l ( λ ) j λ 1 n j A [ j , λ ] s λ ( x i j + 1 , , x i ) ,
and
det M P ( x ) [ i j + 1 , , i | 1 , , j ] = | V x i j + 1 , , x i | l ( λ ) j λ 1 n j A [ j , λ ] s λ ( x i j + 1 , , x i ) .
For polynomial functions, the symmetric functions g P i j , j , j and g P , j , i are
g P i j , j , j ( x 1 , , x j ) = l ( λ ) j λ 1 n j A i j [ j , λ ] s λ ( x 1 , , x j ) , g P , j , i ( x 1 , , x j ) = l ( λ ) j λ 1 n j A [ i , λ ] s λ ( x 1 , , x j )
Using (19), it is possible to find compact formulae for the initial minors of the Wronskian matrix W P ( x ) . They are shown in the next result.
Theorem 4. 
For 1 j i n , we have
det W P ( x ) [ i j + 1 , , i | 1 , , j ] = r = 1 j 1 r ! l ( λ ) j λ 1 n j A i j [ j , λ ] D λ , j x | λ | ,
det W P t ( x ) [ i j + 1 , , i | 1 , , j ] = r = 1 j 1 r ! l ( λ ) j λ 1 n j A [ i , λ ] D λ , j x | λ | ,
where D λ , j are the values defined in (7).
Proof. 
First, we have to evaluate (19) at a single node. The value of a Schur function at a single node can be found by the use of
s λ ( x , , x j ) = x | λ | s λ ( 1 , , 1 j ) = x | λ | D λ , j .
The symmetric functions in (19), evaluated at a single node, take the values
g P i j , j , j ( x , , x ) = l ( λ ) j λ 1 n j A i j [ j , λ ] x | λ | D λ , j , g P , j , i ( x , , x ) = l ( λ ) j λ 1 n j A [ i , λ ] x | λ | D λ , j .
Finally, by inserting the above functions into (14), we obtain the expressions for the minors of the Wronskian matrices.    □
Using (20) and (21), we can find sufficient conditions for the total positivity of the Wronskian matrix W P ( x ) .
Theorem 5. 
Let A be the matrix of change of the basis from the monomial basis of a polynomial basis P and J = ( J i , j ) 1 i , j n , with  J i , j : = ( 1 ) j 1 δ i , j .
(i) 
If A is TP then W P ( x ) is TP for x ( 0 , ) .
(ii) 
If J A J is TP then J W P ( x ) J is TP for x ( , 0 ) .
Proof. 
(i) If A is TP, the matrices A m in (16) are also TP. Now, since all the minors of A m and the quantities D λ , j in (7) are non-negative, from (20) and (21), the initial minors of W P ( x ) are non-negative for x ( 0 , ) . (ii) Let us note that for any matrix M and partition λ = ( λ 1 , , λ j ) , we have
det J M J [ i j + 1 , , i | 1 + λ j , , j + λ 1 ] = ( 1 ) j ( i 1 ) + | λ | det M [ i j + 1 , , i | 1 + λ j , , j + λ 1 ] .
If J A J is TP, then for x ( 0 , ) , we have
0 r = 1 j 1 r ! l ( λ ) j λ 1 n j J A J [ i , λ ] D λ , j x | λ | = r = 1 j 1 r ! l ( λ ) j λ 1 n j ( 1 ) j ( i 1 ) + | λ | A [ i , λ ] D λ , j x | λ | = ( 1 ) j ( i 1 ) det W P t ( x ) [ i j + 1 , , i | 1 , , j ] = det J W P J t ( x ) [ i j + 1 , , i | 1 , , j ] .
Similarly, since J A i j J is TP, we have
0 r = 1 j 1 r ! l ( λ ) j λ 1 n j J A i j J [ j , λ ] D λ , j x | λ | = det J W A J ( x ) [ i j + 1 , , i | 1 , , j ] .
So, J W P ( x ) J is TP at x ( , 0 ) .    □
In Algorithm A1 (see Appendix A), the implementation of the code of Formula (20) is displayed. Note that an algorithm for the implementation of the code of Formula (21) can be obtained similarly.

5. Examples

In order to show the application of Formula (20) for the computation of the initial minors det W P ( x ) [ i j + 1 , , i | 1 , , j ] , let us present two interesting examples.

5.1. Wronskian Matrices of the Bernstein Basis

The Bernstein basis of the space P n 1 [ 0 , 1 ] of polynomials of degree less than or equal to n 1 on the interval [ 0 , 1 ] is B : = ( B 1 n 1 , , B n n 1 ) , with
B i n 1 ( x ) : = n 1 i 1 x i 1 ( 1 x ) n i , x [ 0 , 1 ] , i = 1 , , n .
For n = 4 , ( B 1 3 , B 2 3 , B 3 3 , B 4 3 ) T = A ( 1 , x , x 2 , x 3 ) T , with 
A = 1 3 3 1 0 3 6 3 0 0 3 3 0 0 0 1 ,
and the Wronskian matrix of the Bernstein basis (23) is
W B ( x ) = ( 1 x ) 3 3 x ( 1 x ) 2 3 x 2 ( 1 x ) x 3 3 ( 1 x ) 2 9 x 2 12 x + 3 9 x 2 + 6 x 3 x 2 6 ( 1 x ) 18 x 12 18 x + 6 6 x 6 18 18 6 .
Now, we shall focus on det W B ( x ) [ 2 , 3 , 4 | 1 , 2 , 3 ] to illustrate our proposed procedure for obtaining the initial minors of W B ( x ) .
For computing det W B ( x ) [ 2 , 3 , 4 | 1 , 2 , 3 ] , we have to sum in (20) over all partitions λ , such that l ( λ ) 3 and λ 1 1 . In general, partitions with l ( λ ) j and λ 1 n j are in one-to-one correspondence with integer non-negative decreasing j-tuplas whose first entry does not exceed n j . This observation is taken into account in the present example, and in the implementation of (20) in Algorithm A1 (see Appendix A). In our example, we must consider all non-negative decreasing three-tuplas with initial entry less or equal to 1. Thus, we will be summing over { ( 1 , 1 , 1 ) , ( 1 , 1 , 0 ) , ( 1 , 0 , 0 ) , ( 0 , 0 , 0 ) } . Moreover, from (15) we obtain
A 1 = 3 6 3 0 3 12 9 0 0 6 9 0 0 0 3 0 .
From (18), we can write A 1 [ 3 , λ ] as follows
A 1 [ 3 , ( 1 , 1 , 1 ) ] = det A 1 [ 1 , 1 , 3 | 2 , 3 , 4 ] = 0 , A 1 [ 3 , ( 1 , 1 , 0 ) ] = det A [ 1 , 2 , 3 | 1 , 3 , 4 ] = 0 , A 1 [ 3 , ( 1 , 0 , 0 ) ] = det A 1 [ 1 , 2 , 3 | 1 , 2 , 4 ] = 0 , A 1 [ 3 , ( 0 , 0 , 0 ) ] = det A 1 [ 1 , 2 , 3 | 1 , 2 , 3 ] = 54 .
Taking into account (7), we can derive
D ( 1 , 1 , 1 ) , 3 = 1 , D ( 1 , 1 , 0 ) , 3 = 3 , D ( 1 , 0 , 0 ) , 3 = 3 , D ( 0 , 0 , 0 ) , 3 = 1 ,
and
x | ( 1 , 1 , 1 ) | = x 3 , x | ( 1 , 1 , 0 ) | = x 2 , x | ( 1 , 0 , 0 ) | = x 1 , x | ( 0 , 0 , 0 ) | = x 0 .
Finally, from (20), we obtain
det W B ( x ) [ 2 , 3 , 4 | 1 , 2 , 3 ] = 1 ! 2 ! ( A 1 [ 3 , ( 1 , 1 , 1 ) ] D ( 1 , 1 , 1 ) , 3 x | ( 1 , 1 , 1 ) | + A 1 [ 3 , ( 1 , 1 , 0 ) ] D ( 1 , 1 , 0 ) , 3 x | ( 1 , 1 , 0 ) | + A 1 [ 3 , ( 1 , 0 , 0 ) ] D ( 1 , 0 , 0 ) , 3 x | ( 1 , 0 , 0 ) | + A 1 [ 3 , ( 0 , 0 , 0 ) ] D ( 0 , 0 , 0 ) , 3 x | ( 0 , 0 , 0 ) | ) = 108 .
Following the same reasoning, we can obtain any other initial minor.

5.2. Polynomial Recursive Bases

The example of polynomial recursive bases shows that Equations (20) and (21) can be specially useful for the computation of the initial minors, if the structure of the corresponding matrix of the change of the basis permits a systematic computation of its minors with consecutive rows.
Given values b 1 , , b n , such that b i > 0 , i = 1 , , n , let us define the polynomial recursive basis P : = ( p 1 , , p n ) as
p i = j = 1 i b j x j 1 , i = 1 , , n .
The change of the basis matrix such that ( p 1 , , p n ) T = A ( m 1 , , m n ) T , with  m i ( x ) : = x i 1 , i = 1 , , n , is a nonsingular lower triangular and TP matrix of the following form:
A = b 1 0 0 0 b 1 b 2 0 0 b 1 b 2 b 3 0 0 b 1 b 2 b 3 b n .
So, the basis ( p 1 , , p n ) is TP for x [ 0 , ) . Therefore, by virtue of Theorem 5, the Wronskian matrix W P ( x ) is TP for x [ 0 , ) .
The matrix A 1 is
A 1 = 0 0 0 0 0 b 2 0 0 0 0 b 2 2 b 3 0 0 0 b 2 2 b 3 3 b 4 ( n 1 ) b n 0 ,
and, similarly, by using (16), we obtain A m with m = 1 , , n .
Let us note that the only non-zero minors of A m are
A m [ i j + 1 , , i | l , i j + 2 m , , i m ] = C i , j , m b l + m k = 2 j b i j + k ,
for l = 1 , , i j m + 1 , where A 0 : = A and 
C i , j , m = ( l + m 1 ) m k = 1 j 1 ( i j + k ) m , i j m
with the Pochhammer symbol for descending factorials ( j ) m : = j ( j 1 ) ( j m + 1 ) .
  • Specifically, for the case m = i j , which is relevant in (20), we have
    A i j , [ j , λ ] = det A i j [ 1 , , j | l , 2 , , j ] , l = 1 .
So, the only partition which contributes to (20) is λ = Ø , and 
det W P ( x ) [ i j + 1 , , i | 1 , , j ] = r = 1 j 1 r ! k = 1 j b i j + k ,
which does not depend on x.
In order to compute (21), let us note that, because of (24), the only contribution to the sum in the RHS of (21) comes from the partitions λ = ( λ 1 , , λ j ) with
λ r = i j , r = 1 , , j 1 , λ j = 0 , 1 , , i j ,
Let us call
λ l : = ( i j , , i j j 1 , l 1 ) , l = 1 , , i j + 1 .
Applied on these partitions, Formula (7) reduces to
D λ l , j = i l j 1 .
Thus, finally, we obtain
det W P t ( x ) [ i j + 1 , , i | 1 , , j ] = r = 1 j 1 r ! k = 2 j b i j + k l = 1 i j + 1 b l i l j 1 x ( i j ) ( j 1 ) + l 1 .

6. Conclusions

The elements of the bidiagonal decomposition of a totally positive collocation matrix can be expressed in terms of symmetric functions of the nodes. By examining the relationship between the Wronskian and collocation matrices of a given TP basis of functions, we have developed a method to express the entries of the BD of Wronskian matrices. These entries can be calculated as the values of certain symmetric functions evaluated at a single node. Furthermore, in the case of polynomial bases, we have derived explicit formulas for the entries of the BD of their Wronskian matrices, facilitating a deeper understanding and simpler computation in applications involving TP matrices.

Author Contributions

Conceptualization, P.D., E.M. and B.R.; methodology, P.D., E.M. and B.R.; investigation, P.D., E.M. and B.R.; writing—original draft, P.D., E.M. and B.R.; writing—review and editing, P.D., E.M. and B.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by Spanish research grants PID2022-138569NB-I00 (MCI/AEI) and RED2022-134176-T (MCI/AEI) and by Gobierno de Aragón (E41_23R).

Data Availability Statement

The data and codes used in this work are available under request.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

In Algorithm A1, the implementation of the code of Formula (20) is displayed.
Algorithm A1: Pseudocode for the computation of Formula (20)
   Require: 
i,j,A
   Ensure: 
det W A ( x ) [ i j + 1 , , i | 1 , , j ]
      1:
n = size(A,1);
      2:
m = i − j;
      3:
Am = Am(A,i,j);
      4:
parts = particiones(j,n);
      5:
suma = 0;
      6:
   for k = 1:size(parts,1)
      7:
      lambda = parts(k,:);
      8:
      a = Alambda(Am, j, lambda);
      9:
      d = Dlambdaj(lambda,3);
    10:
      suma = suma + a*d* x s u m ( l a m b d a ) ;
    11:
   end
    12:
prod = 1;
    13:
   for r = 1:j−1
    14:
      prod = prod * factorial(r);
    15:
   end
    16:
total = prod * suma
 
    17:
function part = partitionsR(from, level)
    18:
part = [];
    19:
   for value = from: −1:0
    20:
      if level > 1
    21:
         res = partitionsR(min(from,value),level-1);
    22:
         part = [part; [value .∗ ones(size(res,1),1) res]];
    23:
      else
    24:
         part = [part; value];
    25:
      end
    26:
   end
 
    27:
function partitions = partitions(j, n)
    28:
partitions = partitionsR(n−j,j);
 
    29:
function Alambda = Alambda(A, j, lambda)
    30:
rows = 1:j;
    31:
cols = flip(lambda) + (1:j);
    32:
Alambda = det(A(rows,cols));
 
    33:
function Am = Am(A,i,j)
    34:
m = i-j;
    35:
n = size(A,1);
    36:
r = (m+1):n;
    37:
Am = [A(:,r), zeros(n,m)];
    38:
   for j=1:n-m
    39:
      factorj = factorial(j + m−1)/factorial(j−1);
    40:
      Am(:,j) = Am(:,j).*factorj;
    41:
   end
 
    42:
function d = Dlambdaj(lambda,j)
    43:
d = 1;
    44:
   for kl=nchoosek(1:j,j−1)’
    45:
      k = kl(1);
    46:
      l = kl(2);
    47:
      factor = (lambda(k) − lambda(l) + l − k)/(l − k);
    48:
      d = d * factor;
    49:
   end

References

  1. Ando, T. Totally positive matrices. Linear Algebra Its Appl. 1987, 90, 165–219. [Google Scholar] [CrossRef]
  2. Fallat, S.M.; Johnson, C.R. Totally Nonnegative Matrices; Princeton Series in Applied Mathematics; Princeton University: Princeton, NJ, USA, 2011. [Google Scholar]
  3. Pinkus, A. Totally Positive Matrices; Cambridge Tracts in Mathematics; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  4. Gasca, M.; Peña, J.M. Total positivity and Neville elimination. Linear Algebra Its Appl. 1992, 165, 25–44. [Google Scholar] [CrossRef]
  5. Gasca, M.; Peña, J.M. A matricial description of Neville elimination with applications to total positivity. Linear Algebra Its Appl. 1994, 202, 33–53. [Google Scholar] [CrossRef]
  6. Gasca, M.; Peña, J.M. On factorizations of totally positive matrices. In Total Positivity and Its Applications; Gasca, M., Micchelli, C.A., Eds.; Kluver Academic Publishers: Dordrecht, The Netherlands, 1996; pp. 109–130. [Google Scholar]
  7. Fekete, M. Über ein Problem von Laguerre. Rend. Circ. Mat. Palermo 1913, 34, 89–120. [Google Scholar] [CrossRef]
  8. Karlin, S. Total Positivity; Stanford U.P: Stanford, CA, USA, 1968; Volume 1. [Google Scholar]
  9. Marco, A.; Martínez, J.J. Accurate computation of the Moore–Penrose inverse of strictly totally positive matrices. J. Comput. Appl. Math. 2019, 350, 299–308. [Google Scholar] [CrossRef]
  10. Delgado, J.; Orera, H.; Peña, J.M. Accurate computations with Laguerre matrices. Numer. Linear Algebra Appl. 2019, 26, e2217. [Google Scholar] [CrossRef]
  11. Delgado, J.; Orera, H.; Peña, J.M. Accurate algorithms for Bessel matrices. J. Sci. Comput. 2019, 80, 1264–1278. [Google Scholar] [CrossRef]
  12. Marco, A.; Martínez, J.J.; Peña, J.M. Accurate bidiagonal decomposition of totally positive Cauchy–Vandermonde matrices and applications. Linear Algebra Its Appl. 2017, 517, 63–84. [Google Scholar] [CrossRef]
  13. Demmel, J.; Koev, P. The accurate and efficient solution of a totally positive generalized Vandermonde linear system. SIAM J. Matrix Anal. Appl. 2005, 27, 42–52. [Google Scholar] [CrossRef]
  14. Demmel, J.; Koev, P. Accurate SVDs of polynomial Vandermonde matrices involving orthonormal polynomials. Linear Algebra Its Appl. 2006, 27, 382–396. [Google Scholar] [CrossRef]
  15. Diaz, P.; Mainar, E.; Rubio, B. Polynomial total positivity and high relative accuracy through Schur polynomials. J. Sci. Comput. 2023, 97, 10. [Google Scholar] [CrossRef]
  16. Yang, Z.; Huang, R.; Zhu, W. Accurate computations for eigenvalues of products of Cauchy-polynomial-Vandermonde matrices. Numer. Algorithms 2020, 85, 329–351. [Google Scholar] [CrossRef]
  17. Macdonald, I.G. Symmetric Function and Hall Polynomials, 2nd ed.; Oxford University Press: Oxford, UK, 1998. [Google Scholar]
  18. Peña, J.M. Determinantal criteria for total positivity. Linear Algebra Its Appl. 2001, 332–334, 131–137. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Díaz, P.; Mainar, E.; Rubio, B. Totally Positive Wronskian Matrices and Symmetric Functions. Axioms 2024, 13, 589. https://doi.org/10.3390/axioms13090589

AMA Style

Díaz P, Mainar E, Rubio B. Totally Positive Wronskian Matrices and Symmetric Functions. Axioms. 2024; 13(9):589. https://doi.org/10.3390/axioms13090589

Chicago/Turabian Style

Díaz, Pablo, Esmeralda Mainar, and Beatriz Rubio. 2024. "Totally Positive Wronskian Matrices and Symmetric Functions" Axioms 13, no. 9: 589. https://doi.org/10.3390/axioms13090589

APA Style

Díaz, P., Mainar, E., & Rubio, B. (2024). Totally Positive Wronskian Matrices and Symmetric Functions. Axioms, 13(9), 589. https://doi.org/10.3390/axioms13090589

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop