ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS arxiv: v1 [math.nt] 9 Nov 2013

Transcription

1 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS arxiv:3.2230v [math.nt] 9 Nov 203 Dragan Stankov Abstract. We investigate an infinite sequence of polynomials of the form: a 0 T n (x) + a T n (x) + + a m T n m (x) where (a 0, a,..., a m ) is a fixed m-tuple of real numbers, a 0, a m 0, T i (x) are Chebyshev polynomials of the first kind, n m, m +, m + 2,... Here we analyse the structure of the set of zeros of such polynomial, depending on A and its limit points when n tends to infinity. Also the expression of envelope of the polynomial is given. An application in number theory, more precise, in the theory of Pisot and Salem numbers is presented. 200 Mathematics Subject Classification. Primary B83; Secondary R09, 2D0. Key words and phrases. Chebyshev polynomials, envelope, Pisot numbers, Salem numbers.

2 2 D. STANKOV. Introduction It is well known [6] that the Chebyshev polynomial T n (x) of the first kind is a polynomial in x of degree n, defined by the relation T n (x) cos nθ when x cos θ. Let A (a 0, a,..., a m ) be a (m+)-tuple of real numbers, a 0, a m 0, m. We introduce an infinite sequence of polynomials T n,a (x) a 0 T n (x) + a T n (x) + + a m T n m (x) (n m). We will refer to T n,a (x) as an A-Chebyshev polynomial. We can naturally extend this definition in the case m 0 and A a 0 0: T n,a0 (x) a 0 T n (x). Also, it will be useful to introduce the polynomial P A (x) a 0 x m + a x m + + a m. We will refer to P A (t) as the characteristic polynomial of the A-Chebyshev polynomial. 2. Roots of A-Chebyshev polynomial Let Z n,a denote the set of zeros of the A-Chebyshev polynomial. The aim of this paper is to analyse the structure of the set lim inf Z n,a, when n tends to infinity, depending on A. lim inf Z n,a consists of those elements which are limits of points in Z n,a for all n. That is, x lim inf Z n,a if and only if there exists a sequence of points {x k } such that x k Z k,a and x k x as k.

3 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 3 Example 2.. Using our notation, the Chebyshev polynomial T n (x) is A-Chebyshev polynomial with A. It is well known [5] that the zeros of T n (x) are x n,k cos (nk+ 2 )π n (k, 2,..., n). It is obviously that x n,0 approaches to - and x n,n approaches to when n. Since x n,k are equispaced, it is clear that lim inf Z n, [, ] Example 2.2. What can we say for A-Chebyshev polynomial if A (2, 5, 2)? It is well known [5] that the recurrence relation T n (x) 2xT n (x) T n 2 (x), n 2, 3,... is satisfied. So T n,(2, 5,2) (x) 2T n (x) 5T n (x) + 2T n 2 (x) 4xT n (x) 2T n 2 (x) 5T n (x) + 2T n 2 (x), T n,(2, 5,2) (x) (4x 5)T n (x). Now it is obvious that x 5 is a 4 zero of T n,(2, 5,2) (x), for all n 2, 3,.... So x 5 4 Z n,(2, 5,2) for all n 2, 3,.... Taking into account the previous example we conclude that lim inf Z n,(2, 5,2) [, ] { 5 4 }. Lemma 2.. T n,a (x) 2 (P A(w)w n m +P A (w )w (n m) ) where w x + x 2.

4 4 D. STANKOV Proof. Starting from the definition of A-Chebyshev polynomial and using well known [5] formula T n (x) 2 (wn + w n ) we have: T n,a (x) a i T n i (x) a i 2 (wn i + w n+i ) m 2 ( a i w n i + a i w n+i ) m 2 (wn m a i w m i + w n+m m a i w m+i ) 2 (wn m P A (w) + w n+m P A (w )). One can calculate that if w x + x 2 then x 2 (w + w ). So, from the previous lemma we can deduce next Corollary 2.. If there iz w such as P A (w) P A (w ) 0 then T n,a (x) 0 for x 2 (w + w ) and for all n m. In the previous example we can see that 2, are roots of P 2 A(x) 2x 2 5x + 2, therefore x (2 + ) 5 is a zero of T n,a(x) for all n 2. For the next corollary we need the following definition: the set T of Salem numbers is the set of real algebraic integers τ greater than, such that all its conjugate roots have modulus at most equal to, one at least having a modulus equal to. Corollary 2.2. If τ is a Salem number and P A (x) is its minimal polynomial then T n,a (x) 0 for x 2 (τ + τ ) and for all n m.

