01.12.2014 Views

Variational properties of the zeros of orthogonal polynomials

Variational properties of the zeros of orthogonal polynomials

Variational properties of the zeros of orthogonal polynomials

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Rend. Sem. Mat. Univ. Poi. Torino<br />

Voi. 53,2(1995)<br />

D. Kershaw<br />

VARIATIONAL PROPERTIES OF THE ZEROS<br />

OF ORTHOGONAL POLYNOMIALS<br />

Abstract. The <strong>zeros</strong> <strong>of</strong> P n +i, <strong>the</strong> Legendre polynomial <strong>of</strong> degree n + 1 can be<br />

characterized as being <strong>the</strong> extrema <strong>of</strong> <strong>the</strong> quotient:<br />

J<br />

xp 2 {x) dx<br />

f_..p 2 (x) dx<br />

subject to certain constraints. This result is generalized to include <strong>polynomials</strong><br />

<strong>orthogonal</strong> with respect to a general weight function, and with o<strong>the</strong>r powers <strong>of</strong> <strong>the</strong><br />

variable x. The results are applied to finding ineqùalities for <strong>the</strong> largest zero <strong>of</strong><br />

<strong>orthogonal</strong> <strong>polynomials</strong>.<br />

1. Introduction<br />

Recently, <strong>the</strong> p-version finite element and spectral methods for approximating partial<br />

differential equations have emerged as powerful methods, sometimes yielding exponential<br />

rates <strong>of</strong> convergence. The analysis <strong>of</strong> such methods, and in particular <strong>the</strong> stability, requires<br />

results on <strong>polynomials</strong> that seem to be non-standard. One particular stability result arises in<br />

<strong>the</strong> <strong>the</strong>oretical analysis <strong>of</strong> a posteriori error estimators for p-version finite element method<br />

and was to fìnd an upper bound for <strong>the</strong> quotient:<br />

(1 i) J-i x 2 p 2 (x)dx<br />

S-i P 2 ( X )<br />

where p ranges over ali <strong>polynomials</strong> with real coeffìcients <strong>of</strong> degree n — 1. We shall see<br />

that <strong>the</strong> maximum value is x\ +l where x n+ \ is <strong>the</strong> largest zero <strong>of</strong> P n+ i, <strong>the</strong> Legendre<br />

polynomial <strong>of</strong> degree n + 1, and is attained when<br />

dx<br />

p(x)<br />

rv*J> rp&<br />

The analysis will be used to investigate <strong>the</strong> more general problem <strong>of</strong> finding extrema <strong>of</strong><br />

<strong>the</strong> quotient Q m , given by<br />

, v . f x m p 2 (x) w(x)dx<br />

(1-2) Q m (p) = J \ b : , ' , ' • m = 1,2,<br />

J p 2 (x) w(x)dx


90 D. Kershaw<br />

where w is a nontrivial, non-negative integrable function <strong>of</strong> [a, b] and p ranges over ali real<br />

<strong>polynomials</strong> <strong>of</strong> degree n — 1. The case when m is any positive integer is also amenable<br />

to <strong>the</strong> analysis, but <strong>the</strong> results become progresively more complicated and less interesting<br />

and so will not be investigated in detail.<br />

Some preliminary definitions will be made.<br />

N.<br />

DEFINITION 1.1. a) VM is <strong>the</strong> space <strong>of</strong> real <strong>polynomials</strong> <strong>of</strong> degree not exceeding<br />

b) The orthonormal <strong>polynomials</strong> associated with <strong>the</strong> weight function w will be<br />

denoted by 7T/v, where 7TJV € VM for N = 0,1,... .<br />

e) The <strong>zeros</strong> <strong>of</strong> TTN, ordered in increasing size, are<br />

r(N) IN) (N)<br />

It will always be assumed that <strong>the</strong> coefficient <strong>of</strong> x N in 7T/v is positive. If <strong>the</strong><br />

