matlab中的耿贝尔法,最大降水量多年一遇计算方法及Matlab实现

内蒙古农业科技 2012(1)：75~76 Inner Mongolia Agricultural Science And Technology 最大降水量多年一遇计算方法及 Matlab实现 王 涵． 李 玲 (鄂伦春旗气象局 ，内蒙古 鄂伦春 165450) 中图分类号：Sl61．6 文献标 识码：C 10．3969／j．jssn．1007—0907．2012．O1．042 文章编号：10o7一o9O7(2012)01—0075一O2 当暴雨洪涝 等灾情 出现时 ．“几 年一遇 ”经常作为标志灾害 强度出现在媒体 口中 那么如何计算 和快速计算某一特定气象 要素值究竟是几 年一遇呢?本文通过 Matlab实例对 比分析三种 常见的多年一遇的计算方法。“几年一遇”在数学上实际上就是 概率与数理统计 中极值分布问题。这类 问题 的解决 ，在城市建设 灾害防御规划中 ．有重要的社会 和经济意义。 本文使 用 的 Ma廿ab版 本为 7．10．0 (R2010a)．以 阿里河 镇 1961—2009年 49年最大 日降水量为例 。分别通过指数分布法 、耿 贝尔分布法及皮尔逊 II1分布法讨论及编程计算。具体数据如下 ： 60．8，34,40，36．6，38．8，52．2,48．4，53．2，73．2，51．6，37．3，74．3，43．2，51,40．7， 36．7,49，55．4，70．7，42．1，106．8，72．9，32．2，83．1’47．5，58,46．3，29．6，48．8，7 1．2，66，29．7，53-2，51．9，47．7，43．4，38，26．6，58．2，37．7，58．6，34．8，42，31．8，3 8．2，33．2，63．8，53．6，30．2。 需要说明的是．本文程序使用 的数据为列矩阵 ．所 以如果输 入行矩 阵，需要做转置操作。 极值统计的根本 目的是准确的推断极值序列 的重现期值或 某一极值平均可能在 多少年 内出现一次的重现期 。当 X为其所 在序列 中最大或最小极值 x的事件 ，平均在 T年内出现 1次时 ． 则把这个 T叫做 X 的特定值 X的重现期 ．而在 T年内平均 出现 1次的这个特定值 X叫做重现期值 。如果令 F(x1为 X的分布 函 数 ，x代表 X 的年最大值或最小值 ，则根据概率分 布 ，最 大值或 最小值的再现期T(x)可分布由．T(x) 雨1 和T(x) 南给出。 这就是说 ，重现期就是右侧概率或左侧概率的倒 数。而重现期 内 的极大值或极小值可分别 由 F(x)=1一。 和 F(x) 解 出Xo即 重现期和重现期值互为 函数关系 ．可 以互相解出 ．只要知道函数 Ffx)即 可 。 1 指数分布法 对于极 大值 问题而言 ．很多气候要素 的极值分布都属指数 型的 ，指数分布形式如下 ： x=algT+b 根据求解过程，关键部分 编程如下 ： N=length0i)； laver=-mean(1i)； fori=l：N Ti(i)=(N+1)／i； end Taver=-mean(Ti)； f0r i_ 1：N 收稿 日期：2011-12—28 a1=a1+abs(1i(i)-laver) abs(1og(Ti(i))一log(Taver)； a2=a2+abs(1og(Ti(i))一log(Taver))'2； end a=al／a2； b=laver-a log(Taver)； latT=a log(T)+b； end 通过计算重现期为 1～1000的重现期值得到如下曲线(图 1)