5 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 5 The claim is a direct consequence of a well known property of a Salem number [2] that P A (τ) P A (τ ) 0. Theorem. If there iz a root ω, out of the unit circle, of the polynomial P A, that is P A (ω) 0, ω >, then for every real number ε > 0, there exists a natural number n 0 such that for all n > n 0, there is a root ξ of the A-Chebyshev polynomial T n,a (x) such that ξ (ω + 2 ω ) < ε. Proof. It is convenient to use Lemma 2. to express T n,a (x) P 2 A(w)w n m + P 2 A(w )w (n m) where w x + x 2, or equivalently x x(w) 2 (w + w ). Since x(w) is continuous for w > 0, there is δ > 0 such that if w ω < δ then 2 (w+w ) 2 (ω +ω ) < ε. We can take an δ 2 < ω such that, in the circle {z : z ω δ 2 }, there is no root of P A (w) which is different from ω. Let δ min(δ, δ 2 ) and C {z : z ω δ}. Since C, the boundary of C, is a compact set, P A (w), P A (w ) are continuous on C, there is w min where P A (w) attains its minimum, and w max where P A (w ) attains its maximum on C. Since 2 P A(w max) is constant and ω δ >, there is n 0 such that 2 P A(w min ) ( ω δ) n 0 m > 2 P A(w max). For n n 0, let us denote f(w) 2 P A(w)w n m, g(w) 2 P A(w )w (n m). This notation corresponds to Rouché s theorem which we intend to use. We have to prove

6 6 D. STANKOV that f(w) > g(w) on C. Since w ω δ > we have on C: f(w) 2 P A(w) w n m 2 P A(w min ) ( ω δ) n 0 m > 2 P A(w max) 2 P A(w ) w (n m) g(w). The conditions in Rouché s theorem are thus satisfied. Consequently, since f(w) has root ω, we conclude that f(w) + g(w) has a root, let it be ω, inside the circle C. Clearly, since ω ω < δ, if we denote ξ (ω 2 + ω ), we conclude ξ (ω + 2 ω ) < ε. Finally, we conclude that T n,a (ξ) 2 P A(ω )ω n m + 2 P A(ω )ω (n m) f(ω ) + g(ω ) 0. Theorem 2. If x [, ], then for every real number ε > 0, there exists a natural number n 0 such that for all n > n 0, there is a root ξ of the A-Chebyshev polynomial T n,a (x) such that x ξ < ε. Proof. Directly from the definitions of the Chebyshev polynomial and the A-Chebyshev polynomial we can show that (2.) T n,a (x) a 0 cos nθ+a cos(n )θ+ +a m cos(n m)θ, n m, when x cos θ. Since a k cos(n k)θ a k cos(n m + m k)θ a k (cos(n m)θ cos(m k)θ sin(n m)θ sin(m k)θ) the equation T n,a (x) 0 is equivalent with cos(n m)θ a k cos(m k)θ sin(n m)θ a k sin(m k)θ. k0 k0

