忍者ブログ
日々の研究生活のメモ書きなど
ガンベル分布についての詳しい情報はwikiなどどうぞ↓

■ 参考
ガンベル分布(gumble) wiki



set term png
set output "pin_gumbel_pdf.png"

set title "gumbel pdf(a, b, x)=1/b*exp(-(x-a)/b)*exp(-exp(-(x-a)/b))"
set xlabel "x"
set xlabel "pdf(x)"
set xrange [-5:20]

f(a, b, x)=1/b*exp(-(x-a)/b)*exp(-exp(-(x-a)/b))
plot f(0.5, 2.0, x) w l lw 2 title "a=0.5, b=2.0"
replot f(1.0, 2.0, x) w l lw 2 title "a=1.0, b=2.0"
replot f(1.5, 3.0, x) w l lw 2 title "a=1.5, b=3.0"
replot f(3.0, 4.0, x) w l lw 2 title "a=3.0, b=4.0"
replot f(3.0, 6.0, x) w l lw 2 title "a=3.0, b=6.0"
set term png
set output "pin_gumbel_pdf.png"
replot
こんな感じのスクリプトを書いて
gnupot plot_gumble_pdf.plt
で次のようなグラフができるはず

pin_gumbel_pdf


次の記事でROOTでのグラフの書き方も説明予定







PR
プロフィール
HN:coffee
職業:物理屋(自称)
趣味:映画鑑賞
自己紹介:
#include <stdio.h>
#include "MyProfile.h"

#define TWITTER coffee_pote

#define WISH_LIST
amazonのほしい物リスト
#ifdef RICH_FLAG
// ↑いつも支援いただきありがとうございます m(_ _)m
#endif


int main(void){

printf("\n");
printf("D論・・・? あぁそんな子もいましたね(執筆中)\n");
printf("\n");
printf("猿でもわかるgnuplot を執筆中(こっちの執筆は半年以上何も進んでいない・・・・)\n");
/* 最終更新 2017/07/19 */
return 0;

}
カウンター
ブログ内検索
ツイートするボタン
Flickr

Template "simple02" by Emile*Emilie
忍者ブログ [PR]