广义极值分布线性矩法研究.pdfVIP

广义极值分布线性矩法研究 李淑会，李兴凯，许睢 河海大学，江苏南京（210098 ） 摘 要：在介绍了广义极值分布线性矩与其参数关系后，通过统计试验计算分析了线性矩在 该分布中的统计性能及历史洪水公式的适用性，并与绝对值准则和平方和准则两种不同优化 适线法、矩法作了比较。统计试验表明，线性矩法优于平方和准则适线法及矩法，其参数及 设计值的不偏性很好,与绝对值准则适线法总体上相当。整体上讲，现行句法好与其他传统 方法。 关键词：广义极值分布；线性矩；历史洪水；优化适线法；统计试验 1. 前言 Greenwood 在 1979 年定义概率权重矩之后，Hosking 于 1990 年定义了线性矩 （L-moments ），这引起了国外水文学者的高度兴趣，并作了一系列的研究工作[1][2][3]，主要 研究内容包括在简单样本、不同总体分布下，该估计方法统计性能，并与矩法比较，如何用 于地区综合，线型鉴别等。在 Hosking 研究的基础上，国内对基于 P-Ⅲ型分布、对数正态分 布的线性矩法参数估计进行了研究，并提出了适合于含历史洪水特大值的不连序样本的线性 矩估算公式 [4][5] 。但在国内对基于 GEV 分布的连序样本和非连序样本的线性矩法方面的研 究还很少见，而对评价适线法作为 GEV 分布参数估计性质也不多见，本文主要探讨线性矩 法在 GEV 分布参数估计应用。 2. GEV 的分布函数及其设计值计算公式 ⎧ 1 ⎫ ⎪ ⎡ ⎛x −µ⎞⎤k ⎪ F x P ≥x − − −k k (1) ( ) (ξ ) 1 exp⎨ ⎢1 ⎜ ⎟⎥ ⎬ ≠0 ⎪ ⎣ ⎝ α ⎠⎦ ⎪ ⎩ ⎭ 式中：µ_____位置参数；α_____ 尺度参数（α 0 ）；k_____ 形状参数。 当k 0 时，−∞x ≤µ+α/ k ; 当k ＝0 时，−∞x +∞; 当k 0 时，µ+α/ k ≤x +∞ 。 k ＝0、k0 和 k0 分别代表极值Ⅰ型、Ⅱ型和Ⅲ型分布。当 k ＝0 时，GEV 分布退化为分布 Gumbel 分布。 对任意给定的频率 p ，由p ＝F （xp ）求解与频率p 对应的设计值 xp 为 α { k } xp µ+ 1−[−ln(1−p )] （2 ） k 3. GEV 分布参数与线性矩的关系 线性矩λ 、λ 、τ 、τ 定义见文献[4]，由文献[1]可知线性矩与密度函数中参数或特征 1 2 3 4 参数关系式如下： 若µ 、k 、α 已知，则由文献[1]给出线性矩公式如下（k －1）： + −Γ +k k λ µ α[1 (1 )]/ ⎧1 ⎪ −k − Γ +k k λ α(1 2 ) (1 ) / ⎨⎪ 2