Next Article in Journal
Closeness Centrality of Asymmetric Trees and Triangular Numbers
Next Article in Special Issue
Continuity Equation of Transverse Kähler Metrics on Sasakian Manifolds
Previous Article in Journal
Few-Shot Classification Based on Sparse Dictionary Meta-Learning
Previous Article in Special Issue
Growth Spaces on Circular Domains Taking Values in a Banach Lattice, Embeddings and Composition Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Perturbation Approach to Polynomial Root Estimation and Expected Maximum Modulus of Zeros with Uniform Perturbations

by
Ibrahim A. Nafisah
1,
Sajad A. Sheikh
2,
Mohammed A. Alshahrani
3,
Mohammed M. A. Almazah
4,
Badr Alnssyan
5 and
Javid Gani Dar
6,*
1
Department of Statistics and Operations Research, College of Sciences, King Saud University, Riyadh 11451, Saudi Arabia
2
Department of Mathematics, University of Kashmir, South Campus, Anantnag 192101, India
3
Department of Mathematics, College of Sciences and Humanities, Prince Sattam Bin Abdulaziz University, Alkharj 11942, Saudi Arabia
4
Department of Mathematics, College of Sciences and Arts (Muhyil), King Khalid University, Muhyil 61421, Saudi Arabia
5
Department of Management Information Systems, College of Business and Economics, Qassim University, Buraydah 51452, Saudi Arabia
6
Department of Applied Sciences, Symbiosis Institute of Technology, Symbiosis International (Deemed University), Pune 412115, India
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(19), 2993; https://doi.org/10.3390/math12192993
Submission received: 11 September 2024 / Revised: 24 September 2024 / Accepted: 24 September 2024 / Published: 26 September 2024
(This article belongs to the Special Issue Complex Analysis and Geometric Function Theory, 2nd Edition)

Abstract

:
This paper presents a significant extension of perturbation theory techniques for estimating the roots of polynomials. Building upon foundational results and recent work by Pakdemirli and Yurtsever, as well as taking inspiration from the concept of probabilistic bounds introduced by Sheikh et al., we develop and prove several novel theorems that address a wide range of polynomial structures. These include polynomials with multiple large coefficients, coefficients of different orders, alternating coefficient orders, large linear and constant terms, and exponentially decreasing coefficients. Among the key contributions is a theorem that establishes an upper bound on the expected maximum modulus of the zeros of polynomials with uniformly distributed perturbations in their coefficients. The theorem considers the case where all but the leading coefficient receive a uniformly and independently distributed perturbation in the interval [ 1 , 1 ] . Our approach provides a comprehensive framework for estimating the order of magnitude of polynomial roots based on the structure and magnitude of their coefficients without the need for explicit root-finding algorithms. The results offer valuable insights into the relationship between coefficient patterns and root behavior, extending the applicability of perturbation-based root estimation to a broader class of polynomials. This work has potential applications in various fields, including random polynomials, control systems design, signal processing, and numerical analysis, where quick and reliable estimation of polynomial roots is crucial. Our findings contribute to the theoretical understanding of polynomial properties and provide practical tools for engineers and scientists dealing with polynomial equations in diverse contexts.

1. Introduction

