お子さんが他の何でもなくレゴでまず建築物を設計してしまうもどかしさ

 
すごく・・・よくわかります！
やっぱり慣れたツールが一番ですもんね！！！！１




ちょっと疲れてて目前のやりたいことを処理できませんトゥデイorz



なお、ちょんまげの起源は分からない模様

にほんブログ村
ｂｙミンメイ書房



俺達の戦いは、始まったばかりなんだ！！！(わりとガチで)
暇人大国ニッポン・・・日本からニートが消えたら日本が滅ぶって本当ですか・・・

今日は数式をちょっとお休みしましょうか

 
ここ一週間ぶっ続けで数式三昧でしたからね。読む方もどうもお疲れ様でした。

フィボナッチ数列で書きたかったことがひと通り終わったので、
次何書こうかなぁって思ってたところだったんです。
ちょっと気が抜けました。
ということで気を抜きましょう。


やりたいこと、というか
やりかけでまだ終わってないテーマが2つほどあるので
それに手をつけちゃえばいいのですが
ちょっと英気を養いたいところですね。



僕の英気の養い方は今のところ
「疲れを癒やす」系と「意欲を高める」系の2種類がありまして
「意欲を高める」系にはさらに、「意欲を生み出す」系と「意欲を維持する」系
があるみたいなんですね。


といっても最近気づいたばかりなのでまだ全然形式化されてなくて
色々模索中というか、まだこれしかわかってないのが現状です。



「温泉やプールに入る」という行為も、「疲れを癒やす」系と「意欲を生み出す」系両方にカテゴライズされますし。
まあ本来この手のものはカテゴリがあってないようなものなのかもしれません。
高速ドライブや山登りドライブもいいですね。あと特に行き場とかない自転車とか。


このブログは完全に趣味でやっているものなので収入も支出もないんですが
仕事は別にあって、これがまた趣味も仕事も両方とも肩が凝るんですよ


なので温泉に入るのは最近だとかなりトップにランクインする癒やし効果があります
温泉に入ると、同時にアイデアが沸いてきたり考えが整理されたりもするんですよ。

たとえば波動方程式の境界条件について悩んでいる最中だと、
温泉の水面は水の波に光の波が透過したり反射したり屈折したり干渉したりと
世の中すべて波だらけの状態に見えてきまして
結構楽しくアイデアが湧いてくれるんですよね



意欲を生み出すまではいかなくとも、維持するのには気分転換で漫画読んだりアニメ見たりすると
途中にしてたブログを、じゃあ書き終えようかって気分になります。あんまりケチるもんでもないですねえ



他の方法はなにかないだろうかって思ってるところなんですが
「別の計算を始める」ってのも手ですね
最近Excelで何か始めるのが面倒くさくなってきて困るんですが
始めるとこれ以上落ち着く作業はないってくらいしっくり来るんですねえ


一方で、「食事」に関してはあまり期待できなくなってきました。
逆流性食道炎なんですよ自分・・・
食べたいのに食べたい量食べれないのを目の当たりにするとちょっとキツいですね・・・
まあおかげで脱メタボどころかそのうち質量消えるんじゃないかって思うくらいコツコツ減ってますけど
元々食事にあまり重きをおいてなくてよかった＾＾
といっても、噛むこととかを疎かにしたからこそ逆流性食道炎になったのかもしれませんが。



ここ1週間はホントに、どこの誰だか知らないけれど
一定の高評価を頂いたみたいで嬉しいです。
ありがとうございます。
久々にあっという間でしたよ～この1週間。




それと
フーコーさんおたんじょうびおめでとう！
赤道だと余裕でポールまったく倒れないぜ！

にほんブログ村

おまけ：規格化した固有ベクトルとユニタリとエルミートと、回転行列だったら何度？

 の随伴行列 は規格化してみると面白いですねえ
(たぶん)解除できない二重根号がうっとおしいですけどね・・・

