【技术实现步骤摘要】
基于启发迁移复形竞争进化算法的水文模型参数率定方法
[0001]本专利技术涉及流域水文模型的
,特别是基于启发迁移复形竞争进化算法的水文模型参数率定方法
。
技术介绍
[0002]洗牌复形进化算法是全局寻优算法,是复形竞争进化算法搭载的外部框架
。
相较于其它全局寻优算法,其特殊的洗牌方式极大地提高了各复形内样本点的空间差异性,具有极强的广度搜索能力
。
[0003]水文模型的误差响应面非常复杂,具有高度非线性
、
不连续不可导的特性,存在大量局部最优解,全局寻优算法的深度搜索方案普遍相对简单,寻优效率和解的精度较差
。
大量实验和实际应用证明,复形进化算法采用的单纯形下山搜索法的深度搜索效果优势明显,在参数率定领域应用较广
。
但其搜索过程中仅仅使用了局部最差点和形心点,并未充分利用局部可用信息
。
[0004]因此,原版的洗牌复形进化算法在复形划分时仍有一定的改进空间,可从样本点之间空间距离的角度划分复形,使复形内各样本点之间的空间距离保持足够的差异性,扩大复形内样本点遍布样本空间的范围,增强算法的广度搜索能力,提高水文模型参数率定的寻优效果
。
[0005]申请公开号为
CN106709181A
的专利技术专利申请公开了一种基于并行编程和模块法的分布式水文模型率定方法,所述方法包括:
1)
基于
MPICH
构建的并行编程环境,耦合开源分布 ...
【技术保护点】
【技术特征摘要】 【专利技术属性】
1.
基于启发迁移复形竞争进化算法的水文模型参数率定方法,包括进行数据收集,其特征在于,还包括以下步骤:步骤1:进行复形划分;步骤2:提取单纯形,包括以下子步骤:步骤
21
:根据线性概率分布非重复提取复形样本点填充单纯形,设定单纯形样本数
q
为复形内样本数
n
的一半,根据线性概率分布非重复提取出
n
个复形样本的序列号,依据
n
个复形样本序列号提取复形样本点,填充至单纯形样本空间
X = {x
i
, f
i
, i = 1, ..., q}
;采用所述线性概率分布提取序列号的公式为,其中,
pos
为提取样本的序列号,
npg
为复形成员个数,
genrand_real1
为梅森旋转算法0‑1均匀分布随机数发生器,为向下取整符号;步骤
22
:单纯形内样本点依据目标函数评估值升序排序;步骤
23
:记录单纯形内最优点和最差点及除最差点外其余点的形心点;步骤3:引入所述单纯形的最优信息,生成单纯形搜索点;步骤4:所述单纯形搜索点根据单纯形搜索效果对单纯形最优点的引入比例作出调整;步骤5:替换单纯形最差点为单纯形搜索点,并按原始位置将单纯形样本点重新嵌入复形并按目标函数评估值升序排序复形;步骤6:将各复形填充回总样本,形成新的总样本,完成洗牌复形进化算法的一代进化;步骤7:重复执行步骤1‑
步骤6,直至满足收敛条件,所述收敛条件为总样本中最优点的目标函数评估值满足一定代数下的变化值小于设定阈值;步骤8:输出水文模型最优参数组
。2.
如权利要求1所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法,其特征在于,所述步骤1包括以下子步骤:步骤
11
:在参数可行空间内随机生成总样本;步骤
12
:依据各样本目标函数评估值升序排序;步骤
13
:以复形个数为间隔数间隔提取样本填充复形;步骤
14
:直至所有复形填充完毕,完成复形划分
。3.
如权利要求2所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法,其特征在于,总样本数
s
为复形个数
m
与复形内样本数
n
的乘积,在参数可行空间内,通过梅森旋转算法生成0‑1均匀分布随机数,结合各参数维度允许变化上下限之差的绝对值及各参数维度确定的下限值,随机生成总样本
D = {x
i
, i = 1, ..., s}
,其中,
x
i
为第
i
个样本点
。4.
如权利要求3所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法,其特征在于,所述步骤
12
包括依据各样本点
x
i
代入水文模型生成的模拟径流序列与实测或理想径流序列一同代入目标函数式计算目标函数评估值
f
i
,并对总样本
D = {x
i
, f
i
, i = 1, ..., s}
依据目标函数评估值升序排序
。5.
如权利要求4所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法,其
特征在于,所述步骤
13
包括对于复形
A
k
,从排序后的样本点
x
k
以复形个数
m
作为间隔数间隔提取
n
个样本点填充复形
A
k
,其中,
k=1,
…
,m。6.
如权利要求5所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法,其特征在于,所述步骤
13
还包括将总样本
D
划分为
m
个复形,
D = {A
k
, k = 1, ..., m}
,复形划分表示为,其中,为总样本,为复形1,为复形2,为复形
m。7.
如权利要求6所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法,其特征在于,所述步骤3包括单纯形最差点以单纯形形心点为参考点反射或收缩,并以一定比例向单纯形最优点做线性偏移,生成初步单纯形搜索点:反射点
X
ref_ini
和收缩点
X
con_ini
,公式为
X
ref_ini=
X
w
+2(X
ce
-
X
w
)
,其中,
X
w
为单纯形最差点,
X
ce
为单纯形形心点
。8.
如权利要求7所述的基于启发迁移复形竞争进化算法的水文模型参数率定方法,其特征在于,所述步骤3还包括所述初步单纯形搜索点以比例
技术研发人员:阚光远,侯宇,刘荣华,刘希琛,梁珂,
申请(专利权)人:中国水利水电科学研究院,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。