WANtaroHP (GMT Prob.Paper)

 

Contents

  • プロッティングポジション公式としては,Cunnane 公式が用いられています.
  • このページのグラフのプロットデータは,'F90 Engineering: 水文統計解析' と同じものです(札幌・前橋・京都の雨量データ).
  • 観測値のY軸の値は,X軸の値とY軸の値が線形関係になるよう定めています.Y軸値の変換方法を以下に示します.
確率紙名 クオンタイル 変数変換 Y軸値 X軸値
正規確率紙 xp=μ+σ*zp (xp-μ)/σ=zp Y=zp X=xp
対数正規確率紙 ln(xp)=μ'+σ'*zp {ln(xp)-μ'}/σ'=zp Y=zp X=log(xp)
グンベル確率紙 xp=c-a*ln{-ln(p)} (xp-c)/a=-ln{-ln(p)} Y=-ln{-ln(p)} X=xp
ワイブル確率紙 xp=c+a*{-ln(1-p)}1/k k*ln{(xp-c)/a}=ln{-ln(1-p)} Y=ln{-ln(1-p)} X=log(xp)
指数確率紙 xp=c-a*ln(1-p) (xp-c)/a=-ln(1-p) Y=-ln(1-p) X=xp


スクリプトとバッチファイル

FilenameDescription
awk_prb_paper.txtGMT 作図用データを作成するawkスクリプト
bat_prb_paper.txtGMT実行用バッチファイル(コントロール部)
bat_plot.txtGMT実行用バッチファイル(プロット実行部)
awk クリプトス実行のための引数
gawk -f awk_prb_paper.awk n fnameR > fnameW
n 確率紙の種類を指定する正の整数
fnameR入力ファイル名
fnameW出力ファイル名
n の値確率紙の種類プロット事例
1正規確率紙 fig_M_PRB_ND.png
2対数正規確率紙 fig_M_PRB_LN.png
3グンベル確率紙 fig_M_PRB_GD.png
4ワイブル確率紙 fig_M_PRB_WD.png
5指数確率紙 fig_M_PRB_ED.png


正規確率紙プロット

S_ND M_ND K_ND
札幌前橋京都


対数正規確率紙プロット

S_LN M_LN K_LN
札幌前橋京都


グンベル確率紙プロット

S_GD M_GD K_GD
札幌前橋京都


ワイブル確率紙プロット

S_WD M_WD K_WD
札幌前橋京都


指数確率紙プロット

S_ED M_ED K_ED
札幌前橋京都


pic
inserted by FC2 system