綺麗な形してるだろぉ～？コレ実は｜右上｜＝｜左下｜、｜左上｜＝｜右下｜なんだぜ？


ユニタリ： 
エルミート： 
当たり前っちゃ当たり前なんですが、実数範囲だとどうあがいてもユニタリかエルミートｗｗｗ
特にユニタリだと回転行列なんですよね。
これも当たり前、ユニタリの特性ですが、逆行列が転置行列なんですね。
（3行3列以上は知りませんけど＾＾；狭い！）

  
回転の角度0.5536[rad]=31.7[°]が何を意味するのかわかりませんがｗｗｗ
なんか黄金率とは別のアスペクト比が出てきましたよ


と、そういえば最後に告知：重大な勘違い表記をしていたかもしれません。
先日の、「固有ベクトルを並べた行列P」は随伴行列とは呼ばれないみたいです・・・
 
なんて呼ぶのでしょうか・・・ほら、規格化されてないPとか・・・
規格化されてればユニタリですけどねえ




なんか15年前のことを思い出します
ドヤ顔で提出した制御工学とかのテストがまるっきり0点とかで、なんでや！ヽ｀д´ノって思った思い出・・・

にほんブログ村

対角化の威力でフィボナッチ数列の一般式に迫るｩ～

こんばんは、対角化シリーズ最終回、のつもりです。

初回は、行列表現したフィボナッチ数列を例に、行列の対角化がどんなものなのかをざっくり説明しました。行列同士の掛け算が簡易に行える裏技だと書きました。

つづいて前回は、対角化のための準備である固有値と固有ベクトルについて、フィボナッチ数列の行列表現を例に具体的に計算してみました。

今回は、その対角化を使って、フィボナッチ数列の一般式を導出して締めることにします。


フィボナッチ数列を行列表現した要となる行列A=
を対角化すると、固有値φ1=(1+√5)/2とφ2=(1√-5)/2、そして固有ベクトルからなる随伴行列Pを

を用いて行列Aを対角化すると


こうなることはわかいましたが、だから何なの？と思う方もいるかもしれません。
僕は思ったんですけどね。＾；

この式の両辺に右からP-1を、左からPをかけるとどうなるでしょうか。
Pとその逆行列が約分？されて


このようになることがわかるかと思います。

さて、フィボナッチ数列そのものは、

このように書くことができました。

あとの都合のために、n+1をnに、nをn-1に置き換えて

こうしておきましょう。
この、 が行列Aなので、対角化したものを代入すると


このようになるはずです。

どうも、随伴行列は自分自身の逆行列と打ち消されてべき乗しても形が変わらないようなので

このようにn-1乗は対角化した行列のみにかかってくるようです。

また、対角化していますので行列といえどもべき乗は簡単で、そのまま対角成分のべき乗と書くことができます。

 

これが対角化の威力です。

つまり、フィボナッチ数列の行列表現は

これに帰着でき、もはや手計算で行えてしまうのです。


最終的に行列の上のほうだけわかればいいので、赤で囲った部分にだけ注目し
以下の様な一般式に無事、終着するわけです。＾＾



お疲れ様でしたｍｍ



合体ロボ同人サイト

にほんブログ村

固有値・固有ベクトル・対角化。

前回は、フィボナッチ数列の行列表記を例に、対角化の利用法を見てきました。
対角化は、行列を掛け算する際の裏技のようなもの、といいました。

固有値と固有ベクトルは、その対角化を行うための準備とも書きました。

それではフィボナッチ数列の行列を例に、具体的に固有値と固有ベクトルを計算してみましょう。

まず、固有値の定義はこうです。

Aという行列があったとき
|A-λE|=0となるλが固有値


Eは単位行列で、中身はコレ です。数字でいうところの「1」のようなものです。
何回かけてもEのまま、逆数に相当する逆行列もEそのものです。EじゃなくてIと書いたりもします。