The study of polynomial roots has been a fundamental area of mathematics for centuries, with applications ranging from pure algebra to various fields of science and engineering [1,2,3,4]. In recent years, the use of perturbation theory in estimating polynomial roots has gained significant attention, providing valuable insights into the behavior of roots under small parameter changes [5]. Perturbation theory is also crucial in studying random polynomials, where coefficients are drawn from probability distributions. The analysis of how small random perturbations affect the roots of polynomials provides insights into the stability and sensitivity of polynomial systems under uncertainty [6]. This approach has found applications in various fields, including the study of eigenvalues of random matrices and the behavior of complex systems [7].
Galántai and Hegedüs [8] present perturbation bounds on polynomial zeros that improve the results of Ostrowski, Elsner et al., and Beauzamy while also enhancing the backward stability result of Edelman and Murakami. Pakdemirli and Yurtsever [9] offer theorems to estimate the order of magnitude of roots in polynomials, aiding in a priori estimation before solving the equation. Choque-Rivero and Garza [10] study moment perturbation of matrix orthogonal polynomials, establishing an explicit relation between perturbed and original polynomials. In a significant study, Maiz [11] applies perturbation theory to derive polynomial solutions for real quantum potential systems, providing a useful tool for verifying and improving numerical and approximate methods in quantum theory. Recently, Dmytryshyn [12] developed an algorithm to recover a perturbation of a matrix polynomial from a perturbation of its first companion linearization. The accuracy and stability of numerical algorithms used in perturbation analysis are critical considerations, as highlighted by Higham [13]. These aspects become particularly important when dealing with high-degree polynomials or those with closely spaced roots, where small perturbations can lead to significant changes in root locations. The advancement of computational power has enabled researchers to employ sophisticated statistical models and simulations to study the behavior of polynomial roots under various perturbation scenarios. These computational approaches complement theoretical results and provide valuable insights into the complex dynamics of perturbed polynomial systems.
Ostrowski [14] provided a fundamental result on the continuity of polynomial roots with respect to their coefficients and provided bounds on the distance between the origins of two polynomials based on the differences in their coefficients. His theorem gives bounds on the distance between the roots of two polynomials based on the differences in their coefficients.
Theorem 1 
(Ostrowski [14]). Let p ( z ) = z n + a 1 z n 1 + + a n and q ( z ) = z n + b 1 z n 1
+ + b n . For any root x i of p ( z ) , there exists a root y j of q ( z ) such that
| x i y j | k = 1 n | a k b k | γ n k 1 / n ,
where γ = 2 max 1 k n { | a k | 1 / k , | b k | 1 / k } . Furthermore, the roots of p and q can be enumerated as α 1 , , α n and β 1 , , β n , respectively, in such a way that max 1 i n | α i β i | ( 2 n 1 )   k = 1 n | a k b k | γ n k 1 / n .
Research in this direction has been further developed and refined by subsequent researchers [15,16,17,18], contributing to a more comprehensive understanding of polynomial root perturbation. Bhatia et al. [17] extended the study of root perturbations to matrix eigenvalues. They established an upper bound for the variation of eigenvalues between two matrices.
Definition 1.
Let A , B C n × n . Assume that σ ( A ) = { λ 1 , , λ n } and σ ( B ) = { μ 1 , , μ n } . Let S n be the set of all permutations of { 1 , 2 , , n } . The eigenvalue variation of A and B is defined by v ( A , B ) = min π S n max i | μ π ( i ) λ i | . The quantity v ( A , B ) is also called the (optimal) matching distance between the eigenvalues of A and B .
Bhatia et al. improved Ostrowski’s matrix perturbation theorem in several aspects.
Theorem 2
(Bhatia et al. [17]). Let A , B C n × n . Then v ( A , B ) 4 × 2 1 / n ( A + B ) 1 1 / n A B 1 / n .
Beauzamy further refined the analysis of polynomial root perturbations by considering the multiplicity of roots. For their result, we first need the following definition.
Definition 2.
For the polynomial p ( z ) = j = 0 n a j z n j , the Bombieri norm is defined by
[ p ] B = j = 0 n | a j | 2 n j 1 / 2
The theorem provides a more precise bound that takes into account the derivative of the polynomial at the root.
Theorem 3
(Beauzamy [18]). Let k 1 be an integer and p ( z ) and q ( z ) be two polynomials of degree n, with [ p q ] B ε . If x i is any zero of p ( z ) with multiplicity k, there exists a zero y j of q ( z ) , with | x i y j | n ! ( n k ) ! 1 / k ε 1 / k | q ( k ) ( x i ) | 1 / k ( 1 + | x i | 2 ) n / 2 . Further, if ε ( n k ) ! 2 n ! | p ( k ) ( x i ) | 2 , then
| x i y j | 2 n ! ( n k ) ! 1 / k ε 1 / k | p ( k ) ( x i ) | 1 / k ( 1 + | x i | 2 ) ( n k ) / 2
The work of Pakdemirli and Yurtsever [9] introduced a novel approach to estimating roots of polynomials using perturbation theory. In [9], Pakdemirli and Yurtsever, among other things, established the following theorems.
Theorem 4.
For a polynomial equation of the form
a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0
If all coefficients a i ( i = 0 , , n ) are of the same order of magnitude, then the roots are of O ( 1 ) .
Theorem 5.
For a polynomial equation of the form
a n x n + a n 1 x n 1 + + a m x m + + a 1 x + a 0 = 0
If a m O ( 1 / ϵ k ) ( k > 0 ) with all other coefficients being of O ( 1 ) , then the possible roots may be of either of the following:
  • O ( ϵ k / m ) (for m 0 ).
  • O ( 1 / ϵ k / ( n m ) ) (for m n ).
Their method, based on the order of magnitude of coefficients, provided a framework for predicting the approximate magnitude of roots before solving the equation explicitly. This approach has proven particularly useful in cases where traditional numerical methods may struggle, such as when dealing with ill-conditioned polynomials or when a quick estimate is needed before applying more computationally intensive techniques [19]. These examples illustrate the continued interest and wide interface of perturbation theory, especially in the context of polynomials. The applications of this theory extend to various fields, including control systems [1], signal processing [20], and electrical engineering [3]. Numerical analysis techniques, crucial for implementing perturbation methods, are extensively covered in works such as [2,19,21].
The interdisciplinary nature of polynomial perturbation theory is further evidenced by its applications in nonlinear dynamics [22], computer-aided geometric design [23], and quantum chemistry [24].
In this paper, we aim to further expand the scope of perturbation-based root estimation techniques. We present several new theorems that address previously unexplored polynomial structures, including the following:
  • Polynomials with multiple large coefficients.
  • Polynomials with coefficients of different orders.
  • Polynomials with alternating coefficient orders.
  • Polynomials with large linear and constant terms.
  • Polynomials with exponentially decreasing coefficients.
  • Bounds of zeros of a randomly perturbed polynomial.
