Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Periodic continued fraction

From Wikipedia, the free encyclopedia
This articleneeds additional citations forverification. Please helpimprove this article byadding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Periodic continued fraction" – news ·newspapers ·books ·scholar ·JSTOR
(January 2014) (Learn how and when to remove this message)

Inmathematics, an infiniteperiodic continued fraction is asimple continued fraction that can be placed in the form

x=a0+1a1+1a2+1ak+1ak+1+ak+m1+1ak+m+1ak+1+1ak+2+{\displaystyle x=a_{0}+{\cfrac {1}{a_{1}+{\cfrac {1}{a_{2}+{\cfrac {1}{\quad \ddots \quad a_{k}+{\cfrac {1}{a_{k+1}+{\cfrac {\ddots }{\quad \ddots \quad a_{k+m-1}+{\cfrac {1}{a_{k+m}+{\cfrac {1}{a_{k+1}+{\cfrac {1}{a_{k+2}+{\ddots }}}}}}}}}}}}}}}}}}

where the initial block[a0;a1,,ak]{\displaystyle [a_{0};a_{1},\dots ,a_{k}]} ofk+1 partial denominators is followed by a block[ak+1,ak+2,,ak+m]{\displaystyle [a_{k+1},a_{k+2},\dots ,a_{k+m}]} ofm partial denominators that repeatsad infinitum. For example,2{\displaystyle {\sqrt {2}}} can be expanded to the periodic continued fraction[1;2,2,2,...]{\displaystyle [1;2,2,2,...]}.

This article considers only the case of periodicregular continued fractions. In other words, the remainder of this article assumes that all the partial denominatorsai (i ≥ 1) are positive integers. The general case, where the partial denominatorsai are arbitrary real or complex numbers, is treated in the articleconvergence problem.

Purely periodic and periodic fractions

[edit]

Since all the partial numerators in a regular continued fraction are equal to unity we can adopt a shorthand notation in which the continued fraction shown above is written as

x=[a0;a1,a2,,ak,ak+1,ak+2,,ak+m,ak+1,ak+2,,ak+m,]=[a0;a1,a2,,ak,ak+1,ak+2,,ak+m¯]{\displaystyle {\begin{aligned}x&=[a_{0};a_{1},a_{2},\dots ,a_{k},a_{k+1},a_{k+2},\dots ,a_{k+m},a_{k+1},a_{k+2},\dots ,a_{k+m},\dots ]\\&=[a_{0};a_{1},a_{2},\dots ,a_{k},{\overline {a_{k+1},a_{k+2},\dots ,a_{k+m}}}]\end{aligned}}}

where, in the second line, avinculum marks the repeating block.[1] Some textbooks use the notation

x=[a0;a1,a2,,ak,a˙k+1,ak+2,,a˙k+m]{\displaystyle {\begin{aligned}x&=[a_{0};a_{1},a_{2},\dots ,a_{k},{\dot {a}}_{k+1},a_{k+2},\dots ,{\dot {a}}_{k+m}]\end{aligned}}}

where the repeating block is indicated by dots over its first and last terms.[2]

If the initial non-repeating block is not present – that is, if k = -1, a0 = am and

x=[a0;a1,a2,,am1¯],{\displaystyle x=[{\overline {a_{0};a_{1},a_{2},\dots ,a_{m-1}}}],}

the regular continued fractionx is said to bepurely periodic. For example, the regular continued fraction[1;1,1,1,]{\displaystyle [1;1,1,1,\dots ]} of thegolden ratio φ is purely periodic, while the regular continued fraction[1;2,2,2,]{\displaystyle [1;2,2,2,\dots ]} of2{\displaystyle {\sqrt {2}}} is periodic, but not purely periodic. However, the regular continued fraction[2;2,2,2,]{\displaystyle [2;2,2,2,\dots ]} of thesilver ratioσ=2+1{\displaystyle \sigma ={\sqrt {2}}+1} is purely periodic.

As unimodular matrices

[edit]