ではA-λEの両脇の｜｜は何かと言いますと、これは囲まれた中身の行列式を計算するという記号で
行列式自体はスカラーです。デターミナントと読んで、det(A-λE)とも書きます。

行列式の説明は少し後回しにして、A-λEを実際に計算してみますと

こうなります。

この行列式を求めるのですが、2行×2列の行列式は割りと単純で

右下斜め掛け算－左下斜め掛け算

というスカラー量、つまりただの数になります。


これが=0になるλを求めるわけですから、あとは2次方程式です。

このφ1というのは黄金率です。φ2はオマケのようなもので、2次方程式の解と係数の関係から
φ1φ2=-1(方程式3項目の係数)とか、φ1+φ2=1(方程式2項目係数の符号反転)とかいう関係があります。


固有値に関しては以上です。



＝＝＝＝＝＝＝＝＝＝
次に固有ベクトルの話に入ります。
固有値λ1とλ2を実際にA-λEに入れてみて

となるようなxとyの比を求めます。

この際、連立方程式のようになりますが、この2つは同じ式で、いわゆる永年方程式になるので、具体的にxがいくらいくら、yがいくらいくらとは求まらず、比しか求まりません。
というか、固有値を求める作業自体が永年方程式を作ってるようなものなのです。

λ=φ1を代入してみますと、以下のような行列方程式になり

このような等価な2つの連立方程式(永年方程式)になるので(解と係数の関係を思い出して確認してみてください)

xとyからなる固有ベクトル(列(縦)ベクトル)は１：－φ2から上のようになります。


また、固有値にλ=φ2を代入した際の固有ベクトルは以下のようになります。式変形はほとんど同じなので略します。


この2つの固有ベクトルを横に並べたものを随伴行列と呼んでPという記号で表します。

 

固有ベクトルを求める際に、xとyの比しか求まらないと言いました。
ということは片方の固有ベクトルのxとyの符号を反転して
 ①
こうしたり


あるいは両方とも符号を反転して、こうしても問題はありません。

ですが筋は最後まで通してください。
これらの随伴行列は左側にφ2がきていますが、固有値はφ1が最初(大きい方)で次がφ2(小さい方)
という約束の元で行っています。
それと同様に、勝手にルールを決められるからこそ、同じ約束事で筋を通さないと間違った結果になってしまうのです。

また、これらのベクトルはノルム(ベクトルの向きを無視した大きさや絶対値のようなもの)を1にしていません。(ノルムを1にするとルートの中にルートが入ってしまって面倒なのです＞＜)
このノルムを1にする＝規格化という作業もするかしないかも自由ですが、筋は通してください。


このあと出てくるPの逆行列P-1のところで重要になってくるのです。



＝＝＝＝＝＝＝
今後、印をつけた①にしたがって進めていきます。(都合上、左上プラス、右上がマイナスのほうがいいのです)
①のPの逆行列P-1を求めるのです。

逆行列の求め方はこれを読んでいる方はもうご存知かと思いますが
2行2列の場合は
 
・対角成分(右下がり)は入れ替え、
・そうでないもの(左下がり)は符号反転して
・全体を行列式(デターミナント)で割ります。
 

このPの逆行列を求める際に、①で決めたルールなのかなどというのを守っているのと守っていないのとで結果が違ってくるのです。



＝＝＝＝＝＝＝＝
しまいに
随伴行列のPとP-1が出来上がったら、P-1・A・Pの掛け算をしてみてください。

  
計算結果が2つの固有値で対角化されているのが確認できればOKです。
左上から右下にかけて、φ1(大きい方)→φ2の順番になっていることも確認して下さい。逆はダメです。
※PAP-1ではなくP-1APです。これも順番に気をつけてください。行列の掛け算は順番を気にします。
スカラーは順序気にしません。(√5など)



 次回はいよいよ、フィボナッチ数列の一般式を導出して、対角化の有用性を話してみたいと思います。



まとめ。たのしい！

にほんブログ村
嬉しいです！<´皿｀>

