1

双子な行列の具体例 ~フィボナッチ数とリュカ数~

56
0

概要

この記事は "双子"な行列とその必要十分条件 という記事の続きです。今回は、前回定義した双子という性質を、自分自身の転置行列との間に持つ行列の具体例を紹介します。前回までの記事を読んでいなくてもわかるように書いたつもりです。

復習

この記事から読み始めた人にも伝わるように、前回までの復習をしたいと思います。まず、双子という性質の定義は次のようでした。

双子

A=(aij),B=(bij)はすべての成分が正のN次正方行列とする。
nを正整数とする。
ai1iσ(1)ai2iσ(2)ainiσ(n)=bi1iσ(1)bi2iσ(2)biniσ(n)
が任意のi1,,in{1,2,,N}および任意のσSnで成り立つとき、ABn双子であるという。任意の正整数nに対してn双子であるとき、A,Bは単に双子であるといい、ABと表す。

そして、ABであるかを判定するのに便利な必要十分条件が、次のように示されたのでした。

双子の必要十分条件

任意のp,q,r{1,,N}に対して、
apqaqrarp=bpqbqrbrp
が成り立つことは、ABに必要十分である。

今回はAtAに限るので、この場合の条件に書き直しておきましょう。

命題 1

任意のp,q,r{1,,N}に対して、
apqaqrarp=aqparqapr
が成り立つことは、AtAに必要十分である。

後で使うときのために、この左辺において、pqを入れ替えると右辺になることに注意しておきます。

パスカルの三角形の驚くべき対称性 』では、パスカルの三角形の中に作られる正方行列(すなわちaij=(a+i+jb+j))がAtAを満たすことを見たわけです。

では、他にそのような行列は作れないでしょうか?

具体例

aij=f(i)g(j)

aij=f(i)g(j)のとき、AtAは明らかに成り立ちます。実際、
apqaqrarp=f(p)g(q)f(q)g(r)f(r)g(p)=f(q)g(p)f(r)g(q)f(p)g(r)=aqparqapr
です。これはちょっとつまらないので、スルーしましょう。

aij=f(i)g(j)h(i,j)

この場合は、h(i,j)が対称式のときに成り立ちます(必要条件ではないです)。証明は上とほぼ同じなので省略します。パスカルの三角形のときのaij=(a+i+jb+j)はこの型ですね。

aij=f(i)+g(j)

実際にapqaqrarpを計算してみればわかりますが、これは一筋縄では行きません。
しかし、次のようにすると、うまくAtAが成り立ちます。
f(i)=si, g(j)=sj; aij=si+sj
ただし、s0とします。確かめてみましょう。
apqaqrarp=(sp+sq)(sq+sr)(sr+sp)=sp+q+r+sq+sp+sr+sr+sp+sq+s(p+q+r)=sp+q+r+(sp+sq+sr)+(sp+sq+sr)+s(p+q+r)
これはp,q,rの対称式ですので、当然p,qを入れ替えたaqparqaprに一致します。
さらに、
f(i)=si, g(j)=sj; bij=sisj
とすると、
bpqbqrbrp=(spsq)(sqsr)(srsp)=sp+q+rsqsp+srsr+sp+sqs(p+q+r)=sp+q+r(sp+sq+sr)+(sp+sq+sr)s(p+q+r)
となって、やはりp,qを入れ替えても不変なので、bqpbrqbprに一致します。

例として、s=ϕ=1+52としてみましょう。訳あって、a,Aではなくl,Lを、b,Bではなくf,Fを使うことにします。
L=(lij); lij=ϕi+ϕjF=(fij); fij=ϕiϕj
と定義すれば、
LtL, FtF
が成り立つという訳です。ここで、リュカ数列{Ln}およびフィボナッチ数列{Fn}の一般項を思い出すと、
Ln=ϕn+(ϕ)n, Fn=ϕn(ϕ)n5
でしたから、上で定義した行列L,Fには、次のように対角成分にリュカ数とフィボナッチ数が交互に現れることになります。
L=(5F1L25F3), F=(L15F2L3)

まあ、だからどうしたという話ですが(笑)、面白いなと思ったのでこんなものも作ってみました。

おわりに

AtAを満たす行列として、二項係数、フィボナッチ数、リュカ数が現れたということで、なかなか面白い考察になりました。他にも見つけたよという方はぜひ教えてください。

さて、パスカルの三角形の例を使うと、実はパスカルの三角形の任意の位置の任意の大きさの正方行列の行列式を求める公式が得られるのですが、それは気が向いたら記事にすることにします。

読んでいただきありがとうございました。

投稿日:2021220
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 概要
  2. 復習
  3. 具体例
  4. $a_{ij}=f(i)g(j)$
  5. $a_{ij}=f(i)g(j)h(i,j)$
  6. $a_{ij}=f(i)+g(j)$
  7. おわりに