These theorems provide a comprehensive framework for estimating root magnitudes across a wide range of polynomial types, significantly extending the applicability of perturbation theory in root-finding problems. Moreover, we provide rigorous proofs for each theorem, along with illustrative examples that demonstrate their practical application.
Our work not only contributes to the theoretical understanding of polynomial root behavior but also has potential applications in various fields where quick root estimation is crucial. These include control systems design [1], signal processing [20], and numerical analysis [13], among others.

1.1. Key Concepts

In perturbation theory, the order of a coefficient refers to its magnitude relative to a small parameter ϵ 1 . Several examples follow:
  • O ( 1 ) means the coefficient is of order 1.
  • O 1 ϵ means the coefficient is much larger than 1, approximately 1 ϵ .
  • O ( ϵ ) means the coefficient is much smaller than 1, approximately ϵ .

1.2. Limitations

Despite the advancements presented in this paper, several limitations should be acknowledged. The perturbation-based approach provides estimates of root magnitudes rather than exact root values, which may not be sufficient for all applications. The accuracy of the estimates may decrease for polynomials with very closely spaced roots or in cases where multiple perturbation effects interact in complex ways. It is important to note that the theorems presented are primarily applicable to scalar polynomials and may not directly extend to matrix polynomials or multivariate cases without further development; while the method offers quick estimates, it does not replace the need for rigorous numerical methods in applications requiring high precision. Additionally, the approach may have limited effectiveness for polynomials with highly oscillatory coefficient patterns or those arising from certain transcendental equations. These limitations present opportunities for future research and refinement of the perturbation-based root estimation techniques presented in this paper. Addressing these challenges could further enhance the applicability and reliability of these methods across a broader range of mathematical and engineering problems.
The rest of this paper is organized as follows: Section 2 presents the main theorems and their proofs. Section 3 provides numerical simulations and discussions. Section 4 explores the potential applications and implications of our results. Finally, Section 5 concludes this paper and suggests directions for future research.

1.3. Motivation

The motivation for this research stems from several key factors. Primarily, the ubiquity of polynomial equations in various scientific and engineering fields, including control systems [1], signal processing [20], and electrical engineering [3], underscores the importance of efficient root estimation techniques. There is a pressing need for quick and reliable estimation of polynomial roots, especially in cases where traditional numerical methods may struggle, such as with ill-conditioned polynomials [19].

2. Main Theorems