7 Finally we get ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 7 tan(n m)θ m k0 a k cos(m k)θ m k0 a k sin(m k)θ. The function on the right, let call it R(θ), does not depend on n. The graph of tan(n m)θ consists of parallel equispaced tangents branches. So if we take n : 2n m we double n m and get a new graph which is actually the union of the old one with branches settled in the middle of each pair of neighbouring branches of the old graph. We conclude that all roots of tan(n m)θ R(θ), remain to be the roots of tan 2(n m)θ R(θ), and new roots interlace with old. Finally, changing variables θ arccos x will preserve order and denseness of the roots. 3. Envelope of an A-Chebyshev polynomial Let us observe the Chebishev polynomial T n (x) again. It is well known that, for any n, the graph of the polynomial oscillates between - and when x [, ]. As n increases we have more and more oscillations. Something like that we have in the case an A-Chebyshev polynomial. Example 3.. Let A (, 0, ), so T n,a (x) T n (x) + T n 2 (x) 2xT n (x) T n 2 (x) + T n 2 (x) 2xT n (x). Now it is obvious that T n,a (x) oscillates between lines y ±2x, for x [, ]. We will refer to these lines as an envelope of the A-Chebyshev polynomial. Using the expression (2.) we can study the following Example 3.2. Let A (, 0, ), so T n,a (x) cos nθ cos((n 2)θ) 2 sin((n )θ) sin θ 2 sin((n )θ) cos 2 θ 2 sin((n )θ) x 2. Now it is obvious that T n,a (x) oscillates between upper and

8 8 D. STANKOV lower half of the ellipse y ±2 x 2, for x [, ]. These halves constitute the envelope of the A-Chebyshev polynomial in this case. With the same technique we can find the envelope in the next Example 3.3. Let A (, ), so T n,a (x) cos nθ cos((n )θ) 2 sin((n 2 )θ) sin θ 2 2 sin((n 2 )θ) cos θ 2 2 sin((n )θ) x. Now it is obvious that the envelope of the A-Chebyshev 2 polynomial is a parabola y ± 2 x, for x [, ]. Using previous examples, we can formulate the characteristics that an envelope of the A-Chebyshev polynomial must have. (Env) The Envelope depends only on A. If A is fixed, it is unique for T n,a (x), n N. (Env2) The Envelope is a non negative function. (Env3) The A-Chebyshev polynomial is not greater in modulus than the envelope, x [, ]. (Env4) The Envelope is a smooth function except at its zeros. (Env5) If the envelope and the A-Chebyshev polynomial have equal positive value in x they have also equal the first derivative in x. We shall define the envelope of the A-Chebyshev polynomial as a function which satisfies the characteristics (Env)-(Env5). It is naturally to ask how can we find the envelope for an A-Chebyshev polynomial. The next lemma will be useful. Lemma 3.. Let R(t), I(t) be real differentiable functions of real argument, with R (t), I (t) continuous, E(t) R 2 (t) + I 2 (t), t R. Then following three statements are satisfied:

9 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 9 (i) R(t) E(t), (ii) R(t) E(t) if and only if I(t) 0, (iii) if I(t) 0 and R(t) > 0 then R(t) E(t) and R (t) E (t). Proof. The first and the second statements are straightforward. To demonstrate that R (t) E (t) we need to determinate E (t) 2R(t)R (t)+2i(t)i (t) 2 R 2 (t)+i 2 (t). Using I(t) 0, R(t) > 0 we get the claim. Theorem 3. The envelope E A (x) for an A-Chebyshev polynomial T n,a (x) is the square root of the modulus of m m 2 a 2 i + 2 a i a i+ T (x) + 2 a i a i+2 T 2 (x) + or in more compact form m k + 2 a i a i+k T k (x) + + 2a 0 a m T m (x), E A (x) a i a k T i k (x). k0 Proof. Let z A (t) a 0 cos(nt) + a cos((n )t) + + a m cos((n m)t) + i(a 0 sin(nt) + a sin((n )t) + + a m sin((n m)t)) be an auxiliary function on t R. We can see that T n,a (x) Re(z A (t)) so T n,a (x) 2 z A (t) 2, x cos(t). We will show that E A (x) z(t). ( m 2 ( m z A (t) 2 a k cos((n k)t)) + a k sin((n k)t) k0 k0 ) 2 a i a k cos((n i)t) cos((n k)t)+ a i a k sin((n i)t) sin((n k)t) k0 k0

