ウィローカンパニー(有) | 物理・simulation

サポートベクタマシン実践～その１～

2008.04.10 Thursday 23:44| comments (4) | - |

サポートベクタマシン(SVM)をせっかく習ったので使ってみようと。
今回は
【身長】と【体重】のデータから【性別】を判別するＳＶＭを作りたい。
二次元の直線による分離を実装するのだ。
まぁ、エクセルでいいだろ、エクセル最高！エクセル最高！
（そーいや、エクセルって英語で「優れてる」って意味だよな、上の文はまるで馬から落馬みたいだ）

ＳＴＥＰ①
適当に身長・体重情報を作る。
作るのだるいんでネットから探してきた。

ＳＴＥＰ②
目で見てみる。
「やってみて、触ってみて、確かめて」
だっけ、なんか俳句のような575の文章

おい、物理出のクソ共

2008.03.11 Tuesday 22:56| comments (4) | - |

「単位長さ」とは何か教えてください。

学生の時、フガフガ先生に
「フーガフガフガ！犬君はフーガフガフガ！それに引き換えmona君はフーガフガ！」
と単位長さについて説教されたのを思い出したんだ。
そして、それからもう4年も経っているというのにまだ正しい理解をしていないのは悲しいのでお前らを頼りますよと。

くあああああ、爆笑問題のニッポンの教養を見てるとすげぇなんかこう・・・
大学に帰りたくなりますね、帰らないけど(´∀｀)
とりあえず書籍はコンプリートしよう。
理物は紙とペンとPCさえあればなんとでもなるのだから好きにやらせてもらう。

マルチスレッド・・・(ﾟ∀ﾟ)神のﾖｶｰﾝ

2007.08.31 Friday 00:39| comments (4) | - |

最近コメント欄でkeiがくれた話に関して実験してみた。

お題：マルチスレッドなプログラムを書くだけでCPU二個にきちんと処理を振り分けてくれるかどうか。
計算種目：πの計算
開発環境：Turbo C++ Builder
CPU：Intel Core 2Duo

計算結果
シングルスレッド
・計算回数10^9
・計算時間122秒
・CPU使用率50％くらい

２スレッド
・計算回数10^9
・計算時間61秒
・CPU使用率100％

４スレッド
・計算回数10^9
・計算時間61秒
・CPU使用率100％

結論：
マルチスレッドで計算させた場合、
シングルスレッドの場合にくらべてちょうど半分の時間で済んだ。
ということは完璧に並列化されているということが言える。
また2スレッド以上の場合はCPU使用率100%なので、
ﾃﾞｭｱﾙｺｱをフル活用してるんだろう。
並列化させたかったらスレッドたてるだけで(・∀・)ｲｲﾈ!!

いい加減モデル

2007.08.26 Sunday 21:00| comments (8) | - |

先週から適当に作っている２時限イジングモデル。

プログラミングの能力は学生の時よりそれはもう遥かにパワーアップしているわけで、
同時に二つの系をシミュレーションできるようにしている。
無論N個の系への拡張も容易だ。
これができたので拡張アンサンブル法も簡単に計算できるだろう。

しかし、CPU使用率を見てみると・・・

うは、DualCore CPUの意味が全くねぇｗ
しかし、マルチスレッドプログラムにするのはちょいとめんどくさいな・・・

一様分布を作るのだって楽じゃぁねぇよ

2007.08.24 Friday 21:12| comments (6) | - |

どうも、あっしです。
今日は金曜なのに飲み仲間はサマーバケーションでどっかいってるので一人酒です。
死ぬほど寂しい。
call me please.

～不定期連載・酔っ払いの戯言シリーズ～
-N次元球の表面に一様分布する乱数はどのように生成したらよいか？について-
単純に
角度変数１として０～πの一様乱数をふる
角度変数２として０～πの一様乱数をふる
・・・・
角度変数N-1として０～2πの一様乱数をふる
としてそれらをまとめて一個の変数と思えばよいと思われるかもしれないが、そうではないだろう。
それは早計というものだ。

よくわかる例は我々の知覚できる次元である3次元だ。
よくある物理の教科書に習って角度はθ、φとでもしよう。

φ方向に０～2πの一様乱数をふる（これはよい）
θ方向に０～πの一様乱数を・・・
って、ちょーーーーーーっと待てぇぇいいいいい！！！
そんなことしたらおめぇ、球面上のZ軸上とか下のほうがパンッパンになっちまって、
Z=0のXY平面がスカッスカちなっちまうだろ？

とまぁ、端的にいえばこういうことなんだ。
わかってくれたかい？＞犬

一様だってのは要するに「密度が均一だ」っちゅーことですからね。
上のやりかただとZの値がでかいほうがミッチミチな状態になっちまうのは
図をかけばわかってくれるとおもう。

そこでそいつを補正するのがヤコビアンですよ！（今回は球面内ではなく球面上なので角度部分のみ）
三次元のヤコビアンの角度方向であるsin(θ)ってのはこいつを補正してくれますよねぇ。