Periodic continued fractions are in one-to-one correspondence with the realquadratic irrationals. The correspondence is explicitly provided byMinkowski's question-mark function. That article also reviews tools that make it easy to work with such continued fractions. Consider first the purely periodic part

x=[0;a1,a2,,am¯],{\displaystyle x=[0;{\overline {a_{1},a_{2},\dots ,a_{m}}}],}

This can, in fact, be written as

x=αx+βγx+δ{\displaystyle x={\frac {\alpha x+\beta }{\gamma x+\delta }}}

with theα,β,γ,δ{\displaystyle \alpha ,\beta ,\gamma ,\delta } being integers, and satisfyingαδβγ=1.{\displaystyle \alpha \delta -\beta \gamma =1.} Explicit values can be obtained by writing

S=(1011){\displaystyle S={\begin{pmatrix}1&0\\1&1\end{pmatrix}}}

which is termed a "shift", so that

Sn=(10n1){\displaystyle S^{n}={\begin{pmatrix}1&0\\n&1\end{pmatrix}}}

and similarly a reflection, given by

T(1101){\displaystyle T\mapsto {\begin{pmatrix}-1&1\\0&1\end{pmatrix}}}

so thatT2=I{\displaystyle T^{2}=I}. Both of these matrices areunimodular, arbitrary products remain unimodular. Then, givenx{\displaystyle x} as above, the corresponding matrix is of the form[3]

Sa1TSa2TTSam=(αβγδ){\displaystyle S^{a_{1}}TS^{a_{2}}T\cdots TS^{a_{m}}={\begin{pmatrix}\alpha &\beta \\\gamma &\delta \end{pmatrix}}}

and one has

x=[0;a1,a2,,am¯]=αx+βγx+δ{\displaystyle x=[0;{\overline {a_{1},a_{2},\dots ,a_{m}}}]={\frac {\alpha x+\beta }{\gamma x+\delta }}}

as the explicit form. As all of the matrix entries are integers, this matrix belongs to themodular groupSL(2,Z).{\displaystyle SL(2,\mathbb {Z} ).}

Relation to quadratic irrationals

[edit]

Aquadratic irrational number is anirrational real root of the quadratic equation

ax2+bx+c=0{\displaystyle ax^{2}+bx+c=0}

where the coefficientsa,b, andc are integers, and thediscriminant,b24ac{\displaystyle b^{2}-4ac}, is greater than zero. By thequadratic formula, every quadratic irrational can be written in the form

ζ=P+DQ{\displaystyle \zeta ={\frac {P+{\sqrt {D}}}{Q}}}

whereP,D, andQ are integers,D > 0 is not aperfect square (but not necessarily square-free), andQ divides the quantityP2D{\displaystyle P^{2}-D} (for example(6+8)/4{\displaystyle (6+{\sqrt {8}})/4}). Such a quadratic irrational may also be written in another form with a square-root of a square-free number (for example(3+2)/2{\displaystyle (3+{\sqrt {2}})/2}) as explained forquadratic irrationals.

By considering thecomplete quotients of periodic continued fractions,Euler was able to prove that ifx is a regular periodic continued fraction, thenx is a quadratic irrational number. The proof is straightforward. From the fraction itself, one can construct the quadratic equation with integral coefficients thatx must satisfy.

Lagrange proved the converse of Euler's theorem: ifx is a quadratic irrational, then the regular continued fraction expansion ofx is periodic.[4] Given a quadratic irrationalx one can constructm different quadratic equations, each with the same discriminant, that relate the successive complete quotients of the regular continued fraction expansion ofx to one another. Since there are only finitely many of these equations (the coefficients are bounded), the complete quotients (and also the partial denominators) in the regular continued fraction that representsx must eventually repeat.

Reduced surds

[edit]

