Radius of convergence




In mathematics, the radius of convergence of a power series is the radius of the largest disk in which the series converges. It is either a non-negative real number or {displaystyle infty }infty . When it is positive, the power series converges absolutely and uniformly on compact sets inside the open disk of radius equal to the radius of convergence, and it is the Taylor series of the analytic function to which it converges.




Contents






  • 1 Definition


  • 2 Finding the radius of convergence


    • 2.1 Theoretical radius


    • 2.2 Practical estimation of radius in the case of real coefficients




  • 3 Radius of convergence in complex analysis


    • 3.1 A simple example


    • 3.2 A more complicated example




  • 4 Convergence on the boundary


  • 5 Rate of convergence


  • 6 Abscissa of convergence of a Dirichlet series


  • 7 Notes


  • 8 References


  • 9 External links





Definition


For a power series ƒ defined as:


f(z)=∑n=0∞cn(z−a)n,{displaystyle f(z)=sum _{n=0}^{infty }c_{n}(z-a)^{n},}f(z) =  sum_{n=0}^infty c_n (z-a)^n,

where




a is a complex constant, the center of the disk of convergence,


cn is the nth complex coefficient, and


z is a complex variable.


The radius of convergence r is a nonnegative real number or {displaystyle infty }infty such that the series converges if


|z−a|<r{displaystyle |z-a|<r,}|z-a| < r,

and diverges if


|z−a|>r.{displaystyle |z-a|>r.,}|z-a| > r.,

Some may prefer an alternative definition, as existence is obvious:


r=sup{|z−a| | ∑n=0∞cn(z−a)n  converges }{displaystyle r=sup left{|z-a| left| sum _{n=0}^{infty }c_{n}(z-a)^{n} {text{ converges }}right.right}}{displaystyle r=sup left{|z-a| left| sum _{n=0}^{infty }c_{n}(z-a)^{n} {text{ converges }}right.right}}

On the boundary, that is, where |z − a| = r, the behavior of the power series may be complicated, and the series may converge for some values of z and diverge for others. The radius of convergence is infinite if the series converges for all complex numbers z.[1]



Finding the radius of convergence


Two cases arise. The first case is theoretical: when you know all the coefficients cn{displaystyle c_{n}}c_{n} then you take certain limits and find the precise radius of convergence. The second case is practical: when you construct a power series solution of a difficult problem you typically will only know a finite number of terms in a power series, anywhere from a couple of terms to a hundred terms. In this second case, extrapolating a plot estimates the radius of convergence.



Theoretical radius


The radius of convergence can be found by applying the root test to the terms of the series. The root test uses the number


C=lim supn→|cn(z−a)n|n=lim supn→|cn|n|z−a|{displaystyle C=limsup _{nrightarrow infty }{sqrt[{n}]{|c_{n}(z-a)^{n}|}}=limsup _{nrightarrow infty }{sqrt[{n}]{|c_{n}|}}|z-a|}C = limsup_{nrightarrowinfty}sqrt[n]{|c_n(z-a)^n|} = limsup_{nrightarrowinfty}sqrt[n]{|c_n|}|z-a|

"lim sup" denotes the limit superior. The root test states that the series converges if C < 1 and diverges if C > 1. It follows that the power series converges if the distance from z to the center a is less than


r=1lim supn→|cn|n{displaystyle r={frac {1}{limsup _{nrightarrow infty }{sqrt[{n}]{|c_{n}|}}}}}r = frac{1}{limsup_{nrightarrowinfty}sqrt[n]{|c_n|}}

and diverges if the distance exceeds that number; this statement is the Cauchy–Hadamard theorem. Note that r = 1/0 is interpreted as an infinite radius, meaning that ƒ is an entire function.


The limit involved in the ratio test is usually easier to compute, and when that limit exists, it shows that the radius of convergence is finite.


r=limn→|cncn+1|.{displaystyle r=lim _{nrightarrow infty }left|{frac {c_{n}}{c_{n+1}}}right|.}{displaystyle r=lim _{nrightarrow infty }left|{frac {c_{n}}{c_{n+1}}}right|.}