10 0 D. STANKOV a i a k (cos((n i)t) cos((n k)t) + sin((n i)t) sin((n k)t)) k0 a i a k cos((i k)t). k0 If we substitute x cos(t) in cos((i k)t) we get T i k (x). Since E A (x) does not depend on n (Env) is fulfilled. (Env2), (Env3), (Env4) are straightforward. Using previous lemma if R(t) Re(z(t)), I(t) Im(z(t)), we can easily obtain (Env5). It is useful to calculate the envelope of the A-Chebyshev polynomial for m, 2, 3, 4. Actually, we give the calculation of the square of the envelope, to avoid cumbersome square roots. Using previous formula we have: (m) a a 2 + 2a 0 a x; (m2) a a 2 + a (a 0 a + a a 2 )x + 2a 0 a 2 (2x 2 ) a a 2 + a 2 2 2a 0 a 2 + (2a 0 a + 2a a 2 )x + 4a 0 a 2 x 2 ; (m3) a a 2 + a a (a 0 a + a a 2 + a 2 a 3 )x + 2(a 0 a 2 + a a 3 )(2x 2 ) + 2a 0 a 3 (4x 3 3x) a a 2 + a a 2 3 2a 0 a 2 2a a 3 + (2a 0 a + 2a a 2 + 2a 2 a 3 6a 0 a 3 )x + (4a 0 a 2 + 4a a 3 )x 2 + 8a 0 a 3 x 3 ; (m4) a 2 0+a 2 +a 2 2+a 2 3+a 2 4+2(a 0 a +a a 2 +a 2 a 3 +a 3 a 4 )x+2(a 0 a 2 +a a 3 + a 2 a 4 )(2x 2 )+2(a 0 a 3 +a a 4 )(4x 3 3x)+2a 0 a 4 (8x 4 8x 2 +) a a 2 + a a a 2 4 2a 0 a 2 2a a 3 2a 2 a 4 + (2a 0 a + 2a a 2 + 2a 2 a 3 + a 3 a 4 6a 0 a 3 6a a 4 )x + (4a 0 a 2 + 4a a 3 + 4a 2 a 4 6a 0 a 4 )x 2 + (8a 0 a 3 + 8a a 4 )x 3 + 6a 0 a 4 x 4. Remark 3.. There is a connection with the theory of signal processing. The analytic signal z(t) can be expressed in terms of complex polar

11 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS coordinates, z(t) f(t)+i ˆf(t) A(t)e iφ(t) where A(t) f 2 (t) + ˆf 2 (t), and φ(t) arctan ˆf(t). These functions are respectively called the am- f(t) plitude envelope and instantaneous phase of the signal, ˆf(t) is Hilbert transform of f(t). 4. Connection between the envelope and the characteristic polynomial Until now we used the envelope to describe the graph of the A- Chebyshev polynomial if x is of modulus not greater than. If x > we preferred the characteristic polynomial P A (x). It is natural to ask, is there any connection between the envelope and the characteristic polynomial of the A-Chebyshev polynomial. The next theorem shows that the answer is affirmative. Theorem 4. The envelope of the A-Chebyshev polynomial is the function PA E A (x) (x + x 2 )P A (x x 2 ). Proof. We shall start from the compact form of the envelope given in Theorem 3 and use well known [5] formula T n (x) 2 (wn +w n ) where

12 2 D. STANKOV w x + x 2. E A (x) a i a k T i k (x) k0 a i a k 2 (wi k + w (i k) ) k0 K0 I0 2 a Ka I w K I + 2 a ia k w k i. k0 (Here we renamed i with K and k with I in the first double sum. Now we shall switch the order of summing in the first double sum and apply obvious a K a I a I a K.) I0 K0 m 2 2 a Ia K w K I + 2 a ia k w k i k0 a i a k w m i w m+k k0 a i w m i m k0 a k w m+k 2 a ia k w k i k0 P A (w)p A (w ). PA (x + x 2 )P A (x x 2 ).

13 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 3 Figure. shows graphs of A-Chebyshev polynomials of the first kind T 4,(,0,0,) (x), T 44,(,0,0,) (x) together with their common envelope E(x) 2 + 6x 8x3. Figure. 5. A-Chebyshev polynomial of the second kind It is well known [6] that the Chebyshev polynomial U n (x) of the second kind is a polynomial in x of degree n, defined by the relation U n (x) sin(n+)θ sin θ when x cos θ. Let A (a 0, a,..., a m ) be a (m+)-tuple of real numbers, a 0, a m 0, m. We introduce an infinite sequence of polynomials U n,a (x) a 0 U n (x) + a U n (x) + + a m U n m (x) (n m).

