[改訂新版]C言語による標準アルゴリズム事典 (Software Technology)

作者: 奥村晴彦
出版社/メーカー: 技術評論社
発売日: 2018/04/19
メディア: 単行本（ソフトカバー）
この商品を含むブログ (1件) を見る

C言語による最新アルゴリズム事典という本が1991年に発売されていたのをご存じですか？
専門学校時代に手元にあると良いよと言われ、持ってはいたものの、あまり開いたことがありませんでした。
今日なんとなく開き、アルゴリズムを学びつつ、Kotlinで書いたらどう書くかというのをブログに書くのも面白いのではと思い、やってみました。

改訂新版を先ほど購入したので、内容自体は改訂新版を元にしています。

Weibull分布

p.429

概要

f:id:takattata:20181005225205j:plain
f:id:takattata:20181005225215j:plain

結論: 理屈分かってないけど、そういう分布の乱数の作り方ということは分かった

コード読解

本を元にC++にしたもの

#include <iostream>
#include <cmath>

using namespace std;

double randWeibull(double alpha) {
    double randomValue = rand() / (double)RAND_MAX;
    return pow(-log(1-randomValue), 1/alpha);
}

int main() {
    cout << randWeibull(0.2);
    return 0;
}

https://www.onlinegdb.com/online_c++_compiler

コードのrandomValueが、F(x)のxのところで、購入randomValueヶ月以内に当たる
alphaが〜故障型を決めるパラメータ

引数がalphaだから、
randWeibull関数は任意の故障型の購入経過月数の一様分布(=randomValue)において最初に故障する確率が求められる
という認識で合っているのだろうか🤔

Kotlinで書く

Kotlin v1.3-RC

import kotlin.math.*
import kotlin.random.*

fun randWeibull(alpha: Double): Double {
    val randomValue = Random(10).nextDouble(1.0)
    return -ln(1-randomValue).pow(1/alpha)
}

fun main(args: Array<String>) {
    print(randWeibull(1.3))
}