lectures:乱数
差分
このページの2つのバージョン間の差分を表示します。
両方とも前のリビジョン前のリビジョン次のリビジョン | 前のリビジョン | ||
lectures:乱数 [2022/08/25 14:18] – [Lagged Fibonacci法] kimi | lectures:乱数 [2022/08/25 14:32] (現在) – [線形合同法の特徴] kimi | ||
---|---|---|---|
行 67: | 行 67: | ||
* 結晶構造になる | * 結晶構造になる | ||
+ | === 実数の一様乱数 === | ||
+ | * $[0, M–1]$の乱数$\{x_i\}$を$[0, | ||
+ | * $ u_i = \displaystyle\frac{x_i}{M}$ | ||
+ | |||
+ | 例)$c = 0$, $a = 177$, $M = 2^15$, $X_0 = 1$ | ||
+ | ^ $n$ ^ $x_n$ ^ $u_n$ ^ | ||
+ | | 0| 1| 0.00003051758 | ||
+ | | 1| 177| 0.00540161133 | ||
+ | | 2| 31329| | ||
+ | | 3| 7441| 0.22708129883 | ||
+ | | 4| 6337| 0.19338989258 | ||
+ | |||
+ | === 線形合同法の実装系 === | ||
+ | * $X_n = (aX_{n–1} + c) \pmod M$ | ||
+ | * '' | ||
+ | * $a=1103515245$, | ||
+ | * '' | ||
+ | * $a=25214903917$, | ||
^ $2^{31}$個発生させ、2個ずつペアにして$x$-$y$座標としてプロット | ^ $2^{31}$個発生させ、2個ずつペアにして$x$-$y$座標としてプロット | ||
行 82: | 行 100: | ||
==== Lagged Fibonacci法の実装系 ==== | ==== Lagged Fibonacci法の実装系 ==== | ||
- | * random(Cの現標準乱数) | + | * '' |
* $X_n = (X_{n–r} + X_{n–s}) \pmod M$ | * $X_n = (X_{n–r} + X_{n–s}) \pmod M$ | ||
lectures/乱数.1661404709.txt.gz · 最終更新: 2022/08/25 14:18 by kimi