PⅢ曲线的极大似然估计及应用.pdfVIP

第43卷第21期 人 民 长 江 V01．43．No．21 2012年11月 River Nov．． 2012 Yangtze 文章编号：1001—4179(2012)21一0021—03 P一Ⅲ曲线的极大似然估计及应用 余泱悦1，贺 信2 (1．武汉大学水利水电学院，湖北武汉430072；2．武汉大学经济与管理学院，湖北武汉430072) 摘要：为了使皮尔逊III型曲线(以下简称P—III型曲线)的参数估计精度更高，利用极大似然估计法对其形 状参数、尺度参数及位置参数进行估计。介绍了极大似然估计参数的数学推导过程，提供了查数值表和利用 MATLAB软件这两种方法进行参数估计的计算步骤。以黄家港1968～2008年的年平均径流观测数据为例， 进行了参数估计计算。将该计算结果与矩法估计的结果进行比较表明，极大似然估计具有更高的拟合精度。 关键词：极大似然法；P一Ⅲ型曲线；水文频率计算；MATLAB软件 中图法分类号：P333．9 文献标志码：AA 的信息。另外，极大似然法估计出来的总体分布，能够 1 研究背景 使观测值发生的概率最大化，更符合对水文预测的直 水文频率分析法始于19世纪末期，之后的1921观认识。 年，L．s．HalI设想将皮尔逊曲线族应用于水文频率曲 2极大似然法估计参数的数学推导 线的配线中，H．A．Foster将皮尔逊Ⅲ型曲线(以下简 称P一Ⅲ曲线)应用于曲线拟合中，并制成了离均系数 2．1起始点已知 表，大大简化了频率计算工作。在我国，几代水文工作 对于水文的实际观测量(径流、降雨、洪峰流量 者对我国各流域水文观测数据进行了大量的拟合和分 等)，我们将其视作一个随机变量亭。手服从，(d，卢)， 析，认为P一Ⅲ曲线能较好的拟合暴雨、洪水等系列。 即参数为理和卢的伽玛分布，起始点为n。，则其概率密 水文频率计算的核心就是对P一Ⅲ曲线的参数进行估 度函数为 计。实际水文分析计算中有很多参数估计的方法，其 中，最常用的就是采取矩法估计c，，(变差系数)、c。(偏㈣：』高∽州川e一孚x狐(邶0) 态系数)，并最终得到伽玛分布中参数的估计量。然 【o 戈。。 而R．A．Fisher证明了，伽玛分布参数的矩法估计量是 (1) 无效的。 首先假定起始点(‰)的值已知。例如，已经在多 在实际情况中发现，利用该方法进行估计时，对于 年观测中得出了某个精确的经验值。为了不失一般性， 三阶矩的估算误差较大。因此，本文主要提供一种基 进一步假设。。=0。设亭。，亭：，…，孝。为孝的几个观测值， Likelihood 于MATLAB软件的极大似然法(Maximum 亭为其平均值。极大似然法的估计思想为，找出使观测 Estimation)，直接对P一Ⅲ曲线的形状参数(理)，尺度 事件发生概率最大的参数值。首先得到联合概率密度 参数(口)，位置参数(Ⅱ。)进行估计¨1。尽管，极大似然函数，即似然方程如下： 法对于水文频率的计算结果不太灵敏，但是它与传统