Zeros of Some Self-Reciprocal Polynomials: David Joyner 2010-11-12
Zeros of Some Self-Reciprocal Polynomials: David Joyner 2010-11-12
David Joyner∗
2010-11-12
1 Introduction
This talk is about zeros of a certain family of polynomials which arise nat-
urally in several areas of mathematics. We are especially interested in poly-
nomials which have all their zeros in the unit circle
S 1 = {z ∈ C | |z| = 1}.
Let p be a polynomial
p(z) = a0 + a1 z + . . . + an z n ai ∈ C, (1)
and let p∗ denote1 the reciprocal polynomial or reverse polynomial
1
2 “Smoothness” of roots
A natural question to ask about zeros of polynomials is how “smoothly” do
they vary as a function of the coefficients of the polynomial?
To address this, suppose that the coefficients ai of the polynomial p are
functions of a real parameter t. Abusing notation slightly, identify p(z) =
p(t, z) with a function of two variables (t ∈ R, z ∈ C). Let r = r(t) denote a
root of this polynomial, regarded as a function of t:
p(t, r(t)) = 0.
Example 2 Let
p(z) = 1 + (1 + t) · z + z 2 ,
so we may take
p
−1 − t + (1 + t)2 − 4
r(t) = .
2
Note that r(t) is smooth provided t lies in an interval which does not contain
1 or −3. We can directly verify the lemma holds in this case. Observe (for
later) that if −3 < t < 1 then |r(t)| = 1.
2
Let p(z) = p(t, z) and r = r(t) be as before. Consider the distance
function
d(t) = |r(t)|
of the root r. Another natural question is: How smooth is the distance
function of a root as a function of the coefficients of the polynomial p?
The analog to Lemma 1 holds, with one extra condition.
3 Self-reciprocal polynomials
Let
Rm = {p ∈ R[z]m | p = p∗ }
3
Figure 1: Size of largest root of the polynomial 1 + (1 + t)z + z 2 , −5 < t < 3.
cm (z) = 1 + z + . . . + z m .
If m is even then cm does not have ±1 as roots. Many results in the theory
fall into the following category.
p(z) = a0 + a1 z + . . . + a2n z 2n , ai ∈ R.
Lemma 5 ([DH], §2.1; see also [L2]) The polynomial p ∈ R[z]2n is self-
reciprocal if and only if it can be written
4
n
Y
p(z) = a2n · (1 − αk z + z 2 ), (2)
k=1
Example 6 Note
1 + z + z 2 + z 3 + z 4 = (1 + φ · z + z 2 )(1 + φ · z + z 2 ),
√ √
where φ = 1+2 5 = 1.618 . . . is the “golden ratio,” and φ = 1− 5
2
= −0.618...
is its “conjugate.”
e0 (X1 , X2 , . . . , Xn ) =P
1,
e1 (X1 , X2 , . . . , Xn ) = 1≤j≤n Xj ,
P
e2 (X1 , X2 , . . . , Xn ) = 1≤j<k≤n Xj Xk ,
P
e3 (X1 , X2 , . . . , Xn ) = 1≤j<k<l≤n Xj Xk Xl ,
..
.
en (X1 , X2 , . . . , Xn ) = X1 X2 · · · Xn .
5
Lemma 8 For each n ≥ 1 and αi ∈ C, we have
n
Y 2n
X
2
(z − αk z + 1) = c2n,k z k ,
k=1 k=1
for 0 ≤ k ≤ n.
Here’s another one of those “if it’s near cm then its roots are in S 1 ”-type
results.
has all its roots on S 1 , provided the coefficients satisfy the following condition
m
X
|am | ≥ |aj − am |.
j=0
4 Other characterizations
There are several other characterizations, due to Cohn, Chinen, Chen, and
Fell, which we discuss next.
6
Theorem 12 (Schur-Cohn4 ) Let p ∈ C[z]m be as in (1). Cohn showed that
p has all its zeros on S 1 if and only if
(a) there is a µ ∈ S 1 such that, for all k with 0 ≤ k ≤ n, we have an−k =
µ · ak , and
(b) all the zeros of p0 lie inside or on S 1 .
Corollary 13 p ∈ Rm has all its zeros on S 1 if and only if all the zeros of
p0 lie inside or on S 1 .
This result was proven by Chen [Chen] and later independently by Chi-
nen6 [Ch]. It provides a very large class of self-reciprocal polynomials having
roots on the unit circle.
proof: We shall adapt some ideas from Chinen [Ch] for our argument.
The proof requires recalling the following well-known theorem, discovered
independently by Eneström (in the late 1800’s) and Kakeya (in the early
1900’s).
4
See for example, Chen [Chen], §1.
5
In fact, both are exercises in Marden [Ma].
6
Actually, Chinen’s version is slightly stonger, and it is that version which we are
stating.
7
Theorem 16 (Eneström-Kakeya, version 1) Let f (z) = a0 +a1 z +. . .+ak z k
satisfy a0 > a1 > . . . > ak > 0. Then f (z) has no roots in |z| ≤ 1.
Theorem 17 (Fell) Let P0 (z) and P1 (z) be real monic polynomials of degree
n having zeros in S 1 − {1, −1}. Denote the zeros of P0 (z) by w1 , w2 , . . . , wn
and of P1 (z) by z1 , z2 , . . . , zn . Assume
8
wi 6= zj ,
This theorem is used in the “heuristic argument” given at the end of the
paper.
The autocorrelations of the sequence {ai }ni=0 are the elements of the se-
quence given by
9
n−k
X
ck = ck (p) = aj aj+k , 0 ≤ k ≤ n. (4)
j=0
10
Example 19 The curve
y 2 = x7 − x,
G(z) = N1 z + N2 z 2 /2 + N3 z 3 /3 + · · · .
11
Sage
12
If · denotes the usual inner product,
v · w = v1 w1 + . . . + vn wn ,
C ⊥ = {v ∈ V | v · c = 0, ∀c ∈ C}.
We say C is self-dual if C = C ⊥ .
For each vector v ∈ V , let
supp(v) = {i | vi 6= 0}
denote the support of the vector. The weight of the vector v is wt(v) =
|supp(v)|. The weight distribution vector or spectrum of a code C ⊂ Fn is the
vector
The connection between the weight enumerator of C and that of its dual
is very close, as the following well-known result shows.
13
Definition 24 A polynomial P = PC for which
(xT + (1 − T )y)n AC (x, y) − xn n−d
P (T ) = . . . + T + ... .
(1 − T )(1 − qT ) q−1
is called a Duursma zeta polynomial of C. The Duursma zeta function is
defined in terms of the zeta polynomial by means of
P (T )
ζC (T ) = ,
(1 − T )(1 − qT )
It can be shown that if C1 and C2 are “equivalent” codes then the cor-
responding zeta polynomials are equal. Therefore, these polynomials can be
used to help classify codes.
proof: This is proven in the appendix to Chinen [Ch]. Here is the rough
+y(1−T ))n
idea. Expand (xT(1−T )(1−qT )
in powers of T . The Duursma polynomial is
a polynomial of degree n + 2 − d − d⊥ . Provided d⊥ ≥ 2, the Duursma
polynomial can be written as P (T ) = a0 + a1 T + . . . + an−d T n−d . Now, use
b0,0 b1,0 . . . bn−d,0
an−d An /(q − 1)
0 b1,1 . . . bn−d,1
an−d−1 An−1 /(q − 1)
0 0 b2,2 ...
.. =
.. .
(6)
.. .. .. . .
.
. .
a0 Ad /(q − 1)
0 . . . 0 bn−d,n−d
14
If C is self-dual (i.e., C = C ⊥ ), the Duursma zeta polynomial satisfies a
functional equation of the form
1
P (T ) = q g T 2g P ( ).
qT
√
After a suitable change-of-variable, replacing T by T / q, these polynomials
are self-reciprocal.
In general, the analog of the Riemann hypothesis for curves does not hold
for the Duursma zeta polynomials of self-dual codes. However, when it does
√
hold, after a suitable change-of-variable (replacing T by T / q), we see that
self-dual codes give rise to a class of example of polynomials having roots on
the unit circle.
sage: MS = MatrixSpace(GF(2),4,8)
sage: G = MS([[1,1,1,1,0,0,0,0],[0,0,1,1,1,1,0,0],[0,0,0,0,1,1,1,1],[1,0,1,0,1,0,1,0]])
sage: C = LinearCode(G)
sage: C == C.dual_code()
True
sage: C.zeta_polynomial()
2/5*Tˆ2 + 2/5*T + 1/5
sage: C.<z> = PolynomialRing(CC, "z")
sage: f = (2*zˆ2+2*z+1)/5
sage: rts = f.roots()
sage: [abs(z[0]) for z in rts]
[0.707106781186548, 0.707106781186548]
sage: RR(sqrt(2))
1.41421356237310
sage: RR(1/sqrt(2))
0.707106781186548
PC (T ) = (2T 2 + 2T + 1)/5
15
7.1 Duursma’s conjecture
There is an infinite family of Duursma zeta functions for which Duursma
conjectures that the analog of the RH always holds. The linear codes used
to construct these zeta functions are so-called “extremal self-dual codes.”
Although the construction of these codes is fairly technical (see [JK] for
an expository treatment), we can give some examples.
Example 27 Let
p(z) = a0 + a1 z + . . . + an z n
√
denote the normalized Duursma zeta polynomial, so p(z) = P (z/ q) is self-
reciprocal.
Some examples of the lists of coefficients a0 , a1 , . . ., computed using Sage.
We have normalized the coefficients so that they sum to 10 and represented
the rational coefficients as decimal approximations to give a feeling for their
relative sizes.
• Case Type I:
m = 2: [1.1309, 2.3990, 2.9403, 2.3990, 1.1309]
m = 3: [0.45194, 1.2783, 2.0714, 2.3968, 2.0714, 1.2783, 0.45194]
m = 4: [0.18262, 0.64565, 1.2866, 1.8489, 2.0724, 1.8489, 1.2866,
0.64565, 0.18262]
16
m = 4: [0.26170, 0.75545, 1.3085, 1.7307, 1.8874, 1.7307, 1.3085,
0.75545, 0.26170]
• Case Type IV:
m = 2: [2.8571, 4.2857, 2.8571]
m = 3: [1.6667, 3.3333, 3.3333, 1.6667]
m = 4: [0.97902, 2.4476, 3.1469, 2.4476, 0.97902]
17
then Pt (z) is a multiple of 1 + z + . . . + z d−1 .
Do any of the polynomials Pt (z) have multiple roots (0 < t < 1)? Using
the notation of §2, in the case p(t, z) = Pt (z), we have
P1 (r(t)) − P0 (r(t))
r0 (t) = −pt (t, r(t))/pz (t, r(t)) = .
Pt0 (r(t))
References
[Chen] W. Chen, On the polynomials with all there zeros on the unit circle,
J. Math. Anal. and Appl. 190 (1995)714-724.
[DH] S. DiPippo, E. Howe, Real polynomials with all roots on the unit circle
and abelian varieties over finite fields, J. Number Theory 78(1998)426-
450.
Available: http://arxiv.org/abs/math/9803097
18
[L2] ——, On zeros of reciprocal polynomials, Publ. Math. Debrecen 61,
(2002), 645661.
[S] W. Stein and the Sage group, Sage - Mathematical Software, version
4.5, 2010.
http://www.sagemath.org/.
19