20080511~
13と7と11の倍数の論理積は13と7と11の積の倍数である。
和ァ・・・
パウリ行列はもう書けることがあまりなさそうな気がするので意気消沈してます。
ちょっと意欲が低下気味なので、気分転換がてら新しいテーマに着手しようかなと思いまして。 といっても朝ごはん食べたあとに昼ごはんまで寝ちゃったのが悪いんですがw イモノチベーションかむばーっく! 二度寝やっちゃうとよくないんですよねー 無駄に首は凝るし 時間はなくなるし、昼ごはんまで一向に起きないし。 肉体・精神・時間、3つの意味で余裕を失います>< せめて計算とは関係なくても運動しろっつーの!(大昼寝大会って大2つ絶対入ってるよね) 20時をすぎて、書ける日記なんかないかなんかないか!?って慌てふためいてます。 新しいテーマといっても、Excelで以前シミュレーションした波動方程式(インピーダンスマッチング境界条件とか)のあたりを、2Dから3Dにするだけです。 だけですって自分で言っておきながら、寝る前のこの時刻で適切に処理を進めるのは難しく 朝の爽やかな脳がほしいです><いちおう、意欲を出せるだけの準備はしたつもりです。 明日のー自分をー信じてー!!(信じるなよ・・・)ゲロゲロゲロ・・・|||orz にほんブログ村 間 違 え た!3Dなんて恐れ多い!!2Dto3Dじゃなくて、1Dto2Dの間違いでした!!!!1 強いて言うならgifアニメ化するから、1Dto2D空間だけど、2Dto3D時空ってことっす。 PR
指数関数に行列をぶち込む際に必要になってくる固有値・固有ベクトルや対角化の問題
ネットの海に転がってる、とあるExcelの行列アドインを追加すると、固有値・固有ベクトルを算出してくれる関数が得られるのですが、残念なことに実対称行列だけが対象なのです・・・掛け算は複素に対応してるのになぁ 非対称までとは言わずも、せめてエルミート行列に対応してほしかった~なんて欲を言っても仕方ありませんが 手計算で算出した固有値・固有ベクトルを与えることで 同じアドインで検算なら可能です。 (ソルバーがどうとかいう話もあるようですが僕はまだよくわかりません><) たとえばパウリ行列のy成分σy の固有値は±1で、 固有ベクトルは固有値+1なら 固有値-1ならなので 対角化のための行列Uは になります。 このUの逆行列をExcelで求めてみましょう。 オフィス付属の拡張アドインでも、逆行列を求める関数Minvがあると思いますが 今回は要素が複素数なので、行列アドインに入っているMinvC(エム陰部しぃ)を用います。 Mは行列、Cは複素数の頭文字、invは逆の頭3文字です。 引数が2つあります。2つ目の引数は、複素数をどう扱うのかを決めるオプションです。 このような形であれば1を このような形であれば2を このような形であれば3を選んで末尾にシャープを押してください。 ただ、3の場合は文字列として扱われるのであまり便利とは言えません。 また、1は省略することが出来ます。 インプットもアウトプットも同じ形になります。 今回求めたいのはU=の逆行列なので、 このように記述してMinvCを適当なセルに打ち込み、引数の2行4列を範囲選択し 代表を算出してから、 求めたい範囲である2行4列を範囲選択して、ctrl+shift+enterを押すと、行列全体に反映されます。 結果はこうなりました。 U^-1= つまりこういうことですね。 それではこのσyを対角化してみましょう。 行列の掛け算関数Mmult(えむマルチ)もオフィスにはついているのですが、これも複素数対応ではないので アドインの関数MmultCを使って計算します。使い方は逆行列と同じですが、当然掛け算は2項演算なので、引数が1つ増えます。 U^-1×σy×Uと、3項演算(3重積?)なので、慣れてきたらMmultCのネスト(入れ子)にしてしまいましょう。 MmultC(U-1,MmultC(σy,U)) パウリ行列がエルミートかつユニタリなので、しっかりと絶対値が1で実数の固有値が、+1から-1の順に算出されているのがわかると思います。 あとは先日の日記のとおりです。 第一段階として、パウリ行列σyを、そのままべき乗展開した級数近似を見てみましょう。 0乗(単位行列)から10乗の項までの和を取っています。 4乗から8乗までは略してあります。 左のブロックではべき乗、真ん中のブロックではべき乗の結果を階乗で割り 右のブロックでは累積させています。 角度θ=180°で、しっかりと、マイナスの単位行列-Eになっていることがわかります。(オイラーの等式) 純虚数成分と非対角成分がほぼゼロに収束してることがわかるかと思います。 第二段階として、対角化してからの級数展開の様子を示します。 固有値が±1なので、それに虚数単位と角度iθをかけて計算しています。 左・中・右のブロックは先ほどと同様、べき乗を→階乗で割り→累積 最後に下の段で対角化のためのユニタリ行列で挟んで掛け算しています。 対角化して計算しても、結果は変わらないことがわかると思います。 対角化した際に便利な計算として、第三段階の級数展開を示します。 もはやスカラーです! 級数展開が解析的に計算できてしまいますので!結果が回転行列そのものであることがわかってしまうと思います! そんじゃ任意の角度θで回してみましょうか。 第1段階 行列そのままべき乗 cosやsinが1を超えたら赤くマークしてます。 左が第1段階、右が第3段階でのグラフをモニタしてます。 どうも虚部がちょっと上方修正されて算出されてる感じでしょうか。 第2段階 対角化しつつもべき乗は数値計算する これも、cosやsinが1を超えたら赤マークしてます。 10乗の項までだともう少し収束が足りない感じでしょうか 第3段階 綺麗に円を描いていると思います。 まあよかったらDLでもして遊んでやってください。 注意1:matrixアドインは各自で入れることになると思います 注意2:回転させるのに循環参照を使っているので、動かなかったら設定しといてください。 (第2、第1段階の角度θは第3段階を参照してます:同期) 回し方を100万くらい×(now-today)に変えたこっちも用意しておきます。 循環参照でつまずいたら使ってやってください。(ただしアドインはどうにもなりません) にほんブログ村 しょうもない誇りと引き換えに、読みやすい苗字になれる数少ないチャンスを棒に振りましたことを今も時々航海するのです ちなみに独身ですし、嫁入りする気もされる気もありません。(だが男だ) L・ψ・¬ β宴ー 卍=3 シュタッ (作品が作品なだけに、なるべくカブらない苗字にしよう→法改正→わんさか!) にほんブログ村
文字数が2文字も増えました。
といっても秋か夏の話なんですけどね今年の。 まだ還元されてない濃縮されたココナッツミルクにテキトーな比率でガムシロと牛乳を混ぜていただくのが好きなんですが どういうわけか毎回飲んだ後にアイモテアーがするんですよ まあそれでも飲むんですけどね 「果実型鎧武、ドリアン?」 「そ。素材を活かして調理するとン万円はくだらないそうよ。」 「フン・・・面白い・・・」 「違う!くだらないの!」 「でもこれ生きながら調理されてますよね?」 「だから生け捕りにするのよ。いい?焦がすんじゃないわよ!」 「捕獲レベルは?」 「1円玉100万枚ってところかしらね。」 「うわ・・・この通貨システム持ち歩きづらそう・・・」 「ペンシルisバストっていうでしょ」 「すでに1つの品に2種類の値段がついててシンプルも○ンコもないわ」 にほんブログ村 後藤慎太郎、2代目。ひらがな3つで510 呼ぶときはTOEちゃん。 ====== Which is HONTAI of AMPAMMAMM?
|
カレンダー
カテゴリー
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新記事
(01/01)
(01/03)
(09/23)
(09/23)
(02/11)
(05/30)
(05/28)
(05/28)
(05/27)
(08/04)
(10/24)
(06/08)
(05/22)
(01/13)
(11/04)
最新TB
プロフィール
HN:
量子きのこ
年齢:
43
HP:
性別:
男性
誕生日:
1981/04/04
職業:
WinDOS.N臣T
趣味:
妄想・計算・測定・アニメ
自己紹介:
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます 例:H→7番目、XP→15(P)×26+23(X)=413番目。 A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいズコー
ブログ内検索
アーカイブ
最古記事
(05/11)
(05/11)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/14)
(05/14)
(05/14)
(05/14)
(05/16)
(05/16)
(05/16)
アクセス解析
|