CHEBYSHEV POLYNOMIALS J.C. MASON D.C. HANDSCOMB
CHAPMAN & HALL/CRC
A CRC Press Company Boca Raton London New York Washington, D.C.
© 2003 by CRC Press LLC Library of Congress Cataloging-in-Publication Data
Mason, J.C. Chebyshev polynomials / John C. Mason, David Handscomb. p. cm. Includes bibliographical references and index. ISBN 0-8493-0355-9 (alk. paper) 1. Chebyshev polynomials. I. Handscomb, D. C. (David Christopher) II. Title.. QA404.5 .M37 2002 515′.55—dc21 2002073578
This book contains information obtained from authentic and highly regarded sources. Reprinted material is quoted with permission, and sources are indicated. A wide variety of references are listed. Reasonable efforts have been made to publish reliable data and information, but the author and the publisher cannot assume responsibility for the validity of all materials or for the consequences of their use. Neither this book nor any part may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, microfilming, and recording, or by any information storage or retrieval system, without prior permission in writing from the publisher. The consent of CRC Press LLC does not extend to copying for general distribution, for promotion, for creating new works, or for resale. Specific permission must be obtained in writing from CRC Press LLC for such copying. Direct all inquiries to CRC Press LLC, 2000 N.W. Corporate Blvd., Boca Raton, Florida 33431. Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are used only for identification and explanation, without intent to infringe. Visit the CRC Press Web site at www.crcpress.com
© 2003 by John C. Mason and David Handscomb
No claim to original U.S. Government works International Standard Book Number 0-8493-0355-9 Library of Congress Card Number 2002073578 Printed in the United States of America 1 2 3 4 5 6 7 8 9 0 Printed on acid-free paper
© 2003 by CRC Press LLC In memory of recently departed friends Geoff Hayes, Lev Brutman
© 2003 by CRC Press LLC Preface
Over thirty years have elapsed since the publication of Fox & Parker’s 1968 text Chebyshev Polynomials in Numerical Analysis. This was preceded by Snyder’s brief but interesting 1966 text Chebyshev Methods in Numerical Ap- proximation. The only significant later publication on the subject is that by Rivlin (1974, revised and republished in 1990) — a fine exposition of the theoretical aspects of Chebyshev polynomials but mostly confined to these aspects. An up-to-date but broader treatment of Chebyshev polynomials is consequently long overdue, which we now aim to provide. The idea that there are really four kinds of Chebyshev polynomials, not justtwo,hasstronglyaffectedthecontentofthisvolume.Indeed,theprop- ertiesofthefourkindsofpolynomialsleadtoanextendedrangeofresultsin many areas such as approximation, series expansions, interpolation, quadra- ture and integral equations, providing a spur to developing new methods. We do not claim the third- and fourth-kind polynomials as our own discovery, but we do claim to have followed close on the heels of Walter Gautschi in first adopting this nomenclature. Ordinary and partial differential equations are now major fields of applica- tion for Chebyshev polynomials and, indeed, there are now far more books on ‘spectral methods’ — at least ten major works to our knowledge — than on Chebyshev polynomials per se. This makes it more difficult but less essential to discuss the full range of possible applications in this area, and here we have concentrated on some of the fundamental ideas. We are pleased with the range of topics that we have managed to in- clude. However, partly because each chapter concentrates on one subject area, we have inevitably left a great deal out — for instance: the updating of the Chebyshev–Pad´e table and Chebyshev rational approximation, Chebyshev approximation on small intervals, Faber polynomials on complex contours and Chebyshev (L∞) polynomials on complex domains. For the sake of those meeting this subject for the first time, we have included a number of problems at the end of each chapter. Some of these, in the earlier chapters in particular, are quite elementary; others are invitations to fill in the details of working that we have omitted simply for the sake of brevity; yet others are more advanced problems calling for substantial time and effort. We have dedicated this book to the memory of two recently deceased col- leagues and friends, who have influenced us in the writing of this book. Geoff Hayes wrote (with Charles Clenshaw) the major paper on fitting bivariate polynomials to data lying on a family of parallel lines. Their algorithm retains its place in numerical libraries some thirty-seven years later; it exploits the idea that Chebyshev polynomials form a well-conditioned basis independent
© 2003 by CRC Press LLC of the spacing of data. Lev Brutman specialised in near-minimax approxima- tions and related topics and played a significant role in the development of this field. In conclusion, there are many to whom we owe thanks, of whom we can mention only a few. Among colleagues who helped us in various ways in the writing of this book (but should not be held responsible for it), we must name Graham Elliott, Ezio Venturino, William Smith, David Elliott, Tim Phillips and Nick Trefethen; for getting the book started and keeping it on course, Bill Morton and Elizabeth Johnston in England, Bob Stern, Jamie Sigal and others at CRC Press in the United States; for help with preparing the manuscript, Pam Moore and Andrew Crampton. We must finally thank our wives, Moya and Elizabeth, for the blind faith in which they have encouraged us to bring this work to completion, without evidence that it was ever going to get there. This book was typeset at Oxford University Computing Laboratory, using Lamport’s LATEX2ε package.
John Mason David Handscomb April 2002
© 2003 by CRC Press LLC Contents
1Definitions 1.1Preliminaryremarks 1.2Trigonometricdefinitionsandrecurrences
1.2.1Thefirst-kindpolynomialTn
1.2.2Thesecond-kindpolynomialUn
1.2.3Thethird-andfourth-kindpolynomialsVnandWn (theairfoilpolynomials) 1.2.4Connectionsbetweenthefourkindsofpolynomial 1.3ShiftedChebyshevpolynomials ∗ ∗ ∗ ∗ 1.3.1TheshiftedpolynomialsTn,Un,Vn,Wn 1.3.2Chebyshevpolynomialsforthegeneralrange[a,b] 1.4Chebyshevpolynomialsofacomplexvariable 1.4.1Conformalmappingofacircletoandfromanellipse 1.4.2Chebyshevpolynomialsinz 1.4.3Shabatpolynomials 1.5ProblemsforChapter1
2BasicPropertiesandFormulae 2.1Introduction 2.2Chebyshevpolynomialzerosandextrema 2.3RelationsbetweenChebyshevpolynomialsandpowersofx
2.3.1Powersofxintermsof{Tn(x)}
2.3.2Tn(x)intermsofpowersofx
2.3.3RatiosofcoefficientsinTn(x) 2.4EvaluationofChebyshevsums,products,integralsandderiva- tives 2.4.1EvaluationofaChebyshevsum 2.4.2StabilityoftheevaluationofaChebyshevsum 2.4.3Evaluationofaproduct 2.4.4Evaluationofanintegral 2.4.5Evaluationofaderivative
© 2003 by CRC Press LLC 2.5ProblemsforChapter2
3TheMinimaxPropertyandItsApplications 3.1Approximation—theoryandstructure 3.1.1Theapproximationproblem 3.2Bestandminimaxapproximation 3.3TheminimaxpropertyoftheChebyshevpolynomials 3.3.1WeightedChebyshevpolynomialsofsecond,thirdand fourthkinds 3.4TheChebyshevsemi-iterativemethodforlinearequations 3.5Telescopingproceduresforpowerseries 3.5.1ShiftedChebyshevpolynomialson[0,1] 3.5.2Implementationofefficientalgorithms 3.6Thetaumethodforseriesandrationalfunctions 3.6.1Theextendedtaumethod 3.7ProblemsforChapter3
4OrthogonalityandLeast-SquaresApproximation 4.1Introduction—fromminimaxtoleastsquares 4.2OrthogonalityofChebyshevpolynomials 4.2.1Orthogonalpolynomialsandweightfunctions 4.2.2Chebyshevpolynomialsasorthogonalpolynomials
4.3OrthogonalpolynomialsandbestL2approximations 4.3.1Orthogonalpolynomialexpansions
4.3.2ConvergenceinL2oforthogonalexpansions 4.4Recurrencerelations 4.5Rodrigues’formulaeanddifferentialequations 4.6DiscreteorthogonalityofChebyshevpolynomials 4.6.1First-kindpolynomials 4.6.2Second-kindpolynomials 4.6.3Third-andfourth-kindpolynomials 4.7DiscreteChebyshevtransformsandthefastFouriertransform 4.7.1ThefastFouriertransform
© 2003 by CRC Press LLC 4.8Discretedatafittingbyorthogonalpolynomials:theForsythe– Clenshawmethod 4.8.1Bivariatediscretedatafittingonornearafamilyof linesorcurves 4.9Orthogonalityinthecomplexplane 4.10ProblemsforChapter4
5ChebyshevSeries 5.1Introduction—Chebyshevseriesandotherexpansions 5.2SomeexplicitChebyshevseriesexpansions 5.2.1Generatingfunctions 5.2.2Approximateseriesexpansions 5.3Fourier–ChebyshevseriesandFouriertheory
5.3.1L2-convergence 5.3.2Pointwiseanduniformconvergence 5.3.3BivariateandmultivariateChebyshevseriesexpansions 5.4Projectionsandnear-bestapproximations 5.5Near-minimaxapproximationbyaChebyshevseries
5.5.1Equalityofthenormtoλn 5.6ComparisonofChebyshevandotherorthogonalpolynomial expansions 5.7TheerrorofatruncatedChebyshevexpansion 5.8Seriesofsecond-,third-andfourth-kindpolynomials 5.8.1Seriesofsecond-kindpolynomials 5.8.2Seriesofthird-kindpolynomials 5.8.3MultivariateChebyshevseries 5.9LacunaryChebyshevseries 5.10Chebyshevseriesinthecomplexdomain 5.10.1Chebyshev–Pad´eapproximations 5.11ProblemsforChapter5
6ChebyshevInterpolation 6.1Polynomialinterpolation 6.2Orthogonalinterpolation 6.3Chebyshevinterpolationformulae
© 2003 by CRC Press LLC 6.3.1Aliasing 6.3.2Second-kindinterpolation 6.3.3Third-andfourth-kindinterpolation 6.3.4Conditioning
6.4BestL1approximationbyChebyshevinterpolation 6.5Near-minimaxapproximationbyChebyshevinterpolation 6.6ProblemsforChapter6
7Near-BesLt∞,L1andLpApproximations
7.1Near-bestL∞(near-minimax)approximations
7.1.1Second-kindexpansionsinL∞
7.1.2Third-kindexpansionsinL∞
7.2Near-bestL1approximations
7.3Bestandnear-bestLpapproximations 7.3.1Complexvariableresultsforelliptic-typeregions 7.4ProblemsforChapter7
8IntegrationUsingChebyshevPolynomials 8.1IndefiniteintegrationwithChebyshevseries 8.2Gauss–Chebyshevquadrature 8.3QuadraturemethodsofClenshaw–Curtistype 8.3.1Introduction 8.3.2First-kindformulae 8.3.3Second-kindformulae 8.3.4Third-kindformulae 8.3.5GeneralremarkonmethodsofClenshaw–Curtistype 8.4ErrorestimationforClenshaw–Curtismethods 8.4.1First-kindpolynomials 8.4.2Fittinganexponentialcurve 8.4.3Otherabscissaeandpolynomials 8.5SomeotherworkonClenshaw–Curtismethods 8.6ProblemsforChapter8
© 2003 by CRC Press LLC 9SolutionofIntegralEquations 9.1Introduction 9.2Fredholmequationsofthesecondkind 9.3Fredholmequationsofthethirdkind 9.4Fredholmequationsofthefirstkind 9.5Singularkernels 9.5.1Hilbert-typekernelsandrelatedkernels 9.5.2Symm’sintegralequation 9.6Regularisationofintegralequations 9.6.1Discretedatawithsecondderivativeregularisation 9.6.2Detailsofasmoothingalgorithm(secondderivative regularisation) 9.6.3Asmoothingalgorithmwithweightedfunctionregu- larisation 9.6.4EvaluationofV(λ) 9.6.5Otherbasisfunctions 9.7Partialdifferentialequationsandboundaryintegralequation methods 9.7.1Ahypersingularintegralequationderivedfromamixed boundaryvalueproblemforLaplace’sequation 9.8ProblemsforChapter9
10SolutionofOrdinaryDifferentialEquations 10.1Introduction 10.2Asimpleexample 10.2.1Collocationmethods 10.2.2Errorofthecollocationmethod 10.2.3Projection(tau)methods 10.2.4Erroroftheprecedingprojectionmethod 10.3TheoriginalLanczostau(τ)method 10.4Amoregenerallinearequation 10.4.1Collocationmethod 10.4.2Projectionmethod 10.5Pseudospectralmethods—anotherformofcollocation
© 2003 by CRC Press LLC 10.5.1Differentiationmatrices 10.5.2DifferentiationmatrixforChebyshevpoints 10.5.3Collocationusingdifferentiationmatrices 10.6Nonlinearequations 10.7Eigenvalueproblems 10.7.1Collocationmethods 10.7.2Collocationusingthedifferentiationmatrix 10.8Differentialequationsinonespaceandonetimedimension 10.8.1Collocationmethods 10.8.2Collocationusingthedifferentiationmatrix 10.9ProblemsforChapter10
11ChebyshevandSpectralMethodsforPartialDifferential Equations 11.1Introduction 11.2Interior,boundaryandmixedmethods 11.2.1Interiormethods 11.2.2Boundarymethods 11.2.3Mixedmethods 11.3Differentiationmatricesandnodalrepresentation 11.4Methodofweightedresiduals 11.4.1ContinuousMWR 11.4.2DiscreteMWR—anewnomenclature 11.5ChebyshevseriesandGalerkinmethods 11.6Collocation/interpolationandrelatedmethods 11.7PDEmethods 11.7.1Erroranalysis 11.8SomePDEproblemsandvariousmethods 11.8.1Powerbasis:collocationforPoissonproblem 11.8.2Powerbasis:interiorcollocationfortheL-membrane 11.8.3Chebyshevbasisanddiscreteorthogonalisation 11.8.4Differentiationmatrixapproach:Poissonproblem 11.8.5ExplicitcollocationforthequasilinearDirichletprob- lem:Chebyshevbasis
© 2003 by CRC Press LLC 11.9Computationalfluiddynamics 11.10Particularissuesinspectralmethods 11.11Moreadvancedproblems 11.12ProblemsforChapter11
12Conclusion
Bibliography
Appendices:
ABiographicalNote
BSummaryofNotations,DefinitionsandImportantProper- ties B.1Miscellaneousnotations B.2ThefourkindsofChebyshevpolynomial
CTablesofCoefficients
© 2003 by CRC Press LLC Chapter 1
Definitions
1.1 Preliminary remarks
“Chebyshev polynomials are everywhere dense in numerical anal- ysis.”
This remark has been attributed to a number of distinguished mathematicians and numerical analysts. It may be due to Philip Davis, was certainly spoken by George Forsythe, and it is an appealing and apt remark. There is scarcely any area of numerical analysis where Chebyshev polynomials do not drop in like surprise visitors, and indeed there are now a number of subjects in which these polynomials take a significant position in modern developments — including orthogonal polynomials, polynomial approximation, numerical integration, and spectral methods for partial differential equations. However, there is a different slant that one can give to the quotation above, namely that by studying Chebyshev polynomials one is taken on a journey which leads into all areas of numerical analysis. This has certainly been our personal experience, and it means that the Chebyshev polynomials, far from being an esoteric and narrow subject, provide the student with an opportunity for a broad and unifying introduction to many areas of numerical analysis and mathematics.
1.2 Trigonometric definitions and recurrences
There are several kinds of Chebyshev polynomials. In particular we shall in- troduce the first and second kind polynomials Tn(x)andUn(x), as well as a pair of related (Jacobi) polynomials Vn(x)andWn(x), which we call the ‘Chebyshev polynomials of the third and fourth kinds’; in addition we cover ∗ ∗ ∗ ∗ the shifted polynomials Tn (x), Un(x), Vn (x)andWn (x). We shall, however, only make a passing reference to ‘Chebyshev’s polynomial of a discrete vari- able’, referred to for example in Erd´elyi et al. (1953, Section 10.23), since this last polynomial has somewhat different properties from the polynomials on which our main discussion is based. Some books and many articles use the expression ‘Chebyshev polynomial’ to refer exclusively to the Chebyshev polynomial Tn(x) of the first kind. In- deed this is by far the most important of the Chebyshev polynomials and, when no other qualification is given, the reader should assume that we too are referring to this polynomial.
© 2003 by CRC Press LLC ClearlysomedefinitionofChebyshevpolynomialsisneededrightaway and,asweshallseeasthebookprogresses,wearespoiledforachoiceof definitions.However,whatgivesthevariouspolynomialstheirpowerand relevanceistheircloserelationshipwiththetrigonometricfunctions‘cosine’ and‘sine’.Weareallawareofthepowerofthesefunctionsandoftheir appearanceinthedescriptionofallkindsofnaturalphenomena,andthis mustsurelybethekeytotheversatilityoftheChebyshevpolynomials.We thereforeuseasourprimarydefinitionsthesetrigonometricrelationships.
1.2.1Thefirst-kindpolynomialTn
Definition1.1TheChebyshevpolynomialTn(x)ofthefirstkindisapoly- nomialinxofdegreen,definedbytherelation
Tn(x)=cosnθwhenx=cosθ.(1.1)
Iftherangeofthevariablexistheinterval[−1,1],thentherangeofthe correspondingvariableθcanbetakenas[0,π].Theserangesaretraversedin oppositedirections,sincex=−1correspondstoθ=πandx=1corresponds toθ=0. Itiswellknown(asaconsequenceofdeMoivre’sTheorem)thatcosnθ isapolynomialofdegreenincosθ,andindeedwearefamiliarwiththe elementaryformulae cos0θ=1,cos1θ=cosθ,cos2θ=2cos2θ−1, cos3θ=4cos3θ−3cosθ,cos4θ=8cos4θ−8cos2θ+1,....
Wemayimmediatelydeducefrom(1.1),thatthefirstfewChebyshev polynomialsare 2 T0(x)=1,T1(x)=x,T2(x)=2x −1, 3 4 2 (1.2) T3(x)=4x −3x,T4(x)=8x −8x +1,....
CoefficientsofallpolynomialsTn(x)uptodegreen=21willbefoundin TablesC.2a,C.2binAppendixC.
InpracticeitisneitherconvenientnorefficienttoworkouteachTn(x) from first principles. Rather by combining the trigonometric identity cos nθ +cos(n − 2)θ =2cosθ cos(n − 1)θ with Definition 1.1, we obtain the fundamental recurrence relation
Tn(x)=2xTn−1(x) − Tn−2(x),n=2, 3,..., (1.3a) which together with the initial conditions
T0(x)=1,T1(x)=x (1.3b)
© 2003 by CRC Press LLC recursivelygeneratesallthepolynomials{Tn(x)}veryefficiently. Itiseasytodeducefrom(1.3)thattheleadingcoefficient(thatofxn)in Tn(x)forn>1isdoubletheleadingcoefficientinTn−1(x)andhence,by induction,is2n−1.
Figure1.1:T5(x)onrange[−1,1] Figure1.2:cos5θonrange[0,π]
WhatdoesthepolynomialTn(x)looklike,andhowdoesagraphinthe variablexcomparewithagraphofcosnθinthevariableθ?InFigures1.1 and1.2weshowtherespectivegraphsofT5(x)andcos5θ. It will be noted that the shape of T5(x)on[−1, 1] is very similar to that of cos 5θ on [0,π], and in particular both oscillate between six extrema of equal magnitudes (unity) and alternating signs. However, there are three key differences — firstly the polynomial T5(x) corresponds to cos 5θ reversed (i.e., starting with a value of −1 and finishing with a value of +1); secondly the extrema of T5(x)attheend points x = ±1 do not correspond to zero gradients (as they do for cos 5θ) but rather to rapid changes in the polynomial as a function of x; and thirdly the zeros and extrema of T5(x) are clustered towards the end points ±1, whereas the zeros and extrema of cos 5θ are equally spaced. The reader will recall that an even function f(x) is one for which f(x)=f(−x) for all x and an odd function f(x) is one for which f(x)=−f(−x) for all x. All even powers of x are even functions, and all odd powers of x are odd functions. Equations (1.2) suggest that Tn(x) is an even or odd function, involving only even or odd powers of x, according as n is even or odd. This may be deduced rigorously from (1.3a) by induction, the cases n =0and n = 1 being supplied by the initial conditions (1.3b).
1.2.2 The second-kind polynomial Un
Definition 1.2 The Chebyshev polynomial Un(x) of the second kind is a poly- nomial of degree n in x defined by
Un(x)=sin(n +1)θ/ sin θ when x =cosθ. (1.4)
© 2003 by CRC Press LLC TherangesofxandθarethesameasforTn(x). Elementaryformulaegive sin1θ=sinθ,sin2θ=2sinθcosθ,sin3θ=sinθ(4cos2θ−1), sin4θ=sinθ(8cos3θ−4cosθ),..., sothatweseethattheratioofsinefunctions(1.4)isindeedapolynomialin cosθ,andwemayimmediatelydeducethat 2 U0(x)=1,U1(x)=2x,U2(x)=4x −1, 3 (1.5) U3(x)=8x −4x,....
CoefficientsofallpolynomialsUn(x)uptodegreen=21willbefoundin TablesC.3a,C.3binAppendixC. Bycombiningthetrigonometricidentity sin(n+1)θ+sin(n−1)θ=2cosθsinnθ
withDefinition1.2,wefindthatUn(x)satisfiestherecurrencerelation
Un(x)=2xUn−1(x)−Un−2(x),n=2,3,...,(1.6a) whichtogetherwiththeinitialconditions
U0(x)=1,U1(x)=2x(1.6b) providesanefficientprocedureforgeneratingthepolynomials. Asimilartrigonometricidentity sin(n+1)θ−sin(n−1)θ=2sinθcosnθ leadsustoarelationship
Un(x)−Un−2(x)=2Tn(x),n=2,3,...,(1.7) betweenthepolynomialsofthefirstandsecondkinds. n Itiseasytodeducefrom(1.6)thattheleadingcoefficientofx inUn(x) is2n.
Notethattherecurrence(1.6a)for{Un(x)}isidenticalinformtothe recurrence(1.3a)for{Tn(x)}.Thedifferentinitialconditions[(1.6b)and (1.3b)]yieldthedifferentpolynomialsystems.
InFigure1.3weshowthegraphofU5(x).Itoscillatesbetweensixex- trema,asdoesT5(x)inFigure1.1,butinthepresentcasetheextremahave magnitudes which are not equal, but increase monotonically from the centre towards the ends of the range.
From (1.5) it is clear that the second-kind polynomial Un(x), like the first, is an even or odd function, involving only even or odd powers of x, according as n is even or odd.
© 2003 by CRC Press LLC Figure1.3:U5(x) on range [−1, 1]
1.2.3 The third- and fourth-kind polynomials Vn and Wn (the air- foil polynomials)
Two other families of polynomials Vn and Wn may be constructed, which are related to Tn and Un, but which have trigonometric definitions involving the half angle θ/2(wherex =cosθ as before). These polynomials are sometimes1 referred to as the ‘airfoil polynomials’, but Gautschi (1992) rather appropri- ately named them the ‘third- and fourth-kind Chebyshev polynomials’. First we define these polynomials trigonometrically, by a pair of relations parallel to (1.1) and (1.4) above for Tn and Un. Again the ranges of x and θ are the same as for Tn(x).
Definition 1.3 The Chebyshev polynomials Vn(x) and Wn(x) of the third and fourth kinds are polynomials of degree n in x defined respectively by
1 1 Vn(x)=cos(n + 2 )θ/ cos 2 θ (1.8) and 1 1 Wn(x)=sin(n + 2 )θ/ sin 2 θ, (1.9) when x =cosθ.
1 To justify these definitions, we first observe that cos(n + 2 )θ is an odd 1 polynomial of degree 2n+1 incos 2 θ. Therefore the right-hand side of (1.8) is 1 an even polynomial of degree 2n in cos 2 θ, which is equivalent to a polynomial 2 1 1 of degree n in cos 2 θ = 2 (1 + cos θ) and hence to a polynomial of degree n in cos θ.ThusVn(x) is indeed a polynomial of degree n in x. For example
1 3 1 1 cos(1 + 2 )θ 4cos 2 θ − 3cos 2 θ 2 1 V1(x)= 1 = 1 =4cos 2 θ−3=2cosθ−1=2x−1. cos 2 θ cos 2 θ
We may readily show that
2 V0(x)=1,V1(x)=2x − 1,V2(x)=4x − 2x − 1, 3 2 (1.10) V3(x)=8x − 4x − 4x +1, .... 1See for example, Fromme & Golberg (1981).
© 2003 by CRC Press LLC 1 1 Similarlysin(n+2)θisanoddpolynomialofdegree2n+1insin2θ. Thereforetheright-handsideof(1.9)isanevenpolynomialofdegree2nin 1 21 1 sin2θ,whichisequivalenttoapolynomialofdegreeninsin 2θ=2(1−cosθ) andhenceagaintoapolynomialofdegreenincosθ.Forexample
1 1 31 sin(1+2)θ 3sin2θ−4sin 2θ 21 W1(x)= 1 = 1 =3−4sin 2θ=2cosθ+1=2x+1. sin2θ sin2θ
Wemayreadilyshowthat
2 W0(x)=1,W1(x)=2x+1,W2(x)=4x +2x−1, 3 2 (1.11) W3(x)=8x +4x −4x−1,....
ThepolynomialsVn(x)andWn(x)are,infact,rescalingsoftwoparticular 2 (α,β) 1 1 Jacobi polynomialsPn(x)withα=−2,β=2 andviceversa.Explicitly n 1 1 n 1 1 2 2n (−2,2) 2 2n (2,−2) Vn(x)=2 Pn(x), Wn(x)=2 Pn(x). n n
CoefficientsofallpolynomialsVn(x)andWn(x)uptodegreen=10will befoundinTableC.1inAppendixC. These polynomials too may be efficiently generated by the use of a recur- rence relation. Since
1 1 1 cos(n + 2 )θ +cos(n − 2+ 2 )θ =2cosθ cos(n − 1+ 2 )θ
and 1 1 1 sin(n + 2 )θ +sin(n − 2+ 2 )θ =2cosθ sin(n − 1+ 2 )θ, it immediately follows that
Vn(x)=2xVn−1(x) − Vn−2(x),n=2, 3,..., (1.12a)
and Wn(x)=2xWn−1(x) − Wn−2(x),n=2, 3,..., (1.12b) with V0(x)=1,V1(x)=2x − 1 (1.12c) and W0(x)=1,W1(x)=2x +1. (1.12d)
Thus Vn(x)andWn(x) share precisely the same recurrence relation as Tn(x)andUn(x), and their generation differs only in the prescription of the initial condition for n =1. 2See Chapter 22 of Abramowitz and Stegun’s Handbook of Mathematical Functions (1964).
© 2003 by CRC Press LLC Itisimmediatelyclearfrom(1.12)thatbothVn(x)andWn(x)arepoly- nomialsofdegreeninx,inwhichallpowersofxarepresent,andinwhich theleadingcoefficients(ofxn)areequalto2n.
InFigure1.4weshowgraphsofV5(x)andW5(x). They are exact inverted mirror images of one another, as will be proved in the next section (1.19).
Figure 1.4: V5(x)andW5(x) on range [−1, 1]
1.2.4 Connections between the four kinds of polynomial
We already have a relationship (1.7) between the polynomials Tn and Un.It remains to link Vn and Wn to Tn and Un. This may be done by introducing two auxiliary variables
1 1 1 2 1 1 2 1 u =[2 (1 + x)] =cos2 θ, t =[2 (1 − x)] =sin2 θ. (1.13)
Using (1.8) and (1.9) it immediately follows, from the definitions (1.1) and (1.4) of Tn and Un,that
1 −1 Tn(x)=T2n(u),Un(x)= 2 u U2n+1(u), (1.14) −1 Vn(x)=u T2n+1(u),Wn(x)=U2n(u). (1.15)
Thus Tn(x), Un(x), Vn(x), Wn(x) together form the first- and second-kind polynomials in u,weightedbyu−1 in the case of odd degrees. Also (1.15) shows that Vn(x)andWn(x) are directly related, respectively, to the first- and second-kind Chebyshev polynomials, so that the terminology of ‘Chebyshev polynomials of the third and fourth kind’ is justifiable. From the discussion above it can be seen that, if we wish to establish properties of Vn and Wn, then we have two main options: we can start from the trigonometric definitions (1.8), (1.9) or we can attempt to exploit properties of Tn and Un by using the links (1.14)–(1.15).
Note that Vn and Wn are neither even nor odd (unlike Tn and Un). We n n have seen that the leading coefficient of x is 2 in both Vn and Wn,asitis in Un. This suggests a close link with Un. Indeed if we average the initial conditions (1.12c) and (1.12d) for V1 and W1 we obtain the initial condition
© 2003 by CRC Press LLC (1.6b) for U1, from which we can show that the average of Vn and Wn satisfies the recurrence (1.6a) subject to (1.6b) and therefore that for all n 1 Un(x)= 2 [Vn(x)+Wn(x)]. (1.16)
The last result also follows directly from the trigonometric definitions (1.4), (1.8), (1.9) of Un, Vn, Wn,since sin(n +1)θ sin(n + 1 )θ cos 1 θ +cos(n + 1 )θ sin 1 θ = 2 2 2 2 θ 1 1 sin 2sin 2 θ cos 2 θ 1 1 1 cos(n + 2 )θ sin(n + 2 )θ = 2 1 + 1 . cos 2 θ sin 2 θ
Equation (1.16) is not the only link between the sets {Vn}, {Wn} and {Un}. Indeed, by using the trigonometric relations 1 1 2sin 2 θ cos(n + 2 )θ =sin(n +1)θ − sin nθ, 1 1 2cos 2 θ sin(n + 2 )θ =sin(n +1)θ +sinnθ and dividing through by sin θ, we can deduce that
Vn(x)=Un(x) − Un−1(x), (1.17)
Wn(x)=Un(x)+Un−1(x). (1.18)
Thus Vn and Wn may be very simply determined once {Un} are available. Note that (1.17), (1.18) are confirmed in the formulae (1.5), (1.10), (1.11) and are consistent with (1.16) above.
From the evenness/oddness of Un(x)forn even/odd, we may immediately deduce from (1.17), (1.18) that
Wn(x)=Vn(−x)(n even); (1.19) Wn(x)=−Vn(−x)(n odd). This means that the third- and fourth-kind polynomials essentially transform into each other if the range [−1, 1] of x is reversed, and it is therefore sufficient for us to study only one of these kinds of polynomial. Two further relationships that may be derived from the definitions are
Vn(x)+Vn−1(x)=Wn(x) − Wn−1(x)=2Tn(x). (1.20)
If we were asked for a ‘pecking order’ of these four Chebyshev polynomials Tn, Un, Vn and Wn, then we would say that Tn is clearly the most important and versatile. Moreover Tn generally leads to the simplest formulae, whereas results for the other polynomials may involve slight complications. However, all four polynomials have their role. For example, as we shall see, Un is useful in numerical integration, while Vn and Wn canbeusefulinsituationsinwhich singularities occur at one end point (+1 or −1) but not at the other.
© 2003 by CRC Press LLC 1.3ShiftedChebyshevpolynomials
∗ ∗ ∗ ∗ 1.3.1TheshiftedpolynomialsTn,Un,Vn,Wn
Sincetherange[0,1]isquiteoftenmoreconvenienttousethantherange [−1,1],wesometimesmaptheindependentvariablexin[0,1]tothevariable sin[−1,1]bythetransformation
1 s=2x−1orx=2(1+s), ∗ andthisleadstoashiftedChebyshevpolynomial(ofthefirstkind)Tn(x)of degreeninxon[0,1]givenby
∗ Tn(x)=Tn(s)=Tn(2x−1).(1.21)
Thuswehavethepolynomials
∗ ∗ ∗ 2 T0(x)=1,T1(x)=2x−1,T2(x)=8x −8x+1, ∗ 3 2 (1.22) T3(x)=32x −48x +18x−1,....
∗ From(1.21)and(1.3a),wemaydeducetherecurrencerelationforTn in theform ∗ ∗ ∗ Tn(x)=2(2x−1)Tn−1(x)−Tn−2(x)(1.23a) withinitialconditions
∗ ∗ T0(x)=1,T1(x)=2x−1.(1.23b)
∗ ThepolynomialsTn(x)haveafurtherspecialproperty,whichderivesfrom (1.1)and(1.21):
2 ∗ 2 T2n(x)=cos2nθ=cosn(2θ)=Tn(cos2θ)=Tn(2x −1)=Tn(x ) sothat ∗ 2 T2n(x)=Tn(x ).(1.24) Thispropertymayreadilybeconfirmedforthefirstfewpolynomialsbycom- ∗ √ paringtheformulae(1.2)and(1.22).ThusTn(x)ispreciselyT2n( x),a higherdegreeChebyshevpolynomialinthesquarerootoftheargument,and ∗ relation(1.24)givesanimportantlinkbetween{Tn}and{Tn}whichcom- plementstheshiftrelationship(1.21).Becauseofthisproperty,TableC.2a in Appendix C, which gives coefficients of the polynomials Tn(x)uptode- gree n = 20 for even n, at the same time gives coefficients of the shifted ∗ polynomials Tn (x)uptodegreen = 10. ∗ It is of course possible to define Tn , like Tn and Un, directly by a trigono- metric relation. Indeed, if we combine (1.1) and (1.24) we obtain
∗ 2 Tn (x)=cos2nθ when x =cos θ. (1.25)
© 2003 by CRC Press LLC Thisrelationmightalternativelyberewritten,withθreplacedbyφ/2,inthe form ∗ 2 1 Tn(x)=cosnφwhenx=cos φ/2=2(1+cosφ).(1.26) Indeedthelatterformulacouldbeobtaineddirectlyfrom(1.21),bywriting
Tn(s)=cosnφwhens=cosφ.
∗ NotethattheshiftedChebyshevpolynomialTn(x)isneitherevennorodd, 0 n ∗ andindeedallpowersofxfrom1=x tox appearinTn(x).Theleading n ∗ coefficientofx inTn(x)forn>0maybededucedfrom(1.23a),(1.23b)to be22n−1. ∗ ∗ ∗ ShiftedpolynomialsUn,Vn,Wnofthesecond,thirdandfourthkindsmay bedefinedinpreciselyanalogousways:
∗ ∗ ∗ Un(x)=Un(2x−1),Vn(x)=Vn(2x−1),Wn(x)=Wn(2x−1).(1.27)
∗ ∗ ∗ LinksbetweenUn,Vn,Wnandtheunstarredpolynomials,analogousto(1.24) above,mayreadilybeestablished.Forexample,using(1.4)and(1.27),
sinθU2n−1(x)=sin2nθ=sinn(2θ)=sin2θUn−1(cos2θ) 2 ∗ 2 =2sinθcosθUn−1(2x −1)=sinθ{2xUn−1(x )}
andhence ∗ 2 U2n−1(x)=2xUn−1(x ).(1.28)
∗ ∗ ThecorrespondingrelationsforVn andWn areslightlydifferentinthat theycomplement(1.24)and(1.28)byinvolvingT2n−1andU2n.Firstly,using (1.13),(1.15)and(1.27),
∗ 2 2 −1 Vn−1(u )=Vn−1(2u −1)=Vn−1(x)=u T2n−1(u)
andhence(replacingubyx)
∗ 2 T2n−1(x)=xVn−1(x ).(1.29)
Similarly,
∗ 2 2 Wn−1(u )=Wn−1(2u −1)=Wn−1(x)=U2n(u)
andhence(replacingubyx)
∗ 2 U2n(x)=Wn(x ).(1.30)
Becauseoftherelationships(1.28)–(1.30),TablesC.3b,C.2b,C.3ainAp- pendix C, which give coefficients of Tn(x)andUn(x)uptodegreen = 20, at ∗ ∗ thesametimegivethecoefficientsoftheshiftedpolynomialsUn(x), Vn (x), ∗ Wn (x), respectively, up to degree n = 10.
© 2003 by CRC Press LLC 1.3.2 Chebyshev polynomials for the general range [a, b]
In the last section, the range [−1, 1] was adjusted to the range [0, 1] for conve- nience, and this corresponded to the use of the shifted Chebyshev polynomials ∗ ∗ ∗ ∗ Tn , Un, Vn , Wn in place of Tn, Un, Vn, Wn respectively. More generally we may define Chebyshev polynomials appropriate to any given finite range [a, b] of x, by making this range correspond to the range [−1, 1] of a new variable s under the linear transformation 2x − (a + b) s = . (1.31) b − a
The Chebyshev polynomials of the first kind appropriate to [a, b]arethus Tn(s), where s is given by (1.31), and similarly the second-, third- and fourth- kind polynomials appropriate to [a, b]areUn(s), Vn(s), and Wn(s).
Example 1.1: The first-kind Chebyshev polynomial of degree three appropriate to the range [1, 4] of x is
3 2x − 5 2x − 5 2x − 5 1 3 2 T3 =4 − 3 = (32x − 240x + 546x − 365). 3 3 3 27
Note that in the special case [a, b] ≡ [0, 1], the transformation (1.31) be- comes s =2x−1, and we obtain the shifted Chebyshev polynomials discussed in Section 1.3.1.
Incidentally, the ‘Chebyshev Polynomials Sn(x)andCn(x)’ tabulated by the National Bureau of Standards (NBS 1952) are no more than mappings of Un and 2Tn to the range [a, b] ≡ [−2, 2]. Except for C0, these polynomials all have unit leading coefficient, but this appears to be their only recommending feature for practical purposes, and they have never caught on.
1.4 Chebyshev polynomials of a complex variable
We have chosen to define the polynomials Tn(x), Un(x), Vn(x)andWn(x) with reference to the interval [−1, 1]. However, their expressions as sums of powers of x can of course be evaluated for any real x, even though the substitution x =cosθ is not possible outside this interval. For x in the range [1, ∞), we can make the alternative substitution
x =coshΘ, (1.32)
© 2003 by CRC Press LLC with Θ in the range [0, ∞), and it is easily verified that precisely the same polynomials (1.2), (1.5), (1.10) and (1.11) are generated by the relations
Tn(x)=coshnΘ, (1.33a) sinh(n +1)Θ Un(x)= , (1.33b) sinh Θ 1 cosh(n + 2 )Θ Vn(x)= 1 , (1.33c) cosh 2 Θ 1 sinh(n + 2 )Θ Wn(x)= 1 . (1.33d) sinh 2 Θ
For x in the range (−∞, −1] we can make use of the odd or even parity of the Chebyshev polynomials to deduce from (1.33) that, for instance,
n Tn(x)=(−1) cosh nΘ
where x = − cosh Θ.
It is easily shown from (1.33) that none of the four kinds of Chebyshev polynomials can have any zeros or turning points in the range [1, ∞). The same applies to the range (−∞, −1]. This will later become obvious, since we shall show in Section 2.2 that Tn, Un, Vn and Wn each have n real zeros in the interval [−1, 1], and a polynomial of degree n can have at most n zeros in all.
The Chebyshev polynomial Tn(x) can be further extended into (or initially defined as) a polynomial Tn(z) of a complex variable z. Indeed Snyder (1966) and Trefethen (2000) both start from a complex variable in developing their expositions.
1.4.1 Conformal mapping of a circle to and from an ellipse
For convenience, we consider not only the variable z but a related complex variable w such that 1 −1 z = 2 (w + w ). (1.34)
Then, if w moves on the circle |w| = r (for r>1) centred at the origin, we have w = reiθ = r cos θ +ir sin θ, (1.35) w−1 = r−1e−iθ = r−1 cos θ − ir−1 sin θ, (1.36) and so, from (1.34), z = a cos θ +ib sin θ (1.37)
© 2003 by CRC Press LLC where 1 −1 1 −1 a=2(r+r ),b=2(r−r ).(1.38) Hencezmovesonthestandardellipse x2 y2 + =1(1.39) a2 b2 centredattheorigin,withmajorandminorsemi-axesa,bgivenby(1.38).It iseasytoverifyfrom(1.38)thattheeccentricityeofthisellipseissuchthat ae= a2−b2=1, andhencetheellipsehasfociatz=±1. Inthecaser=1,wherewmovesontheunitcircle,wehaveb=0and theellipsecollapsesintotherealinterval[−1,1].However,ztraversesthe intervaltwiceaswmovesroundthecircle:from−1to1asθmovesfrom−π to0,andfrom1to−1asθmovesfrom0toπ.
Figure1.5:Thecircle|w|=r=1.5anditsimageinthezplane
Figure1.6:Thecircle|w|=1anditsimageinthezplane
Thestandardcircle(1.35)andellipse(1.39)areshowninFigure1.5,and thespecialcaser=1isshowninFigure1.6.SeeHenrici(1974–1986)for further discussions of this mapping.
© 2003 by CRC Press LLC From (1.34) we readily deduce that w satisfies
w2 − 2wz +1=0, (1.40)
a quadratic equation with two solutions 2 w = w1,w2 = z ± z − 1. (1.41)
This means that the mapping from w to z is2to1,withbranchpointsat√ z = ±1. It is convenient to define the complex square root z2 − 1sothatit lies in the same quadrant as z (except for z on the real interval [−1, 1], along which the plane must be cut), and to choose the solution 2 w = w1 = z + z − 1, (1.42)
so that |w| = |w1|≥1. Then w depends continuously on z along any path that does not intersect the interval [−1, 1], and it is easy to verify that −1 2 w2 = w1 = z − z − 1, (1.43)
with |w2|≤1. −1 If w1 moves on |w1| = r,forr>1, then w2 moves on |w2| = w1 = r−1 < 1. Hence both of the concentric circles −1 Cr := {w : |w| = r} ,C1/r := w : |w| = r
transform into the same ellipse defined by (1.37) or (1.39), namely 2 Er := z : z + z − 1 = r . (1.44)
1.4.2 Chebyshev polynomials in z
Defining z by (1.34), we note that if w lies on the unit circle |w| = 1 (i.e. C1), then (1.37) gives z =cosθ (1.45) and hence, from (1.42). w = z + z2 − 1=eiθ. (1.46)
Thus Tn(z) is now a Chebyshev polynomial in a real variable and so by our standard definition (1.1), and (1.45), (1.46),
1 inθ −inθ 1 n −n Tn(z)=cosnθ = 2 (e +e )= 2 (w + w ).
This leads us immediately to our general definition, for all complex z, namely 1 n −n Tn(z)= 2 (w + w ) (1.47)
© 2003 by CRC Press LLC where 1 −1 z=2(w+w ).(1.48)
AlternativelywemaywriteTn(z)directlyintermsofz,using(1.42)and (1.43),as 1 2 n 2 n Tn(z)=2{(z+ z −1) +(z− z −1) }.(1.49)
IfzliesontheellipseEr,thelocusof(1.48)when|w|=r>1,thenit followsfrom(1.47)thatwehavetheinequality
1 n −n 1 n −n 2(r −r )≤|Tn(z)|≤2(r +r ).(1.50)
InFig.1.7weshowthelevelcurvesoftheabsolutevalueofT5(z),anditcan easilybeseenhowtheseapproachanellipticalshapeasthevalueincreases.
Figure1.7:Contoursof|T5(z)| in the complex plane
We may similarly extend polynomials of the second kind. If |w| =1,so that z =cosθ, we have from (1.4), sin nθ Un−1(z)= . sin θ Hence, from (1.45) and (1.46), we deduce the general definition wn − w−n Un−1(z)= (1.51) w − w−1
1 −1 where again z = 2 (w + w ). Alternatively, writing directly in terms of z, √ √ 2 n 2 n 1 (z + z − 1) − (z − z − 1) Un−1(z)= √ . (1.52) 2 z2 − 1
If z lies on the ellipse (1.44), then it follows directly from (1.51) that
rn − r−n rn + r−n ≤|Un−1(z)|≤ ; (1.53) r + r−1 r − r−1
© 2003 by CRC Press LLC however,whereasthebounds(1.50)on|Tn(z)|areattainedontheellipse,the bounds(1.53)on|Un−1(z)|areslightlypessimistic.Forasharpupperbound, wemayexpand(1.51)into n−1 n−3 3−n 1−n Un−1(z)=w +w +···+w +w (1.54) givingus n−1 n−3 3−n 1−n |Un−1(z)|≤ w + w +···+ w + w =rn−1+rn−3+···+r3−n+r1−n rn−r−n = ,(1.55) r−r−1 whichliesbetweenthetwoboundsgivenin(1.53).InFig.1.8weshowthe level curves of the absolute value of U5(z).
Figure 1.8: Contours of |U5(z)| in the complex plane
The third- and fourth-kind polynomials of degree n in z may readily be defined in similar fashion (compare (1.51)) by
n+ 1 −n− 1 w 2 w 2 V z + , n( )= 1 − 1 (1.56) w 2 + w 2 n+ 1 −n− 1 w 2 − w 2 W z n( )= 1 − 1 (1.57) w 2 − w 2
1 where w 2 is consistently defined from w. More precisely, to get round the ambiguities inherent in taking square roots, we may define them by wn+1 + w−n Vn(z)= , (1.58) w +1 wn+1 − w−n Wn(z)= (1.59) w − 1 It is easily shown, by dividing denominators into numerators, that these give 1 −1 polynomials of degree n in z = 2 (w + w ).
© 2003 by CRC Press LLC 1.4.3 Shabat polynomials
Shabat & Voevodskii (1990) introduced the concept of ‘generalised Chebyshev polynomials’ (or Shabat polynomials), in the context of trees and number theory. The most recent survey paper in this area is that of Shabat & Zvonkin (1994). They are defined as polynomials P (z) with complex coefficients having two critical values A and B such that
P (z)=0=⇒ P (z)=A or P (z)=B.
The prime example of such a polynomial is Tn(z), a first-kind Chebyshev polynomial, for which A = −1andB = +1 are the critical values.
1.5 Problems for Chapter 1
1. The equation x =cosθ defines infinitely many values of θ corresponding to a given value of x in the range [−1, 1]. Show that, whichever value is chosen, the values of Tn(x), Un(x), Vn(x)andWn(x) as defined by (1.1), (1.4), (1.8) and (1.9) remain the same.
2. Determine explicitly the Chebyshev polynomials of first and second kinds of degrees 0, 1, 2, 3, 4 appropriate to the range [−4, 6] of x.
3. Prove that Tm(Tn(x)) = Tmn(x) and that
Um−1(Tn(x))Un−1(x)=Un−1(Tm(x))Um−1(x)=Umn−1(x).
4. Verify that equations (1.33) yield the same polynomials for x>1asthe trigonometric definitions of the Chebyshev polynomials give for |x|≤1.
5. Using the formula
1 −1 1 −1 z = 2 (r + r )cosθ + 2 i(r − r )sinθ, (r>1)
which defines a point on an ellipse centred at 0 with foci z = ±1,
(a) verify that 2 1 −1 1 −1 z − 1= 2 (r − r )cosθ + 2 i(r + r )sinθ
and hence √ (b) verify that z + z2 − 1 = r.
© 2003 by CRC Press LLC 6. By expanding by the first row and using the standard three-term recur- rence for Tr(x), show that 2x −100··· 000 −12x −10··· 000 0 −12x −1 ··· 000 Tn(x)= . ...... ...... 0000··· −12x −1 ··· − x 0000 0 1 (n × n determinant)
Write down similar expressions for Un(x), Vn(x)andWn(x). 7. Given that the four kinds of Chebyshev polynomial each satisfy the same recurrence relation Xn =2xXn−1 − Xn−2,
with X0 = 1 in each case and X1 = x, 2x, 2x +1, 2x − 1 for the four respective families, use these relations only to establish that
(a) Vi(x)+Wi(x)=2Ui(x),
(b) Vi(x) − Wi(x)=2Ui−1(x),
(c) Ui(x) − 2Ti(x)=Ui−2(x),
(d) Ui(x) − Ti(x)=xUi−1(x). 8. Derive the same four formulae of Problem 7, this time using only the trigonometric definitions of the Chebyshev polynomials. 9. From the last two results in Problem 7, show that
(a) Ti(x)=xUi−1(x) − Ui−2(x),
(b) Ui(x)=2xUi−1(x) − Ui−2(x).
© 2003 by CRC Press LLC Chapter 2
Basic Properties and Formulae
2.1 Introduction
The aim of this chapter is to provide some elementary formulae for the manip- ulation of Chebyshev polynomials and to summarise the key properties which will be developed in the book. Areas of application will be introduced and discussed in the chapters devoted to them.
2.2 Chebyshev polynomial zeros and extrema
The Chebyshev polynomials of degree n>0 of all four kinds have precisely n zeros and n + 1 local extrema in the interval [−1, 1]. In the case n =5,this is evident in Figures 1.1, 1.3 and 1.4. Note that n − 1 of these extrema are interior to [−1, 1], and are ‘true’ alternate maxima and minima (in the sense that the gradient vanishes), the other two extrema being at the end points ±1 (where the gradient is non-zero).
From formula (1.1), the zeros for x in [−1, 1] of Tn(x)mustcorrespondto the zeros for θ in [0,π]ofcosnθ,sothat
1 nθ =(k − 2 )π, (k =1, 2,...,n).
Hence, the zeros of Tn(x)are
1 (k − 2 )π x = xk =cos , (k =1, 2,...,n). (2.1) n
Example 2.1: For n = 3, the zeros are √ √ π 3 3π 5π 3 x = x1 =cos = ,x2 =cos =0,x3 =cos = − . 6 2 6 6 2
Note that these zeros are in decreasing order in x (corresponding to in- creasing θ), and it is sometimes preferable to list them in their natural order as (n − k + 1 )π x =cos 2 , (k =1, 2,...,n). (2.2) n
© 2003 by CRC Press LLC Note, too, that x = 0 is a zero of Tn(x) for all odd n, but not for even n,and that zeros are symmetrically placed in pairs on either side of x =0.
The zeros of Un(x) (defined by (1.4)) are readily determined in a similar way from the zeros of sin(n +1)θ as kπ x = yk =cos , (k =1, 2,...,n) (2.3) (n +1) or in their natural order (n − k +1)π x =cos , (k =1, 2,...,n). (2.4) n +1
One is naturally tempted to extend the set of points (2.3) by including the further values y0 =1andyn+1 = −1, giving the set kπ x = yk =cos , (k =0, 1,...,n+1). (2.5) (n +1)
These are zeros not of Un(x), but of the polynomial 2 (1 − x )Un(x). (2.6) However, we shall see that these points are popular as nodes in applications to integration.
The zeros of Vn(x)andWn(x) (defined by (1.8), (1.9)) correspond to zeros 1 1 of cos(n + 2 )θ and sin(n + 2 )θ, respectively. Hence, the zeros of Vn(x) occur at 1 (k − 2 )π x =cos 1 , (k =1, 2,...,n) (2.7) n + 2 or in their natural order 1 (n − k + 2 )π x =cos 1 , (k =1, 2,...,n), (2.8) n + 2
while the zeros of Wn(x) occur at kπ x =cos 1 , (k =1, 2,...,n) (2.9) n + 2 or in their natural order (n − k +1)π x =cos 1 , (k =1, 2,...,n). (2.10) n + 2
Note that there are natural extensions of these point sets, by including the value k = n + 1 and hence x = −1 in (2.7) and the value k = 0 and hence x = 1 in (2.9). Thus the polynomials
(1 + x)Vn(x)and(1− x)Wn(x)
© 2003 by CRC Press LLC haveaszerostheirnaturalsets(2.7)fork=1,...,n+1and(2.9)fork= 0,1,...,n,respectively.
TheinternalextremaofTn(x)correspondtotheextremaofcosnθ,namely thezerosofsinnθ,since d d d dx −nsinnθ Tn(x)= cosnθ= cosnθ = . dx dx dθ dθ −sinθ
Hence,includingthoseatx=±1,theextremaofTn(x)on[−1,1]are kπ x=cos ,(k=0,1,...,n)(2.11) n orintheirnaturalorder (n−k)π x=cos ,(k=0,1,...,n).(2.12) n
2 Thesearepreciselythezerosof(1−x )Un−1(x),namelythepoints(2.5) above(withnreplacedbyn−1).Notethattheextremaareallofequal magnitude(unity)andalternateinsignatthepoints(2.12)between−1and +1,asshowninFigure1.1.
The extrema of Un(x), Vn(x), Wn(x) are not in general as readily deter- mined; indeed finding them involves the solution of transcendental equations. For example, d d sin(n +1)θ −(n +1)sinθ cos(n +1)θ +cosθ sin(n +1)θ Un(x)= = dx dx sin θ sin3 θ and the extrema therefore correspond to values of θ satisfying the equation tan(n +1)θ =(n +1)tanθ =0.
All that we can say for certain is that the extreme values of Un(x)have magnitudes which increase monotonically as |x| increases away from 0, until the largest magnitude of n + 1 is achieved at x = ±1. On the other hand, from the definitions (1.4), (1.8), (1.9), we can show that 2 1 − x Un(x)=sin(n +1)θ, √ √ 1 1+xVn(x)= 2cos(n + 2 )θ, √ √ 1 1 − xWn(x)= 2sin(n + 2 )θ; √ √ 2 Hence√ the extrema of the weighted polynomials 1 − x Un(x), 1+xVn(x), 1 − xWn(x) are explicitly determined and occur, respectively, at (2k +1)π 2kπ (2k +1)π x =cos ,x=cos ,x=cos (k =0, 1,...,n). 2(n +1) 2n +1 2n +1
© 2003 by CRC Press LLC 2.3 Relations between Chebyshev polynomials and powers of x
It is useful and convenient in various applications to be able to express Cheb- yshev polynomials explicitly in terms of powers of x, and vice versa. Such formulae are simplest and easiest to derive in the case of the first kind poly- nomials Tn(x), and so we concentrate on these.
2.3.1 Powers of x in terms of {Tn(x)}
The power xn can be expressed in terms of the Chebyshev polynomials of degrees up to n, but, since these are alternately even and odd, we see at once that we need only include polynomials of alternate degrees, namely Tn(x), n Tn−2(x), Tn−4(x), . . . . Writing x =cosθ, we therefore need to express cos θ in terms of cos nθ,cos(n − 2)θ,cos(n − 4)θ, ..., andthisisreadilyachieved by using the binomial theorem as follows: n n (eiθ +e−iθ)n =einθ + ei(n−2)θ + ···+ e−i(n−2)θ +e−inθ 1 n − 1 n =(einθ +e−inθ)+ (ei(n−2)θ +e−i(n−2)θ)+ 1 n + (ei(n−4)θ +e−i(n−4)θ)+···. (2.13) 2 Here we have paired in brackets the first and last terms, the second and second-to-last terms, and so on. The number of such brackets will be
n/2 +1
where m denotes the integer part of m.Whenn is even, the last bracket in (2.13) will contain only the one (middle) term e0θ [= 1]. Now using the fact that
(eiθ +e−iθ)n =(2cosθ)n =2n cosn θ
we deduce from (2.13) that