The quadraticsurdζ=P+DQ{\displaystyle \zeta ={\frac {P+{\sqrt {D}}}{Q}}} is said to bereduced ifζ>1{\displaystyle \zeta >1} and itsconjugateη=PDQ{\displaystyle \eta ={\frac {P-{\sqrt {D}}}{Q}}}satisfies the inequalities1<η<0{\displaystyle -1<\eta <0}. For instance, the golden ratioϕ=(1+5)/2=1.618033...{\displaystyle \phi =(1+{\sqrt {5}})/2=1.618033...} is a reduced surd because it is greater than one and its conjugate(15)/2=0.618033...{\displaystyle (1-{\sqrt {5}})/2=-0.618033...} is greater than −1 and less than zero. On the other hand, the square root of two2=(0+8)/2{\displaystyle {\sqrt {2}}=(0+{\sqrt {8}})/2} is greater than one but is not a reduced surd because its conjugate2=(08)/2{\displaystyle -{\sqrt {2}}=(0-{\sqrt {8}})/2} is less than −1.

Galois proved that the regular continued fraction which represents a quadratic surd ζ is purely periodic if and only if ζ is a reduced surd. In fact, Galois showed more than this. He also proved that if ζ is a reduced quadratic surd and η is its conjugate, then the continued fractions for ζ and for (−1/η) are both purely periodic, and the repeating block in one of those continued fractions is the mirror image of the repeating block in the other. In symbols we have

ζ=[a0;a1,a2,,am1¯]1η=[am1;am2,am3,,a0¯]{\displaystyle {\begin{aligned}\zeta &=[{\overline {a_{0};a_{1},a_{2},\dots ,a_{m-1}}}]\\[3pt]{\frac {-1}{\eta }}&=[{\overline {a_{m-1};a_{m-2},a_{m-3},\dots ,a_{0}}}]\,\end{aligned}}}

where ζ is any reduced quadratic surd, and η is its conjugate.

From these two theorems of Galois a result already known to Lagrange can be deduced. Ifr > 1 is a rational number that is not a perfect square, then

r=[a0;a1,a2,,a2,a1,2a0¯].{\displaystyle {\sqrt {r}}=[a_{0};{\overline {a_{1},a_{2},\dots ,a_{2},a_{1},2a_{0}}}].}

In particular, ifn is any non-square positive integer, the regular continued fraction expansion ofn contains a repeating block of lengthm, in which the firstm − 1 partial denominators form apalindromic string.

Length of the repeating block

[edit]

By analyzing the sequence of combinations

Pn+DQn{\displaystyle {\frac {P_{n}+{\sqrt {D}}}{Q_{n}}}}

that can possibly arise whenζ=P+DQ{\displaystyle \zeta ={\frac {P+{\sqrt {D}}}{Q}}} is expanded as a regular continued fraction,Lagrange showed that the largest partial denominatorai in the expansion is less than2D{\displaystyle 2{\sqrt {D}}}, and that the length of the repeating block is less than 2D.

More recently, sharper arguments[5][6][7] based on thedivisor function have shown that the length of the repeating block for a quadratic surd of discriminantD ison the order ofO(DlnD).{\displaystyle {\mathcal {O}}({\sqrt {D}}\ln {D}).}

Canonical form and repetend

[edit]

The following iterative algorithm[8] can be used to obtain the continued fraction expansion in canonical form (S is anynatural number that is not aperfect square):

m0=0{\displaystyle m_{0}=0\,\!}
d0=1{\displaystyle d_{0}=1\,\!}
a0=S{\displaystyle a_{0}=\left\lfloor {\sqrt {S}}\right\rfloor \,\!}
mn+1=dnanmn{\displaystyle m_{n+1}=d_{n}a_{n}-m_{n}\,\!}
dn+1=Smn+12dn{\displaystyle d_{n+1}={\frac {S-m_{n+1}^{2}}{d_{n}}}\,\!}
an+1=S+mn+1dn+1=a0+mn+1dn+1.{\displaystyle a_{n+1}=\left\lfloor {\frac {{\sqrt {S}}+m_{n+1}}{d_{n+1}}}\right\rfloor =\left\lfloor {\frac {a_{0}+m_{n+1}}{d_{n+1}}}\right\rfloor \!.}