フィボナッチ数列を、フィボナッチ行列として、捕まえる。

フィボナッチ数列はよく、「漸化式は単純なのに一般式を探そうとすると複雑」な、いい例として使われます。

しかしながら行列を用いれば一般式っぽく書くことが可能です。
これはただのこじつけではなく
行列という数学の道具を用いることで、スカラー(ただの数)での一般式を導き出す手順にもつながります。


n番目とn-1番目のフィボナッチ数列をFn、Fn-1とし、列ベクトルで書くと 
になりますが、1つ次のフィボナッチ数列にするには、以下のようなルールを持たせた行列を掛け算することで変換が可能です。
 

また、これを拡張すると、最初のフィボナッチ数列つまり1番目(上)と0番目(下)のセットからn+1番目(上)とn番目(下)のセットを作り出したいときは行列A＝ をn回掛け算(n乗)すればよいことがわかります。
つまりこうです。
 


しかし、行列の掛け算は手計算でやるにはやや複雑です。
そのため、掛け算を簡易化するために固有値や固有ベクトル、対角化という行列の道具を使うと便利です。


固有値や固有ベクトルというのはなんとも説明しがたい概念なのですが
とりあえず今回は対角化(掛け算の簡略化)のための準備と認識してもらえばOKかと思います。

対角化すると何が便利かといいますと
Aという行列があって、その正体が実はB＝ という対角成分だけの行列だったらn乗するとき
 
になってさぞかし楽だろうなーという希望的観測に基づいた理論といっても過言ではないかもしれません。


これが実はできてしまうのです。

このようになる行列Pとスカラーa、bが割りとしょっちゅう存在するのです。
ここでいうaやbというスカラー(数)が固有値、行列Pを構成するベクトルが固有ベクトルと呼ばれているものです。
フィボナッチ数列の場合は2行2列の行列を扱いますが
一般にn行×n列の行列にはn個の固有値と、n列の固有ベクトルがnセットあることが多いです。


 次回は具体的に行列の固有値と固有ベクトルを計算してみます。


にほんブログ村
ｼｬｴｴｰ

バインバインふぃぼなっち(一般式の威力)

 
一番上、左辺が2つの黄金率(x^2-x-1=0の解2つ)なんですけど
ルート(「平方」根)だけが入った多項式なので、2乗しても2項で済むんですね
そうやって2乗をさらに2乗・・・としていくと128乗、つまり128番目のフィボナッチ数列の種までがたった6回の手計算で行えちゃいまして

しかも2乗する際のルールが、
 
・1項目：元の2項目の2乗に5をかけて1項目の2乗に足して2で割る
・2項目：1項目×2項目

 
と、至極単純であり、
この「種」からフィボナッチ数列を取り出す際も、2項目の分数やら√5やらをバッサバッサちぎっては取っ払った整数がそのまんまフィボナッチ数列というなんとも単純明快なものなので

0番目：0

1番目：1
2番目：1
4番目：3
8番目：21
16番目：987
32番目：2,178,309
64番目：10,610,209,857,723
128番目：251,728,825,683,549,488,150,424,261

と、あっさり計算できてしまうわけです。

最後のあたりでExcelの計算限界に達したので電卓の力を借りました＾＾
電卓の限界がいまいちよくわかりませんｗ

にほんブログ村

PC電卓がExcelを超える日

Excel：ふぃぼなちってない　　電卓：ちゃんとふぃぼなちってる 

超える日っちゅうか一時的に超える場合があるよってだけですが
ところどころExcelは電卓の上位互換だったり下位互換だったりって当たり前なんですけどね。
120番目くらいまでいけるみたいです。
5,358,359,254,990,966,640,871,840

ちなみに、このくらい数が大きくなると
黄金率：φ＝(1+√5)/2として
Fn=[φ^n-(-1/φ)^n]/√5の2項目がほとんど意味をなさなくなって
Fn≒φ^n/√5になるんですよね。電卓で打つのが楽になりますｗ


