テトラナッチ数列で類推する、ボナッチ系数列の一般式の一般化～その２～

前回の続きです。

テトラボナッチ数列の一般式を出すために計算していて、|P|までは計算が終わりました。

|P|Q11を計算してみましょう。


2列目から3列目を引いて2列目に代入しますと

1列目から3列目を引いて1列目に代入しますと

行列式を2次に縮めます。

1列目と2列目をそれぞれ、r2-r4とr3-r4で約分すると

1列目から2列目を引いて1列目に代入


|P|で割ると




同様に、|P|Q21も計算しますと


ここまでは同じ方法なので、よかったらやってみてください。
それから、また同じように|P|で割ります。









あとは、これのQ11～Q41に代入すればいいだけなので




こういう風になるわけです。
ちゃんとwikiとも一致していますね。つづく

テトラナッチ数列で類推する、ボナッチ系数列の一般式の一般化

こう書かれるわけですが(あ、rは固有値です)


Pというのが以下のようにあらわされ

Pの逆行列をQとすると


このように計算されますが
 
n番目のFしか求める必要がないのと、初期値がF3以外全部ゼロなことから

必要な計算は青の網掛け以外のところだけとなり

ここまで省略することができます。

また、Q11～Q41は

このうちの

こんな風になります。

まず|P|を計算してみましょう。
 
2列目から1列目を引いたものを、2列目に代入しますと

さらに、3列目から4列目を引いたものを3列目に代入しますと
 
4列目から1列目を引いたものを4列目に代入しますと

ここで、4行目のほとんどがゼロになったのでくくりだして3次の行列に縮めたいのですが
1行1列から数えて偶数番目に1があるので、縮める際にマイナスを付け忘れないでください。

1列目、2列目、3列目がそれぞれ、r2-r1、r3-r4、r4-r1で約分できるので

2列目から3列目を引いて2列目に代入しますと

1列目から3列目を引いて1列目に代入しますと

こうなるので3行目ががら空きになりました。
そこでまた行列式を2次に縮めます。
今回の1は1行1列から数えて奇数番目にあるので符号はプラスです。


ここでも1列目と2列目で、それぞれr2-r4とr3-r1で約分できるので、外に出してしまいますと

こうなるので、1列目から2列目を引いて1列目に代入しますと

あとはサラスの方法を使っても使わなくても迷わないと思いますので

いちおう、添え字の小さいほうから大きいほうを引く約束にして、符号を調整しておきましょう

|P|はこのようになりました。


とりあえず今回はこの辺で。つづく。

仮面ライダー王蛇「

サヴァ・フィナーレ！」

なんとかナッチ数列系の一般式が導けそうな気がしてきた

これな

トリボナッチの分子が分母の整数倍かどうかは任意番目について確認できる見通しは立っていないのに

こっちの見通しのほうが早いんかい

行列式の変形ギミックすごーい！たーのしー！


しかしなんだ、サーバダウンか。まじ綯エル・プサイ・コングルゥ。

仮面ライダーレーザー

生きているとき→バイクになる→上田麗奈さんがはねになって東山奈央さんが鈴菌になる
死んでいるとき→自転車になる→上田麗奈さんがひろみになって東山奈央さんがあみになる


変身して自転車になったら別に動力必要ないから、ただの無機物化の状態変化になってしまう

あるいは、回転する手足・・・？
ペダル回したらプーリーで後ろ足に動力が伝わって、手足をクルクル回しながら・・・なにそれキモい。御堂筋くんでも乗らんで。

乗るのはゴーストくらいだろう。それゴーストライカーや。耳が聞こえる



シャカリキにペダル回したあとは大好きー大好きー(歌：へご)



ああ＾～学外で同人タイムマシン開発してコミケで売りたいんじゃ＾～
そろそろ学外で部活もの流行らそうぜ
っていうかもうそろそろ大学生を主人公にしようぜよー

将来的には大学院博士課程まで卒業しないと就職は難しい時代がきて
そのあとは人工知能によって生涯学生制度が確立されるんだから
もっと大人げない大人の主人公が出てきてもいいじゃねえかよお

パンツの研究しかしない大人が、数学にしか興味がなかった人と、学外部活を通じてギルドを組むアニメやろうよ～

回りながら増えるダンジョン

完成図見て思ったのは、「昔のスクリーンセーバーみたいだなー」って。
癒し？残念ながら僕は迷路系のスクリーンセーバーに不安しか覚えなかったんだわ
続きどうなんの！？って。有限のコンピュータに無限を求めちゃいけない。
あえていうなら、再現性のある迷路生成をせいぜいお願いするんだな。


