忍者ブログ
20080511~ 13と7と11の倍数の論理積は13と7と11の積の倍数である。 和ァ・・・
[3774] [3773] [3772] [3771] [3770] [3769] [3768] [3767] [3766] [3765] [3764]
 先日までの日記参照



これを規格化ってのはなかなか死ぬる・・・。

まあ、各変数を一様乱数使って割り当ててるから、めったに固有値がカブったり
固有ベクトルの一部が本質的にゼロになったりはしないので
ジョルダン標準形やグラムシュミットの正規直交化法を使う必要はほとんどないんだけども
それでもかなりめんどい。


なにか近道はできないものかと考えて、ちょっと思いついたのは
この行列Pの「逆行列がPのエルミート行列そのもの」であることを使うこと。

つまりEを3次の単位行列として、PP†=Eのことなんだけども

この内訳を考えると

こういうことなので総和記号を用いて、以下のように、l行目k列目の要素は

といえる。
また、l=kのときは1、それ以外はゼロなので、クロネッカーのデルタを用いてこのようなこともいうことができるし

複素内積(u,v)を

このように定義すると

たったこれだけで記述できてしまう。




ここで、今回は行列指数関数を求めたいんだった。

間にこのような、対角行列の指数関数が入ってしまうので、さっきのように簡単にはいかないが、楽をするヒントは得られそうだ。

つまり行列指数関数expの各要素elkは

とできそうなので、
少し計算を省けそうなのだ。

 
この2つの式を用いると

l=kの場合は

なのだし、vl3はlが1から3まで実数なのだから
|vl2|^2が逆算できて、

少~しだけ手間が省けるし、δ=0になる非対角要素も


これを逆算して

さっきよりももうちょっとだけ大掛かりに手間が省けそうだ。

ところで、Σ(λm)=0なのだが、3つではなく実質2つのλで構成された行列指数関数になりはしないかと若干心配だったりする。

拍手[0回]

PR

コメント


コメントフォーム
お名前
タイトル
文字色
メールアドレス
URL
コメント
パスワード
  Vodafone絵文字 i-mode絵文字 Ezweb絵文字


忍者ブログ [PR]
カレンダー
10 2024/11 12
S M T W T F S
1 2
3 4 5 6 7 8 9
10 11 12 13 14 15 16
17 18 19 20 21 22 23
24 25 26 27 28 29 30
ブログランキング
ブログランキング参戦中
にほんブログ村 アニメブログ 深夜アニメへ
にほんブログ村 漫画ブログ SF・ファンタジー漫画へ
にほんブログ村 科学ブログ 自然科学へ
よかったらポチッとお願いします^^
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新TB
プロフィール
HN:
量子きのこ
年齢:
43
性別:
男性
誕生日:
1981/04/04
職業:
WinDOS.N臣T
趣味:
妄想・計算・測定・アニメ
自己紹介:
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます
例:H→7番目、XP→15(P)×26+23(X)=413番目。
A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいズコー
バーコード
ブログ内検索
アクセス解析