weight function is distinguished by a suffix <strong>the</strong>n this suffix will be attached as a superfix<br />

to 7r. For example ir N is <strong>the</strong> normalized polynomial <strong>of</strong> degree N which is <strong>orthogonal</strong> to<br />

1, x,..., x^"" 1 with respect to <strong>the</strong> weight function w\.<br />

DEFINITION 1.2. The weights <strong>of</strong> <strong>the</strong> N-point Gaussian quadrature formula<br />

associated with <strong>the</strong> weight function w will be denoted by<br />

r(N) (N) JN) , rr(N) ff(N) rr(N)<br />

(The quadrature points are<br />

T(N) r(N) AN),<br />

The following will be needed, <strong>the</strong> first is a standard result, and <strong>the</strong> second is easily<br />

proved.<br />

LEMMA 1.3. a) H^N) > 0, i = 1,2,...,A^,<br />

b) If a 4- b = 0 and w(x) = w(—x) <strong>the</strong>n<br />

x i w +xW j+1 =0,'i = l,2,...,JV.<br />

2. Rayleigh's Principle<br />

The solution <strong>of</strong> <strong>the</strong> problem is essentially contained in <strong>the</strong> following <strong>the</strong>orem which<br />

follows from Rayleigh's principle, (see for example [2]), however some special cases will<br />

be investigated.


<strong>Variational</strong> Properties <strong>of</strong> <strong>the</strong> Zeros <strong>of</strong> Orthogonal Polynomials 91<br />

THEOREM 2.1. a) The maximum (minimum) value <strong>of</strong><strong>the</strong> quotient<br />

(2.i) Qm{P)=£ *y w ;


92 D. Kersliaw<br />

3. Case m = 1<br />

In this section we restrict our attention to <strong>the</strong> problem when m — I. The following<br />

sections will be devoted to <strong>the</strong> case m = 2. The results in <strong>the</strong> general case become<br />

progressively less interesting with increasing ra, and we shall not dwell on <strong>the</strong>m.<br />

THEOREM 3.1. Let Qi(p) be defined as in (2.1). The we have<br />

(3.1) x^1'<br />

= max Qi(p) = min Qi{q), r = 1,2.... ,n,<br />

p<br />

q<br />

where <strong>the</strong> maximum is to be taken over ali <strong>polynomials</strong> p e V n -i such that<br />

(3.2.) / p(x)<br />

Ja<br />

X — Xs<br />

X<br />

(n) w(x)dx = 0, s = r -1-1, r + 2,..., n,<br />

and thè minimum is taken over <strong>polynomials</strong> q e V n -i such that<br />

(3.3) / q(x) ^ ^ w(x)dx = 0, s = 1,2,... ,r - 1.<br />

«/ a X X g<br />

In each case <strong>the</strong> extremum is attained when<br />

(3.4) p(x) = q(x) = - n(x)<br />

(n)<br />

Pro<strong>of</strong>. We shall give <strong>the</strong> pro<strong>of</strong> in <strong>the</strong> case <strong>of</strong> <strong>the</strong> maximum, <strong>the</strong> one for <strong>the</strong><br />

minimum is similar.<br />

Let p be such that (3.2) is true. The evaluation <strong>of</strong> <strong>the</strong> integrai in (3.2) by <strong>the</strong><br />

appropriate n-point Gaussian quadrature formula gives<br />

ffJ»V(*S">) } = o,<br />

(<strong>the</strong> prime denotes <strong>the</strong> derivative) and so, since Ha and ir' n {xs ) do not vanish, we have<br />

(3.5) p(a4 n) )=0, s = r + l,r + 2,...,n.-<br />

With <strong>the</strong> same Gaussian quadrature formula, <strong>the</strong> quotient Q\(p) now reduces to<br />

which, since H^1<br />