This is shown as follows. The ratio test says the series converges if


limn→|cn+1(z−a)n+1||cn(z−a)n|<1.{displaystyle lim _{nto infty }{frac {|c_{n+1}(z-a)^{n+1}|}{|c_{n}(z-a)^{n}|}}<1.} lim_{ntoinfty} frac{|c_{n+1}(z-a)^{n+1}|}{|c_n(z-a)^n|} < 1.

That is equivalent to


|z−a|<1limn→|cn+1||cn|=limn→|cncn+1|.{displaystyle |z-a|<{frac {1}{lim _{nto infty }{frac {|c_{n+1}|}{|c_{n}|}}}}=lim _{nto infty }left|{frac {c_{n}}{c_{n+1}}}right|.} |z - a| < frac{1}{lim_{ntoinfty} frac{|c_{n+1}|}{|c_n|}} = lim_{ntoinfty} left|frac{c_n}{c_{n+1}}right|.


Practical estimation of radius in the case of real coefficients




Domb–Sykes plot of the function f(ε)=ε(1+ε3)1+2ε.{displaystyle f(varepsilon )={frac {varepsilon (1+varepsilon ^{3})}{sqrt {1+2varepsilon }}}.}f(varepsilon)=frac{varepsilon(1+varepsilon^3)}{sqrt{1+2varepsilon}}.[2] On the left (a) is a straightforward plot of the ratio of the power-series coefficients cn−1/cn{displaystyle c_{n-1}/c_{n}}{displaystyle c_{n-1}/c_{n}} as a function of index n{displaystyle n}n; on the right, (b) is the Domb–Sykes plot of cn/cn−1{displaystyle c_{n}/c_{n-1}}c_{n}/c_{{n-1}} as a function of 1/n{displaystyle 1/n}1/n. The solid green line is the straight-line asymptote in the Domb–Sykes plot, which intercepts the vertical axis at −2 and has a slope +1. Thus there is a singularity at ε=−1/2{displaystyle varepsilon =-1/2}{displaystyle varepsilon =-1/2} and so the radius of convergence is r=1/2.{displaystyle r=1/2.}{displaystyle r=1/2.}


Usually, in scientific applications, only a finite number of coefficients cn{displaystyle c_{n}}c_{n} are known. Typically,[vague] as n{displaystyle n}n increases, these coefficients settle into a regular behavior determined by the nearest radius-limiting singularity. In this case, two main techniques have been developed, based on the fact that the coefficients of a Taylor series are roughly exponential with ratio 1/r{displaystyle 1/r}1/r where r is the radius of convergence.



  • The basic case is when the coefficients ultimately share a common sign or alternate in sign. As pointed out earlier in the article, in many cases the limit limn→cn/cn−1{displaystyle lim _{nto infty }{{c_{n}}/{c_{n-1}}}}{displaystyle lim _{nto infty }{{c_{n}}/{c_{n-1}}}} exists, and in this case 1/r=limn→cn/cn−1.{displaystyle {1}/{r}=lim _{nto infty }{{c_{n}}/{c_{n-1}}}.}{displaystyle {1}/{r}=lim _{nto infty }{{c_{n}}/{c_{n-1}}}.} Negative r{displaystyle r}r means the convergence-limiting singularity is on the negative axis. Estimate this limit, by plotting the cn/cn−1{displaystyle c_{n}/c_{n-1}}c_{n}/c_{{n-1}} versus 1/n{displaystyle 1/n}1/n, and graphically extrapolate to 1/n=0{displaystyle 1/n=0}1/n=0 (effectively n=∞{displaystyle n=infty }n=infty ) via a linear fit. The intercept with 1/n=0{displaystyle 1/n=0}1/n=0 estimates the reciprocal of the radius of convergence, 1/r{displaystyle 1/r}1/r. This plot is called a Domb–Sykes plot.

  • The more complicated case is when the signs of the coefficients have a more complex pattern. Mercer and Roberts proposed the following procedure.[3] Define the associated sequence


