■マチアセビッチとフィボナッチ数生成関数(その95)

 ここでは,√5に収束する数列を考えることにします.

2^2-5・1^2=-1

9^2-5・4^2=+1

===================================

 まず

  (2+√5)^n=an+bn√5

  (2−√5)^n=an−bn√5

を満足させるような整数列{an},{bn}を考えます.これらの数列は

  an^2−5n^2=(−1)^n

となる関係式で結ばれていて,

  an/bn→ √5

ですから,√5に最も近い分数を与えることがわかります(最良近似).

  an+1+bn+1√5=(2+√5)(an+bn√5)

          =(2an+5bn)+(an+2bn)√5

より

  an+1=2an+5bn,bn+1=an+2bn

  an+1=2an+5bn=2an+5(an-1+2bn-1)

 =2an+an-1+2(2an-1+5bn-1)=4an+an-1

  bn+1=an+2bn=(2an-1+5bn-1)+2bn

 =2(an-1+2bn-1)+2bn+bn-1=4bn+bn-1

より

  an+1=4an+an-1,bn+1=4bn+bn-1

 α,βを2次方程式x^2−4x−1=0の根{2+/-√5}として,

  an+1−αan=β(an−αan-1)=β^2(an-1−αan-2)=・・・=β^(n-1)(a2−αa1)

α,βを入れ替えると

  an+1−βan=α^(n-1)(a2−βa1)

  an+1−αan=β^(n-1)(a2−αa1)

 したがって,整数列{an}の一般項は

  an={α^(n-1)(a2−βa1)−β^(n-1)(a2−αa1)}/(α−β)

α=(2+√5),β=(2-√5),初期値をa1=2,a2=9とすると

  an=1/2√5{{(2+√5}^(n-1)(5+2√5)-{(2-√5)}^(n-1)(5-2√5)}

 整数列{bn}でも同じ漸化式ですから,同じ一般項になります.

  bn={α^(n-1)(b2−βb1)−β^(n-1)(b2−αb1)}/(α−β)

初期値をb1=1,b2=4とすると

  bn=1/2√5{(2+√5)}^(n-1)(2+√5)-{(2-√5}^(n-1)(2-√5)}

 ここで,n→∞のとき(2−√5)^n→0ですから

  an/bn→ √5

となるのを確かめることができます.

===================================

  an=1/2{{(2+√5}^(n-1)(√5+2)-{(2-√5)}^(n-1)(√5-2)}

  an=1/2{{(2+√5}^(n)+{(2-√5)}^(n)}

  bn=1/2√5{(2+√5)}^(n)-{(2-√5}^(n)}

===================================