Theorem 6
(multiple large coefficients). Consider a polynomial equation of the form
a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0
Assume there are multiple coefficients of order O 1 ϵ k for some k > 0 , where these coefficients correspond to indices p and q satisfying 1 p q n 1 . The remaining coefficients, including the constant term a 0 and the leading term a n , are of order O ( 1 ) . Then, the roots are of order O ϵ k p or O 1 ϵ k n q , where p is the smallest index and q is the largest index among the large coefficients.
Proof. 
Let p be the smallest index and q be the largest index such that a p , a q O 1 ϵ k . We will consider two cases for the possible order of magnitude of the roots.
Case 1: Assume | x | O ( ϵ r ) , where r > 0 .
In this case, the term a p x p will be the dominant term among those with large coefficients, as it has the smallest power of x. We need to balance this term with the constant term a 0 :
a p x p a 0
Substituting the orders of magnitude:
1 ϵ k · ( ϵ r ) p O ( 1 )
ϵ p r ϵ k O ( 1 )
ϵ p r k O ( 1 )
For relation (2) to hold, the exponent must satisfy
p r k = 0
r = k p
Therefore, one possible order for the roots is O ϵ k p .
Case 2: Assume | x | O 1 ϵ r , where r > 0 .
In this case, the term a n x n will be the dominant term among those with O ( 1 ) coefficients. We need to balance this term with a q x q , which is the dominant term among those with large coefficients:
a n x n a q x q
Substituting the orders of magnitude:
O ( 1 ) · 1 ϵ r n 1 ϵ k · 1 ϵ r q
1 ϵ n r 1 ϵ q r + k
For Equation (3) to hold, the exponents must satisfy
n r = q r + k
r ( n q ) = k
r = k n q
Therefore, another possible order for the roots is O 1 ϵ k n q .
Thus, we have shown that the roots are of order O ϵ k p or O 1 ϵ k n q , where p is the smallest index and q is the largest index among the large coefficients. □
Theorem 7
(coefficients of different orders). Consider a polynomial equation where the coefficients have different orders of magnitude:
a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0
If a i O 1 ϵ k i for i = 0 , 1 , , n , where some k i may be zero or positive, then the possible orders of the roots are given by
O ϵ k i k j i j
for any pair of distinct indices i j .
Proof. 
We will prove this theorem by considering the balance between any two terms in the polynomial and showing that this leads to the stated order of magnitude for the roots.
Let i and j be any two distinct indices from the set { 0 , 1 , , n } . Consider balancing the i-th and j-th terms of the polynomial:
a i x i a j x j
Substituting the orders of magnitude for the coefficients:
1 ϵ k i · x i 1 ϵ k j · x j
Assume that x ϵ r for some real number r. Substituting this into the above equation:
1 ϵ k i · ( ϵ r ) i 1 ϵ k j · ( ϵ r ) j
ϵ r i ϵ k i ϵ r j ϵ k j
ϵ r i k i ϵ r j k j
For the relation (4) balance to hold, the exponents must be equal:
r i k i = r j k j
Solving for r:
r i r j = k i k j
r ( i j ) = k i k j
r = k i k j i j
Therefore, one possible order for the roots is O ϵ k i k j i j .
Since this argument holds for any pair of distinct indices i and j, we have multiple possible orders for the roots, each corresponding to a different pair of terms being balanced.
Thus, we have proven that the possible orders of the roots are given by O ϵ k i k j j i for any pair of distinct indices i j . □
Theorem 8
(alternating coefficient orders). Consider a polynomial equation of the form
a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0
where the coefficients alternate between two orders of magnitude:
a i O ( 1 / ϵ k ) i f i i s e v e n O ( 1 ) i f i i s o d d
with k > 0 and ϵ 1 . Then the roots of this polynomial are of order O ( 1 / ϵ k ) .
Proof. 
Without loss of generality, let n be even. We consider two cases:
Case 1: Assume | x | O ( ϵ r ) , where r > 0 .
Balancing the constant term with the first-order term:
a 0 a 1 x
Substituting the orders of magnitude:
1 ϵ k O ( 1 ) · ϵ r
The expression (6):
ϵ k ϵ r r = k
Case 2: Assume | x | O ( 1 / ϵ r ) , where r > 0 .
Balancing the n-th term with the ( n 1 ) -th term:
a n x n a n 1 x n 1
Substituting the orders of magnitude:
1 ϵ k · 1 ϵ r n O ( 1 ) · 1 ϵ r n 1
Expression (7) simplifies to
1 ϵ k + n r 1 ϵ ( n 1 ) r
For Equation (8) to hold, the exponents must be equal:
k + n r = ( n 1 ) r k = r
giving us
r = k
In both cases, we find r = k . Therefore, the roots are of order O ( 1 / ϵ k ) . □
Theorem 9
(polynomial with large linear and constant terms). Consider a polynomial equation of the form
x n + a n 1 x n 1 + + a 2 x 2 + a 1 x + a 0 = 0
where a 1 O 1 ϵ k and a 0 O 1 ϵ m with m > k > 0 and ϵ 1 . All other coefficients are O ( 1 ) . Then, the roots of this polynomial are of order O ( ϵ m k ) , O ϵ k n 1 , or O ( 1 ) .
Proof. 
We consider three cases based on the assumed order of the roots.
Case 1: Assume | x | O ( ϵ r ) , where r > 0 .
Balancing the constant term with the linear term:
a 0 a 1 x
Substituting the orders of magnitude:
1 ϵ m 1 ϵ k · ϵ r
Expression (9) gives
ϵ m ϵ r k m = r k
r = m k
For r to be positive, we require m > k . Therefore, the roots are of order O ϵ m k .
Case 2: Assume | x | O 1 ϵ r , where r > 0 .
Balancing the n-th term with the linear term:
x n a 1 x
Substituting the orders of magnitude:
1 ϵ r n 1 ϵ k · 1 ϵ r
The expression (10) simplifies to
1 ϵ n r 1 ϵ k + r
For Equation (11) to hold, the exponents must satisfy
n r = k + r
r ( n 1 ) = k
r = k n 1
Thus, the roots are of order O ϵ k n 1 .
Case 3: Assume | x | O ( 1 ) .
In this case, all terms are O ( 1 ) , and thus, the roots can be of order O ( 1 ) .
Therefore, the roots are of order O ϵ m k (from Case 1), O ϵ k n 1 (from Case 2), or O ( 1 ) (from Case 3). □
Theorem 10
(exponentially decreasing coefficients). Consider a polynomial equation of the form
a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0
where the coefficients decrease exponentially:
a i O 1 ϵ k · 2 i
with k > 0 , ϵ 1 , and i = 0 , 1 , , n . Then, the roots of this polynomial are of order O ϵ k · 2 i for each i = 0 , 1 , , n 1 .
Proof. 
To determine the order of the roots, we consider the balance between consecutive terms in the polynomial.
Consider balancing the i-th and ( i + 1 ) -th terms of the polynomial:
a i x i a i + 1 x i + 1
Substituting the orders of magnitude:
1 ϵ k · 2 i · x i 1 ϵ k · 2 i + 1 · x i + 1
Assume x ϵ r for some real number r. Substituting this into the above equation:
1 ϵ k · 2 i · ( ϵ r ) i 1 ϵ k · 2 i + 1 · ( ϵ r ) i + 1
Simplifying:
ϵ r i k · 2 i ϵ r ( i + 1 ) k · 2 i + 1
Equate the exponents:
r i k · 2 i = r ( i + 1 ) k · 2 i + 1
Simplifying:
r i k · 2 i = r i + r k · 2 i + 1
Subtract r i from both sides:
k · 2 i = r k · 2 i + 1
Rearranging terms:
r = k · 2 i + 1 k · 2 i
r = k · 2 i ( 2 1 )
r = k · 2 i
Therefore, the order of the roots corresponding to the balance between the i-th and ( i + 1 ) -th terms is O ϵ k · 2 i .
Since this argument holds for each pair of consecutive indices i and i + 1 where i = 0 , 1 , , n 1 , we conclude that the roots of the polynomial are of order O ϵ k · 2 i for each i = 0 , 1 , , n 1 . That proves the theorem. □