にほんブログ村

Excelの限界～ver.フィボナッチ数列～

 
囲った部分、途中からフィボナッチになってないの図

Excelの限界を示したつもりが、逆に浮動小数点の化け物じみた正確性を示すことになってしまったでござる
フィボナッチ数列の漸化式と一般式、間違えるところも、間違え方もまったく一緒ｗこりゃすげえ

にほんブログ村

エクセルの限界～フィボナッチ数列版～

  
これ以上桁が大きいと精度を維持できないとかいう・・・
πでやったけど、有効数字15桁同士の浮動小数点を重ねるシステムの構築とかもうやだ(´・ω・｀)
当時は非オブジェクティブだったけど、仮にオブジェクティブにできるようになったとしてプログラミングそのものがｲﾔﾀﾞ (´・ω・｀)
早く恐怖症を治したい・・・かも？

にほんブログ村

フィボナッチ数列と母関数と行列の固有値・固有ベクトル・対角化

大学付属の図書館に行ってきました。
「数学ガール」が５巻揃っているのです。貸し出されずに。
貸出手続きは面倒臭いのでその場で読みます。ノート持参です。一人ぼっちでしたー。

無印の数学ガールです。コミカライズ版が手元にあるのでだいたいの流れは知ってるんですが
やっぱり数式部分はコミックだけでは弱いので、いつか読もうと思っていました。
コミック版と違って高いんですよねー


やっぱり無印だけあって、最初の数十ページはちょい退屈でしたが
フィボナッチ数列のくだりは見たことない方法(コミック版でちょい見ましたが)で一般項を導出してました。


一旦「母関数」とかいうのに世界を移して、それから元の世界に戻って導出終わり。
という流れなんですが、ようやく理解はしましたが、いまいち納得がいかないというか・・・寄り道した必然的な意味がよくわからないというか・・・

本書でも「回りくどいやり方」とは書いているのですが本当に回りくどいです。どう見てもありがとうございましたｗｗｗｗ
でもこれ、なんか「数列の一般式」に必須の知識みたいで・・・なんなんですかもう。。。。



僕の知ってる方法、といっても本当に知っているのか疑問ですが
知ってる方法と違うんですよね。
コレジャナイロボですよ。


僕の知ってる方法は行列を使った方法なんです。
でも最後までつながってなかったので、いい機会なのでネットでそれ系のサイトを見てきました。

実は「対角化」の利点がわからなかったのです。
意味は知ってました。でも「だからなんなのよ」だったのです。

Aっちゅう行列があってそいつを対角化したいとき、
Aの「固有値」を求めるんですね
それで、固有値から「固有ベクトル」を求めてそいつを規格化すると、随伴行列でしたっけ？
みたいなのができるんです。この行列の性質がユニタリなんすわ。「大きさが１」みたいな感じです。これをPと呼ぶことにすると
Pの逆行列はPの転置行列そのものなんですね。あ、いや複素だとエルミート共役でしたっけ？今は実数なんで転置でいいですけど。

で、PAP^-1だかtPAPだかすると、見事に固有値だけが対角成分に現れた対角行列が出来上がるんですよ。


そこで僕は止まっていたのです。「だから何？」状態だったのです。


でもこの真骨頂はtPAP=[対角化]じゃなくてP[対角化]P^1=A
こちらの方にあったんですわ。まさにコペルニクス逆回転。


何に使えるかって言うと、Aのn乗がスンナリ出せるんです。
A^n=P×[対角化]^n×P^-1なんで、ただ単に行列の中身の固有値をn乗するだけでAのn乗が求まるんですよ。まさに手計算レベルの簡易さっすよ。


これの一連の流れにフィボナッチ数列のシステムをぶち込むと、
黄金率その１の何乗引く黄金率その２の何乗
とかいう一般項の式が求められるんですわ。

