Next Article in Journal
On the Periodicity of General Class of Difference Equations
Next Article in Special Issue
Trapezium-Type Inequalities for an Extension of Riemann–Liouville Fractional Integrals Using Raina’s Special Function and Generalized Coordinate Convex Functions
Previous Article in Journal
Natural Paracontact Magnetic Trajectories on Unit Tangent Bundles
Previous Article in Special Issue
Quantum Trapezium-Type Inequalities Using Generalized ϕ-Convex Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximating Functions of Positive Compact Operators by Using Bell Polynomials

by
Diego Caratelli
1,2,
Pierpaolo Natalini
3 and
Paolo Emilio Ricci
4,*
1
Department of Research and Development, The Antenna Company, High Tech Campus 29, 5656 AE Eindhoven, The Netherlands
2
Department of Electrical Engineering, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven, The Netherlands
3
Department of Mathematics and Physics, Roma Tre University, Largo San Leonardo Murialdo 1, 00146 Rome, Italy
4
Department of Mathematics, International Telematic University UniNettuno, Corso Vittorio Emanuele II 39, 00186 Rome, Italy
*
Author to whom correspondence should be addressed.
Axioms 2020, 9(3), 73; https://doi.org/10.3390/axioms9030073
Submission received: 10 June 2020 / Revised: 23 June 2020 / Accepted: 24 June 2020 / Published: 30 June 2020
(This article belongs to the Special Issue Special Functions and Their Applications)

Abstract

:
After recalling the most important properties of the Bell polynomials, we show how to approximate a positive compact operator by a suitable matrix. Then, we derive a representation formula for functions of the obtained matrix, which can be considered as an approximate value for the functions of the corresponding operator.

Graphical Abstract

1. Introduction

The Bell polynomials [1] are a standard mathematical tool for computing the nth derivative of a composite function. They have also been applied in order to solve different problems as the Blissard problem (see [2], p. 46), the representation of Lucas polynomials of the first and second kind [3], the representation formulas for the elementary symmetric functions of a countable set of numbers [4]. In this article, after recalling the most important properties of the Bell polynomials, we will show another possible application that is their connection with the orthogonal invariants of a positive compact operator (shortly PCO).
The theory of the orthogonal invariants for a PCO was introduced by Fichera in order to approximate the eigenvalues of linear elliptic differential problems satisfying suitable conditions [5,6]. An important tool in this framework is given by the so-called Robert’s formulas which permit to reduce the order (or the degree) of orthogonal invariants. In [4], it was proven that the Robert’s formulas are nothing but the recurrence relation and the Faà di Bruno formula for the Bell polynomials. Further relations among orthogonal invariants, proven in the same article, are reported in Section 5. All equations recalled there generalize the algebraic Newton–Girard formulas to the elementary symmetric functions of a countable set of numbers, which seems to be a useful result, since, in Section 7, the results on matrix functions derived by using the Dunford–Taylor (also called Riesz–Fantappiè) integral [7] are extended for approximating matrix functions of strictly positive compact operators.
The first part of this article contains a survey about the use of Bell polynomials in connection with the above-mentioned theory of orthogonal invariants. We think that the results achieved in [4] are necessary for the understanding of what is subsequently exposed. It should also be noted that the Journal on which this paper has been published is not accessible on the web.
It is well known that positive compact operators, according to Fredholm’s theory, are the limit in the norm of finite dimension operators, which can be represented by positive definite matrices. This fact has been used several times in literature (see, for example, [8,9,10] and the references therein), but referring to the eigenvalues of the matrices involved. In this article, instead, exploiting the theory of the orthogonal invariants, we have used matrix entries (avoiding the knowledge of eigenvalues), Bell’s polynomials, and Robert’s formulas. This is, to our knowledge, an unusual approach.
A few examples of the effectiveness of this procedure are shown in Section 8, by using the Mathematica © computer algebra program.

2. Recalling the Bell Polynomials