and p 2 {x^1')<br />

are positive, gives<br />

QX(P) < 4 n \<br />

and equality is attained when p satisfies ei<strong>the</strong>r (3.2) or (3.5).<br />

It is worth emphasising <strong>the</strong> result for <strong>the</strong> extreme <strong>zeros</strong>.


<strong>Variational</strong> Properties <strong>of</strong> <strong>the</strong> Zeros <strong>of</strong> Orthogonal Polynomìals 93<br />

COROLLARY 3.2. We have<br />

(3.6) max Q 1 (p) = x^\<br />

pev„.-i<br />

and <strong>the</strong> maximum is attained when<br />

(3.7) p(x) = K n (x)/(x~xM),<br />

(3.8). min Qi(p) = x ( \<br />

P€Vn-l<br />

and <strong>the</strong> minimum is attained when<br />

(3.9) p(x) = ir n (x)/(x-x


94 D. Kerskaw<br />

If <strong>the</strong>re is sufficient symmetry <strong>the</strong> bounds can be attained, as will be seen from <strong>the</strong><br />

following result.<br />

Only <strong>the</strong> pro<strong>of</strong> for <strong>the</strong> maximum will be given.<br />

THEOREM 4.3. Let a — -b, and w(x) = w(-x), <strong>the</strong>n<br />

(4.4) max Q 2 (p) = (x^Y,<br />

Witti equality when<br />

(4.5) p(x)= ' r " +l(l)<br />

Pro<strong>of</strong>. This follows from <strong>the</strong>orem 4.1 since <strong>the</strong> largest zero <strong>of</strong> F2 n +i will be <strong>the</strong><br />

largest zero <strong>of</strong> 7r n+ i. But more directly by lemma 1.3 (b) we have<br />

(xr i) ) 2 =(-a i) ) 2 .<br />

and <strong>the</strong> polynomial given by (4.5) when inserted into Qi will give <strong>the</strong> desired result.<br />

The section closes by giving an upper bound for <strong>the</strong> quotient Q2.<br />

THEOREM 4.4.<br />

Ifa>0<strong>the</strong>n<br />

(4.6) (^ti 1 ') 2 > ggLQM > ( x " n) ) 2 -<br />

Pro<strong>of</strong>. The leading coefficients <strong>of</strong> 7r n and 7r n+ i are positive, and so<br />

TT n (x) > 0 for x > x^\<br />

VlW > ° for x > x( n+l\<br />

and because <strong>of</strong> <strong>the</strong> interlacing <strong>of</strong> <strong>the</strong> <strong>zeros</strong> <strong>of</strong> <strong>the</strong> <strong>polynomials</strong> ir n and 7r n +i we can conclude<br />

that<br />

7r n+1 (s£°) < 0-<br />

The <strong>zeros</strong> <strong>of</strong> each polynomial He in (a, 6), consequently n n {x) and 7r„ + i(x) have opposite<br />

signs for x < 0.<br />

The sign <strong>of</strong> F2 n +i(xn ) is <strong>the</strong> same as that <strong>of</strong> 7r n +i (^n) )7r(-»n n) ), which is<br />

positive.<br />

The sign <strong>of</strong> F 2n+ i(x (^' 1) 1 ) is opposite to that <strong>of</strong> TT n+ i{-x^^ )^n{x^\ )» and<br />

so is negative.<br />

It follows that F2 n +i has a zero in (xn jX^+i )• There are no more <strong>zeros</strong> <strong>of</strong> F2 n +i<br />

above x^^ ' since it is positive <strong>the</strong>re; and so because ali its <strong>zeros</strong> are positive <strong>the</strong> largest<br />

zero will provide <strong>the</strong> upper bound.


<strong>Variational</strong> Propertìes <strong>of</strong> <strong>the</strong> Zeros <strong>of</strong> Orthogonal Polynomials 95<br />

5. Comparison Theorems<br />

In this section we shall apply some <strong>of</strong> <strong>the</strong> above results to fìnding some inequalitìes<br />

satisfied by <strong>the</strong> largest zero <strong>of</strong> <strong>orthogonal</strong> <strong>polynomials</strong>. There will be corresponding results<br />

for <strong>the</strong> smallest zero, <strong>the</strong> statements and pro<strong>of</strong>s for those will not be given since it will be<br />

clear what <strong>the</strong>y should be.<br />

THEO REM 5.1. Let<br />

(5.1) ^(r)<br />

«"(?) =<br />

if<br />

f xp 2 (t) w r (t)dt<br />

$ h a p>(t) w r (t)dt , r = l,2,<br />

(5.2) Q Ì Ì\P)-Q Ì Ì\P) <<br />

for ali <strong>polynomials</strong> p <strong>of</strong> degree not exceeding n — 1, <strong>the</strong>n<br />

(5.3)<br />

.(»»)<br />

(»),<br />

x^(w 1 )-x^(w 2 )<br />

< e.<br />

Pro<strong>of</strong>. The upper inequality gives for p G V<br />

QÌ 1) (p)


96 D. Kershaw<br />

It follows that<br />

(5.7)<br />

x n (w 2 ) x n [w 1 )><br />

The right hand side <strong>of</strong> (5.7) is positive if<br />

which can be written out as<br />

This is equal to<br />

(xq ì q)2{q ì q)i ><br />

(q ì q)2/(yq ì q)i ì<br />

PO '6 nb no<br />

/ / {x-y)q 2 (x)q 2 (y) w 2 (x)wi(y)dijdx.<br />

Ja Ja<br />

fa<br />

nx<br />

(x - y)q 2 {x)q 2 W2{x) W2(y)<br />

{y) dydx,<br />

wi{x) wi(y)<br />

which is positive because <strong>of</strong> (5.4). Consequently vve have <strong>the</strong> required result. •<br />

An inclusion <strong>the</strong>orem is possible when <strong>the</strong> range <strong>of</strong> integration is in [0, oo). In<br />

particular it is applicable to Laguerre <strong>polynomials</strong>.<br />

THEOREM 5.3. Let. a> 0, and W2(x) = xwi(x), <strong>the</strong>n<br />

(5.8) v&RXwx) > 4 n) (-W2) > a£ n) (fi)-<br />