bn2=cn+1cn−1−cn2cncn−2−cn−12n=3,4,5,….{displaystyle b_{n}^{2}={frac {c_{n+1}c_{n-1}-c_{n}^{2}}{c_{n}c_{n-2}-c_{n-1}^{2}}}quad n=3,4,5,ldots .}{displaystyle b_{n}^{2}={frac {c_{n+1}c_{n-1}-c_{n}^{2}}{c_{n}c_{n-2}-c_{n-1}^{2}}}quad n=3,4,5,ldots .}

Plot the finitely many known bn{displaystyle b_{n}}b_{n} versus 1/n{displaystyle 1/n}1/n, and graphically extrapolate to 1/n=0{displaystyle 1/n=0}1/n=0 via a linear fit.

The intercept with 1/n=0{displaystyle 1/n=0}1/n=0 estimates the reciprocal of the radius of convergence, 1/r{displaystyle 1/r}1/r.


This procedure also estimates two other characteristics of the convergence limiting singularity. Suppose the nearest singularity is of degree p{displaystyle p}p and has angle ±θ{displaystyle pm theta }{displaystyle pm theta } to the real axis. Then the slope of the linear fit given above is (p+1)/r{displaystyle -(p+1)/r}{displaystyle -(p+1)/r}. Further, plot 12(cn−1bncn+cn+1cnbn){displaystyle {frac {1}{2}}left({frac {c_{n-1}b_{n}}{c_{n}}}+{frac {c_{n+1}}{c_{n}b_{n}}}right)}{displaystyle {frac {1}{2}}left({frac {c_{n-1}b_{n}}{c_{n}}}+{frac {c_{n+1}}{c_{n}b_{n}}}right)} versus 1/n2{displaystyle 1/n^{2}}1/n^{2}, then a linear fit extrapolated to 1/n2=0{displaystyle 1/n^{2}=0}{displaystyle 1/n^{2}=0} has intercept at cos⁡θ{displaystyle cos theta }cos theta .


Radius of convergence in complex analysis


A power series with a positive radius of convergence can be made into a holomorphic function by taking its argument to be a complex variable. The radius of convergence can be characterized by the following theorem:


The radius of convergence of a power series ƒ centered on a point a is equal to the distance from a to the nearest point where ƒ cannot be defined in a way that makes it holomorphic.

The set of all points whose distance to a is strictly less than the radius of convergence is called the disk of convergence.




A graph of the functions explained in the text: Approximations in blue, circle of convergence in white


The nearest point means the nearest point in the complex plane, not necessarily on the real line, even if the center and all coefficients are real. For example, the function


f(z)=11+z2{displaystyle f(z)={frac {1}{1+z^{2}}}}{displaystyle f(z)={frac {1}{1+z^{2}}}}

has no singularities on the real line, since 1+z2{displaystyle 1+z^{2}}1+z^2 has no real roots. Its Taylor series about 0 is given by


n=0∞(−1)nz2n.{displaystyle sum _{n=0}^{infty }(-1)^{n}z^{2n}.}sum_{n=0}^infty (-1)^n z^{2n}.

The root test shows that its radius of convergence is 1. In accordance with this, the function ƒ(z) has singularities at ±i, which are at a distance 1 from 0.


For a proof of this theorem, see analyticity of holomorphic functions.



A simple example


The arctangent function of trigonometry can be expanded in a power series:


arctan⁡(z)=z−z33+z55−z77+⋯.{displaystyle arctan(z)=z-{frac {z^{3}}{3}}+{frac {z^{5}}{5}}-{frac {z^{7}}{7}}+cdots .}{displaystyle arctan(z)=z-{frac {z^{3}}{3}}+{frac {z^{5}}{5}}-{frac {z^{7}}{7}}+cdots .}

It is easy to apply the root test in this case to find that the radius of convergence is 1.



A more complicated example


Consider this power series:


zez−1=∑n=0∞Bnn!zn{displaystyle {frac {z}{e^{z}-1}}=sum _{n=0}^{infty }{frac {B_{n}}{n!}}z^{n}}{displaystyle {frac {z}{e^{z}-1}}=sum _{n=0}^{infty }{frac {B_{n}}{n!}}z^{n}}

