風力タービンによるエネルギー生産を分析する
WeibullDistributionを使って風速をモデル化し，GE 1.5MW風力タービンについて測定された効率曲線を使って予想される1年間のエネルギー生産量を求める．
dist = WeibullDistribution[2, 14 / Sqrt[π]];
turbine = Interpolation[{{0., 0.}, {0.5, 0.}, {1., 0.}, {1.5, 0.}, {2., 0.}, {2.5, 0.}, {3., 0.}, {3.5, 0.}, {4., 36.}, {4.5, 66.}, {5., 104.}, {5.5, 150.}, {6., 205.}, {6.5, 269.}, {7., 344.}, {7.5, 428.}, {8., 528.}, {8.5, 644.}, {9., 774.}, {9.5, 926.5}, {10., 1079.}, {10.5, 1211.}, {11., 1342.}, {11.5, 1401.}, {12., 1460.}, {12.5, 1477.}, {13., 1494.}, {13.5, 1500.}, {14., 1500.}, {14.5, 1500.}, {15., 1500.}, {15.5, 1500.}, {16., 1500.}, {16.5, 1500.}, {17., 1500.}, {17.5, 1500.}, {18., 1500.}, {18.5, 1500.}, {19., 1500.}, {19.5, 1500.}, {20., 1500.}, {20.5, 1500.}, {21., 1500.}, {21.5, 1500.}, {22., 1500.}, {22.5, 1500.}, {23., 1500.}, {23.5, 1500.}, {24., 1500.}, {24.5, 1500.}, {25., 1500.}, {25.5, 0.}, {26., 0.}, {26.5, 0.}, {27., 0.}, {27.5, 0.}, {28., 0.}, {28.5, 0.}, {29., 0.}, {29.5, 0.}, {30., 0.}}];Row[{Plot[24 365 PDF[dist, x], {x, 0, 24}, AxesLabel -> (Style[#, Bold]&) /@ {"m/s", "時間"}, Filling -> Axis, BaseStyle -> {FontFamily -> "Verdana"}, PlotLabel -> (Style[#, 12, Bold]&)@"風速分布", ImageSize -> 250], Plot[turbine[x], {x, 0, 24}, AxesLabel -> (Style[#, Bold]&) /@ {"m/s", "kW"}, BaseStyle -> {FontFamily -> "Verdana"}, PlotLabel -> (Style[#, 12, Bold]&)@"風力タービンの効率", ImageSize -> 250]}]NExpectation[24 365 turbine[x], xdist]