14 4 D. STANKOV We will refer to U n,a (x) as an A-Chebyshev polynomial of the second kind. We can naturally extend this definition in the case m 0 and A a 0 0: U n,a0 (x) a 0 U n (x). We will refer to the polynomial P A (x) a 0 x m + a x m + + a m as the characteristic polynomial of the A-Chebyshev polynomial. Lemma 5.. U n,a (x) w w (w n+ m P A (w) w n +m P A (w )) where w x + x 2. Proof. Starting from the definition of A-Chebyshev polynomial and using well known [5] formula U n (x) wn+ w n w w we have: U n,a (x) a i U n i (x) a i w n+ i w n +i w w w w ( m a i w n+ i (wn+ m w w a i w n +i ) a i w m i w n +m m w w (wn+ m P A (w) w n +m P A (w )). a i w m+i ) Theorem 5. If there is a root ω, out of the unit circle, of the polynomial P A, that is P A (ω) 0, ω >, then for every real number ε > 0,

15 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 5 there exists a natural number n 0 such that for all n > n 0, there is a root ξ of the A-Chebyshev polynomial of the second kind U n,a (x) such that ξ (ω + 2 ω ) < ε. Proof. It is convenient to use the previous lemma to express U n,a (x) w w (w n+ m P A (w) w n +m P A (w )) where w x + x 2 or equivalently x x(w) 2 (w +w ). Since x(w) is continuous for w > 0, there is δ > 0 such that if w ω < δ then 2 (w+w ) 2 (ω+ω ) < ε. We can take an δ 2 < ω such that, in the circle {z : z ω δ 2 }, there is no root of P A (w) which is different from ω. Let δ min(δ, δ 2 ) and C {z : z ω δ}. Since C, the boundary of C, is a compact set, P A (w), P A (w ) are continuous on C, there is w min where P A (w) gets its minimum and w max where P A (w ) gets its maximum on C. Since w w P A (w max) is constant and ω δ >, there is n 0 such that w w P A (w min ) ( ω δ) n 0+ m > w w P A (w max). us denote f(w) w w w n+ m P A (w), g(w) For n n 0 let w n +m P w w A (w ). This notation corresponds with Rouché s theorem which we intend to use. We have to prove that f(w) > g(w) on C. Since w ω δ > we have on C: f(w) > w w P A(w) w n+ m w w P A(w min ) ( ω δ) n 0+ m w w P A(wmax) w w P A(w ) w (n+ m) g(w).

16 6 D. STANKOV The conditions in Rouché s theorem are thus satisfied. Consequently, since f(w) has root ω, we conclude that f(w) + g(w) has a root, let it be ω, inside the circle C. Clearly, since ω ω < δ, if we denote ξ (ω 2 + ω ), we conclude ξ (ω + 2 ω ) < ε. Finally, we conclude that U n,a (ξ) w w P A (ω )ω n+ m w w P A (ω )ω (n+ m) f(ω ) + g(ω ) An application in number theory Recall that q > is a Pisot number if q is an algebraic integer, whose other conjugates are of modulus strictly less than. Salem proved that every Pisot number is a limit point of the set T of Salem numbers. Let Q(x) x m P ( x ) a 0 + a x + + a m x m be the reciprocal polynomial of the polynomial P (x) a 0 x m + a x m + + a m. Salem showed that if P (x) is the minimal polynomial of a Pisot number q then R k (x) x k P (x) + Q(x) is polynomial with a root τ k that is a Salem number, and the limit of the sequence τ k is q, k. There is a connection between Salem sequence R k (x) with A-Chebyshev polynomials T n,a (x) the characteristic polynomial of which is P (x). We have seen that T n,a ( 2 (w + w )) T n,a (x) m a it n i (x) m a i 2 (wn i + w n+i ). Now we can show that 2w n T n,a ( 2 (w + w )) R 2n m (w). Really, we obtain 2w n T n,a ( 2 (w + w )) m a i(w 2n i + w i ) w 2n m m a i(w m i ) + m a i(w i ) w 2n m P (w) + Q(w). The question is what is going on if we use A-Chebyshev polynomial of the second kind instead of T n,a (x). We will demonstrate that one more sequence of Salem numbers, which converges to the Pisot number q, appears. It is obvious that U n,a ( 2 (w+w )) U n,a (x) m a iu n i (x)

