kogelab.memo

2011-05-07 MetropolisサンプラーをHTML5のcanvasで描いてみる

ついでにGibbs samplerもcanvasで書いてみる

algorithm, Programming, javascript

MetropolisがGibbsと並んで有名って言うなら、Gibbsもやってみなければいけない。

というわけでやってみた。

ここであそべる

おおー。たのしい。

ただ、無相関の2変量ガウスでこれをやっても楽しさ半減ですね…

Permalink | コメント(0) | トラックバック(0) | 22:45

MetropolisサンプラーをHTML5のcanvasで描いてみる

algorithm, javascript, Programming

MCMCの中でもGibbs samplerと並んで有名なMetropolis samplingをjsでやってみる。

HTML5のcanvasを使って視覚化もしてみる。

このへんであそべる

赤いところがAcceptされたサンプル、青いのがRejectされた動き

moreで1ステップとか2ステップづつ進むとニョロニョロ動いて楽しい。

対象は相関のない2変量ガウス（正規）分布。

乱数生成器はMath.random()なので、本当はサンプラーに使うべきではないですが、どの道モックなので気にしない。

Permalink | コメント(0) | トラックバック(0) | 17:21

2011-05-01

切符の番号とかナンバープレートの4つの数字と＋－×÷を使って10を作るアレを解いてみる

algorithm, Programming, javascript

　切符にプリントされている４個の数字とか、ナンバープレートに書かれている４桁の数字を、四則演算のみを使って10にするという遊びあるじゃないですか。

そう、あれ。あれです。誰か名前を付けていてもよさそうなんだけど、特別名前が付いている訳でもなさそうなあれ。

　どうも海外では、これと似たようなゲームを出題するTV番組があったらしく、そこからCountdown numbers gameという名前で呼ばれているらしいです。

ルールはもうちょっと複雑で、数字が6個だったり、10にするんじゃなくて、その都度違ったり、時間制限があるとか。

　これを計算機で解く場合、おそらく一番単純な方法は、葉ノードが数字、接点が演算子の木を全通り出力　→　評価という流れかと。

そういうわけで、実装してみた。例によってJavaScriptです。

このへんであそべる

　遊び方は見たまんま。

　計算したい条件（ナンバープレートの数字と、作りたい10という数字）を入力すると、可能な計算式が全て出ます。4つの数字以上にも対応、途中計算で分数が出てくるものも解として出力します。

　大昔に、auの携帯向けのフリーアプリで、まさにこの問題をプレイヤーに解かせるというゲームがあったのですが、困ったことに、自力で解けないときに回答が出てこない（すごいストレス）＆なんか解けそうに無い問題が出てくる（理不尽）という代物で、非常に不満だったのですが、これがあれば万事解決。

　これを使えば暇なドライブでどうしても解けないナンバープレートを見つけたときに、「解なし」と自信たっぷりにいえる、その快感に酔いしれること請け合いです。

　GWの単調なドライブのお供に是非。

追記

　10 puzzle とか Make 10とかいうらしい。解く前にググれ、か。

Permalink | コメント(0) | トラックバック(0) | 02:33

2011-04-29

モンテカルロ法とJavaScriptでπ

javascript, algorithm, Programming

　これはもう説明が全く不要な、コード書く人なら間違いなく通る、モンテカルロシミュレーションでπの値を計算しようという遊び。

これはとても楽しい。視覚化されると楽しさ10倍。というわけでわりと投げやりにJavaScriptで書いた。

ここであそべる

画像にあるサンプル点は全て1個1個がdivなので、1万サンプルとかやると真面目に死にますので気をつけて。

あと、計算誤差については何も考慮せず作ったので、ある程度の精度以上は出ません。というか、非常に精度が悪い。

電車内で暇に任せてiPod Touchで書いたコードに枝葉を死ぬほどつけたらなんかできたのでとりあえず公開。

Permalink | コメント(0) | トラックバック(0) | 04:10

2011-04-27 括弧対応をとる

括弧対応をとる

雑記, Programming

一瞬必要だったので。corr_tableをいじればどんな組でも可。括弧の交差はダメ。

var corr_table = { '(' : ')', '[' : ']' };
function checkBrackets( t ){
	var stack = [];
	for( var i = 0; i < t.length; ++i ){
		for( var j in corr_table ){
			if( t[i] == j ){
				stack.push( j );
				break;
			} else if( t[i] == corr_table[j] ){
				if( stack.pop() != j ){
					return false;
				} else {
					break;
				}
			}
		}
	}
	return stack.length == 0;
}