Applications in Finding Bounds of Zeros

As an extension of our results to the study of bounds of zeros of random polynomials, taking inspiration from the work of Sheikh and Mir [25] about the zero bounds of random polynomials, we present the following theorem that concerns the roots of a randomly perturbed monic polynomial with uniformly perturbed coefficients. In the study of polynomials with perturbed coefficients, understanding the impact of these perturbations on the roots is crucial for both theoretical and practical applications. The following theorem addresses this by considering two classes of perturbed polynomials by employing the Cauchy bound for roots. Specifically, we consider the case when all but the leading coefficient receives a uniformly distributed perturbation in the interval [ 1 , 1 ] , each modeled as independent uniform random variables. By deriving an upper bound on the expected maximum modulus of the polynomial’s zeros, the theorem provides insights into how these specific perturbations influence the distribution of the roots.
Theorem 11
(expected maximum modulus of zeros with uniform perturbations). Let q ( z ) = z n + ( a 1 + δ 1 ) z n 1 + + ( a n + δ n ) be a polynomial where each δ k is a random variable representing a small perturbation, modeled as δ k U n i f o r m ( 1 , 1 ) where U n i f o r m ( a , b ) means being uniformly distributed over the interval [ a , b ] . The expected maximum modulus of the zeros of q ( z ) is bounded above by
E max i | z i | R + n 1 n + 1
where z i are the zeros of q ( z ) , where R = 1 + max 0 i n 1 | a i | is the Cauchy bound of the zeros of polynomial
p ( z ) = z n + a n 1 z n 1 + + a 0
Before proceeding to the proof, we need a lemma which we state and prove as follows:
Lemma 1 (expected maximum of uniform variables).
For δ k U n i f o r m ( 1 , 1 ) , the expected maximum of n such variables is
E [ M ] = n 1 n + 1
Proof of Lemma 1.
Let M = max ( δ 1 , δ 2 , , δ n ) , where each δ k is independently and identically distributed as Uniform ( 1 , 1 ) .
The cumulative distribution function (CDF) of a single uniform variable δ k is given by
F ( x ) = x + 1 2 for 1 x 1
The CDF of the maximum M of n such independent variables is
F M ( x ) = x + 1 2 n
The probability density function (PDF) of M is the derivative of the CDF:
f M ( x ) = d d x x + 1 2 n = n x + 1 2 n 1 · 1 2
Thus, the PDF is
f M ( x ) = n 2 x + 1 2 n 1
The expected value of M is
E [ M ] = 1 1 x · f M ( x ) d x
Substituting the PDF:
E [ M ] = 1 1 x · n 2 x + 1 2 n 1 d x
Change variables: let u = x + 1 2 , then d u = 1 2 d x and d x = 2 d u . The limits of integration change from x = 1 to x = 1 to u = 0 to u = 1 .
E [ M ] = 0 1 ( 2 u 1 ) · n · u n 1 d u = n 0 1 ( 2 u n u n 1 ) d u = n 2 u n + 1 n + 1 u n n 0 1 = n 2 n + 1 1 n = n 1 n + 1
Thus, the expected maximum of n uniformly distributed variables is
E [ M ] = n 1 n + 1
Proof of Theorem 11.
Let M = max 1 k n | δ k | . We claim that the maximum modulus of the zeros of q ( z ) is bounded by R + M , where R is the bound on the maximum modulus of the zeros of the unperturbed polynomial z n + a 1 z n 1 + + a n .
Consider the perturbed polynomial:
q ( z ) = z n + ( a n 1 + δ n 1 ) z n 1 + + ( a 0 + δ 0 )
where each δ k is a random variable representing a small perturbation, modeled as δ k Uniform ( 1 , 1 ) .
For the unperturbed polynomial:
p ( z ) = z n + a n 1 z n 1 + + a 0
Cauchy’s bound provides an upper limit on the magnitude of the roots:
| z i | 1 + max ( | a n 1 | , | a n 2 | , , | a 0 | )
Define this bound as R:
R = 1 + max ( | a n 1 | , | a n 2 | , , | a 0 | )
For the perturbed polynomial q ( z ) , the coefficients are a k + δ k . The bound on the magnitude of the roots becomes
| z i | 1 + max ( | a n 1 + δ n 1 | , | a n 2 + δ n 2 | , , | a 0 + δ 0 | )
The maximum can be decomposed as
max ( | a n 1 + δ n 1 | , | a n 2 + δ n 2 | , , | a 0 + δ 0 | ) max ( | a n 1 | , | a n 2 | , , | a 0 | ) + max ( | δ n 1 | , | δ n 2 | , , | δ 0 | )
Define the maximum perturbation as M:
M = max ( | δ n 1 | , | δ n 2 | , , | δ 0 | )
Substituting the decomposed maximum into the bound for the perturbed polynomial, we have
| z i | 1 + max ( | a n 1 | , | a n 2 | , , | a 0 | ) + M
Thus, the bound for the perturbed polynomial is
| z i | R + M
This derivation shows that the bound R + M accounts for both the original bound R and the maximum effect of the perturbations M, ensuring that the magnitude of the roots is appropriately bounded.
Thus, the expected maximum modulus of the zeros is
E max i | z i | E [ R + M ]
Since R is a constant, we can separate the expectation:
E [ R + M ] = R + E [ M ]
Substituting the bound for E [ M ] using Lemma 1, we obtain:
E max i | z i | R + n 1 n + 1
E max i | z i | E [ R + M ]
That establishes the theorem. □

