Weibull分布参数的粒子群算法估计.pdfVIP

第42卷第6期 中国海洋大学学报 42(6)：120～125 2012年6月 PERl0DICALOFOCEANUNIVERSITYOFCHINA Jum，2012 Weibull分布参数的粒子群算法估计’ 董 胜1，韩 意1，陶山山1，樊敦秋2 (1．中围海洋大学工程学院，山东青岛266100；2．胜利石油管理局钻井T22研究院，山东东营，257017) 摘要： 基于广西涠洲岛海洋监测站3个方向的年极值波高观测资料，在假设其服从Weibull分布的基础上，运用最小 中4种方法得出的拟合结果及运算效率进行比较分析，说明了粒子群算法在估计极值分布参数中的优势。 关键词：Weibull分布；参数估计；粒子群算法；波高 中图法分类号： U697．1 文献标志码：A 文章编号： 1672—5174(2012)06—120-06 在海洋技术、水利工程和城市建筑设计中， 求出初始解，再运用基于粒子群算法的极大似然估计 Weibull分布以其较强的适应性和灵活性被广泛地运法迭代求解。 用于多年一遇极值风速、波高、水位、降水量等的长期 上述成果说明粒子群算法正在被广泛应用于实际 预测中[1]。目前，Weibull分布的参数估计方法有矩工程中。本文将粒子群算法引入对Weibull分布的参 法、极大似然法、概率权重矩法、相关系数法及Bayes 估计法等fl伽。严晓东等人[41对三参数Weibull分布的 法，即最小二乘法、矩估计法和最速下降法，与新引入 多种估计方法进行比较研究，并指出各种方法的适用 的粒子群算法．以广西洞洲岛海洋监测站3个方向的 范围。极大似然法需要迭代求解联立的3个超越方 年极值波高观测数据为例。对Weibull分布的参数进 程，计算相当复杂is3，而且对初值的要求比较高心]；概 行了估计，并对各种方法所得结果进行了比较分析。 率权重矩法(PWM)不适合于子样较少的情况，当子样 少于10个时，计算精度很差口1；相关系数法的计算比 1 Weibull分布参数的估计方法 较复杂[5]；Bayes方法将先验分布和后验分布结合起 三参数Weibull分布的分布函数可表示如下： 来f5]，但后验分布比较复杂【3]，并且一般需要和其他方 法联合使用[5]。2001年郑红星和李丽娟将最速下降法 (1) 应用于水质模型参数的优化中[6]，并将结果与模拟退 式中：a，b和c分别表示位置参数、尺度参数和形状参 火算法和遗传算法进行比较，得出最速下降法的计算 数，且6o，c0。其对应的密度函数为 结果易受初始搜索条件的影响，寻优过程容易陷入局 部最优，而且精度低于非数值随机优化算法，但是最速 下降法的计算速度明显高于非数值随机优化算法。 Weibull分布的参数估计方法有很多，比如最小二 近几年，粒子群算法的研究引起了许多学者的关 乘法、矩估计法、最速下降法等，在引入Weibull分布 注。该算法简洁、收敛速度快等优势被广泛地应用于 的粒子群算法前，本文将简单介绍一下这些方法。以 参数的非线性拟合中。江燕掣70将其应用于新安江模 型参数的优选中。郭建青等一8]在求解河流水质参数的 布的独立同分布样本。 函数优化问题时引入了该算法，取得了良好的效果。 1．1最小二乘法 路志强等f93运用该法对暴雨强度公式参数进行了优 在式(1)中令F(z)一p，则式(1)变形为： 化。郭建青等u0]研究了含水层参数的函数优化问题， (3) 认为粒子群算法给出的结果优良。罗航等。1¨结合双线 p_1-exp[