where the rational numbers Bn are the Bernoulli numbers. It may be cumbersome to try to apply the ratio test to find the radius of convergence of this series. But the theorem of complex analysis stated above quickly solves the problem. At z = 0, there is in effect no singularity since the singularity is removable. The only non-removable singularities are therefore located at the other points where the denominator is zero. We solve


ez−1=0{displaystyle e^{z}-1=0,}e^z-1=0,

by recalling that if z = x + iy and e iy = cos(y) + i sin(y) then


ez=exeiy=ex(cos⁡(y)+isin⁡(y)),{displaystyle e^{z}=e^{x}e^{iy}=e^{x}(cos(y)+isin(y)),,}e^z = e^x e^{iy} = e^x(cos(y)+isin(y)),,

and then take x and y to be real. Since y is real, the absolute value of cos(y) + i sin(y) is necessarily 1. Therefore, the absolute value of e z can be 1 only if e x is 1; since x is real, that happens only if x = 0. Therefore z is pure imaginary and cos(y) + i sin(y) = 1. Since y is real, that happens only if cos(y) = 1 and sin(y) = 0, so that y is an integer multiple of 2π. Consequently the singular points of this function occur at



z = a nonzero integer multiple of 2πi.

The singularities nearest 0, which is the center of the power series expansion, are at ±2πi. The distance from the center to either of those points is 2π, so the radius of convergence is 2π.



Convergence on the boundary


If the power series is expanded around the point a and the radius of convergence is r, then the set of all points z such that |za| = r is a circle called the boundary of the disk of convergence. A power series may diverge at every point on the boundary, or diverge on some points and converge at other points, or converge at all the points on the boundary. Furthermore, even if the series converges everywhere on the boundary (even uniformly), it does not necessarily converge absolutely.


Example 1: The power series for the function ƒ(z) = 1/(1 − z), expanded around z = 0, which is simply


n=0∞zn,{displaystyle sum _{n=0}^{infty }z^{n},} sum_{n=0}^infty z^n,

has radius of convergence 1 and diverges at every point on the boundary.


Example 2: The power series for g(z) = −ln(1 − z), expanded around z = 0, which is


n=1∞1nzn,{displaystyle sum _{n=1}^{infty }{frac {1}{n}}z^{n},} sum_{n=1}^infty frac{1}{n} z^n,

has radius of convergence 1, and diverges for z = 1 but converges for all other points on the boundary. The function ƒ(z) of Example 1 is the derivative of g(z).


Example 3: The power series


n=1∞1n2zn{displaystyle sum _{n=1}^{infty }{frac {1}{n^{2}}}z^{n}}{displaystyle sum _{n=1}^{infty }{frac {1}{n^{2}}}z^{n}}

has radius of convergence 1 and converges everywhere on the boundary absolutely. If h is the function represented by this series on the unit disk, then the derivative of h(z) is equal to g(z)/z with g of Example 2. It turns out that h(z) is the dilogarithm function.


Example 4: The power series


i=1∞aizi where ai=(−1)n−12nn for n=⌊log2⁡(i)⌋+1, the unique integer with 2n−1≤i<2n,{displaystyle sum _{i=1}^{infty }a_{i}z^{i}{text{ where }}a_{i}={frac {(-1)^{n-1}}{2^{n}n}}{text{ for }}n=lfloor log _{2}(i)rfloor +1{text{, the unique integer with }}2^{n-1}leq i<2^{n},}{displaystyle sum _{i=1}^{infty }a_{i}z^{i}{text{ where }}a_{i}={frac {(-1)^{n-1}}{2^{n}n}}{text{ for }}n=lfloor log _{2}(i)rfloor +1{text{, the unique integer with }}2^{n-1}leq i<2^{n},}

has radius of convergence 1 and converges uniformly on the entire boundary {|z| = 1}, but does not converge absolutely on the boundary.[4]



Rate of convergence


If we expand the function


