|
20080511～
13と7と11の倍数の論理積は13と7と11の積の倍数である。
和ァ・・・
これ↓ができあがるまで組み上げてみましょう。
まずは増幅器(定数倍)から。 数値計算のところのGAIN_fを貼ります。ほかにも色々GAINがあるのですが、とりあえず_fでいいみたいです。 次に、ステップ関数を貼ります。 「信号源」から選んでやります。 次に、時間シミュレーションのための時計を貼り付けます。これも信号源からです。 赤矢印と黒矢印の時計がありますが、赤矢印のほうを使うようです。 クロック_cと書いてあるやつのようです。 信号を混ぜるための、アナログ加算器を貼ります。入力が3つあるようですが、ここでは2つしか使いません。数値計算のパレットにあります。 積分器を貼ります。「連続時間システム」のところにあります。 何の略かはわかりませんが、_mとついた積分マークを貼ります。 観測者を用意します。 「出力・表示」のcスコープとかいうやつですね。 赤矢印1つ、黒矢印1つのやつです。 増幅器に、逆に向いてもらいたいので、増幅器右クリックで、「反転」させます。 黒矢印・赤矢印同士をつなげて、システムを構築します。 線の途中から線を引こうとすると、普通に線が分岐します。 斜めになっても、つながった時点で大概まっすぐに直してくれます。 増幅器マークをダブルクリックして、のGAINをマイナス1.25倍にします(負帰還) システムができあがったので、細かい設定に入ります。 まずはステップ関数の設定をしましょう。 次に、スコープの設定をしましょう。 最小マイナス２、最大２になるようなグラフにしたいので、このように記載します。 また、先ほどは５までの時間でステップが発生する仕組みでしたので、その2倍の１０までの時間スパンで見てみることにします。 最後に、「シミュレーション→設定」から 積分終了時間を設定して、実行すれば「➡」(こんな感じの再生ボタンがあるはずです) 以下の過渡応答ようなグラフが現れるはずです。 僕も習いたてなのでわからないことだらけですが どの部品がどのジャンルのパレットにあるのか、ちょっと紛らわしい気がしますね。 あと、貼るべき部品もちょっと紛らわしいです。 積分や時計、スコープや増幅器の部品がいくつかあるので、どれを使ったらいいのか迷いそうです。
PR ステップ関数を積分した出力を、1.25倍して負帰還したときの応答の例のようです どこぞの過渡応用のような時定数が目立つ曲線になってました。 これをオペアンプとかでやりたいんだけどなあ 困ったことに、直流電池と抵抗だけのシンプルな回路においても何かがうまくいきません 陰解法がどうとかいいやがります。 電流計を設置しようとしたらエラーが出るし、いまんとこわけわからんです (本が10年前のもののようで、当時の「scicos」記載だから、なんか不安なんすよね・・・)
この状態遷移図に相当する有向重みなしグラフの隣接行列は以下のようになります。 つまりこういうことです。 この2乗である は何を意味するのかというと このように、2フレーズ連続で歌う歌い方がも6通りしかないことを意味しています。(行列の中に1が6つしかない) 3乗や4乗、5乗も同様で３，４，５フレーズ連続で歌う歌い方も6通りしかないことを意味しています。 そして、6乗は単位行列です。 これが意味するのは 起点のことなる自分自身への一周が、6パターンだけあるよということです。 無向グラフや、重みつきになると、もう少し複雑になると思います。 (重みつきはまだ計算したことがありません。我々は習いたてなので)
トレースがゼロなので、うっかり固有値を求めてみたくなりました。 まず行列式は、 このようなマトリョーシカ展開の繰り返しなので、－１になります。 この トレース＝ゼロ、デターミナント＝マイナス１をヒントにしておきましょう。 これも、同じような掃き出し法の繰り返しで、マトリョーシカ展開ができ、 簡単に解ける6次方程式になりました。 ガウス平面の上下左右対称、ど真ん中の六角形の状態に固有値が並んでいるわけです。 つまり、非負行列に対応するペロンフロベニウスの定理は＞(大なり)ではなく≧(大なりイコール)であって 横長・右寄りだけでなく、上下左右対称のど真ん中も含んでいたということになります。 では、固有ベクトルはどうなるかといいますと おそらくこうなるはずです。 規格化もしてやると、ユニタリ行列の性質「逆行列がエルミート共役」が使えて便利でしょう。 しかしながら、このような煩雑な固有値・固有ベクトルを求める以前に この行列は6乗すると単位行列に戻るのです。A^6=A^0=Eということです。 つまり、固有値・固有ベクトル、これらは少しも大事ではなかったのです！！！！！１
こいつの特性方程式が6次方程式なのも忘れて無我夢中で解いてしまった。 λ^6=1になったからいいものの、ガロア理論の範疇超えてたの忘れてたから 一歩間違えたら死んでいたかもしれない かといってまったくあてがなかったわけでもなく トレースがゼロで行列式がマイナス１になるから、 つい、なんとかなると思ってしまったんだ。 いや、なんとかなると思ったかどうかも定かではない。 固有値がカブるかどうか気になって仕方なかったんだ λ=exp(jnπ/3) やぁ、ゴールデンウィークははしゃいでしまうね。
重みなし有向グラフ理論 隣接行列 それが大事 アルゴリズム 状態遷移図 行列のできる力学相談所 量子力学？ ターミネーター不在につき、インピーダンスマッチングできませーん！
図1のような、向きが決まっていて重みが平等な五角形のグラフがあるとします。 ある点が「反時計回りにだけ」隣に移動するための行列Aなので、Aは図1のように定義されます。 このAのべき乗は何を意味するのかといいますと たとえば3乗だったら「3つの辺で成り立っているパーツはどこにいくつあるのか」 といった情報を教えてくれます。(図3) 向きがない、すなわち矢印が両方の向きについた「無向グラフ」や、 「重みを考慮したグラフ」だともう少し複雑になるのですが このようにモデルを非常に単純化することで、 モデルと行列の意味をリンクさせてシンプルに考えることができます。 興味深いのは、5乗つまり1周すると5次の単位行列になって、元に戻ることです。 この組み合わせは、自分自身のところに戻る、たった1パターンしかないことがわかります。(図5) また、この行列には部分的に単位行列やゼロ行列が含まれているため 適切に分割することで、要素数が大きな素数でも、 2次行列に帰着させて手計算で計算が行えるのがメリットです。
歩道がいくつであっても、1通りしかない 当然ながら、4乗で単位行列にリセットされる。 ちなみに{}の記号は反交換関係で、 4次行列の掛け算がめんどいので、ディラック行列みたいに2次ずつに分けた。 要素数が大きな素数だったりすると、行列の次数も連動して大きな素数になるだろうし こういう風な分け方ができるとすごくありがたいんだけども (とはいっても実数行列限定になるので、Excelとかに計算させれば全然負担はないんだけども あーでも、固有ベクトルがな、複素数になりえるんだよな まあそうなったらscilabやwolframαの力を借りればいいだけの話かもしれない) それにしても歯が痛くてやってられん
状態遷移図にできるなら、グラフ理論にしたっていいじゃない
0から1、1から宇宙の果てまで～照らし出す光は
回路ショートくん 鬼畜こけし ほら君も目をつぶって太陽拳～ 昨日あたりから頭痛がひどいので、とりあえずできた素材を時々ばらまく日々が続いてすみません リニア方眼紙がペンローズダイヤグラムの方眼紙に遷移する図です。 とりあえず今は光中心版ができたばかりなので、ミンコフスキー時空に入れる際は 左右どちらかに斜め45度傾けてイメージしてください ちなみに僕、とても不器用なので、光の振動数がGHzオーダーの回路相手だと 検査するだけで回路壊します。 例： 抵抗値を測る ↓ 測ったはいいが不安になってもう一度測る(強迫性障害) ↓ ショートしてはいけない部分をショートする ↓ テスタの内部電源が測定対象の許容電圧を超えている ↓ きな臭い ↓ うっすらと視界が濁る ↓ 測定に再現性がなくなったので上司に質問 ↓ 怒られる シノの目、誰の目？
ヘキサボナッチ数列自体は無理数にも複素数にも拡張されてくれます 6つの白丸を全部足すと黒丸になります。格子同士の間でも、複素数でもなります。 相変わらず青丸が実整数の範疇に収まっているのがわかるかと思います  にほんブログ村
マイナス番目のフィボナッチ数列は
定義 F(n+2)=F(n+1)+F(n) を移項して F(n)=F(n+2)-F(n+1) 整数mを定義して、整数nにm-2を代入すれば F(m-2)=F(m)-F(m-1) これで済むんだけど じゃあ実数番目のフィボナッチ数列はどうすんのってなると・・・積分・・・！？ でもそんな記事は見当たらない 単に離散数nを連続数xに置き換えて F(x+2)=F(x+1)+F(x) にすればいいだけ・・・？ それでいいなら任意のボナッチに拡張はすぐなんだけども・・・ 複素数だとしても、実数は実数同士、純虚数は純虚数同士で演算すればなんの問題もないかもね・・・
さっきの日記見直してて
フィボナッチの虚部はゼロだよなぁ～って思いながらぐぐってたら、 ゼロじゃなかった！ えええええー！？ だったら素数かどうかってのが根底から揺らぐじゃん！！！！ ってことで、おとなしくscilabに計算させました。 for L=1:1:20 A=[ones(1,L);eye(L-1,L-1),zeros(L-1,1)] for n=1:1:200 B=A^((n-100)/10); m(n,L)=imag(B(L,1)); end end これはっつけたら全自動でやってくれます・・・ 奇数ボナッチが多少値が出てたのはやっぱり計算誤差だったっぽいですすみませんorz はぁー・・・素数関係ないのか・・・ちょっとがっかり。
ちなみに偶数ボナッチおよびフィボナッチ数列の虚部は、正負実数番目に拡張しても全部ゼロらしいです。
小さい数列(テトラ)ほど、大きな数列の値から、急峻な時定数で収束していることがわかる。 もちろん周期は整数の2だ。 exp(-)*sin の形だとは思うけども。 偶数の中で2だけがのけものなのがなんか引っかかる。まるで素数じゃないか 奇数ボナッチの正負実数番目の虚部はホントにゼロでいいんだろうか 小さな値しか出てないのもあれば、変に膨張してるのもある。ヘプタボナッチとか・・・ でもヘプタボナッチは整数番目でも虚部がゼロにならないのがおかしい。 n+=0.1なんて雑なことをしてるからいかんのだろうか 整数型m++しつつn=m/10とかせにゃならんかなあ だいたい、「2が素数だから虚部がゼロ」って論理だと、素数7のヘプタボナッチでおかしな値が出るのとは別の理由がいるだろうが。 「1が素数ですらないから数列がモノボナッチがまったく変化しない」ってのは、まあ理に適ってるとしてもだよ？ 計算誤差で片付けていいのかなあ？あんまり不規則すぎてよくわかんね サンプルすうが10までってのもデータ数が足りないんだよなあ これ以上やるにはExcelでは荷が重いし、scilabに頼ったほうがいいとは思うけど、あんまり乗りこなしてないのよね～アイツ ちなみに、10までで唯一合成数の奇数ノナボナッチ うーん・・・やはり偶奇ではなく素数かどうかがカギか・・・ 数列の値が小さすぎて計算しきれてないっぽいな 15のペンタデカボナッチでは計算不可能だろうなぁ。間空きすぎなんだよもう ノナボナッチを-10じゃなくて-100から始めても、計算誤差が邪魔してダメだった～ 倍精度の型同士をつなげて4倍精度みたいにするオブジェクトがあればいけるだろうけどねえ 20桁以上の円周率を計算するときとかに使うアレよアレ あったとしてもまだ僕には使いこなせない サイラボにもそんなのねえかなー あーそうか。なんとなくわかったかも。 量子数が固有値になって、打ち消さなくなった的なやつじゃね？
元データの行列式を用意しなさい！
そうしたらできる範囲で因数分解してあげるから！ パズルゲームって行列式に似てるよね。僕パズルゲームあんまやらないんだけどね ルールが・・・覚えられないんだ。 ルービックキューブってあるだろ あれってノギスに似たトラップがあるやんか ギミックに夢中になってそのあとの説明がまるで頭に入らないアレ あの色のついたシールが内側にくることはあるの！？とか思うやん？ 分解したくなるやん？ぼく絶対直せない自信あるけど！ ノギスはモアレ現象の原理が気になって、測定器であることを忘れてしまうんや！
|
カレンダー
カテゴリー
ブログランキング
ブログランキング参戦中
よかったらポチッとお願いします＾＾
最新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
プロフィール
ＨＮ：
量子きのこ
年齢：
43
ＨＰ：
性別：
男性
誕生日：
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)
アクセス解析
|