Notice thatmn,dn, andan are always integers.The algorithm terminates when this triplet is the same as one encountered before.The algorithm can also terminate on ai when ai = 2 a0,[9] which is easier to implement.

The expansion will repeat from then on. The sequence[a0;a1,a2,a3,]{\displaystyle [a_{0};a_{1},a_{2},a_{3},\dots ]} is the continued fraction expansion:

S=a0+1a1+1a2+1a3+{\displaystyle {\sqrt {S}}=a_{0}+{\cfrac {1}{a_{1}+{\cfrac {1}{a_{2}+{\cfrac {1}{a_{3}+\,\ddots }}}}}}}

Example

[edit]

To obtain114 as a continued fraction, begin withm0 = 0;d0 = 1; anda0 = 10 (102 = 100 and 112 = 121 > 114 so 10 chosen).

114=114+01=10+114101=10+(11410)(114+10)114+10=10+114100114+10=10+1114+1014.{\displaystyle {\begin{aligned}{\sqrt {114}}&={\frac {{\sqrt {114}}+0}{1}}=10+{\frac {{\sqrt {114}}-10}{1}}=10+{\frac {({\sqrt {114}}-10)({\sqrt {114}}+10)}{{\sqrt {114}}+10}}\\&=10+{\frac {114-100}{{\sqrt {114}}+10}}=10+{\frac {1}{\frac {{\sqrt {114}}+10}{14}}}.\end{aligned}}}
m1=d0a0m0=1100=10.{\displaystyle m_{1}=d_{0}\cdot a_{0}-m_{0}=1\cdot 10-0=10\,.}
d1=Sm12d0=1141021=14.{\displaystyle d_{1}={\frac {S-m_{1}^{2}}{d_{0}}}={\frac {114-10^{2}}{1}}=14\,.}
a1=a0+m1d1=10+1014=2014=1.{\displaystyle a_{1}=\left\lfloor {\frac {a_{0}+m_{1}}{d_{1}}}\right\rfloor =\left\lfloor {\frac {10+10}{14}}\right\rfloor =\left\lfloor {\frac {20}{14}}\right\rfloor =1\,.}

So,m1 = 10;d1 = 14; anda1 = 1.

114+1014=1+114414=1+1141614(114+4)=1+1114+47.{\displaystyle {\frac {{\sqrt {114}}+10}{14}}=1+{\frac {{\sqrt {114}}-4}{14}}=1+{\frac {114-16}{14({\sqrt {114}}+4)}}=1+{\frac {1}{\frac {{\sqrt {114}}+4}{7}}}.}

Next,m2 = 4;d2 = 7; anda2 = 2.

114+47=2+114107=2+147(114+10)=2+1114+102.{\displaystyle {\frac {{\sqrt {114}}+4}{7}}=2+{\frac {{\sqrt {114}}-10}{7}}=2+{\frac {14}{7({\sqrt {114}}+10)}}=2+{\frac {1}{\frac {{\sqrt {114}}+10}{2}}}.}
114+102=10+114102=10+142(114+10)=10+1114+107.{\displaystyle {\frac {{\sqrt {114}}+10}{2}}=10+{\frac {{\sqrt {114}}-10}{2}}=10+{\frac {14}{2({\sqrt {114}}+10)}}=10+{\frac {1}{\frac {{\sqrt {114}}+10}{7}}}.}
114+107=2+11447=2+987(114+4)=2+1114+414.{\displaystyle {\frac {{\sqrt {114}}+10}{7}}=2+{\frac {{\sqrt {114}}-4}{7}}=2+{\frac {98}{7({\sqrt {114}}+4)}}=2+{\frac {1}{\frac {{\sqrt {114}}+4}{14}}}.}
114+414=1+1141014=1+1414(114+10)=1+1114+101.{\displaystyle {\frac {{\sqrt {114}}+4}{14}}=1+{\frac {{\sqrt {114}}-10}{14}}=1+{\frac {14}{14({\sqrt {114}}+10)}}=1+{\frac {1}{\frac {{\sqrt {114}}+10}{1}}}.}
114+101=20+114101=20+14114+10=20+1114+1014.{\displaystyle {\frac {{\sqrt {114}}+10}{1}}=20+{\frac {{\sqrt {114}}-10}{1}}=20+{\frac {14}{{\sqrt {114}}+10}}=20+{\frac {1}{\frac {{\sqrt {114}}+10}{14}}}.}