f(x)=sin⁡x=∑n=0∞(−1)n(2n+1)!x2n+1=x−x33!+x55!− for all x{displaystyle f(x)=sin x=sum _{n=0}^{infty }{frac {(-1)^{n}}{(2n+1)!}}x^{2n+1}=x-{frac {x^{3}}{3!}}+{frac {x^{5}}{5!}}-cdots {text{ for all }}x}f(x)=sin x = sum^{infin}_{n=0} frac{(-1)^n}{(2n+1)!} x^{2n+1} =  x - frac{x^3}{3!} + frac{x^5}{5!} - cdotstext{ for all } x

around the point x = 0, we find out that the radius of convergence of this series is {displaystyle scriptstyle infty }scriptstyleinfty meaning that this series converges for all complex numbers. However, in applications, one is often interested in the precision of a numerical answer. Both the number of terms and the value at which the series is to be evaluated affect the accuracy of the answer. For example, if we want to calculate f(0.1) = sin(0.1) accurate up to five decimal places, we only need the first two terms of the series. However, if we want the same precision for x = 1 we must evaluate and sum the first five terms of the series. For f(10), one requires the first 18 terms of the series, and for f(100) we need to evaluate the first 141 terms.


So the fastest convergence of a power series expansion is at the center, and as one moves away from the center of convergence, the rate of convergence slows down until you reach the boundary (if it exists) and cross over, in which case the series will diverge.



Abscissa of convergence of a Dirichlet series


An analogous concept is the abscissa of convergence of a Dirichlet series


n=1∞anns.{displaystyle sum _{n=1}^{infty }{a_{n} over n^{s}}.}sum_{n=1}^infty {a_n over n^s}.

Such a series converges if the real part of s is greater than a particular number depending on the coefficients an: the abscissa of convergence.



Notes





  1. ^ Mathematical Analysis-II. Krishna Prakashan Media..mw-parser-output cite.citation{font-style:inherit}.mw-parser-output .citation q{quotes:"""""""'""'"}.mw-parser-output .citation .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/6/65/Lock-green.svg/9px-Lock-green.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-limited a,.mw-parser-output .citation .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/d/d6/Lock-gray-alt-2.svg/9px-Lock-gray-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/a/aa/Lock-red-alt-2.svg/9px-Lock-red-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration{color:#555}.mw-parser-output .cs1-subscription span,.mw-parser-output .cs1-registration span{border-bottom:1px dotted;cursor:help}.mw-parser-output .cs1-ws-icon a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/4/4c/Wikisource-logo.svg/12px-Wikisource-logo.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output code.cs1-code{color:inherit;background:inherit;border:inherit;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;font-size:100%}.mw-parser-output .cs1-visible-error{font-size:100%}.mw-parser-output .cs1-maint{display:none;color:#33aa33;margin-left:0.3em}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration,.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left,.mw-parser-output .cs1-kern-wl-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right,.mw-parser-output .cs1-kern-wl-right{padding-right:0.2em}


  2. ^ See Figure 8.1 in: Hinch, E.J. (1991), Perturbation Methods, Cambridge Texts in Applied Mathematics, 6, Cambridge University Press, p. 146, ISBN 0-521-37897-4


  3. ^ Mercer, G.N.; Roberts, A.J. (1990), "A centre manifold description of contaminant dispersion in channels with varying flow properties", SIAM J. Appl. Math., 50 (6): 1547–1565, doi:10.1137/0150091


  4. ^ Sierpiński, Wacław (1918), "O szeregu potęgowym który jest zbieżny na całem swem kole zbieżności jednostajnie ale nie bezwzględnie", Prace matematyka-fizyka, 29, pp. 263–266




References




  • Brown, James; Churchill, Ruel (1989), Complex variables and applications, New York: McGraw-Hill, ISBN 978-0-07-010905-6


  • Stein, Elias; Shakarchi, Rami (2003), Complex Analysis, Princeton, New Jersey: Princeton University Press, ISBN 0-691-11385-8



External links


  • What is radius of convergence?



這個網誌中的熱門文章

12.7 cm/40 Type 89 naval gun

Shark

Wiciokrzew