まあ、僕にとっては虐待でも、きっとこれをご褒美と思ってくれる業界もあるだろうから載せる
これは再現性あるよ。だってgifだもん

トリボナッチ数列の一般式の、分子が実数であることは確かめる

続きです



分母が実数であることはわかりました。

じゃあ分子も実数でないと、この数列は整数はおろか、実数にはなりませんよね。


そういうことで、分子も実数であることを確かめます。


分子はこちらになります。

第一項はBが実数なので、その2乗ももちろん実数。

2項目と3項目を変形してみましょう。


はい、これで、2項目と3項目の和も実数であることがわかりましたね。はい。

しかし、分母の整数倍かどうかを確かめるには程遠い様子・・・

任意のnで確かめるのはちょっと荷が重そうです。

nが3か4ぐらいでしたら、まだいけそうなんですけどねー

二項定理使って直交座標系でいくべきか、ド・モアブル使って極座標系でいくべきか
おそらく前者・・・あわよくば数学的帰納法みたいな感じで任意のnに対応できればうれしい


無理数の無理性から、約分するほかに無理数を有理数にする道はないはずですからねえ

ヘキサボナッチ数列の一般式など

よくまとめられてるなおい・・・

フィボナッチ数の類似の数列を参照すると、類推でどこまでもいけるようにできてる。


たとえばL=6ならヘキサボナッチ数列の一般式が出せる。


この6次方程式を解いた6つの複素解がa1～a6に当たるわけだけど
これを使って

こういう風に定義してやれば、たちまち、
nが0～4までのFnがゼロで、n=5になったところからポッと有限の値が出て、
それまでの数列6つの和を取る、まるで規則性のあるようなないようなヘキサボナッチ数列

を、いきなりn=40とかやっても即座に出せる。

ほんともうびっくりだよ。理論上はL=500とかでもできるし
そしたら500番目くらいまで数列がゼロで、501番目とかでいきなり有限の値が現れるんだ
そういうのが無数に作れる
狂気か！

トリボナッチ数列の一般式の分母を整理する

昨日の続きです。


トリボナッチ数列の一般式はこれでした。


この分母を整理したいと思います。

まずはBの2乗を求めましょう。

だったので、その2乗は

こうなります。

続いて、|b|^2を求めましょう。bは

このような複素数だったので、その絶対値の2乗は

こうなります。

最後に
-2BRe(b)を求めましょう。Bとbは先ほど出したのを使って

こうなります。

ではこの3つを合わせて分母にしますと

このようになることがわかりました。無理数ですが、実数ですね。
ということは、分子も実数で、nがいくつであっても、この分母の整数倍であることが予想されます。

すごーい！たーのしー！

君はトリボナッチ数列の一般式の分母のフレンズなんだね！

トリボナッチ数列の一般式を行列で求める

 
トリボナッチ数列はこのようになっていて、
フィボナッチ数列では「前の2つの合計」だったのが、「前の3つの合計」に変わっています。

これを行列で表すとこのような漸化式から一般式に変わります。



ここで、行列のべき乗が現れます。これを効率的に求める方法として、対角化がありますが、
対角化の準備段階として、固有値と固有ベクトルを求める必要があります。

ここに、以下に示すAという行列があって、Aのべき乗を求めたいとしましょう。


Aの固有値が3つバラバラで、対角成分だけが固有値のみでできているJという行列があり


行列JとAとの間に

AP=PJ

という変換法則があったとしたら 

Aのn乗は次のように簡単に求められます。


このPは、固有ベクトルである縦ベクトルを、横に並べたものとなります。


まず、固有値を求めてみましょう。

固有値はこのように求めることができます。

この3次方程式の解は、カルダノの方法で


このように求まります。
ここで、λ1は実数、λ2とλ3は複素数で、λ3はλ2と複素共役であることがわかるので
B=λ1、b=λ2とおくことにします。

また、固有ベクトルは、この行列の場合
固有値がBのとき、bのとき、b*(複素共役)のときでそれぞれ

とわかっているので、AとJの間の変換のための行列Pは

であることがわかりますし、Pの逆行列は
行列式|P|と余因子adjPを用いて

こう書けるので

行列式は

こうなり

余因子は

こうなります。


一般式は

この行列の、一番下の列だけあればよいですし、縦ベクトル(1,0,0)の最初の1以外はゼロなので

青く塗った部分の計算はいりません。

結局、一般式Tnは
 
このようになります。

複素数でなおかつ、無理数なのに、
不思議なことに、トリボナッチ数列は自然数として算出されるのです。

