北国ぺんぎん

この正規分布乱数って工業的に使われてるみたいっす。
・・・自分、使ったことないっすよ？（工場勤務）
どうもぺんぎんっす( ◎v◎ )


イン・ワールドで統計解析はしないっすけどね・・・
欲しいと言われたもんっすから調べたっすけど。

方法は「ボックス・ミューラー法」ってやつっす。
1958年に・・・って今年50周年っすね。
おめでとうございますっす(?)

一様分布乱数2つから正規分布乱数2つを作るっす。
そのうちの1つを出したっす。

平均0で、標準偏差1っす。

---

///////////////////////////////
// The MIT Lisenses
// Copyright (c) 2008 Naoya Bellic

// Permission is hereby granted, free of charge, to any person obtaining a copy
// of this software and associated documentation files (the "Software"), to deal
// in the Software without restriction, including without limitation the rights
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
// copies of the Software, and to permit persons to whom the Software is
// furnished to do so, subject to the following conditions:

// The above copyright notice and this permission notice shall be included in
// all copies or substantial portions of the Software.

// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
// THE SOFTWARE.
///////////////////////////////
float BoxMuller()
{
　　return llSqrt(-2 * llLog(llFrand(1.0))) * llCos(TWO_PI * llFrand(1.0));
}

---

もうひとつのほうはcosをsinにしたやつっす。

端は厳密には違うっす。
使う一様分布乱数は(0,1]でないといけないっす。
0だと対数が定義できないっすね。(真数条件ってやつっす)
llFrand(1.0)は[0,1)を返してくるっす。

とはいえ、統計解析はインワールドではしないっすよね？
他にコレって何に使うんっすかね・・・


09/02/12
端の問題を改善したコードがリンク先にあるっす。