Consider the composite function Φ ( t ) : = f ( g ( t ) ) , where x = g ( t ) and y = f ( x ) , are differentiable functions (up to the considered order), defined in suitable intervals of the real axis, so that Φ ( t ) can be differentiated n times with respect to t, by using the chain rule.
Here, and in what follows, we use the notations:
Φ m : = D t m Φ ( t ) , f h : = D x h f ( x ) | x = g ( t ) , g k : = D t k g ( t ) .
Then, the n-th derivative of Φ ( t ) is represented by
Φ n = Y n ( f 1 , g 1 ; f 2 , g 2 ; ; f n , g n ) ,
where the Y n denote the Bell polynomials.
For example, one has:
Y 1 ( f 1 , g 1 ) = f 1 g 1 Y 2 ( f 1 , g 1 ; f 2 , g 2 ) = f 1 g 2 + f 2 g 1 2 Y 3 ( f 1 , g 1 ; f 2 , g 2 ; f 3 , g 3 ) = f 1 g 3 + f 2 ( 3 g 2 g 1 ) + f 3 g 1 3 .
Further examples can be found in [2], p. 49.
Proposition 1.
The Bell polynomials satisfy the recurrence relation:
{ Y 0 : = f 1 ; Y n + 1 ( f 1 , g 1 ; ; f n , g n ; f n + 1 , g n + 1 ) = = k = 0 n ( n k ) Y n k ( f 2 , g 1 ; f 3 , g 2 ; ; f n k + 1 , g n k ) g k + 1 .
An explicit expression for the Bell polynomials is given by the Faà di Bruno formula:
Φ n = Y n ( f 1 , g 1 ; f 2 , g 2 ; ; f n , g n ) = = π ( n ) n ! r 1 ! r 2 ! r n ! f r [ g 1 1 ! ] r 1 [ g 2 2 ! ] r 2 [ g n n ! ] r n ,
where the sum runs over all partitions π ( n ) of the integer n, r i denotes the number of parts of size i, and r = r 1 + r 2 + + r n denotes the number of parts of the considered partition.
The proof of the Faà di Bruno formula can be found in the Riordan book [2]. See also [11,12], where the proof is based on the umbral calculus.
However, the Faà di Bruno is not suitable by the computational point of view, owing the higher complexity with respect to the recursion formula.
The traditional form of the Bell polynomials [13] is given by:
Y n ( f 1 , g 1 ; f 2 , g 2 ; ; f n , g n ) = h = 0 n k B n , k ( g 1 , g 2 , , g n k + 1 ) f k
and the B n , k coefficients satisfy the recursion:
B n , k ( g 1 , g 2 , , g n k + 1 ) = h = 0 n k ( n 1 h 1 ) B n h 1 , k ( g 1 , g 2 , , g n h k + 1 ) g h + 1 .

3. An Extension of the Newton–Girard Formulas

In our opinion, one of the most important applications of the Bell polynomials is the following extension of the Newton–Girard formulas (see e.g., [14]):
Consider the (finite or infinite) sequence of real or complex numbers μ 1 , μ 2 , μ 3 , , and denote by
σ 1 = i μ i , σ 2 = i < j μ i μ j ,     ,     σ k = i 1 < i 2 < < i k μ i 1 μ i 2 μ i k ,     ,
the relevant elementary symmetric functions, and by
s 1 = i μ i ,         s 2 = i μ i 2 ,     ,     s k = i μ i k ,     ,
the symmetric functions power sums.
Of course, in case of infinite sets of indices, we assume that all of the above considered expansions are convergent.
In [4], the following result is proven:
Proposition 2.
For any integer k, the following representation formulas hold true:
σ k = ( 1 ) k k ! Y k ( 1 , s 1 ; 1 , s 2 ; 1 , 2 ! s 3 ; ; 1 , ( k 1 ) ! s k )
s k = 1 ( k 1 ) ! Y k ( 1 , σ 1 ; 1 , 2 ! σ 2 ; ; ( 1 ) k 1 ( k 1 ) ! , ( 1 ) k k ! σ k )
The above Formulas (6)–(7) constitute an extension of the Newton–Girard formulas and their inverse, since we have, in particular:
σ 1 = s 1 σ 2 = 1 2 ( s 1 2 s 2 ) σ 3 = 1 6 ( s 1 3 3 s 1 s 3 + 2 s 3 ) σ 4 = 1 24 ( s 1 4 6 s 1 2 s 2 + 8 s 1 s 3 + s 2 2 6 s 4 ) s 1 = σ 1 s 2 = σ 1 2 2 σ 2 s 3 = σ 1 3 3 σ 1 σ 2 + 3 σ 3 s 4 = σ 1 4 4 σ 1 2 σ 2 + 4 σ 1 σ 3 + 2 σ 2 2 4 σ 4
Remark 1.
Note that the above equation implies the following well known recursion:
σ 1 = s 1 σ 2 = 1 2 ( σ 1 s 1 s 2 ) σ 3 = 1 3 ( σ 2 s 1 σ 1 s 2 + s 3 ) σ 4 = 1 4 ( σ 3 s 1 σ 2 s 2 + σ 1 s 3 s 4 )
which permits recovering the elementary symmetric functions.

4. Orthogonal Invariants of PCO and Robert’s Formulas

The eigenvalues μ k of a positive compact operator T in a complex Hilbert space H satisfy
μ 1 μ 2 μ 3 0
and, when infinite many eigenvalues exist, zero is an accumulation point.
If the PCO is strictly positive, the above condition becomes:
μ 1 μ 2 μ 3 > 0 .
A classical example of eigenvalue problem for a strictly positive operator is given by
T ϕ = K ϕ : = A K ( x , y ) ϕ ( y ) d y = μ ϕ ( x ) ,
where the kernel K ( x , y ) of the second kind Fredholm operator K belongs to L 2 ( A × A ) , and is such that K ( x , y ) = K ( y , x ) ¯ , ( K ϕ , ϕ ) > 0 if ϕ 0 L 2 ( A ) (see Mikhlin [15] and Tricomi [16]).
The orthogonal invariants, introduced by Fichera [5,6], are, by definition, symmetric functions of the eigenvalues of T:
I s n ( T ) = k 1 < k 2 < < k s [ μ k 1 μ k 2 μ k s ] n .
The index s is called the order, and the exponent n the degree of the invariant.
In the infinite dimensional case, a connection with the basic invariants I s 1 ( T ) ( s = 1 , 2 , , n ) , or I 1 n ( T ) ( n = 1 , 2 , , s ) have been derived by Robert [17].
The relative formulas are as follows:
I s n ( T ) = 1 s q = 1 s ( 1 ) q 1 I 1 q n ( T ) I s q n ( T )
I s n ( T ) = ( 1 ) s k = 1 s ( 1 ) k k ! 1 q i s q 1 + + q k = s I 1 n q 1 ( T ) I 1 n q k ( T ) q 1 q k ,
which allow for reducing the orthogonal invariant I s n ( T ) to I 1 h ( T ) ( h = 1 , 2 , , n s ) .
Since the eigenvalues of T n are given by μ i n , in the following, denoting by n the smallest integer such that I s n ( T ) < , we will put T : = T n , so that I s 1 ( T ) = I s 1 ( T n ) = I s n ( T ) , and the above Equations (13) and (14) become:
I s 1 ( T ) = 1 s q = 1 s ( 1 ) q 1 I 1 q ( T ) I s q 1 ( T )
I s 1 ( T ) = ( 1 ) s k = 1 s ( 1 ) k k ! 1 q i s q 1 + + q k = s I 1 q 1 ( T ) I 1 q k ( T ) q 1 q k .
Under the above condition, the operator T is called to belong to the Trace-class [18], p. 521.
Remark 2.
It is worth noting that there exist PCO not satisfying the above-mentioned condition that requires the existence of an integer n such that I s n ( T ) = I s 1 ( T ) < [5,6]; however, this condition is satisfied by all the PCO occurring in applications.
It is well known that, in the Hilbert space H = L 2 ( A ) , the following results hold:
n , s N , I s n ( T ) < i f f T n φ = A K ( x , y ) φ ( y ) d y
where
K ( x , y ) = A H ( x , z ) H ( z , y ) d z , H ( x , y ) = H ( y , x ) ¯ L 2 [ A × A ] .
Furthermore, the orthogonal invariants can be represented by the multiple integrals [6]:
I s n ( T ) = 1 s ! A A f ( x 1 , , x s ) d x 1 d x s ,
where f ( x 1 , , x s ) denotes the Fredholm determinant:
f ( x 1 , x 2 , , x s ) : = | K ( x 1 , x 1 ) K ( x 1 , x 2 ) K ( x 1 , x s ) K ( x 2 , x 1 ) K ( x 2 , x 2 ) K ( x 2 , x s ) K ( x s , x 1 ) K ( x s , x 2 ) K ( x s , x s ) | .
In particular, for s = 1 we have:
I 1 n ( T ) = A K ( x , x ) d x = A × A | H ( x , y ) | d x d y ,
and, for s = 2 :
I 2 n ( T ) = 1 2 A × A K ( x , x ) K ( y , y ) | K ( x , y ) | 2 d x d y .
In what follows, we consider only the operators T that belong to the Trace class. We put T = T n , so that it results in:
I 1 n ( T ) = I 1 1 ( T n ) = I 1 1 ( T ) = A K ( x , x ) d x < 1 p t 26 p t I 2 n ( T ) = I 2 1 ( T ) = 1 2 A × A K ( x , x ) K ( y , y ) | K ( x , y ) | 2 d x d y .
As a consequence, all the invariants are bounded.

5. Orthogonal Invariants’ Reduction Formulas

Writing the representation formulas of Proposition 3.1 in terms of orthogonal invariants, we obtain:
I k 1 ( T ) = ( 1 ) k k ! Y k ( 1 , I 1 1 ( T ) ; 1 , I 1 2 ( T ) ; 1 , 2 ! I 1 3 ( T ) ; ; 1 , ( k 1 ) ! I 1 k ( T ) )
I 1 k ( T ) = = 1 ( k 1 ) ! Y k ( 1 , I 1 1 ( T ) ; 1 , 2 ! I 2 1 ( T ) ; ; ( 1 ) k 1 ( k 1 ) ! , ( 1 ) k k ! I k 1 ( T ) ) .
In [4], the following results have been proven:
Proposition 3.
The first Robert formula is equivalent to the recurrence relation of the Bell polynomials.
Proposition 4.
The second Robert formula is equivalent to the Faà di Bruno representation formula for the Bell polynomials.
In particular, from the above Equation (17), we find:
I 2 1 ( T ) = 1 2 [ ( I 1 1 ( T ) ) 2 I 1 2 ( T ) ] I 3 1 ( T ) = 1 6 [ ( I 1 1 ( T ) ) 3 3 I 1 1 ( T ) I 1 2 ( T ) + 2 I 1 3 ( T ) ] I 4 1 ( T ) = 1 24 [ ( I 1 1 ( T ) ) 4 6 ( I 1 1 ( T ) ) 2 ( I 1 2 ( T ) ) + 8 ( I 1 1 ( T ) ) ( I 1 3 ( T ) ) + 3 ( ( I 1 2 ( T ) ) 2 6 ( I 1 4 ( T ) ) ] .
Remark 3.
Note that, recalling the representation in terms of seconk kind Bell polynomials, Equation (17) becomes:
I k 1 ( T ) = ( 1 ) k k ! B n , k ( I 1 1 ( T ) , I 1 2 ( T ) , 2 ! I 1 3 ( T ) , , ( k 1 ) ! I 1 k ( T ) )
so that the I k 1 ( T ) can be computed by using the recursion (3).
Proposition 5.
For any integer k 1 , the orthogonal invariant I 1 k ( T ) , is expressed in terms of I h 1 ( T ) , k = 1 , 2 , , s , by
I 1 k ( T ) = π ( k ) ( 1 ) k + s k ( s 1 ) ! r 1 ! r 2 ! r k ! ( I 1 1 ( T ) ) r 1 ( I 2 1 ( T ) ) r 2 ( I k 1 ( T ) ) r k ,
where π ( k ) denotes the sum running on all partitions of k = r 1 + 2 r 2 + + k r k and s = r 1 + r 2 + + r k .
In particular, we have:
I 1 2 ( T ) = ( I 1 1 ( T ) ) 2 2 I 2 1 ( T ) I 1 3 ( T ) = ( I 1 1 ( T ) ) 3 3 I 1 1 ( T ) I 2 1 ( T ) + 3 I 3 1 ( T ) I 1 4 ( T ) = ( I 1 1 ( T ) ) 4 4 ( I 1 1 ( T ) ) 2 I 2 1 ( T ) + 4 I 1 1 ( T ) I 3 1 ( T ) + 2 ( I 2 1 ( T ) ) 2 4 I 4 1 ( T ) .
Note that Equation (19), making use of Bell polynomials, is more convenient by the computational point of view with respect to Equation (20).
The recursion (7) becomes:
I 1 1 ( T ) = I 1 1 ( T ) I 2 1 ( T ) = 1 2 [ ( I 1 1 ( T ) ) 2 I 1 2 ( T ) ] I 3 1 ( T ) = 1 3 [ I 2 1 ( T ) I 1 1 ( T ) I 1 1 ( T ) I 1 2 ( T ) + I 1 3 ( T ) ] I 4 1 ( T ) = 1 4 [ I 3 1 ( T ) I 1 1 ( T ) + I 2 1 ( T ) I 1 2 ( T ) + I 1 1 ( T ) I 1 3 ( T ) I 1 4 ( T ) ]

6. Matrix Functions

We consider a r × r matrix A = { a h , k } , whose invariants
{ u 1 : = t r A = a 11 + a 22 + + a r r u 2 : = i < j 1 , r | a i i a i j a j i a j j | u r : = d e t A
are known, so that its characteristic polynomial is given by
P ( λ ) = λ r u 1 λ r 1 + u 2 λ r 2 + + ( 1 ) r u r .
Let f be a holomorphic in an open set Δ , containing all the eigenvaues of A . Then, the matrix functions f ( A ) are given by the Dunford–Taylor (also called Riesz–Fantappiè) integral [18,19]:
f ( A ) = 1 2 π i [ k = 1 r γ f ( λ ) h = 0 k 1 ( 1 ) h u h λ k h 1 P ( λ ) d λ A r k ]
where γ denotes a simple contour enclosing all the zeros of P ( λ ) .

7. Approximation of a Strictly PCO

Equation (24) is used in what follows, since we approximate a strictly positive compact operator T by a r × r matrix A , where the integer r is such that the rth eigenvalue μ r of T is negligible, i.e.,
μ r < ε
where ε is a prescribed small quantity.
We assume, for the characteristic polynomial, the following form:
P ( λ ) = λ r I 1 1 ( T ) λ r 1 + + ( 1 ) r 1 I r 1 1 ( T ) λ + ( 1 ) r I r 1 ( T ) .
Therefore, the companion matrix of A is given by
T A = ( I 1 1 ( T ) 1 0 0 0 I 2 1 ( T ) 0 1 0 0 I 3 1 ( T ) 0 0 1 0 ( 1 ) r 2 I r 1 1 ( T ) 0 0 0 1 ( 1 ) r 1 I r 1 ( T ) 0 0 0 0 ) .
and Equation (24) becomes:
f ( A ) = 1 2 π i [ k = 1 r γ f ( λ ) h = 0 k 1 ( 1 ) h I h 1 ( T ) λ k h 1 P ( λ ) d λ A r k ]
Note that the powers of the companion matrix A can be computed by using the results in [20], Theorem 2.

8. Numerical Examples

8.1. Kernels of Positive Compact Operators

In the literature, there are a few classical kernels of PCO reported, such as:
  • A kernel considered by S.G. Miklin (see [21], p. 364):
    K ( x , y ) = { x ( 1 y ) if 0 x y 1 y ( 1 x ) if 0 y x 1
  • A kernel considered by Ostrowski (see [22,23]):
    K ( x , y ) = { 2 log y ( 1 x ) x y if 0 x y 1 2 log x ( 1 y ) x y if 0 y x 1
  • A kernel considered by Tricomi (see [16]:
    K ( x , y ) = cos n π ( x y ) ,   n = 1 , 2 ,
    or, more generally, for a fixed N:
    K ( x , y ) = n = 0 N a n cos n π ( x y ) , a n > 0 .
    The kernel # 2 is considered as an example in what follows.

8.1.1. Example 1

Let us consider the Ostrowski’s kernel K ( x , y ) in (30). As demonstrated by means of the inverse iteration method in [22,23], the relevant operator T displays the eigenvalues:
λ 1 1 . 205965762856735 λ 2 0 . 734006078303625 λ 3 0 . 528546963103265 λ 4 0 . 413767155184915 λ 5 0 . 342859803772745
and, therefore, can be approximated by a matrix of order r = 5 with a characteristic polynomial:
P ( λ ) = k = 1 r λ λ k .
Starting from (34), the invariants are readily evaluated as:
u 1 = I 1 1 ( T ) 3 . 225145763221285 u 2 = I 2 1 ( T ) 3 . 920164426948163 u 3 = I 3 1 ( T ) 2 . 263632211150296 u 4 = I 4 1 ( T ) 0 . 625036213142174 u 5 = I 5 1 ( T ) 0 . 066372901341124
so that the companion matrix A can be written as:
Axioms 09 00073 i001
By making use of the Dunford–Taylor integral Formula (24), the matrix function:
f ( A ) = W ( A )
with W ( · ) denoting the product logarithm, can be computed as:
Axioms 09 00073 i002

8.1.2. Example 2

Let us consider the kernel K ( x , y ) relevant to the positive compact operator T which describes the transverse vibrations of a wedge-shaped beam [24]. By using the inverse iteration method, one can prove that T features the following eigenvalues:
λ 1 1 . 9717553675 × 10 2 λ 2 1 . 1930198795 × 10 3 λ 3 2 . 3684105195 × 10 4 λ 4 7 . 5156070990 × 10 5 λ 5 3 . 0836121615 × 10 5 λ 6 1 . 4889262790 × 10 5 λ 7 8 . 0483648760 × 10 6 λ 8 4 . 7300057890 × 10 6 λ 9 2 . 9688104105 × 10 6 λ 10 1 . 9692639370 × 10 6 λ 11 1 . 3727098565 × 10 6 λ 12 1 . 0037480315 × 10 6 λ 13 7 . 6076195215 × 10 7
and, therefore, can be approximated by a matrix of order r = 13 with characteristic polynomial:
P ( λ ) = k = 1 r ( λ λ k ) .
Starting from (38), the invariants are readily evaluated as:
u 1 = I 1 1 ( T ) 2 . 128914972669765 × 10 2 u 2 = I 2 1 ( T ) 3 . 147982936967773 × 10 5 u 3 = I 3 1 ( T ) 9 . 746686526750313 × 10 9 u 4 = I 4 1 ( T ) 9 . 799739677748166 × 10 13 u 5 = I 5 1 ( T ) 4 . 058777349187134 × 10 17 u 6 = I 6 1 ( T ) 8 . 015136142635608 × 10 22 u 7 = I 7 1 ( T ) 8 . 297854829897979 × 10 27 u 8 = I 8 1 ( T ) 4 . 787022852309753 × 10 32 u 9 = I 9 1 ( T ) 1 . 593001334376288 × 10 37 u 10 = I 1 1 0 ( T ) 3 . 084794075371912 × 10 43 u 11 = I 1 1 1 ( T ) 3 . 391992870896023 × 10 49 u 12 = I 1 1 2 ( T ) 1 . 945386222042645 × 10 55 u 13 = I 1 1 3 ( T ) 4 . 485002305839361 × 10 62
so that the companion matrix A can be written as:
Axioms 09 00073 i003
By making use of the Dunford–Taylor integral Formula (24), the matrix function:
f ( A ) = e 10 A A 10 + 1
can be computed as:
Axioms 09 00073 i004

8.2. Kernels of Inverse of Differential PCO

In [5], many problems relevant to strictly positive definite differential operators have been considered. In particular, for the one connected with the free vibrations of a clamped plate, large numbers of eigenvalues were determined by using the orthogonal invariants method. It is also recalled there that the Fredholm integer function associated with a given kernel K ( x , y ) , where T ϕ = A K ( x , y ) ϕ ( y ) d y = μ ϕ ( x ) is given by:
G ( μ , K ) = 1 + k = 1 ( 1 ) k I k 1 ( T ) μ k .
Putting λ = μ 1 , and considering the partial sum, of order r, of the series (30), we recover the polynomial P ( λ ) in Equation (26) as an approximation of the Fredholm integer function (41).
In [5], Fichera, considering the biharmonic operator
Δ 4 u ( x , y ) = Δ 2 Δ 2 u = 4 u x 4 + 2 4 u x 2 y 2 + 4 u y 4
was able to find lower and upper approximations of the eigenvalues of the biharmonic problem
{ Δ 4 u ( x , y ) = λ u ( x , y ) ( in D ) u = u x = u y = 0 ( on D ) .
both for the circular and the square clamped plate.
In what follows, we choose the mean value, as the most probable one, of the very close approximation of the eigenvalues in [5]. Then, considering the reciprocal of these numbers, we apply the results of Section 7 to the inverse operator of the problem (43), which is the relevant Green function, finding, as examples, a few elementary functions of this operator.

8.2.1. Example 3

In the case of a square clamped plate with D { 0 x π 0 y π } , the operator T which describes the relevant free vibrations features the eigenvalues:
λ 1 0 . 07522320605817612 λ 2 = λ 3 0 . 018083584860929095 λ 4 0 . 008318170002599427 λ 5 0 . 005626648256273572 λ 6 0 . 005573530190419659 λ 7 = λ 8 0 . 0035787778473651247 λ 9 = λ 10 0 . 0021993489926981614 λ 11 0 . 0020129229654381126 λ 12 0 . 0016618597373763058
and, therefore, can be approximated by a matrix of order r = 12 with characteristic polynomial:
P ( λ ) = k = 1 r λ λ k .
Starting from (45), the invariants are readily evaluated as:
u 1 = I 1 1 ( T ) 0 . 146139760612268 u 2 = I 2 1 ( T ) 7 . 435124123661688 × 10 3 u 3 = I 3 1 ( T ) 1 . 923841312742190 × 10 4 u 4 = I 4 1 ( T ) 2 . 933369354737111 × 10 6 u 5 = I 5 1 ( T ) 2 . 845078672719340 × 10 8 u 6 = I 6 1 ( T ) 1 . 830841553974511 × 10 10 u 7 = I 7 1 ( T ) 7 . 975992951293202 × 10 13 u 8 = I 8 1 ( T ) 2 . 357390178302073 × 10 15 u 9 = I 9 1 ( T ) 4 . 646069994825235 × 10 18 u 10 = I 1 1 0 ( T ) 5 . 833548132502803 × 10 21 u 11 = I 1 1 1 ( T ) 4 . 213005422976779 × 10 24 u 12 = I 1 1 2 ( T ) 1 . 329866701687217 × 10 27
so that the companion matrix A can be written as:
Axioms 09 00073 i005
By making use of the Dunford–Taylor integral Formula (24), the matrix function:
f ( A ) = log ( A )
can be computed as:
Axioms 09 00073 i006

8.2.2. Example 4

In the case of a circular clamped plate with D x 2 + y 2 1 , the operator T which describes the relevant free vibrations features the eigenvalues:
λ 1 0 . 009581894052997456 λ 2 = λ 3 0 . 0022123673547497866 λ 4 0 . 0008220946148692253 λ 5 0 . 005626648256273572 λ 6 0 . 0006322131569880101 λ 7 = λ 8 0 . 00038401499194528553 λ 9 = λ 10 0 . 00027026150503226925
and, therefore, can be approximated by a matrix of order r = 10 with characteristic polynomial:
P ( λ ) = k = 1 r λ λ k .
Starting from (49), the invariants are readily evaluated as:
u 1 = I 1 1 ( T ) 2 . 239613778458294 × 10 2 u 2 = I 2 1 ( T ) 1 . 834047170883283 × 10 4 u 3 = I 3 1 ( T ) 7 . 231844608387320 × 10 7 u 4 = I 4 1 ( T ) 1 . 549026411656578 × 10 9 u 5 = I 5 1 ( T ) 1 . 916470543643900 × 10 12 u 6 = I 6 1 ( T ) 1 . 421413780191233 × 10 15 u 7 = I 7 1 ( T ) 6 . 388300810511680 × 10 19 u 8 = I 8 1 ( T ) 1 . 698427046738193 × 10 22 u 9 = I 9 1 ( T ) 2 . 451210499521412 × 10 26 u 10 = I 1 1 0 ( T ) 1 . 477289762215714 × 10 30
so that the companion matrix A can be written as:
Axioms 09 00073 i007
By making use of the Dunford–Taylor integral Formula (24), the matrix function:
f ( A ) = tan 1 ( A )
can be computed as:
Axioms 09 00073 i008

9. Conclusions

The use of Bell polynomials allows for extending to a countable set of numbers the Newton–Girard formulas originally stated in the framework of polynomials.
Then, the Robert’s formulas for the orthogonal invariants of a positive compact operator are recognized as the recursion or the Faà di Bruno [25] formula for the Bell polynomials. Therefore, the reduction of the order (or the grade) of the orthogonal invariants is easily established.
By using the results on matrix functions, based on the Dunford–Taylor integral, an approximation to holomorphic function of positive compact operators is obtained considering a truncation of Fredholm’s integral function whose coefficients are precisely the orthogonal invariants of the PCO. This is possible since the eigenvalues of such operators often tend to zero in a very fast way so that a finite set of eigenvalues gives a sufficient idea of the behavior of the considered function. Numerical examples obtained by using the computer algebra program Mathematica © are included. The results matches with the theoretical approach.

Author Contributions

Methodology, P.E.R.; software, D.C.; validation, D.C.; formal analysis, P.E.R. and P.N.; investigation, D.C.; data curation, D.C.; writing—original draft preparation, P.E.R.; writing—review and editing, D.C.; visualization, D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bell, E.T. Exponential polynomials. Ann. Mathematics 1934, 35, 258–277. [Google Scholar] [CrossRef]
  2. Riordan, J. An Introduction to Combinatorial Analysis; J. Wiley & Sons: Chichester, UK, 1958. [Google Scholar]
  3. Bruschi, M.; Ricci, P.E. I polinomi di Lucas e di Tchebycheff in più variabili. Rend. Mat. 1980, 13, 507–530. [Google Scholar]
  4. Cassisa, C.; Ricci, P.E. Orthogonal invariants and the Bell polynomials. Rend. Mat. Appl. 2000, 20, 293–303. [Google Scholar]
  5. Fichera, G. Abstract and Numerical Aspects of Eigenvalue Theory; Lecture Notes; Department of Mathematics, The University of Alberta: Edmonton, AB, Canada, 1973. [Google Scholar]
  6. Fichera, G. Metodi e risultati concernenti l’analisi numerica e quantitativa. Atti Accad. Nazion. Lincei Mem. 1974, 7, 1–202. [Google Scholar]
  7. Caratelli, D.; Ricci, P.E. A Numerical Method for Computing the Roots of Non-Singular Complex-Valued Matrices. Symmetry 2020, 12, 966. [Google Scholar] [CrossRef]
  8. Dodds, P.G. Positive compact operators. Quaest. Math. 1995, 18, 21–45. [Google Scholar] [CrossRef]
  9. Osborn, J.E. Spectral Approximation for Compact Operators. Math. Comput. 1975, 29, 712–715. [Google Scholar] [CrossRef]
  10. Zerner, M. Quelques propriétés spectrales des opérateurs positifs. J. Funct. Anal. 1987, 72, 381–417. [Google Scholar] [CrossRef] [Green Version]
  11. Roman, S.M. The Faà di Bruno Formula. Am. Math. Mon. 1980, 87, 805–809. [Google Scholar] [CrossRef]
  12. Roman, S.M.; Rota, G.C. The umbral calculus. Adv. Math. 1978, 27, 95–188. [Google Scholar] [CrossRef] [Green Version]
  13. Comtet, L. Advanced Combinatorics: The Art of Finite and Infinite Expansions; D. Reidel Publishing Co.: Dordrecht, The Netherlands, 1974; Available online: https://doi.org/10.1007/978-94-010-2196-8 (accessed on 12 February 2020).
  14. Kurosh, A. Cours d’Algèbre Supérieure; Éditions Mir: Moscou, Russia, 1971. [Google Scholar]
  15. Mikhlin, S.G. Integral Equations and Their Applications, 2nd ed.; Pergamon Press: Oxford, UK, 1964. [Google Scholar]
  16. Tricomi, F.G. Su di una classe di nuclei definiti positivi. Boll. Dell’Unione Mat. Ital. 1967, 22, 1–3. [Google Scholar]
  17. Robert, D. Invariants orthogonaux pour certaines classes d’operateurs. Ann. Mathém. Pures Appl. 1973, 52, 81–114. [Google Scholar]
  18. Kato, T. Perturbation Theory for Linear Operators; Springer: Berlin/Heidelberg, Germany, 1966. [Google Scholar]
  19. Bruschi, M.; Ricci, P.E. An explicit formula for f(A) and the generating function of the generalized Lucas polynomials. SIAM J. Math. Anal. 1982, 13, 162–165. [Google Scholar] [CrossRef]
  20. Caratelli, D.; Ricci, P.E. A note on integer powers of a Companion matrix and Applications. Tbilisi Math. J. 2020. to appear. [Google Scholar]
  21. Krall, G. Meccanica Tecnica Delle Vibrazioni; Veschi: Roma, Italy, 1970; Volume I. [Google Scholar]
  22. Fichera, G.; Sneider, M.A. Un problema di autovalori proposto da Alexander M. Ostrowski. Rend. Mat. 1975, 8, 201–224. [Google Scholar]
  23. Natalini, P.; Noschese, S.; Ricci, P.E. An iterative method for computing the eigenvalues of second kind Fredoholm operators and applications. Electron. Trans. Numer. Anal. 1999, 9, 128–136. [Google Scholar]
  24. Belingeri, C.; Germano, B. Numerical approximation of eigenvalues for transverse vibrations of a wedge-shaped beam. Appl. Math. Inf. Tbil. 1999, 4, 1–10. [Google Scholar]
  25. Faà di Bruno, F. Théorie des Formes Binaires; Brero: Turin, Italy, 1876. [Google Scholar]

Share and Cite

MDPI and ACS Style

Caratelli, D.; Natalini, P.; Ricci, P.E. Approximating Functions of Positive Compact Operators by Using Bell Polynomials. Axioms 2020, 9, 73. https://doi.org/10.3390/axioms9030073

AMA Style

Caratelli D, Natalini P, Ricci PE. Approximating Functions of Positive Compact Operators by Using Bell Polynomials. Axioms. 2020; 9(3):73. https://doi.org/10.3390/axioms9030073

Chicago/Turabian Style

Caratelli, Diego, Pierpaolo Natalini, and Paolo Emilio Ricci. 2020. "Approximating Functions of Positive Compact Operators by Using Bell Polynomials" Axioms 9, no. 3: 73. https://doi.org/10.3390/axioms9030073

APA Style

Caratelli, D., Natalini, P., & Ricci, P. E. (2020). Approximating Functions of Positive Compact Operators by Using Bell Polynomials. Axioms, 9(3), 73. https://doi.org/10.3390/axioms9030073

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