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 , where and , are differentiable functions (up to the considered order), defined in suitable intervals of the real axis, so that can be differentiated n times with respect to t, by using the chain rule.
Here, and in what follows, we use the notations:
Then, the
n-th derivative of
is represented by
where the
denote the Bell polynomials.
Further examples can be found in [
2], p. 49.
Proposition 1. The Bell polynomials satisfy the recurrence relation: An explicit expression for the Bell polynomials is given by the Faà di Bruno formula:
where the sum runs over all partitions
of the integer
n,
denotes the number of parts of size
i, and
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:
and the
coefficients satisfy the recursion:
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
, and denote by
the relevant elementary symmetric functions, and by
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: The above Formulas (
6)–(
7) constitute an extension of the Newton–Girard formulas and their inverse, since we have, in particular:
Remark 1. Note that the above equation implies the following well known recursion:which permits recovering the elementary symmetric functions. 4. Orthogonal Invariants of PCO and Robert’s Formulas
The eigenvalues
of a positive compact operator
T in a complex Hilbert space
satisfy
and, when infinite many eigenvalues exist, zero is an accumulation point.
If the PCO is
strictly positive, the above condition becomes:
A classical example of eigenvalue problem for a strictly positive operator is given by
where the kernel
of the second kind Fredholm operator
belongs to
, and is such that
if
(see Mikhlin [
15] and Tricomi [
16]).
The orthogonal invariants, introduced by Fichera [
5,
6], are, by definition, symmetric functions of the eigenvalues of
T:
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
, or
have been derived by Robert [
17].
The relative formulas are as follows:
which allow for reducing the orthogonal invariant
to
.
Since the eigenvalues of
are given by
, in the following, denoting by
the smallest integer such that
, we will put
, so that
, and the above Equations (
13) and (
14) become:
Under the above condition, the operator
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 [5,6]; however, this condition is satisfied by all the PCO occurring in applications. It is well known that, in the Hilbert space
, the following results hold:
where
Furthermore, the orthogonal invariants can be represented by the multiple integrals [
6]:
where
denotes the Fredholm determinant:
In particular, for
we have:
and, for
:
In what follows, we consider only the operators
T that belong to the Trace class. We put
, so that it results in:
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:
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:
Remark 3. Note that, recalling the representation in terms of seconk kind Bell polynomials, Equation (17) becomes: so that the can be computed by using the recursion (3). Proposition 5. For any integer , the orthogonal invariant , is expressed in terms of , by where denotes the sum running on all partitions of and .
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:
6. Matrix Functions
We consider a
matrix
, whose invariants
are known, so that its characteristic polynomial is given by
Let
f be a holomorphic in an open set
, containing all the eigenvaues of
. Then, the matrix functions
are given by the Dunford–Taylor (also called Riesz–Fantappiè) integral [
18,
19]:
where
denotes a simple contour enclosing all the zeros of
.
7. Approximation of a Strictly PCO
Equation (
24) is used in what follows, since we approximate a strictly positive compact operator
by a
matrix
, where the integer
r is such that the
rth eigenvalue
of
is negligible, i.e.,
where
is a prescribed small quantity.
We assume, for the characteristic polynomial, the following form:
Therefore, the companion matrix of
is given by
and Equation (
24) becomes:
Note that the powers of the companion matrix
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:
8.1.1. Example 1
Let us consider the Ostrowski’s kernel
in (
30). As demonstrated by means of the inverse iteration method in [
22,
23], the relevant operator
displays the eigenvalues:
and, therefore, can be approximated by a matrix of order
with a characteristic polynomial:
Starting from (
34), the invariants are readily evaluated as:
so that the companion matrix
can be written as:
By making use of the Dunford–Taylor integral Formula (
24), the matrix function:
with
denoting the product logarithm, can be computed as:
8.1.2. Example 2
Let us consider the kernel
relevant to the positive compact operator
which describes the transverse vibrations of a wedge-shaped beam [
24]. By using the inverse iteration method, one can prove that
features the following eigenvalues:
and, therefore, can be approximated by a matrix of order
with characteristic polynomial:
Starting from (
38), the invariants are readily evaluated as:
so that the companion matrix
can be written as:
By making use of the Dunford–Taylor integral Formula (
24), the matrix function:
can be computed as:
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
, where
is given by:
Putting
, and considering the partial sum, of order
r, of the series (
30), we recover the polynomial
in Equation (
26) as an approximation of the Fredholm integer function (
41).
In [
5], Fichera, considering the biharmonic operator
was able to find lower and upper approximations of the eigenvalues of the biharmonic problem
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
, the operator
which describes the relevant free vibrations features the eigenvalues:
and, therefore, can be approximated by a matrix of order
with characteristic polynomial:
Starting from (
45), the invariants are readily evaluated as:
so that the companion matrix
can be written as:
By making use of the Dunford–Taylor integral Formula (
24), the matrix function:
can be computed as:
8.2.2. Example 4
In the case of a circular clamped plate with
, the operator
which describes the relevant free vibrations features the eigenvalues:
and, therefore, can be approximated by a matrix of order
with characteristic polynomial:
Starting from (
49), the invariants are readily evaluated as:
so that the companion matrix
can be written as:
By making use of the Dunford–Taylor integral Formula (
24), the matrix function:
can be computed as:
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
- Bell, E.T. Exponential polynomials. Ann. Mathematics 1934, 35, 258–277. [Google Scholar] [CrossRef]
- Riordan, J. An Introduction to Combinatorial Analysis; J. Wiley & Sons: Chichester, UK, 1958. [Google Scholar]
- Bruschi, M.; Ricci, P.E. I polinomi di Lucas e di Tchebycheff in più variabili. Rend. Mat. 1980, 13, 507–530. [Google Scholar]
- Cassisa, C.; Ricci, P.E. Orthogonal invariants and the Bell polynomials. Rend. Mat. Appl. 2000, 20, 293–303. [Google Scholar]
- Fichera, G. Abstract and Numerical Aspects of Eigenvalue Theory; Lecture Notes; Department of Mathematics, The University of Alberta: Edmonton, AB, Canada, 1973. [Google Scholar]
- Fichera, G. Metodi e risultati concernenti l’analisi numerica e quantitativa. Atti Accad. Nazion. Lincei Mem. 1974, 7, 1–202. [Google Scholar]
- 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]
- Dodds, P.G. Positive compact operators. Quaest. Math. 1995, 18, 21–45. [Google Scholar] [CrossRef]
- Osborn, J.E. Spectral Approximation for Compact Operators. Math. Comput. 1975, 29, 712–715. [Google Scholar] [CrossRef]
- Zerner, M. Quelques propriétés spectrales des opérateurs positifs. J. Funct. Anal. 1987, 72, 381–417. [Google Scholar] [CrossRef] [Green Version]
- Roman, S.M. The Faà di Bruno Formula. Am. Math. Mon. 1980, 87, 805–809. [Google Scholar] [CrossRef]
- Roman, S.M.; Rota, G.C. The umbral calculus. Adv. Math. 1978, 27, 95–188. [Google Scholar] [CrossRef] [Green Version]
- 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).
- Kurosh, A. Cours d’Algèbre Supérieure; Éditions Mir: Moscou, Russia, 1971. [Google Scholar]
- Mikhlin, S.G. Integral Equations and Their Applications, 2nd ed.; Pergamon Press: Oxford, UK, 1964. [Google Scholar]
- Tricomi, F.G. Su di una classe di nuclei definiti positivi. Boll. Dell’Unione Mat. Ital. 1967, 22, 1–3. [Google Scholar]
- Robert, D. Invariants orthogonaux pour certaines classes d’operateurs. Ann. Mathém. Pures Appl. 1973, 52, 81–114. [Google Scholar]
- Kato, T. Perturbation Theory for Linear Operators; Springer: Berlin/Heidelberg, Germany, 1966. [Google Scholar]
- 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]
- Caratelli, D.; Ricci, P.E. A note on integer powers of a Companion matrix and Applications. Tbilisi Math. J. 2020. to appear. [Google Scholar]
- Krall, G. Meccanica Tecnica Delle Vibrazioni; Veschi: Roma, Italy, 1970; Volume I. [Google Scholar]
- Fichera, G.; Sneider, M.A. Un problema di autovalori proposto da Alexander M. Ostrowski. Rend. Mat. 1975, 8, 201–224. [Google Scholar]
- 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]
- 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]
- Faà di Bruno, F. Théorie des Formes Binaires; Brero: Turin, Italy, 1876. [Google Scholar]
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).