» , .(2)<br />

Pro<strong>of</strong>. Now Xn (^2) is <strong>the</strong> largest zero <strong>of</strong> ir n , which is such that<br />

*6<br />

/ n^(x)x r xwi(x)dx = 0, r = 0,1,... ,71 - 1.<br />

J a<br />

.(2)<br />

Since xir n e V n +i and is <strong>orthogonal</strong> to ali <strong>polynomials</strong> <strong>of</strong> degree n — 1, it must be that<br />

linear combination <strong>of</strong> ir n and ir n +i which vanishes at x = 0. Hence we can write<br />

As in <strong>the</strong> pro<strong>of</strong> <strong>of</strong> <strong>the</strong>orem (5.2) we see that 7Tn (0) and ^^(0) have opposite signs, and<br />

<strong>the</strong> pro<strong>of</strong> continues as before. •<br />

THEOREM 5.4. Let a — —b and wi,W2 be weight functions sudi that<br />

wi(-x) =wi(x), w 2 {-x) = w 2 (x).<br />

If<br />

(5.9)<br />

<strong>the</strong>n<br />

(5.10)<br />

w 2 (x)<br />

wi{x)<br />

w 2 {y)<br />

> 0, 0 < y < x < ò,<br />

wi(y)<br />

x(w 2 ) > 4 n) K)-


<strong>Variational</strong> Propertìes <strong>of</strong> <strong>the</strong> Zeros <strong>of</strong> Orthogonal Polynomials 97<br />

Pro<strong>of</strong>. Just as in <strong>the</strong> pro<strong>of</strong> <strong>of</strong> <strong>the</strong>orem 5.2, in order that<br />

should be true we must have<br />