3. Discussion of Simulation Results of Uniform Perturbations

Our numerical simulations and visualizations provide valuable insights into the behavior of polynomial zeros under different perturbation scenarios, as described in Theorem 11. In our simulations depicted in Figure 1 for Theorem 11, we considered polynomials of degree n = 10 with coefficients perturbed with the Cauchy bound R = 2 for the underlying polynomial p ( z ) . We conducted a total of 10,000 trials to ensure statistical significance. The results demonstrate the effectiveness of the derived upper bound for the expected maximum modulus of zeros under uniform perturbations. The average maximum modulus of zeros obtained from our simulations was approximately 1.7011, which is consistently below the theoretical upper bound of 2.8182, validating the theorem’s prediction. The simulation results for Theorem 11 demonstrate the effectiveness of the derived upper bound for the expected maximum modulus of zeros under uniform perturbations. We observed that the average maximum modulus of zeros obtained from our simulations is consistently below the theoretical upper bound, validating the theorem’s prediction. Additionally, the distribution of maximum moduli shows a right-skewed shape, indicating that while most perturbed polynomials have zeros with relatively small moduli, there are occasional instances of larger maximum moduli. Furthermore, the theoretical bound provides a reliable estimate of the worst-case scenario, as it consistently lies to the right of the bulk of the distribution.

4. Conclusions

In this paper, we have significantly extended the application of perturbation theory to the estimation of polynomial roots. Building upon the foundational work of Pakdemirli and Yurtsever [9], we developed and proved several novel theorems that address a wide range of polynomial structures while developing a theorem illustrating how the work can be extended to randomly perturbed polynomials. The results we study include polynomials with multiple large coefficients, coefficients of different orders, alternating coefficient orders, large linear and constant terms, and exponentially decreasing coefficients. These theorems provide a comprehensive framework for estimating the order of magnitude of polynomial roots based on the structure and magnitude of their coefficients. Our results offer several key contributions to the field. First, they extend the applicability of perturbation-based root estimation to a broader class of polynomials, including those with complex coefficient patterns. Second, they provide quick and reliable estimates of root magnitudes without the need for explicit root-finding algorithms, which can be particularly valuable in time-sensitive applications or as a precursor to more detailed numerical analysis. Third, they offer insights into the relationship between coefficient patterns and root behavior, deepening our theoretical understanding of polynomial properties.
The numerical examples presented for the randomly perturbed polynomials with uniformly and identically distributed perturbation errors in this paper demonstrate the utility of our theorems across various scenarios that might arise in practical situations. These results suggest that our approach could be particularly valuable in fields such as control systems, signal processing, and numerical analysis, where quick estimation of polynomial roots is often necessary.
However, it is important to note the limitations of our approach. While these theorems provide excellent estimates for the order of magnitude of roots, they do not give exact root values. Additionally, in some cases with very close root magnitudes, the method may not distinguish between individual roots.
Future research directions could include extending these techniques to multivariate polynomials, investigating the application of these methods to transcendental equations, developing algorithms that combine these estimation techniques with traditional root-finding methods for improved efficiency, and exploring the implications of these results in specific application areas, such as filter design or stability analysis of dynamical systems.
In conclusion, this work represents a significant advancement in the field of perturbation-based polynomial root estimation. By providing a more comprehensive set of tools for analyzing polynomial roots, we hope to facilitate advancements in both theoretical mathematics and practical applications across various scientific and engineering disciplines.

