ワイブル分布の積分

　風力発電に関する本を読んでいると、風速の度数分布はワイブル分布で表されることが多いと書いてあった。ワイブル分布とは、品質管理などでも使われる分布関数である。

　風速の度数分布の確率密度関数f(V)と、風速がVa以下の確率F(Va)はそれぞれ以下のようにあらわされる。

　f(V)=(k/c)(V/c)^(k-1)exp〔-(V/c)^k〕（V：風速、k：形状係数、c：尺度係数）

　F(Va)=1-exp〔-(Va/c)^k〕

　f(V)からF(Va)を導出する過程は書いてないが、これはf(V)を0～Vaまで積分すればよいはず。理工書を読むのに一番大切なのは数式を理解するということだ。出てきた数式は、自分の手で追ってみる。これが理工書を読みこなす秘訣である。やってみるとミスプリが見つかったりして意外と楽しい。

　導出過程の書いてあるものはもちろんのこと、導出過程が省略されているものについても自分でその行間を埋めていくとよい。そうやって初めて理工書に書いてあることが理解できるようになるのだ。

　最初は生真面目に部分積分法などを使ってやってみたが、私は計算自体はあまり得意でないので、どんどん泥沼に落ち込んでいく（笑）。

　ここでふと気が付いたのは、(V/c)の形。係数のほうは(k-1)乗、exp関数の中はｋ乗という形になっている。

　ここで次の積分を考える。

　∫x^(k-1)exp(-x^k)dx

　今、t=-x^k と変数変換すると、dt=-kx^(k-1)dx

　∴　与式=-(1/k)∫exp(t)dt=-(1/k)exp(t)　

　変数を元に戻すと、与式=-(1/k)exp(-x^k)　・・・(1)

　となり、この式は簡単に積分できることが分かる。

　そこで最初のf(V)とF(Va)の式に戻る。f(V)の積分においては、(1)においてx=V/cとおいて0～Vaまで積分するとF(Va)が簡単に導出できる。

　関数の形は複雑だが、それに惑わされないで自分の手で計算して確かめてほしい。それにしてもブログって数式を表すときになんとも不便だなあ。

　またワイブル分布において、k=2の場合はレイリー分布と呼ばれ、こちらも風速の度数分布を表すのに使われるようだ。