皮尔逊曲线计算原理.docVIP

3　理论P - Ⅲ曲线的绘制 311　用GAMMA INV函数计算对应频率P的xp 值 GAMMAD IST为γ分布函数,其反函数GAMMA INV 可返回具有给定频率的γ累积分布的区间点。但不能直接应用GAMMA INV函数,因为GAMMAD IST与P - Ⅲ曲线的γ分布函数关系式不尽相同。 P - Ⅲ曲线的频率密度函数为: f ( x) =（βα/Γ (α)）( x - α0 )α- 1 e-β( x -α0) (2) 式中: α = 4 /C2s ; β = 2 / ( x·Cs ·Cv ) ; α0 = x (1 - 2Cv /Cs ) x为随机变量; f ( x)为频率密度函数;Γ(α)为α的γ分布函数。 P与xp 之间的关系为: P = P ( x ≥ xp ) =βα/Γ (α) ∫∞xp( x - α0 )α- 1 e-β( x - a0) dx (3) 　　令tp =β( x -α0 ) ,有: P = P ( x ≥ xp ) =1/Γ (α) ∫∞tpt α- 1 e- t dt = 1 -1/Γ (α) ∫tp0t α- 1 e- t dt (4) 　　水文计算中,随机变量xp 可用下式表示: xp = Kp ·x = (1 +Φp Cv ) x 离均系数Φp =（xp – x）/（Cv ·x） =（Cs/2）tp -（2/Cs） (5) 　　式(4)变换得: 1 - P =1/Γ (α) ∫tp0t α- 1 e- t dt (6) 式(6)右边的积分式称为标准Γ分布函数; tp 为标准Γ 分布分位数, excel内置函数GAMMA INV为返回Γ累积分布 函数的反函数,故tp = GAMMA INV (1 - P,4/C2s， 1) (7) 由式(5) 、式(7)得: Kp = 1 +CV (Cs/2·GAMMA INV (1 -P,4/C2s, 1) -2/Cs) (8) 312　理论P - Ⅲ曲线的绘制 根据公式(8)计算各频率的模比系数Kp 值,同时按照公 式(1)计算各频率对应的横向坐标值,如表3所示。 在“流量频率曲线”工作表中设置源数据。选择“理论频率 曲线”工作表中的A列值为x值, E列值为y值,图表类型选择 “无数据点平滑线散点图”绘制理论P - Ⅲ曲线,如图1所示。