基于蒙特·卡洛算法的新安江模型参数优化方法与流程

技术特征：

1.一种基于蒙特·卡洛算法的新安江模型参数优化方法，其特征在于，包括构造概率模型与目标函数：

1)概率模型

公式(1)是为使用蒙特·卡洛算法对新安江模型进行参数优化而构建的概率模型；

φ_i＝f(max_i,min_i)＝(max_i-min_i)*Rnd+min_i (1)

其中：

φ_i为第i个参数，Rnd是0到1之间的随机数；公式(1)的含义是为参数生成max_i到min_i之间的随机数组。

2)目标函数

定义了相对误差目标函数与日效率系数目标函数来评价实测流量与模拟流量的拟合程度；

相对误差：由单位年实测流量与模拟流量决定，是实测雨量年流量和与模拟流量年流量和的绝对差占实测流量年流量和的比例；相对误差公式见(2)。

$f(p,q)=\frac{\left|\underset{j=1}{\overset{N}{\Σ}}{\left(\underset{i=1}{\overset{M}{\Σ}}{p}_i\right)}_j-\underset{j=1}{\overset{N}{\Σ}}{\left(\underset{i=1}{\overset{M}{\Σ}}{q}_i\right)}_j\right|}{\underset{j=1}{\overset{N}{\Σ}}{\left(\underset{i=1}{\overset{M}{\Σ}}{p}_i\right)}_j}*100---\left(2\right)$

其中：

p_i为第i日实测流量，q_i为第i日模拟流量，N为率定数据的总年数，M为率定数据的每月天数；

日效率系数：由每天实测流量与模拟流量的值决定，日效率系数公式见(3)。

$f(p,q)=1-\frac{\underset{i=1}{\overset{N}{\Σ}}{(p_i-q_i)}^2}{\underset{i=1}{\overset{N}{\Σ}}{(p_i-k)}^2}---\left(3\right)$

其中：

p_i为第i日实测流量，q_i为第i日模拟流量，k为常数；k常数固定时，f(p,q)值越大，则参数优化的值越优。