5. Potential Applications in Various Scientific Fields

The novel theorems and techniques presented in this paper for estimating polynomial roots using perturbation theory have potential applications across a wide range of scientific and engineering disciplines. These applications stem from the ubiquity of polynomial equations in modeling complex systems and phenomena. In the field of random polynomials, where coefficients are drawn from probability distributions, perturbation analysis helps in understanding the stability and sensitivity of roots to small random variations [6]. This approach is particularly useful in fields such as statistical physics, where it aids in studying phase transitions and critical phenomena. Furthermore, perturbation techniques applied to random polynomials contribute to the analysis of eigenvalue distributions in random matrix theory, which has far-reaching implications in quantum chaos, wireless communications, and financial mathematics In the field of control systems engineering, our results can significantly impact the analysis and design of feedback control systems. The characteristic equation of a linear time-invariant system is often represented as a polynomial, and its roots determine the system’s stability and transient response [1]. Our quick estimation methods can provide engineers with rapid insights into system behavior during the design phase, potentially streamlining the iterative process of controller tuning.
Signal processing is another area that could benefit from these advancements. Many digital filters are designed based on polynomial equations, where the roots of these polynomials correspond to the filter’s poles and zeros [20]. Our theorems could offer a more intuitive understanding of how coefficient changes affect filter characteristics, possibly leading to more efficient filter design algorithms.
In the realm of numerical analysis, our work has implications for improving root-finding algorithms. By providing good initial estimates of root magnitudes, our methods could enhance the convergence speed and reliability of iterative root-finding techniques such as the Newton–Raphson method or Muller’s method [21]. This could be particularly valuable when dealing with high-degree polynomials or ill-conditioned systems.
The field of computational physics often involves solving polynomial equations that arise from discretized differential equations. Our perturbation-based approach could offer physicists quick insights into the behavior of these discretized systems, potentially guiding the choice of numerical methods or helping to identify potential instabilities in simulations [2].
In electrical engineering, particularly in circuit analysis and design, polynomial equations frequently appear in transfer functions and network analysis. Our methods could provide engineers with rapid estimates of circuit behavior, assisting in the early stages of design or in troubleshooting complex systems [3].
The study of dynamical systems, which spans multiple disciplines including mathematics, physics, and biology, often involves analyzing polynomial equations. Our theorems could offer new tools for understanding bifurcations and stability changes in these systems, potentially leading to new insights in areas such as population dynamics, chemical kinetics, and climate modeling [22].
In the field of computer graphics and geometric modeling, polynomial equations are used extensively in curve and surface representation. Our quick estimation techniques could potentially optimize algorithms for intersection detection or shape manipulation, leading to more efficient rendering and modeling software [23].
Quantum chemistry calculations often involve solving polynomial equations that arise from approximations to the Schrödinger equation. Our methods could potentially assist in the rapid estimation of molecular orbital energies or in the analysis of potential energy surfaces [24].
While these potential applications are diverse, they all share a common thread: the need for quick, reliable estimates of polynomial roots. Our work provides a theoretical foundation that could be adapted and refined for these specific contexts, potentially leading to advancements across multiple scientific disciplines. As researchers and practitioners in these fields begin to apply and extend our results, we anticipate a rich interplay between theory and application, driving further innovations in both perturbation theory and its practical implementations.

Author Contributions