[φ^n-(-1/φ)^n]/[φ-(-1/φ)]
でしたっけ。φ：黄金率その１。x^2-x-1=0の解の大きい方、プラスのほうです。




ところで、ミルカさんがテトラさんを蹴飛ばしたくだりの「ぼく」のフォローは本当になにもないんですね(笑)らしいっちゃらしいっちゅうか



角な機体はしないでくださいねー

にほんブログ村

だいすう！

京アニのfree見て思ったんですが、なにかと五人組のアニメが成功しやすい感じなんですかね？
男だらけでも女だらけでも、男女混合でも。

そしてその五人組は、４人と最終兵器１人という構成・・・

リンというサメみたいな歯の人はまるであずにゃんポジじゃないですか。



そういえば以前、けいおん！のwikipediaのパロで、けいさん！ってのをやったことがあったんですが
別に人を人に対応させる必要ないじゃないですか。
我々には擬人化という強力な武器がある。



じゃあもう、数学の概念を擬人化しちゃっていいんじゃないでしょうか！？


遅れてきた５人目といえば・・・そうだ代数方程式！


そのまんま１次式、２次式、３次式、４次式の４人組がいて
最終兵器が５次式なんですよ！こいつが気難しいやつでしてね！なかなか解けない。


といっても僕にはガロア理論がいまいちわかっていないので、もしわからずじまいでしたら
誰か書いてくれると嬉しいなぁ（チラッ




にほんブログ村

同じ数を見てる

①ωは1の3乗根の1ではない2つのうちどちらか1つである。

②ωは1の3乗根の実数ではない2つのうちどちらか1つである。


③ω=(-1+i√3)/2
または
ω =(-1-i√3)/2 


④ω=-cos(π/3)+isin(π/3)
または
ω=-cos(π/3)-isin(π/3) 


⑤ω=cos(2π/3)+isin(2π/3)
または
ω=cos(2π/3)-isin(2π/3) 


⑥ω=cos(2π/3)+isin(2π/3)
または
ω=cos(-2π/3)+isin(-2π/3) 



⑦ω=exp(i2π/3)
または
ω=exp(-i2π/3) 


⑧ω=exp(i2π/3)
またはその複素共役

にほんブログ村

あずまんがサイクラ

春日「電波受信したらfaxできるやんかー」

よみ「まあそうだけど」

春日「じゃあ重力波受信したら」

よみ「なんでや！テレポーテーション関係ないやろ！いいかげんにしろ！」





とも「それでさーこの電動コケシの充電池がさーまたコケシみたいで」

春日「おおー、なんかマトリョーシカみたいや」

ちよ「ねじまきの動力を充電池で蓄えるんですか・・・？」

とも「うん。こいつ、電磁波浴びても充電できるんよ。太陽電池ってやつ？電子と陽電ホールが分極してね」

よみ「おい今物質中でモノポールくらいのことサラッと言ったぞ！？お前それ何eVの電磁波なんだ！？テラじゃなかろうな！？」

にほんブログ村

mod(○,9)

78を9で割ったあまりとか考えてる時に
7と8を足して一旦15にしてから1と5足して6にすんじゃなくて
7から1引いて6にすると、なんかやたら落ち着くんだよなぁ。

無駄に大きい数とかあんま好きじゃないんだ
数値は好きだけど番号は嫌いなほうで
桁が上がれば数値ってより番号とみなしたくなるんだよ


え、それ数値じゃなくて番号ですよね？社会科ですか？あーダメダメ、暗記科目アウトなんすわー
って、かんじ。


物理学者さんたちがどうして統一理論を求めてるのか知ってるかい
彼らはあれで暗記がすごく苦手なんだよ
だから、必要最小限の暗記しかしなくていいように、なるべくコンパクトな自然界の記述方法を模索してるんだ
PCでいうところの、CPU(と内部記憶装置)だけスペック高い人って感じだね
メモリは外部記憶装置に入りますかって？そこはまあ・・・察してくれ。

にほんブログ村