(5.11) f (\x-y)q 2 {x)q 2 {y)<br />

J-bJ-b<br />

The left hand side <strong>of</strong> (5.11) can be written<br />

«S"W > sW(toi)<br />

w 2 (x)<br />

wi(x)<br />

w 2 (y)<br />

wi(y)<br />

dydx > 0.<br />

-0 px<br />

x PÒ PU PÒ PX<br />

rb p0 pb px<br />

w 2 {x) w 2 (y)<br />

dydx.<br />

-b /<br />

wi(x) wi(y)<br />

+ + (x-y)q 2 (x)q 2 (y)<br />

The middle term -b <strong>of</strong> JO (5.12) J-b is always JO Jo positive because <strong>of</strong> (5.9). The outer terms can be<br />

combined to give <strong>the</strong> expression<br />

'6 px<br />

0 JO<br />

(x - y) [q 2 {x)q 2 (y) - q 2 {-x)q 2 (-y)}<br />

w 2 (x)<br />

wi(x)<br />

The determinant here is positive, and so it remains to show that<br />

q 2 (x)q 2 (y)-q 2 (-x)q 2 (-y)>0.<br />

Now<br />

w 2 (y)<br />

wi(y)<br />

dydx.<br />

7T (1)<br />

q[x) = r^<br />

(n) / \<br />

X-Xh, J {Wi)<br />

where Xn (wi) is <strong>the</strong> largest zero <strong>of</strong> nii .<br />

Suppose first that n is even, <strong>the</strong>n since <strong>the</strong> <strong>zeros</strong> are symmetrically placed about <strong>the</strong><br />

origin (lemma 1.3 (b)) we can write (superfìxes will now be omitted)<br />

q{x) = (x 2 - x\){x 2 -x 2 2)... {x 2 - X^_!)(.T + x m )<br />

=<br />

Jm\X)\X -\- X m ))<br />

where f m is an even function and x m > 0. Then<br />

q(x)q(y) + q(-x)q(-y) = 2f m (x)f m (y) (xy -f xty ,<br />

and<br />

It follows that<br />

q{x)q{y) - q(-x)q(-y) = 2f m {x)f m {y)(x + y)x rn .<br />

q 2 (x)q 2 {y) - q 2 {-x)q 2 (-y) = if^{x)f^{y){xy + x 2 m)(x + y)x mì<br />

which is non-negative for 0 < x,y < b. Hence for n even <strong>the</strong> sum <strong>of</strong> integrals in (5.12) is<br />

positive, and so<br />

x^(w»)<br />

> 4*»(«n).


98 D. Kershaw<br />

The pro<strong>of</strong> when n is odd is similar, <strong>the</strong> only difference being that we have to use<br />

<strong>the</strong> representation <strong>of</strong> q given by<br />

where n ~ 2m + 1.<br />

q(x) = x{x 2 - x\){x 2 -xl)... (x 2 - a4_i)(x + x m ),<br />

•<br />

Acknowledgments. My thanks are due to Dr M. Ainsworth <strong>of</strong> Leicester University for<br />

drawing my attention to <strong>the</strong> problem, and also to <strong>the</strong> referee whose helpful comments have<br />

improved <strong>the</strong> presentation.<br />

REFERENCES<br />

[1] A. ERDÉLYI (et al.) Higher Transcendental Functions, Volume IL McGraw Hill (1953).<br />

[2] J.N. FRANKLIN, Matrix Theory. Prentice-Hall (1968).<br />

[3] G.SZEGÒ, Orthogonal Polynomials. A.M.S. Colloquium Publications, second edition, volume<br />

XXIII (1959).<br />

Donald KERSHAW<br />

Department <strong>of</strong> Ma<strong>the</strong>matics and Statistics, Fylde College,<br />

Lancaster University, Lancaster LAI 4YF, United Kingdom.<br />

Lavoro pervenuto in redazione il 20.10.1994 e, informa definitiva, il 10.4.1995.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!