[发明专利]一种绘制P‑III型分布频率曲线的优化方法

摘要

本发明公开了一种绘制P‑III型分布频率曲线的优化方法，通过建立P‑III型分布频率曲线优化模型，对理论频率和经验频率的相对误差赋予不同的权重，以相对误差权重的平方和最小为目标函数，并根据流域概况和水库的实际情况添加参数的约束条件，以矩法估计的参数为初值，采用改进的非常快速模拟退火算法求解模型，得到参数的最优解，从而绘制出P‑III型分布频率曲线。该方法参数估计结果唯一，避免了人为因素产生的任意性，还能够照顾到历史洪水和特大洪水，拓展了P‑III型分布曲线的适用性，为绘制P‑III型分布频率曲线提供了重要的参考依据。

主权项

一种绘制P‑III型分布频率曲线的优化方法，其特征在于包括如下步骤：步骤（1），收集历史洪水和实测洪水资料，对于收集到的历史洪水，设其最大调查期为N，调查的历史洪水个数为m，实测系列的年数为n，实测系列中含特大洪水个数为l；是样本容量为n‑l+m的样本系列，并按从小到大的顺序排列，用表示系列中的第j个值；步骤（2），计算洪水序列的经验频率，并用矩法初步估计分布函数的参数；其中，所采取的经验频率公式为不连序样本的经验频率公式；对于步骤（1）中得到的样本系列，将它的经验频率记作；步骤（3），建立P‑III型分布频率曲线优化模型；将步骤（1）中的洪水数据代入步骤（2）得到的频率分布函数，求得不同的洪水点据对应的理论频率，记作，从而得到理论频率与经验频率的相对误差，对不同的洪水点据，赋予相应的权重；采用相对误差权重的平方和最小为目标函数，并根据流域概况和水库的实际情况添加参数的约束条件，建立P‑III型分布频率曲线优化模型；步骤（4），以步骤（2）中估计的参数为初值，采用改进的非常快速模拟退火算法求解步骤（3）中的模型，得到P‑III型分布函数的各项参数，从而绘制频率曲线。