用matlab绘制P三曲线,科学网—水文频率曲线及MATLAB绘制 - 张凌的博文

1. 水文频率曲线

水文分析计算中使用的概率分布曲线俗称水文频率曲线，习惯上把由实测资料(样本)绘制的频率曲线称为经验频率曲线，而把由数学方程式所表示的频率曲线称为理论频率曲线。所谓水文频率分布线型是指所采用的理论频率曲线(频率函数)的型式(水文中常用线型为正态分布型、极值分布型、皮尔逊Ⅲ型分布型等)，它的选择主要取决于与大多数水文资料的经验频率点据的配合情况。分布线型的选择与统计参数的估算，一起构成了频率计算的两大内容。

2. Pearson Type III 频率曲线2.1 皮尔逊Ⅲ型曲线的概率密度函数    皮尔逊Ⅲ型曲线是一条一端有限一端无限的不对称单峰、正偏曲线(见图4-4-3)，数学上常称伽玛分布，其概率密度函数为:

式中：Γ(α)―α的伽玛函数；α、β、a0―分别为皮尔逊Ⅲ型分布的形状尺度和位置未知参数，α﹥0，β﹥0 。

显然，三个参数确定以后，该密度函数随之可以确定。可以推论，这三个参数与总体三个参数

、Cv、CS具有如下关系：(4-4-9)

2.2 皮尔逊Ⅲ型频率曲线及其绘制

水文计算中，一般需要求出指定频率P所相应的随机变量取值xp，也就是通过对密度曲线进行积分，即:(4-4-10)求出等于及大于xp的累积频率P值。直接由式(4-4-10)计算P值，非常麻烦，实际做法是通过变量转换，变换成下面的积分形式:(4-4-11)式(4-4-11)中被积函数只含有一个待定参数CS，其它两个参数

、Cv都包含在

中。 ，x是标准化变量，

称为离均系数。的均值为0，标准差为1。因此，只需要假定一个CS值，便可从式(4-4-11)通过积分求出

与之间的关系。对于若干个给定的CS值，

的对应数值表，已先后由美国福斯特和前苏联雷布京制作出来，见附表1″皮尔逊Ⅲ型频率曲线的离均系数

值表”。由就可以求出相应频率的x值：(4-4-12)

附表1 皮尔逊Ⅲ型频率曲线的离均系数