Now, loop back to the second equation above.

Consequently, the simple continued fraction for the square root of 114 is

114=[10;1,2,10,2,1,20¯].{\displaystyle {\sqrt {114}}=[10;{\overline {1,2,10,2,1,20}}].\,} (sequenceA010179 in theOEIS)

114 is approximately 10.67707 82520. After one expansion of the repetend, the continued fraction yields the rational fraction211941985{\displaystyle {\frac {21194}{1985}}} whose decimal value is approx. 10.67707 80856, a relative error of0.0000016% or 1.6 parts in 100,000,000.

Generalized continued fraction

[edit]

A more rapid method is to evaluate itsgeneralized continued fraction. From the formula derivedthere:

z=x2+y=x+y2x+y2x+y2x+=x+2xy2(2zy)yy22(2zy)y22(2zy){\displaystyle {\begin{aligned}{\sqrt {z}}={\sqrt {x^{2}+y}}&=x+{\cfrac {y}{2x+{\cfrac {y}{2x+{\cfrac {y}{2x+\ddots }}}}}}\\&=x+{\cfrac {2x\cdot y}{2(2z-y)-y-{\cfrac {y^{2}}{2(2z-y)-{\cfrac {y^{2}}{2(2z-y)-\ddots }}}}}}\end{aligned}}}

and the fact that 114 is 2/3 of the way between 102=100 and 112=121 results in

114=10263=322+23=323+2/364+264+264+264+=323+2192+18192+18192+,{\displaystyle {\begin{aligned}{\sqrt {114}}={\cfrac {\sqrt {1026}}{3}}={\cfrac {\sqrt {32^{2}+2}}{3}}&={\cfrac {32}{3}}+{\cfrac {2/3}{64+{\cfrac {2}{64+{\cfrac {2}{64+{\cfrac {2}{64+\ddots }}}}}}}}\\&={\cfrac {32}{3}}+{\cfrac {2}{192+{\cfrac {18}{192+{\cfrac {18}{192+\ddots }}}}}},\end{aligned}}}

which is simply the aforementioned[10;1,2,10,2,1,20,1,2]{\displaystyle [10;1,2,\,10,2,1,\,20,1,2]} evaluated at every third term. Combining pairs of fractions produces

114=322+23=323+64/3205011205012050=323+64615039615096150,{\displaystyle {\begin{aligned}{\sqrt {114}}={\cfrac {\sqrt {32^{2}+2}}{3}}&={\cfrac {32}{3}}+{\cfrac {64/3}{2050-1-{\cfrac {1}{2050-{\cfrac {1}{2050-\ddots }}}}}}\\&={\cfrac {32}{3}}+{\cfrac {64}{6150-3-{\cfrac {9}{6150-{\cfrac {9}{6150-\ddots }}}}}},\end{aligned}}}

which is now[10;1,2,10,2,1,20,1,2¯]{\displaystyle [10;1,2,{\overline {10,2,1,20,1,2}}]} evaluated at the third term and every six terms thereafter.

See also

[edit]

Notes

[edit]
  1. ^Pettofrezzo & Byrkit 1970, p. 158.
  2. ^Long 1972, p. 187.
  3. ^Khinchin 1964.
  4. ^Davenport 1982, p. 104.
  5. ^Hickerson 1973.
  6. ^Cohn 1977.
  7. ^Podsypanin 1982.
  8. ^Beceanu 2003.
  9. ^Gliga 2006.

References

[edit]
  • Long, Calvin T. (1972).Elementary Introduction to Number Theory (3 Sub ed.). Waveland Pr Inc.LCCN 77-171950.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Periodic_continued_fraction&oldid=1283498694"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp