新年、あけましておめでとうございます。 今年も「りょうとのITブログ」をよろしくお願いします。 さて、新年1回目のエントリは、「プログラミングについて」です。 久々ですね。 しかも言語はR! 果たしてどれだけの需要があるのか?そんなものはガン無視です。 能書きはこれくらいにして、本題に入ります。 やることは、タイトルにありますように、 「モンテカルロ法で円周率を計算」 です。 「モンテカルロ法とは?」「どうやって円周率を計算するのか?」 といった事にも触れます。 本エントリの大筋は、 1. モンテカルロ法とは 2. モンテカルロ法で円周率を計算するアルゴリズムについて 3. Rで円を描画 4. Rによる実装及び計算結果 5.
文部科学省発行「高等学校情報科『情報Ⅰ』教員研修用教材」の「学習16」にある「確定モデルと確率モデル」では確率モデルを使ったシミュレーション手法としてモンテカルロ法による円周率の計算が紹介されています。こちらの内容をJavaScriptとグラフライブラリのPlotly. jsで学習する方法を紹介いたします。 サンプルプロジェクト モンテカルロ法による円周率計算(グラフなし) (zip版) モンテカルロ法による円周率計算(グラフあり) (zip版) その前に、まず、円周率の復習から説明いたします。 円周率とはなんぞや? モンテカルロ法で円周率を求めるのをPythonで実装｜shimakaze_soft｜note. 円の面積や円の円周の長さを求めるときに使う、3. 14…の数字です、π(パイ)のことです。 πは数学定数の一つだそうです。JavaScriptではMathオブジェクトのPIプロパティで円周率を取ることができます。 alert() 正方形の四角形の面積と円の面積 正方形の四角形の面積は縦と横の長さが分かれば求められます。 上記の図は縦横100pxの正方形です。 正方形の面積 = 縦 * 横 100 * 100 = 10000です。 次に円の面積を求めてみましょう。 こちらの円は直径100pxの円です、半径は50です。半径のことを「r」と呼びますね。 円の面積 = 半径 * 半径 * π πの近似値を「3」とした場合 50 * 50 * π = 2500π ≒ 7500 です。 当たり前ですが正方形の方が円よりも面積が大きいことが分かります。図で表してみましょう。 どうやって円周率を求めるか? まず、円の中心から円周に向かって線を何本か引いてみます。 この線は中心から見た場合、半径の長さであり、今回の場合は「50」です。 次に、中心から90度分、四角と円を切り出した次の図形を見て下さい。 モンテカルロ法による円周率の計算では、この図に乱数で点を打つ 上記の図に対して沢山の点をランダムに打ちます、そして円の面積に落ちた点の数を数えることで円周率が求まります!
5なので、 (0. 5)^2π = 0. 25π この値を、4倍すればπになります。 以上が、戦略となります。 実はこれがちょっと面倒くさかったりするので、章立てしました。 円の関数は x^2 + y^2 = r^2 (ピタゴラスの定理より) これをyについて変形すると、 y^2 = r^2 - x^2 y = ±√(r^2 - x^2) となります。 直径は1とする、と2. で述べました。 ですので、半径は0. 5です。 つまり、上式は y = ±√(0. 25 - x^2) これをRで書くと myCircleFuncPlus <- function(x) return(sqrt(0. 25 - x^2)) myCircleFuncMinus <- function(x) return(-sqrt(0. 25 - x^2)) という2つの関数になります。 論より証拠、実際に走らせてみます。 実際のコードは、まず x <- c(-0. 5, -0. 4, -0. 3, -0. モンテカルロ法による円周率の計算 | 共通教科情報科「情報Ⅰ」「情報Ⅱ」に向けた研修資料 | あんこエデュケーション. 2, -0. 1, 0. 0, 0. 2, 0. 3, 0. 4, 0. 5) yP <- myCircleFuncPlus(x) yM <- myCircleFuncMinus(x) plot(x, yP, xlim=c(-0. 5, 0. 5), ylim=c(-0. 5)); par(new=T); plot(x, yM, xlim=c(-0. 5)) とやってみます。結果は以下のようになります。 …まあ、11点程度じゃあこんなもんですね。 そこで、点数を増やします。 単に、xの要素数を増やすだけです。以下のようなベクトルにします。 x <- seq(-0. 5, length=10000) 大分円らしくなってきましたね。 (つなぎ目が気になる、という方は、plot関数のオプションに、type="l" を加えて下さい) これで、円が描けたもの、とします。 4. Rによる実装 さて、次はモンテカルロ法を実装します。 実装に当たって、細かいコーディングの話もしていきます。 まず、乱数を発生させます。 といっても、何でも良い、という訳ではなく、 ・一様分布であること ・0. 5 > |x, y| であること この2つの条件を満たさなければなりません。 (絶対値については、剰余を取れば良いでしょう) そのために、 xRect <- rnorm(1000, 0, 0.
モンテカルロ法の具体例として,円周率の近似値を計算する方法,およびその精度について考察します。 目次 モンテカルロ法とは 円周率の近似値を計算する方法 精度の評価 モンテカルロ法とは 乱数を用いて何らかの値を見積もる方法をモンテカルロ法と言います。 乱数を用いるため「解を正しく出力することもあれば,大きく外れることもある」というランダムなアルゴリズムになります。 そのため「どれくらいの確率でどのくらいの精度で計算できるのか」という精度の評価が重要です。そこで確率論が活躍します。 モンテカルロ法の具体例として有名なのが円周率の近似値を計算するアルゴリズムです。 1 × 1 1\times 1 の正方形内にランダムに点を打つ(→注) 原点(左下の頂点)から距離が 1 1 以下なら ポイント, 1 1 より大きいなら 0 0 ポイント追加 以上の操作を N N 回繰り返す,総獲得ポイントを X X とするとき, 4 X N \dfrac{4X}{N} が円周率の近似値になる 注: [ 0, 1] [0, 1] 上の 一様分布 に独立に従う二つの乱数 ( U 1, U 2) (U_1, U_2) を生成してこれを座標とすれば正方形内にランダムな点が打てます。 図の場合, 4 ⋅ 8 11 = 32 11 ≒ 2. 91 \dfrac{4\cdot 8}{11}=\dfrac{32}{11}\fallingdotseq 2. 91 が π \pi の近似値として得られます。 大雑把な説明 各試行で ポイント獲得する確率は π 4 \dfrac{\pi}{4} 試行回数を増やすと「当たった割合」は に近づく( →大数の法則 ) つまり, X N ≒ π 4 \dfrac{X}{N}\fallingdotseq \dfrac{\pi}{4} となるので 4 X N \dfrac{4X}{N} を の近似値とすればよい。 試行回数 を大きくすれば,円周率の近似の精度が上がりそうです。以下では数学を使ってもう少し定量的に評価します。 目標は 試行回数を◯◯回くらいにすれば,十分高い確率で,円周率として見積もった値の誤差が△△以下である という主張を得ることです。 Chernoffの不等式という飛び道具を使って解析します!
024\)である。 つまり、円周率の近似値は以下のようにして求めることができる。 N <- 500 count <- sum(x*x + y*y < 1) 4 * count / N ## [1] 3. モンテカルロ法 円周率 考え方. 24 円周率の計算を複数回行う 上で紹介した、円周率の計算を複数回行ってみよう。以下のプログラムでは一回の計算においてN個の点を用いて円周率を計算し、それを\(K\)回繰り返している。それぞれの試行の結果を に貯めておき、最終的にはその平均値とヒストグラムを表示している。 なお、上記の計算とは異なり、第1象限の1/4円のみを用いている。 K <- 1000 N <- 100000 <- rep(0, times=K) for (k in seq(1, K)) { x <- runif(N, min=0, max=1) y <- runif(N, min=0, max=1) [k] <- 4*(count / N)} cat(sprintf("K=%d N=%d ==> pi=%f\n", K, N, mean())) ## K=1000 N=100000 ==> pi=3. 141609 hist(, breaks=50) rug() 中心極限定理により、結果が正規分布に従っている。 モンテカルロ法を用いた計算例 モンティ・ホール問題 あるクイズゲームの優勝者に提示される最終問題。3つのドアがあり、うち1つの後ろには宝が、残り2つにはゴミが置いてあるとする。優勝者は3つのドアから1つを選択するが、そのドアを開ける前にクイズゲームの司会者が残り2つのドアのうち1つを開け、扉の後ろのゴミを見せてくれる。ここで優勝者は自分がすでに選んだドアか、それとも残っているもう1つのドアを改めて選ぶことができる。 さて、ドアの選択を変更することは宝が得られる確率にどの程度影響があるのだろうか。 N <- 10000 <- floor(runif(N) * 3) + 1 # 宝があるドア (1, 2, or 3) <- floor(runif(N) * 3) + 1 # 最初の選択 (1, 2, or 3) <- floor(runif(N) * 2) # ドアを変えるか (1:yes or 0:no) # ドアを変更して宝が手に入る場合の数を計算 <- (! =) & () # ドアを変更せずに宝が手に入る場合の数を計算 <- ( ==) & () # それぞれの確率を求める sum() / sum() ## [1] 0.
参考文献: [1] 河西朝雄, 改訂C言語によるはじめてのアルゴリズム入門, 技術評論社, 1992.
?露出癖に目覚めた女達の過激な野外露出画像 7/24 (土) 21:30 【盗撮】「週1回、秘書が社長と絶対セックスする日」とあるオフィスから流出したヤバすぎる映像がこちら エロマーゾフの兄弟 7/25 (日) 5:05 【正常位】平野綾さん、舞台のセックスシーンの稽古風景がこちら… 7/30 (金) 5:33 おっぱいを見られてるのに気づいてない…ゆるゆるの胸元がセクシー過ぎる胸チラエロ画像 エロ画像CLUB 7/24 (土) 7:00 すっぽんぽんなのに堂々と野外を闊歩できる! ?全裸の外国人が楽しそうなボディーペイントエロ画像 7/21 (水) 7:00 【鬼畜ぶっかけ】街で見かけた女の子に背後から近づいてブッカケて逃げる極悪非道な"ぶっかけ痴漢"のエロ画像 キャプエロ画像ず 7/24 (土) 12:00 【マンスジ画像】割れ目に喰い込んだパンツが見せる奇跡の縦線"マンスジ"のエロ画像(33枚) なにもかもあけっぴろげに…全裸セックス画像100枚 えっちなお姉さん。 7/26 (月) 18:00 今夜のオナネタ。エロ画像詰め合わせ Vol. 474 7/31 (土) 14:00 鈴村あいり 『スズムラ無双』 10P乱交がエロ過ぎる【100作品記念特別企画】 おっぱいの楽園♪ 7/18 (日) 19:06 鷲尾めい×クリムゾン 漫画のようなJcup神乳美女一年間痴かんされ続けた 7/19 (月) 22:04 【体験談・小説】好みの同僚M(25)をレイプしてしてまった。 エチケン - エッチな体験談と写メ 7/26 (月) 21:00 【体験談・小説】ご近所の6人家族のお話。その後の後のごの5 7/29 (木) 22:23 これがリアル…ごく一般的なカップルの性交!強烈な生々しさが興奮するネットに流出したハメ撮り画像 エロ写メサークル 7/26 (月) 9:00 激しすぎるセックスにドン引き…マン汁飛び交う集団FUCK! 【画像】美人なのにガキみたいに小さいおっぱいの女ｗｗｗ | エロ年表アンテナ. ?エロ過ぎる乱交GIF画像 7/23 (金) 9:00 天宮花南 美巨乳グラドルのAVデビュー画像 エロ画像 PinkLine 7/24 (土) 15:02 堀中未来 経験人数500人の生徒会長セックス画像 7/25 (日) 9:02 【今井夏帆】パイパン黒ギャルの巨乳美女が強制連続中出しセックス エロ画像が見たいんだ! 7/22 (木) 21:20 【八木奈々】透明感あふれる清純派美少女のイチャラブ密着セックス 7/19 (月) 21:12 休みの日は公園に行け!ぱんちらパンチラまみれで最高だな!
エロ画像が見たいんだ!
名前 タイトル メールアドレス URL パスワード 非公開コメント 管理者にだけ表示を許可する
盗撮カメラlog 無料エロ動画AV倶楽部 メガえろ動画 サイトについて 当サイトは管理人が個人的にネット上で集めた画像や動画などを記事にしております。 権利を侵害する目的は一切御座いませんので、掲載について問題がある場合はお手数ですが、 削除フォーム よりご連絡下さい。確認後、迅速に対処または削除いたします。 当サイトは相互リンクを随時募集しています。リンクフリーですので、アクセスを頂ければ確認とれ次第リンク致します。 多くのアクセスを頂いているサイト様は、画像ピックアップ等でアクセスをお返しします! スマートフォン版 いつでも暇つぶし程度にエロアングルの画像でも見て下さいw お手持ちのスマホでQRコードを読み込んで下さい。 メニュー 上に戻る 投稿のRSS サイトアップ お問い合わせ
3. 2」 7/21 (水) 16:00 【人妻エロ漫画】 近所の巨乳美人妻の変態趣味を目撃した口止めに熟れた身体に中出しする! エロマンガ熟女 7/19 (月) 11:51 【熟女エロ漫画】 「お母さんとエッチしたい」とマザコン息子が本気で母親に告白してみたら…! 声優の松岡禎丞さん、キリトみたいな主人公のキャラばっかりやってるイメージあるけど、けっこういろんな役を演じてるよね - カチカチ通信. 7/20 (火) 11:31 【いのうえみたん】おっぱい10回も揉んだらエッチに発展するに決まってる件 エロマンガ番長 7/22 (木) 1:01 【あきは@】水泳部の女子たちがエロすぎてコーチどころじゃない件 7/21 (水) 0:01 最近お腹のお肉が気になる若妻さんがマッサージにいったら旦那さんが横の部屋にいるのにエッチな性感マッサージをされてアへっちゃう【無料 エロ同人】 エロ同人カフェ~エロ漫画・エロ同人・エロアニメ~ 7/19 (月) 8:13 巨乳美女が緊縛拘束されて複数ちんこに囲まれ乱暴に陵辱レイプされちゃってる【無料 エロ同人】 7/19 (月) 8:11 新人OL泥酔レイプ!お酒を飲ませてヤリタイ放題犯しまくる泥酔強姦 メス狩り 7/19 (月) 21:50 麻酔を吸引させられ意識のない無抵抗の美人患者を弄ぶ歯科医師 7/24 (土) 21:28