nが自然数ならトリボナッチ数列も自然数で合ってるのですが、nは実は負の整数まで拡張でき
トリボナッチ数列自体も、負の整数まで拡張されます。


もうちょっと綺麗にできそうな気がしますね。


テトラナッチ数列を行列の固有値から求めることもできなくはないんですが
4次方程式の解である固有値そのものが、解析的に化け物じみた狂気をはらんでいるため
テトラナッチは諦めましたｗｗｗｗ
狂気の様を見たいなら、ウルフラムαに「行列の固有値として」代入してみると見れます。
なぜかx^4-x^3-x^2-x-1=0という高次方程式にすると近似値だけが返ってくるんですよね。
3次方程式x^3-x^2-x-1=0も同様です。

僕がやりたいのは、解析的に整数に落ち着くまでの姿を見たいのです；ω；

バーチャルゲージ粒子「バレなきゃ犯罪じゃないんですよ」

バレバレだよ！お前らが犯してなかったらこんなに宇宙豊かじゃないよってくらい恩恵被ってるよ！！！

方程式メーカーベンチマーク＠scilab

やる気が出ないのでぼーっと、行列2重ループに頼らない方程式メーカーを考えてました。
2378次方程式の時点で、スタックが足りないとか言い出して計算できなくなりました。
ハードの問題なのか、それともコードが要改良なのか、知るか。

m=2577
x=poly(0,'x')
y=zeros(m+1,1);
for n=0:1:m
　y(n+1)=x^n;
end

z=y(1);
for n=1:1:m
　z=z+y(n+1);
end

z=z-2*x^m;

w=roots(z);
re=real(w);im=imag(w);
w2=[re,im];


出てきたw2をExcelのcomplex関数で複素数(文字型)c=a+ibにしてやって、imargument(c)でソートしてやれば
線のグラフにしても編み物にならずシンプルにグラフ化されますよっと

ゼーガペイン×このすば

めっちゃ混ざってた。

一緒にログアウトしたパーティーメンバーの長が覆面で、更衣室で素顔覗いたらウォーズマン。

イェルさん！？作りかけなら無理してログアウトしなくていいんですよ？

あと、ゼーガペインアルティールがドローンだった。やだ、なにこれダサカッコいい。なんかブラックタイガーみたいのもいる。え、巨大ロボットじゃなくてパワードスーツだったの？等身大のくせに走らずに滑るし、opでは謎スカイダイビングは基本。みんな滑って落ちろってか。季節感あるよな

胸はね、さすがゲーム。ブラしなくても垂れない。変なとこだけゲームなんだから。

微分方程式は生涯スポーツだったんだよ！！！ΩΩΩ＜な、なんだってー

いまだに中学生くらいのときのスキー板をつけてスキー場に行ったりします。
体型が変わらないので、ウェアもそのまんま、全身ダッサダサの
1年1回行ければ十分ってな感じで、やースキーってのは金がかかるから
年1回くらいしか行けないんですわ。


幸い、僕の場合は体が覚えててくれるのでブランクが5年程度だとまだ余裕で思い出せます。


微分方程式はそれに似てるんですね。
いつだったか、帰省中に船の中で暇だったから簡単な1F微分方程式を解いたことがありました。
ノートとペンがあれば、ほかに何も参考になるものが必要ないため、気軽にできる遊びでした。
ケータイはあったでしょうが、スマホはまだない時代、ケータイもネットにつなげる時代ではあったと思いますが、海の上ですからねえ、圏外だったかもしれません



20歳未満の記憶では有効なのかと思いきや
これがそうでもなく、紙媒体のみに許された技能でもない。

Excelでも割りと同じ作業を繰り返しやったりすることが多くて
というか部品化があまり好きではないからっていうのもあるんですが
何度も何度も同じことをするんですね。

プログラミングについても同様の癖があることが最近わかりました。
scilabに特殊ユニタリSU(3)とかSU(4)とか入れるときに
scipad(今のscinotes)を使わない癖があるんですよ


そうすると手が勝手に動きだすんですねえ


文法を覚える際には便利な癖だと思ってます。






追記
ところで、エピソード記憶についてなんですが
たとえば右に72、左に56があったとしましょう。
どっちも8の倍数ですよね。で、両方を割ったら9と7になりますね
平均すると8じゃないですか。
じゃあ72+56っていくつよ？8×8=64の2倍で128じゃね？

こうやって遊んだことを覚えると、この4桁の数字をなかなか忘れないんですわ。
実際、2軒√5軒斜め隣の家まで記憶が保てましたからね！