■ディオファントス近似(その6)
連分数に対しては殆ど計算規則は知られていない.
[a0:a1,a2,・・・]+n=[a0+n:a1,a2,・・・]
[a0:a1,a2,・・・]^-1=[0:a0,a1,a2,・・・]
以外,たとえば,
[a0:a1,a2,・・・]+[b0:b1,b2,・・・]=[c0:c1,c2,・・・]
のciをaiとbiで表す公式はない.
===================================
【1】周期1の整数
最も素朴な循環連分数は
√m=[q0;2q0,2q0,2q0,・・・]
で表されるものと考えられます.→コラム「整数の平方根の連分数」,「もうひとつの白銀比」参照
このとき,
P=2q0^2+1,Q=2q0
より,mは
(2q0^2+1)^2−m・4q0^2=±1
を満たす整数となるのですが,結局,このようなmは
m=q0^2+1=2,5,10,・・・
となることが導き出されます.
√2=[1;2,2,2,・・・]
√5=[2;4,4,4,・・・]
√10=[3;6,6,6,・・・]
√101=[10;20,20,20,・・・]
[Q]長方形から正方形を2つ切り取った後に残る長方形がもとの長方形と相似になるのは?
[A]1:x=x−2:1 → x=1+√2
[Q]長方形から正方形をn個切り取った後に残る長方形がもとの長方形と相似になるのは?
[A]1:x=x−n:1 → x=(n+√(n^2+4))/2
これは黄金比のもうひとつの一般化であるが,この操作は無限連分数
(n+√(n^2+4))/2=[n:n,n,n,,n,・・・]
で表されることと同義である.
φ=[1:1,1,1,,1,・・・]
1+√2=[2:2,2,2,2,・・・]
(3+√13)/2=[3:3,3,3,3,・・・]
===================================
【2】周期1から周期2へ
[Q]大長方形から正方形をa個切り取った後に残る中長方形から,正方形をb個切り取った後に残る小長方形が中長方形と相似になるのは?
[A]単純循環連分数
L=[a:b,b,b,b,・・・]
で表される数Lを求めてみることにしましょう.
L−a=R=[0:b,b,b,b,・・・]=1/(b+R)
R^2+bR−1=0 → R=(−b+(b^2+4)^(1/2))/2
L=a+R=a−b/2+(b^2/4+1)^(1/2)
a=1,b=1 → L=(1+√5)/2=[1;1,1,1,1,1,・・・]
a=1,b=2 → L=√2=[1;2,2,2,2,2,・・・]
a=2,b=4 → L=√5=[2;4,4,4,・・・]
[Q]大長方形から正方形をa個切り取った後に残る中長方形から,正方形をb個切り取り,さらに正方形をc個切り取った後に残る小長方形が中長方形と相似になるのは?
[A]同様に,2項が循環する連分数は
L=[a:b,c,b,c,・・・]
L−a=R=[0:b,c,b,c,・・・]=1/(b+1/(c+R))
bR^2+bcR−c=0 → R=(−c+(c^2+4c/b)^(1/2))/2
L=a−c/2+(c^2/4+c/b)^(1/2)
a=1,b=1,c=2 → L=√3=[1;1,2,1,2,1,2,・・・]
a=2,b=2,c=4 → L=√6=[2;2,4,2,4,2,・・・]
実数xの整数部分,小数部分を
x=[x]+{x}
で表します.同様に,連分数表示を整数部分と小数部分に分けます.
[a:b,b,b,b,・・・]=a+[0:b,b,b,b,・・・]=a+<b,b,b,b,・・・>
[a:b,c,b,c,・・・]=a+[0:b,c,b,c,・・・]=a+<b,c,b,c,・・・>
たとえば,小数部分<b,b,b,・・・>となる数をxとおくと,
x=1/(b+x) → x=(√(n^2+4)−n)/2
このように循環連分数は整数係数の2次方程式の解となります(ラグランジュの定理).
(√5−1)/2=[0:1,1,1,,1,・・・]
√2−1=[0:2,2,2,2,・・・]
(√13−3)/2=[0:3,3,3,3,・・・]
<1,1,1,・・・>=xはx=1/(1+x)の正の解(√5−1)/2に収束します.この近似分数は,フィボナッチ数列
a0=0,a1=1,an+2=an+an+1
1,2,3,5,8,13,・・・
の相隣る2項の比となります.
また,<2,2,2,・・・>=xはx=1/(2+x)の正の解√2−1に収束します.この近似分数は,数列
a0=1,a1=1,an+2=an+2an+1
1,2,6,12,29,70,・・・
の相隣る2項の比となります.この数列は√2−1の最良近似数列です.
===================================
【3】黄金比のもうひとつの一般化
最も素朴な循環連分数は
√m=[q0;2q0,2q0,2q0,・・・]
で表されます.このとき,
P=2q0^2+1,Q=2q0
より,mは
(2q0^2+1)^2−m・4q0^2=±1
を満たす整数となるのですが,結局,このようなmは
m=q0^2+1=2,5,10,・・・
となることが導き出されます.
√2=[1;2,2,2,・・・]
√5=[2;4,4,4,・・・]
それに対して,高貴な数とは黄金比のようにその連分数展開が無限に多くの1で終わる数
[a0:a1,・・・,an,1,1,1,・・・]
として定義される.→コラム「高貴な数」参照
たとえば,
1/τ^2=1/(1+τ)=[0:2,1,1,1,・・・]
であり,黄金比
τ=[1:1,1,1,1,・・・]
は無理数の中の最も無理な数であり,高貴な数の中で最も高貴なものである.
なお,w=[0:a1,a2,・・・]とすれば
1−w=[0:1,a1−1,a2,・・・]
wが高貴であるならば1−wもそうである.また,もし,a1−1=0ならば,1−(1−w)=wより
[・・・,am,0,am+2,・・・]=[・・・,am+am+2,・・・]
1/(1+τ)=[0:2,1,1,1,・・・]
−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
√2=[1:2,2,2,2,・・・]
を得るのに,
Ak=2Ak-1+Ak-2,A0=1,A1=1
Bk=2Bk-1+Bk-2,B0=0,B1=1
を用いて
Ak=1,1,3,7,17,41,99,・・・
Bk=0,1,2,5,12,29,70,・・・
を得る.Ak,Bkはx^2−2y^2=+1とx^2−2y^2=−1の解と交互に対応する.そして,Ak/Bkは上下から交互に√2に近づく.
ν=[a0:a1,・・・,an,1,1,1,・・・]
によって定義される高貴な数νは,
ν=(τAn+An-1)/(τBn+Bn-1)
のように表現できる.ここでAk,Bkは[a0:a1,・・・,an]のk番目の近似分数の分子と分母である.
===================================