17 m ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 7 a i (w n+ i w n +i ) m w w a i(w n i + w n i 2 + w n i 4 + w n+i+2 + w n+i )). We claim that S 2n (w) w n U n,a ( 2 (w + w )) is polynomial of degree 2n with a root τ 2n that is a Salem number, and the limit of the sequence τ 2n is q, n. Using the previous theorem it is clear that there is a root τ 2n of S 2n (w) such as τ 2n q. It is obvious that S 2n (w) is a reciprocal polynomial. It remains to be proved that all other roots of S 2n (w) are in the unit circle. We shall apply the method Salem (communicated by Hirschman) used to prove the same property of his sequence R k (x) [7]. Using Lemma 5. we have S 2n (w) w n U n,a ( 2 (w + w )) w n w w (wn+ m P A (w) w n +m P A (w )) wn+ w 2 (wn+ m P A (w) w n +m P A (w )) w 2 (w2n+2 m P A (w) w m P A (w )) w 2 (w2n+2 m P A (w) Q(w)). We denote by ɛ a positive number and consider the equation ( + ɛ)w 2n+2 m P A (w) Q(w) 0. Since for w we have P (w) Q(w), it follows by Rouché s theorem that inside the circle w the number of roots of the last equation is equal to the number of roots of w 2n+2 m P (w), that is, (2n+2 m)+m. As ɛ 0, these roots vary continuously. Hence, for ɛ 0 we have 2n + roots with modulus. It is obvious that two roots are, so the

18 8 D. STANKOV fraction can be reduced with w 2. Finally we conclude that at most one root of S 2n (w) is outside the unit circle. Example 6.. Let q be the golden ratio, the greater root of P (x) x 2 x. Then R k (w) w k+2 w k+ w k w 2 w + is the Salem s sequence. Our sequence of Salem numbers τ 2m which converge to q consists of the greatest in modulus roots of the polynomials S 2n (w) w 2 (w2n+2 m P A (w) Q(w)) w 2 (w2n (w 2 w ) + w 2 + w ) (w 2n + w 2n 2 + w 2n 4 + w 2 + ) (w 2n + w 2n 3 + w 2n 5 + w 3 + w) (w 2n 2 + w 2n 4 + w 2n 6 + w 4 + w 2 ). Finally S 2n (w) w 2n (w 2n + w 2n 3 + w 2n 5 + w 3 + w) +. References [] H. Bavinck, On the zeros of certain linear combinations of Chebyshev polynomials, J. Comput. Appl. Math., 65, 9-26 (995). [2] M.-J. Bertin, A. Decomps-Guilloux, M. Grandet-Hugot, M. Pathiaux- Delefosse, and J.-P. Schreiber, Pisot and Salem numbers, Birkhäuser, Basel, (992). [3] J. Cigler, A simple approach to q-chebyshev polynomials, arxiv: [4] A. Dubickas, Totally real algebraic integers in small intervals, Lith. Math. J., 40, (2000). [5] J.C. Mason and D.C. Handscomb, Chebyshev Polynomials, Rhapman & Hall CRC, Boca Raton London New York Washington D.C. (2003). [6] T.J. Rivlin, The Chebyshev Polynomials, John Wiley, New York (974).

19 ON LINEAR COMBINATIONS OF CHEBYSHEV POLYNOMIALS 9 [7] R. Salem, Algebraic Numbers and Fourier Analysis, Heath, Boston (963). [8] D. Stankov, A class of discrete spectra of non-pisot numbers, Publ. Inst. Math., 83, 9-4 (2008). [9] D. Stankov, On spectra of neither Pisot nor Salem algebraic integers, Monatsh. Math., 59, 5-3 (200). Katedra Matematike RGF-a, Universitet u Beogradu, 000 Beograd, Dušina 7, Serbia address: dstankov@rgf.bg.ac.rs