忍者ブログ
20080511~ 13と7と11の倍数の論理積は13と7と11の積の倍数である。 和ァ・・・
[3814] [3813] [3812] [3811] [3810] [3809] [3808] [3807] [3806] [3805] [3804]
a1=ones(1,500)
a2=eye(500,500)←2重ループなので多少時間かかります
a3=zeros(500,1)
A=[a1,1;a2,a3]←2重
x=poly(0,'x')
q=det(x*eye(A)-A)←ここで一番時間食いそうです
p=roots(q)
re=real(p)
im=imag(p)

これを慣れたExcelに張り付けて、
re
im
complex(re,im)
imargument←これでソートして、re対imでグラフ化


さすがにちょっと時間かかりましたね。パソコンが

ベンチマークとしていかがでしょうwwwwww


もうね、解が、ほぼほぼ線です。吹き出しみたくなってます

1001次は終わる気しなかったので止めました。


実数だけども、非対称行列だとspecとbdiagの結果が食い違いますねえ。
specのほうがもっともらしく見えるけど、
もしかしてどっちも合ってるなんてことあるかな??
今度、実際に対角化して確かめてみましょうか。

拍手[0回]

PR

コメント


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


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