Methodology, I.A.N. and M.M.A.A.; Software, J.G.D.; Validation, B.A.; Formal analysis, M.A.A.; Writing—original draft, S.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through the Large Research Project under grant number “RGP2/41/45”, and this study is supported via funding from Prince Sattam bin Abdulaziz University project number (PSAU/2024/R/1445).

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author. The authors declare that no third-party data or copyrighted data has been used.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ogata, K. Modern Control Engineering, 5th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2010. [Google Scholar]
  2. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes: The Art of Scientific Computing, 3rd ed.; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  3. Nilsson, J.W.; Riedel, S.A. Electric Circuits, 10th ed.; Pearson: Upper Saddle River, NJ, USA, 2015. [Google Scholar]
  4. Wilkinson, J.H. Rounding Errors in Algebraic Processes; Courier Corporation: North Chelmsford, MA, USA, 1994. [Google Scholar]
  5. Nayfeh, A.H. Perturbation Methods; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  6. Edelman, A.; Kostlan, E. How many zeros of a random polynomial are real? Bull. Am. Math. Soc. 1995, 32, 1–37. [Google Scholar] [CrossRef]
  7. Tao, T.; Vu, V. Random matrices: Universality of local eigenvalue statistics. Acta Math. 2012, 206, 127–204. [Google Scholar]
  8. Galántai, A.; Hegedüs, C. On the limits of the roots of an algebraic equation with positive coefficients. Numer. Math. 2008, 109, 77–100. [Google Scholar]
  9. Pakdemirli, M.; Yurtsever, H. Estimating roots of polynomials using perturbation theory. Appl. Math. Comput. 2007, 118, 2025–2028. [Google Scholar]
  10. Choque-Rivero, A.E.; Garza, L. Moment perturbation of matrix polynomials. Integral Transform. Spec. Funct. 2015, 26, 177–191. [Google Scholar] [CrossRef]
  11. Maiz, F. Development of the perturbation theory using polynomial solutions. J. Maths. Phys. 2019, 60, 013505. [Google Scholar] [CrossRef]
  12. Dmytryshyn, A. Recovering a perturbation of a matrix polynomial from a perturbation of its first companion linearization. BIT Numer. Math. 2022, 62, 69–88. [Google Scholar] [CrossRef]
  13. Higham, N.J. Accuracy and Stability of Numerical Algorithms; SIAM: Philadelphia, PA, USA, 2002. [Google Scholar]
  14. Ostrowski, A. Recherches sur la méthode de Gräffe et les zeros des polynômes et des series de Laurent. Acta Math. 1940, 99–257. [Google Scholar] [CrossRef]
  15. Chu, E.K. Generalization of the Bauer-Fike Theorem. Numer. Math. 1986, 49, 685–691. [Google Scholar] [CrossRef]
  16. Song, Y. A note on the variation of the spectrum of an arbitrary matrix. Linear Algebra Appl. 2002, 342, 41–46. [Google Scholar] [CrossRef]
  17. Bhatia, R.; Elsner, L.; Krause, G. Bounds for the variation of the roots of a polynomial and the eigenvalues of a matrix. Linear Algebra Appl. 1990, 142, 195–209. [Google Scholar] [CrossRef]
  18. Beauzamy, B. How the roots of a polynomial vary with its coefficients: A local quantitative result. Can. Math. Bull. 1999, 42, 3–12. [Google Scholar] [CrossRef]
  19. Trefethen, L.N. Approximation Theory and Approximation Practice; SIAM: Philadelphia, PA, USA, 2019. [Google Scholar]
  20. Oppenheim, A.V.; Schafer, R.W. Discrete-Time Signal Processing, 3rd ed.; Pearson: Upper Saddle River, NJ, USA, 2015. [Google Scholar]
  21. Burden, R.L.; Faires, J.D.; Burden, A.M. Numerical Analysis, 10th ed.; Cengage Learning: Boston, MA, USA, 2015. [Google Scholar]
  22. Strogatz, S.H. Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  23. Farin, G.; Hoschek, J.; Kim, M.S. (Eds.) Handbook of Computer Aided Geometric Design; Elsevier: Amsterdam, The Netherlands, 2002. [Google Scholar]
  24. Szabo, A.; Ostlund, N.S. Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory; Dover Publications: Mineola, NY, USA, 1996. [Google Scholar]
  25. Sheikh, S.A.; Mir, M.I. Probabilistic zero bounds of certain random polynomials. J. Appl. Math. Stat. Inform. 2024, 20, 53–66. [Google Scholar] [CrossRef]
Figure 1. Distribution of maximum modulus of zeros.
Figure 1. Distribution of maximum modulus of zeros.
Mathematics 12 02993 g001
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

Nafisah, I.A.; Sheikh, S.A.; Alshahrani, M.A.; Almazah, M.M.A.; Alnssyan, B.; Dar, J.G. Perturbation Approach to Polynomial Root Estimation and Expected Maximum Modulus of Zeros with Uniform Perturbations. Mathematics 2024, 12, 2993. https://doi.org/10.3390/math12192993

AMA Style

Nafisah IA, Sheikh SA, Alshahrani MA, Almazah MMA, Alnssyan B, Dar JG. Perturbation Approach to Polynomial Root Estimation and Expected Maximum Modulus of Zeros with Uniform Perturbations. Mathematics. 2024; 12(19):2993. https://doi.org/10.3390/math12192993

Chicago/Turabian Style

Nafisah, Ibrahim A., Sajad A. Sheikh, Mohammed A. Alshahrani, Mohammed M. A. Almazah, Badr Alnssyan, and Javid Gani Dar. 2024. "Perturbation Approach to Polynomial Root Estimation and Expected Maximum Modulus of Zeros with Uniform Perturbations" Mathematics 12, no. 19: 2993. https://doi.org/10.3390/math12192993

APA Style

Nafisah, I. A., Sheikh, S. A., Alshahrani, M. A., Almazah, M. M. A., Alnssyan, B., & Dar, J. G. (2024). Perturbation Approach to Polynomial Root Estimation and Expected Maximum Modulus of Zeros with Uniform Perturbations. Mathematics, 12(19), 2993. https://doi.org/10.3390/math12192993

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