･･･30分後・・・
OK BOY、やっとスマートな説明を思いついた。
球面上の面積要素ってr^2sin(θ)dθdφだよね？
（ここでrは一定なんだからそこは無視してＯＫだ）
じゃぁ、そもそもの問題を読み替えると
【この面積要素の上に一様に点をまくためにはどうしたらよいか】
という問題になるわけだ。
そうすると・・・OH,MY GOD!!
sin(θ)のウェイトでθは分布を作らないとだめだよね？
ってことですわ。
一定　＝　乱数の分布/sin(θ)dφdθ　
ってかんじかな。だからsin(θ）に比例するように乱数を作りなさいよと。

あぁ、実際に紙の上に球面上に等間隔で点を打った場合を創造するといいよ。
やっぱＺ軸の上とか下のほうがミッチミチになるだろ？
だからＺ＝０のＸＹ平面はsin(90°）=１の一番大きなウェイトをつけて
点を打つことでそのなんだ、等間隔に打ってしまったためにできた密度差を埋めるようにするわけだ。

・・・・えらそうにウンチク書いてるけど、間違ってたらｽﾏｿwww
酔っ払いなんで。
こうやってweb上に誤った情報が蓄積されていくのは悲しいね。
（↑じゃぁ消せよってわけだが、酔っ払いなんで許してくれ。飲酒執筆だ飲酒執筆）

ポアソン分布について

2006.09.20 Wednesday 21:47| comments (3) | trackbacks (0) |

Z=X+Y
(X,Yはポアソン分布の確率変数)
のとき
Zの平均＝Xの平均+Yの平均
Zの分散＝Xの分散+Yの分散
となり、Zの分布関数自体もポアソン分布になるが
Z＝XーY
の時、Zの分布関数はポアソン分布

ではない

今日は仕事後にらーめん行ってきた。
週末は一応大学の同期飲みらしいので週末の予定を先取りだ!
そこでビール一杯飲んだのだが…いいね。

音波

2006.03.04 Saturday 22:28| comments (0) | trackbacks (0) |

頭の悪そうな人間の車から『ドフッ!ドフッ!』と重低音だけがやたら響いてきたり、妹の部屋からこれまた『ドフッ!ドフッ!』とただの雑音にしか聞こえない音楽が漏れてきて不快なことこの上ないのですが、注意したところで原因を理解し、状況を改善するほどの頭を持ってないのはもうこれ自明なので最近は放置です。

んでまぁ、このやり場のない怒りを「なんで低音ばかりがよく聞こえるんかいな?」という方向へ向けてみた。

高木君に聞いたら「知りません」とか言われて
(´A｀).oO(好奇心のかけらもない奴だな、逝ってよし)
と思ったので自分で調べました。

したらば、あっさりみつかった。ここ、もういっちょ、どっせぃ
音響学会なんてあんのね。

結局回折の波長依存性の問題だったようだ。
遮蔽物での吸収も関係しようだけどまぁいいやー
あとはこいつを数式で書くことを考えYO!

弾性エラスティカ

2006.02.13 Monday 21:08| comments (4) | trackbacks (0) |

連続体力学を勉強中なあっしなわけだが、
「弾性エラスティカ」なるものが登場して来てはぁはぁだ。

かの天才オイラーが命名したらしいが…

し、椎名林檎の曲のタイトルっぽい…

そんだけ。

これでぼかぁ棒の曲がりやら紐を垂らしたときの曲がり方を計算できるようになったわけだ。
本当に計算できるかはしらんがな。
こたつ

「がな」つながりです…(二回目)

独占型企業における利益の値段依存性について・その2

2006.02.12 Sunday 23:10| comments (0) | trackbacks (0) |

人間の所得格差をポアソン分布で導入しようと思ったが、
離散分布なのでグラフがガタガタでムッカームカムカなのでやめ。
毎度おなじみ正規分布にした。
あと「金持は心がゆたかパラーメータ」も導入した。
金持喧嘩せず。

状況

人間は平均2000円分散2000円の所持金格差を持つ

4000円で所持金をスケールした「金持は心がゆたかパラーメータ」により物を購入するかしないかを乱数で降ってみた

「縦軸：売り上げ、横軸：物の値段」でsimulation

結果

なんか僕ののぞんだ形(多価関数)ではなくてがっかりだ。
今度はライバル企業でもいれるか…

書き忘れてたけど、これをやってる目的は
「ほぉら、売り上げは多価関数だっただろ?
だから値下げしてみろ!意外と売り上げがあがるかもしれないぞ!
(そして、あっしはその商品を安く買えてうれしいなっと…)」
ということです。

独占型企業における利益の値段依存性について・その1

2006.02.08 Wednesday 00:28| comments (7) | trackbacks (0) |

ある一社が独占的に販売している物の値段をいくらに設定すれば
利益がもっとも大きくなるのかを考えてみた。
暇だから。

状況と条件
・人は一様ランダムで所持金1～1000円を持つ
・所持金が物の値段以上であればそれを買う
・物の値段を1～1000円の範囲で変化させる

「y:売り上げ,x:物の値段」のグラフ
人は1000人。

どうみても一価関数&極大値のポイントが所持金の平均値に一致しとる…
所持金の分布を現実に即してポアソン分布あたりに変えてみようかねぇ

1/1

<< September 2010 >>
Sun	Mon	Tue	Wed	Thu	Fri	Sat
			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