一种数据采集点部署方法及系统

1.一种数据采集点部署方法，其特征在于，所述方法包括：
确定待检测水域的网格覆盖率和水质特征非均匀分布环境场的标量场重构误差；
结合所述网格覆盖率和所述标量场重构误差，构建适应度函数；
根据适应度函数，利用粒子群优化算法与引力搜索算法对采样点位置集合进行寻优处
理，确定采样点位置集合的最优解及所述最优解对应的适应度值；
按照所述采样点位置集合的最优解对应的采样点位置在所述待检测水域部署数据采
样点；
所述确定水质特征非均匀分布环境场的标量场重构误差，具体包括：
确定误差分析的位置点(x,y)个数为I；
获取所述水质特征非均匀分布环境场中任意的所述位置点的温度真实值Z和与所述位
置点的温度真实值Z对应的估计值Z'；
利用公式 计算所述标量场重构误差；
所述构建适应度函数具体包括：
以所述网格覆盖率的最大化和所述标量场重构误差的最小化为目标，构建适应度函
数；所述适应度函数为：fitness＝a*(1‑f1)+b*f2；
其中，a，b为常数因子，f1为所述网格覆盖率，f2为所述标量场重构误差。
2.根据权利要求1所述的数据采集点部署方法，其特征在于，所述确定待检测水域的网
格覆盖率具体包括：
将所述待检测水域中二维平面S区域离散化为等单位边长的网格；所述网格的数量为
C；
在所述网格上随机部署D个采样点；每个所述采样点的位置为Xd＝(xd,yd),d＝1,2,…,
D；
采用布尔感知模型判断所述采样点所在的网格是否被覆盖，得到判断结果；所述布尔
感知模型为：
其中，Dis(Xd,p)≤RS时，P(Xd,p)＝1，表示所述采样点
所在的网格被覆盖；Dis(Xd,p)＞RS时，P(Xd,p)＝0，表示所述采样点所在的网格未被覆盖；
Dis(Xd,p)为所述采样点d与任意所述网格的中心点p(x,y)的欧几里得距离：
Rs为监测点的有效感知半径；
根据所述判断结果确定被覆盖的网格数目为Cs；
按照公式 计算所述网格覆盖率。
3.根据权利要求1所述的数据采集点部署方法，其特征在于，所述利用所述粒子群优化
算法与引力搜索算法对采样点位置集合进行寻优处理，确定采样点位置集合的最优解及所
述最优解对应的适应度值，具体包括：
(1)初始化粒子群中每个粒子的位置和速度，所述位置满足公式
所述速度满足公式
所述粒子代表一组采样点位置集合，即 其中Xid＝(xid,yid)，i＝1,
2,…,N，d＝1,2,…,D，其中，Xid表示采样点d的位置， 表示粒子i的位置，Vid表示采样点d
的速度，表示粒子i的速度；
(2)确定最大迭代次数Gmax；
(3)根据所述适应度函数计算所述粒子的适应度值
(4)比较所述粒子群中的所有所述粒子的适应度值的大小，确定所述适应度值的最小
值，以所述适应度值的最小值对应的所述粒子作为所述采样点位置集合的最优解；以所述
适应度值的最小值做为所述最优解对应的适应度值
(5)判断迭代次数是否达到最大迭代次数，若是，输出所述采样点位置集合的最优解及
所述最优解对应的适应度值 若否，利用公式
更新每一所述粒子的所述速
度和所述位置，并返回步骤(3)；
其中， 表示粒子i的速度矩阵的转置， 表示粒子i的更新后的速度矩阵的转置，
表示粒子i的位置矩阵的转置， 表示粒子i的更新后的位置矩阵的转置，c表示学习因
子，randi表示[0,1]的均匀随机数，Lbest为寻优过程中的全局最优解，w为惯性权重；
所述惯性权重的计算公式为：
其中wmax为所述惯性权重的最大值、wmin为所
述惯性权重的最小值；f表示粒子的适应度函数值，favg表示所有粒子的平均适应度函数值，
fmin表示所有粒子的最小适应度函数值。
4.一种数据采集点部署系统，其特征在于，所述系统包括：
网格覆盖率和标量场重构误差确定单元，用于确定待检测水域的网格覆盖率和水质特
征非均匀分布环境场的标量场重构误差；
适应度函数构建单元，用于结合所述网格覆盖率和所述标量场重构误差，构建适应度
函数；
寻优单元，用于根据适应度函数，利用粒子群优化算法与引力搜索算法对采样点位置
集合进行寻优处理，确定采样点位置集合的最优解及所述最优解对应的适应度值；
部署单元，用于按照所述采样点位置集合的最优解对应的采样点位置在所述待检测水
域部署数据采样点；
所述网格覆盖率和标量场重构误差确定单元还包括标量场重构误差确定子单元，所述
标量场重构误差确定子单元用于确定水质特征非均匀分布环境场的标量场重构误差，所述
标量场重构误差确定子单元具体包括：
误差分析点确定模块，用于确定误差分析的位置点(x,y)个数为I；
温度获取模块，用于获取所述水质特征非均匀分布环境场中任意的所述位置点的温度
真实值Z和与所述位置点的温度真实值Z对应的估计值Z'；
标量场重构误差计算模块，用于利用公式 计算所述标量场重构误差；
所述适应度函数构建单元用于构建适应度函数，所述适应度函数构建单元具体包括：
适应度函数构建子单元，用于以所述网格覆盖率的最大化和所述标量场重构误差的最
小化为目标，构建适应度函数；所述适应度函数为：fitness＝a*(1‑f1)+b*f2；
其中，a，b为常数因子，f1为所述网格覆盖率，f2为所述标量场重构误差。
5.根据权利要求4所述的数据采集点部署系统，其特征在于，所述网格覆盖率和标量场
重构误差确定单元包括网格覆盖率确定子单元，所述网格覆盖率确定子单元，用于确定待
检测水域的网格覆盖率；所述网格覆盖率确定子单元具体包括：
网格划分模块，用于将所述待检测水域中二维平面S区域离散化为等单位边长的网格；
所述网格的数量为C；
随机部署模块，用于在所述网格上随机部署D个采样点；每个所述采样点的位置为Xd＝
(xd,yd),d＝1,2,…,D；
布尔感知模型判断模块，用于采用布尔感知模型判断所述采样点所在的网格是否被覆
盖，得到判断结果；所述布尔感知模型为：
其中，Dis(Xd,p)≤RS时，P(Xd,p)＝1，表示所述采样点
所在的网格被覆盖；Dis(Xd,p)＞RS时，P(Xd,p)＝0，表示所述采样点所在的网格未被覆盖；
Dis(Xd,p)为所述采样点d与任意所述网格的中心点p＝(x,y)的欧几里得距离：
Rs为监测点的有效感知半径；
被覆盖网格数目确定模块，用于根据所述判断结果确定被覆盖的网格数目为Cs；
网格覆盖率计算模块，用于按照公式 计算所述网格覆盖率。
6.根据权利要求4所述的数据采集点部署系统，其特征在于，所述寻优单元利用所述粒
子群优化算法与引力搜索算法对采样点位置集合进行寻优处理，确定采样点位置集合的最
优解及所述最优解对应的适应度值，所述寻优单元具体包括：
粒子初始化子单元，用于初始化粒子群中每个粒子的位置和速度，所述位置满足公式
所述速度满足公式
所述粒子代表一组采样点位置集合，即
其中Xid＝(xid,yid)，i＝1,2,…,N，d＝1,2,…,D，其中，Xid表示采样点d
的位置， 表示粒子i的位置，Vid表示采样点d的速度， 表示粒子i的速度；
迭代次数确定子单元，用于确定最大迭代次数Gmax；
适应度函数计子单元，用于根据所述适应度函数计算所述粒子的适应度值
适应度值比较子单元，用于比较所述粒子群中的所有所述粒子的适应度值的大小，确
定所述适应度值的最小值，以所述适应度值的最小值对应的所述粒子作为所述采样点位置
集合的最优解；以所述适应度值的最小值做为所述最优解对应的适应度值
迭代次数判断子单元，用于判断迭代次数是否达到最大迭代次数，若是，输出所述采样
点位置集合的最优解及所述最优解对应的适应度值 若否，利用公式
更新每一所述粒子的所述速
度和所述位置，并返回所述适应度函数计子单元；
其中， 表示粒子i的速度矩阵的转置， 表示粒子i的更新后的速度矩阵的转置，
表示粒子i的位置矩阵的转置， 表示粒子i的更新后的位置矩阵的转置，c表示学习因
子，randi表示[0,1]的均匀随机数，Lbest为寻优过程中的全局最优解，w为惯性权重；
所述惯性权重的计算公式为：
其中wmax为所述惯性权重的最大值、wmin为所
述惯性权重的最小值；f表示粒子的适应度函数值，favg表示所有粒子的平均适应度函数值，
fmin表示所有粒子的最小适应度函数值。

一种数据采集点部署方法及系统\n技术领域\n[0001] 本发明涉及水质监测领域，特别是涉及一种数据采集点部署方法及系统。\n背景技术\n[0002] 在对水环境场进行采样点部署的时候，通常以考虑监测覆盖率的方式进行采样点\n的均匀部署，即构建基于覆盖率的代价函数，利用粒子群优化算法寻优方式进行迭代搜索，\n获得最优的采样点部署结构。这种只考虑覆盖率而均匀部署采样点的方式仅仅适用于水环\n境特征均匀分布的情况，而且粒子群算法具有收敛速度慢、易陷入局部最优等缺点，对于水\n环境特征为高斯分布的情况下，通过这种方法获得的数据信息不足以反映整个区域的水质\n状态。\n发明内容\n[0003] 本发明的目的是提供一种数据采样点部署方法及系统，使得依据其采样值能够较\n好地反映整个水环境监测区域的水质状态。\n[0004] 为实现上述目的，本发明提供了如下方案：\n[0005] 一种数据采集点部署方法，其特征在于，所述方法包括：\n[0006] 确定待检测水域的网格覆盖率和水质特征非均匀分布环境场的标量场重构误差；\n[0007] 结合所述网格覆盖率和所述标量场重构误差，构建适应度函数；\n[0008] 根据适应度函数，利用粒子群优化算法与引力搜索算法对采样点位置集合进行寻\n优处理，确定采样点位置集合的最优解及所述最优解对应的适应度值；\n[0009] 按照所述采样点位置集合的最优解对应的采样点位置在所述待检测水域部署数\n据采样点。\n[0010] 可选的，所述确定待检测水域的网格覆盖率具体包括：\n[0011] 将所述待监测水域中二维平面S区域离散化为等单位边长的网格；所述网格的数\n量为C；\n[0012] 在所述网格上随机部署D个采样点；每个所述采样点的位置为Xd＝(xd,yd),d＝1,\n2,…,D；\n[0013] 采用布尔感知模型判断所述采样点所在的网格是否被覆盖，得到判断结果；所述\n布尔感知模型为：\n[0014] 其中，Dis(Xd,p)≤RS时，P(Xd,p)＝1，表示所述采\n样点所在的网格被覆盖；Dis(Xd,p)＞RS时，P(Xd,p)＝0，表示所述采样点所在的网格未被覆\n盖；Dis(Xd,p)为所述采样点d与任意所述网格的中心点p(x,y)的欧几里得距离：\nRs为监测点的有效感知半径；\n[0015] 根据所述判断结果确定被覆盖的网格数目为Cs；\n[0016] 按照公式 计算所述网格覆盖率。\n[0017] 可选的，所述确定水质特征非均匀分布环境场的标量场重构误差，具体包括：\n[0018] 确定误差分析的位置点(x,y)个数为I；\n[0019] 获取所述水质特征非均匀分布环境场中任意的所述位置点的温度真实值Z和与所\n述位置点的温度真实值Z对应的估计值Z'；\n[0020] 利用公式 计算所述标量场重构误差。\n[0021] 可选的，所述构建适应度函数具体包括：\n[0022] 以所述网格覆盖率的最大化和所述标量场重构误差的最小化为目标，构建适应度\n函数；所述适应度函数为：fitness＝a*(1‑f1)+b*f2；\n[0023] 其中，a，b为常数因子，f1为所述网格覆盖率，f2为所述标量场重构误差。\n[0024] 可选的，利用所述粒子群优化算法与引力搜索算法对采样点位置集合进行寻优处\n理，确定采样点位置集合的最优解及所述最优解对应的适应度值，具体包括：\n[0025] (1)初始化粒子群中每个粒子的位置和速度，所述位置满足公式\n[0026] 所述速度满足公式\n所述粒子代表一组采样点位置集合，即\n其中Xid＝(xid,yid)，i＝1,2,…,N，d＝1,2,…,D，其中，Xid表示采样点d的位置， 表示粒子\ni的位置，Vid表示采样点d的速度， 表示粒子i的速度；\n[0027] (2)确定最大迭代次数Gmax；\n[0028] (3)根据所述适应度函数计算所述粒子的适应度值\n[0029] (4)比较所述粒子群中的所有所述粒子的适应度值的大小，确定所述适应度值的\n最小值，以所述适应度值的最小值对应的所述粒子作为所述采样点位置集合的最优解；以\n所述适应度值的最小值做为所述最优解对应的适应度值\n[0030] (5)判断迭代次数是否达到最大迭代次数，若是，输出所述采样点位置集合的最优\n解及所述最优解对应的适应度值 若否，利用公式\n更新每一所述粒子的所述速\n度和所述位置，并返回步骤(3)；\n[0031] 其中， 表示粒子i的速度矩阵的转置， 表示粒子i的更新后的速度矩阵的转\n置， 表示粒子i的位置矩阵的转置， 表示粒子i的更新后的位置矩阵的转置，c表示学\n习因子，randi表示[0,1]的均匀随机数，Lbest为寻优过程中的全局最优解，w为惯性权重；\n[0032] 所述惯性权重的计算公式为：\n为所述惯性权重的最大值、wmin为所述惯性权重的最小值；f表示粒子的适应度函数值，favg\n表示所有粒子的平均适应度函数值，fmin表示所有粒子的最小适应度函数值。\n[0033] 一种数据采集点部署系统，其特征在于，所述系统包括：\n[0034] 网格覆盖率和标量场重构误差确定单元，用于确定待检测水域的网格覆盖率和水\n质特征非均匀分布环境场的标量场重构误差；\n[0035] 适应度函数构建单元，用于结合所述网格覆盖率和所述标量场重构误差，构建适\n应度函数；\n[0036] 寻优单元，用于根据适应度函数，利用粒子群优化算法与引力搜索算法对采样点\n位置集合进行寻优处理，确定采样点位置集合的最优解及所述最优解对应的适应度值；\n[0037] 部署单元，用于按照所述采样点位置集合的最优解对应的采样点位置在所述待检\n测水域部署数据采样点。\n[0038] 可选的，所述网格覆盖率和标量场重构误差确定单元包括网格覆盖率确定子单\n元，所述网格覆盖率确定子单元，用于确定待检测水域的网格覆盖率；所述网格覆盖率确定\n子单元具体包括：\n[0039] 网格划分模块，用于将所述待监测水域中二维平面S区域离散化为等单位边长的\n网格；所述网格的数量为C；\n[0040] 随机部署模块，用于在所述网格上随机部署D个采样点；每个所述采样点的位置为\nXd＝(xd,yd),d＝1,2,…,D；\n[0041] 布尔感知模型判断模块，用于采用布尔感知模型判断所述采样点所在的网格是否\n被覆盖，得到判断结果；所述布尔感知模型为：\n[0042] 其中，Dis(Xd,p)≤RS时，P(Xd,p)＝1，表示所述采\n样点所在的网格被覆盖；Dis(Xd,p)＞RS时，P(Xd,p)＝0，表示所述采样点所在的网格未被覆\n盖；Dis(Xd,p)为所述采样点d与任意所述网格的中心点p＝(x,y)的欧几里得距离：\nRs为监测点的有效感知半径；\n[0043] 被覆盖网格数目确定模块，用于根据所述判断结果确定被覆盖的网格数目为Cs；\n[0044] 网格覆盖率计算模块，用于按照公式 计算所述网格覆盖率。\n[0045] 可选的，所述网格覆盖率和标量场重构误差确定单元还包括标量场重构误差确定\n子单元，所述标量场重构误差确定子单元用于确定水质特征非均匀分布环境场的标量场重\n构误差，所述标量场重构误差确定子单元具体包括：\n[0046] 误差分析点确定模块，用于确定误差分析的位置点(x,y)个数为I；\n[0047] 温度获取模块，用于获取所述水质特征非均匀分布环境场中任意的所述位置点的\n温度真实值Z和与所述位置点的温度真实值Z对应的估计值Z'；\n[0048] 标量场重构误差计算模块，用于利用公式 计算所述标量场重构误\n差。\n[0049] 可选的，所述适应度函数构建单元用于构建适应度函数，所述适应度函数构建单\n元具体包括：\n[0050] 适应度函数构建子单元，用于以所述网格覆盖率的最大化和所述标量场重构误差\n的最小化为目标，构建适应度函数；所述适应度函数为：fitness＝a*(1‑f1)+b*f2；\n[0051] 其中，a，b为常数因子，f1为所述网格覆盖率，f2为所述标量场重构误差。\n[0052] 可选的，所述寻优单元利用所述粒子群优化算法与引力搜索算法对采样点位置集\n合进行寻优处理，确定采样点位置集合的最优解及所述最优解对应的适应度值，所述寻优\n单元具体包括：\n[0053] 粒子初始化子单元，用于初始化粒子群中每个粒子的位置和速度，所述位置满足\n公式 所述速度满足公式\n[0054] 所述粒子代表一组采样点位置集合，即\n其中Xid＝(xid,yid)，i＝1,2,…,N，d＝1,2,…,D，其中，Xid表示采样点\nd的位置， 表示粒子i的位置，Vid表示采样点d的速度，表示粒子i的速度；\n[0055] 迭代次数确定子单元，用于确定最大迭代次数Gmax；\n[0056] 适应度函数计算子单元，用于根据所述适应度函数计算所述粒子的适应度值\n[0057] 适应度值比较子单元，用于比较所述粒子群中的所有所述粒子的适应度值的大\n小，确定所述适应度值的最小值，以所述适应度值的最小值对应的所述粒子作为所述采样\n点位置集合的最优解；以所述适应度值的最小值做为所述最优解对应的适应度值\n[0058] 迭代次数判断子单元，用于判断迭代次数是否达到最大迭代次数，若是，输出所述\n采样点位置集合的最优解及所述最优解对应的适应度值 若否，利用公式\n更新每一所述粒子的所述速\n度和所述位置，并返回步骤(3)；\n[0059] 其中， 表示粒子i的速度矩阵的转置， 表示粒子i的更新后的速度矩阵的转\n置， 表示粒子i的位置矩阵的转置， 表示粒子i的更新后的位置矩阵的转置，c表示学\n习因子，randi表示[0,1]的均匀随机数，Lbest为寻优过程中的全局最优解，w为惯性权重；\n[0060] 所述惯性权重的计算公式为： 其中wmax为\n所述惯性权重的最大值、wmin为所述惯性权重的最小值；f表示粒子的适应度函数值，favg表\n示所有粒子的平均适应度函数值，fmin表示所有粒子的最小适应度函数值。\n[0061] 根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明在构建适应\n度函数时，将同时考虑覆盖性能和能反映采样数据优劣的重构性能，使得采样点达到对环\n境标量场的最大化覆盖，在搜索最优采样点时，将传统粒子群优化算法与引力搜索算法结\n合，对采样点位置集合进行寻优处理，确定采样点位置集合的最优解及所述最优解对应的\n适应度值，以提高算法收敛速度和精度，使得依据其采样值能够较好地反映整个水环境监\n测区域的水质状态。\n附图说明\n[0062] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所\n需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施\n例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获\n得其他的附图。\n[0063] 图1为本发明的实施例1提供的数据采集点部署方法的流程图；\n[0064] 图2为本发明的实施例2提供的数据采集点部署系统的结构框图。\n具体实施方式\n[0065] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完\n整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于\n本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他\n实施例，都属于本发明保护的范围。\n[0066] 本发明的目的是提供一种数据采集点部署方法及系统，依据其采样值能够较好地\n反映整个水环境监测区域的水质状态。\n[0067] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实\n施方式对本发明作进一步详细的说明。\n[0068] 实施例1\n[0069] 如图1所示，本实施例提供的数据采集点部署方法包括：\n[0070] 步骤101：确定待检测水域的网格覆盖率和水质特征非均匀分布环境场的标量场\n重构误差；\n[0071] 在本实施例中，具体采用以下步骤确定待检测水域的网格覆盖率和水质特征非均\n匀分布环境场的标量场重构误差。\n[0072] 所述确定待检测水域的网格覆盖率具体包括：\n[0073] 将所述待监测水域中二维平面S区域离散化为等单位边长的网格；所述网格的数\n量为C；\n[0074] 在所述网格上随机部署D个采样点；每个所述采样点的位置为Xd＝(xd,yd),d＝1,\n2,…,D；\n[0075] 采用布尔感知模型判断所述采样点所在的网格是否被覆盖，得到判断结果；所述\n布尔感知模型为：\n[0076] 其中，Dis(Xd,p)≤RS时，P(Xd,p)＝1，表示所述采\n样点所在的网格被覆盖；Dis(Xd,p)＞RS时，P(Xd,p)＝0，表示所述采样点所在的网格未被覆\n盖；Dis(Xd,p)为所述采样点d与任意所述网格的中心点p(x,y)的欧几里得距离：\nRs为监测点的有效感知半径；\n[0077] 根据所述判断结果确定被覆盖的网格数目为Cs；\n[0078] 按照公式 计算所述网格覆盖率。\n[0079] 根据布尔感知模型确定的网格覆盖率，能够更加准确地反映实际覆盖情况，可以\n直接作为水环境检测覆盖率的度量指标。\n[0080] 所述确定水质特征非均匀分布环境场的标量场重构误差，具体包括：\n[0081] 确定误差分析的位置点(x,y)个数为I；\n[0082] 获取所述水质特征非均匀分布环境场中任意的所述位置点的温度真实值Z和与所\n述位置点的温度真实值Z对应的估计值Z'；\n[0083] 利用公式 计算所述标量场重构误差。\n[0084] 针对水质特征非均匀分布环境场，将较多的采样位置部署在变化较大、特征明显\n的区域中，而在变化平缓的区域中部署相对少的采样点，这相比于完全均匀部署的方式能\n够获得更好的环境标量场重构性能。\n[0085] 步骤102：结合所述网格覆盖率和所述标量场重构误差，构建适应度函数；\n[0086] 所述构建适应度函数具体包括：\n[0087] 以所述网格覆盖率的最大化和所述标量场重构误差的最小化为目标，构建适应度\n函数；所述适应度函数为：fitness＝a*(1‑f1)+b*f2；\n[0088] 其中，a，b为常数因子，f1为所述网格覆盖率，f2为所述标量场重构误差，1‑f1将覆\n盖率的最大化转化为覆盖空洞的最小化，使得寻优目标变成求解函数最小值的问题。\n[0089] 由适应度函数建立了基于最大化覆盖率和最小化标量场的重构误差的水环境检\n测性能衡量准则，使得采样点部署能达到对环境标量场的最大化覆盖。\n[0090] 步骤103：根据适应度函数，利用粒子群优化算法与引力搜索算法对采样点位置集\n合进行寻优处理，确定采样点位置集合的最优解及所述最优解对应的适应度值；\n[0091] 所述利用所述粒子群优化算法与引力搜索算法对采样点位置集合进行寻优处理，\n确定采样点位置集合的最优解及所述最优解对应的适应度值，具体包括：\n[0092] (1)初始化粒子群中每个粒子的位置和速度，所述位置满足公式\n[0093] 所述速度满足公式\n[0094] 所述粒子代表一组采样点位置集合，即 其中Xid＝(xid,yid)，i\n＝1,2,…,N，d＝1,2,…,D，其中，Xid表示采样点d的位置， 表示粒子i的位置，Vid表示采样\n点d的速度， 表示粒子i的速度；\n[0095] (2)确定最大迭代次数Gmax；\n[0096] (3)根据所述适应度函数计算所述粒子的适应度值\n[0097] (4)比较所述粒子群中的所有所述粒子的适应度值的大小，确定所述适应度值的\n最小值，以所述适应度值的最小值对应的所述粒子作为所述采样点位置集合的最优解；以\n所述适应度值的最小值做为所述最优解对应的适应度值\n[0098] (5)判断迭代次数是否达到最大迭代次数，若是，输出所述采样点位置集合的最优\n解及所述最优解对应的适应度值 若否，利用公式\n更新每一所述粒子的所述速\n度和所述位置，并返回步骤(3)；\n[0099] 其中， 表示粒子i的速度矩阵的转置， 表示粒子i的更新后的速度矩阵的转\n置， 表示粒子i的位置矩阵的转置， 表示粒子i的更新后的位置矩阵的转置，c表示学\n习因子，randi表示[0,1]的均匀随机数，Lbest为寻优过程中的全局最优解，w为惯性权重；\n[0100] 所述惯性权重，能够体现粒子继承当前速度的能力，为了平衡算法的全局搜索能\n力和局部改良能力，所述惯性权重采用自适应权重，计算公式为：\n[0101] 其中wmax为所述惯性权重的最大值、wmin\n为所述惯性权重的最小值；f表示粒子的适应度函数值，favg表示所有粒子的平均适应度函\n数值，fmin表示所有粒子的最小适应度函数值。\n[0102] 利用粒子群优化算法与引力搜索算法对采样点位置集合进行寻优处理，确定采样\n点位置集合的最优解及所述最优解对应的适应度值，实现对水域的最大化覆盖监测的同\n时，能够较好的重构整个水域的水质分布特征，从而较好的反映整个水环境监测区域的水\n质状态。\n[0103] 步骤104：按照所述采样点位置集合的最优解对应的采样点位置在所述待检测水\n域部署数据采样点。\n[0104] 本申请提出的数据采集点部署方法，通过网格覆盖率和标量场重构误差两者相结\n合，作为水环境测试性能衡量标准，建立与衡量标准相关的适应度函数，利用寻优，确定采\n样点位置集合的最优解及所述最优解对应的适应度值，实现对水域的最大化覆盖监测的同\n时，能够较好的重构整个水域的水质分布特征，从而较好的反映整个水环境监测区域的水\n质状态。\n[0105] 实施例2\n[0106] 如图2所示，本实施例提供的数据采集点部署系统包括：网格覆盖率和标量场重构\n误差确定单元201、适应度函数构建单元202、寻优单元203和部署单元204。\n[0107] 网格覆盖率和标量场重构误差确定单元201，用于确定待检测水域的网格覆盖率\n和水质特征非均匀分布环境场的标量场重构误差。\n[0108] 所述网格覆盖率和标量场重构误差确定单元201包括网格覆盖率确定子单元和标\n量场重构误差确定子单元；\n[0109] 所述网格覆盖率确定子单元具体包括：\n[0110] 网格划分模块，用于将所述待监测水域中二维平面S区域离散化为等单位边长的\n网格；所述网格的数量为C；\n[0111] 随机部署模块，用于在所述网格上随机部署D个采样点；每个所述采样点的位置为\nXd＝(xd,yd),d＝1,2,…,D；\n[0112] 布尔感知模型判断模块，用于采用布尔感知模型判断所述采样点所在的网格是否\n被覆盖，得到判断结果；所述布尔感知模型为：\n[0113] 其中，Dis(Xd,p)≤RS时，P(Xd,p)＝1，表示所述采\n样点所在的网格被覆盖；Dis(Xd,p)＞RS时，P(Xd,p)＝0，表示所述采样点所在的网格未被覆\n盖；Dis(Xd,p)为所述采样点d与任意所述网格的中心点p＝(x,y)的欧几里得距离：\nRs为监测点的有效感知半径；\n[0114] 被覆盖网格数目确定模块，用于根据所述判断结果确定被覆盖的网格数目为Cs；\n[0115] 网格覆盖率计算模块，用于按照公式 计算所述网格覆盖率。\n[0116] 根据布尔感知模型确定的网格覆盖率，能够更加准确地反映实际覆盖情况，可以\n直接作为水环境检测覆盖率的度量指标。\n[0117] 所述标量场重构误差确定子单元具体包括：\n[0118] 误差分析点确定模块，用于确定误差分析的位置点(x,y)个数为I；\n[0119] 温度获取模块，用于获取所述水质特征非均匀分布环境场中任意的所述位置点的\n温度真实值Z和与所述位置点的温度真实值Z对应的估计值Z'；\n[0120] 标量场重构误差计算模块，用于利用公式 计算所述标量场重构误\n差。\n[0121] 针对水质特征非均匀分布环境场，将较多的采样位置部署在变化较大、特征明显\n的区域中，而在变化平缓的区域中部署相对少的采样点，这相比于完全均匀部署的方式能\n够获得更好的环境标量场重构性能。\n[0122] 适应度函数构建单元202，用于结合所述网格覆盖率和所述标量场重构误差，构建\n适应度函数。所述适应度函数构建单元202具体包括：\n[0123] 适应度函数构建子单元，用于以所述网格覆盖率的最大化和所述标量场重构误差\n的最小化为目标，构建适应度函数；所述适应度函数为：fitness＝a*(1‑f1)+b*f2；\n[0124] 其中，a，b为常数因子，f1为所述网格覆盖率，f2为所述标量场重构误差，1‑f1将覆\n盖率的最大化转化为覆盖空洞的最小化，使得寻优目标变成求解函数最小值的问题。\n[0125] 由适应度函数建立了基于最大化覆盖率和最小化标量场的重构误差的水环境检\n测性能衡量准则，使得采样点部署能达到对环境标量场的最大化覆盖。\n[0126] 寻优单元203，用于根据适应度函数，利用粒子群优化算法与引力搜索算法对采样\n点位置集合进行寻优处理，确定采样点位置集合的最优解及所述最优解对应的适应度值。\n所述寻优单元203具体包括：\n[0127] 粒子初始化子单元，用于初始化粒子群中每个粒子的位置和速度，所述位置满足\n公式 所述速度满足公式\n[0128]\n[0129] 所述粒子代表一组采样点位置集合，即 其中Xid＝(xid,yid)，i\n＝1,2,…,N，d＝1,2,…,D，其中，Xid表示采样点d的位置， 表示粒子i的位置，Vid表示采样\n点d的速度， 表示粒子i的速度；\n[0130] 迭代次数确定子单元，用于确定最大迭代次数Gmax；\n[0131] 适应度函数计算子单元，用于根据所述适应度函数计算所述粒子的适应度值\n[0132] 适应度值比较子单元，用于比较所述粒子群中的所有所述粒子的适应度值的大\n小，确定所述适应度值的最小值，以所述适应度值的最小值对应的所述粒子作为所述采样\n点位置集合的最优解；以所述适应度值的最小值做为所述最优解对应的适应度值\n[0133] 迭代次数判断子单元，用于判断迭代次数是否达到最大迭代次数，若是，输出所述\n采样点位置集合的最优解及所述最优解对应的适应度值 若否，利用公式\n更新每一所述粒子的所述速\n度和所述位置，并返回所述适应度函数计算子单元；\n[0134] 其中， 表示粒子i的速度矩阵的转置， 表示粒子i的更新后的速度矩阵的转\n置， 表示粒子i的位置矩阵的转置， 表示粒子i的更新后的位置矩阵的转置，c表示学\n习因子，randi表示[0,1]的均匀随机数，Lbest为寻优过程中的全局最优解，w为惯性权重；\n[0135] 所述惯性权重，能够体现粒子继承当前速度的能力，为了平衡算法的全局搜索能\n力和局部改良能力，所述惯性权重采用自适应权重，计算公式为：\n[0136] 其中wmax为所述惯性权重的最大值、wmin\n为所述惯性权重的最小值；f表示粒子的适应度函数值，favg表示所有粒子的平均适应度函\n数值，fmin表示所有粒子的最小适应度函数值。\n[0137] 利用粒子群优化算法与引力搜索算法对采样点位置集合进行寻优处理，确定采样\n点位置集合的最优解及所述最优解对应的适应度值，实现对水域的最大化覆盖监测的同\n时，能够较好的重构整个水域的水质分布特征，从而较好的反映整个水环境监测区域的水\n质状态。\n[0138] 部署单元204，用于按照所述采样点位置集合的最优解对应的采样点位置在所述\n待检测水域部署数据采样点。\n[0139] 本申请提出的数据采集点部署系统，通过网格覆盖率确定子单元和标量场重构误\n差确定子单元两者相结合，作为水环境测试性能衡量标准，建立与衡量标准相关的适应度\n函数构建单元，利用寻优单元，确定采样点位置集合的最优解及所述最优解对应的适应度\n值，实现对水域的最大化覆盖监测的同时，能够较好的重构整个水域的水质分布特征，从而\n较好的反映整个水环境监测区域的水质状态。\n[0140] 对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比\n较简单，相关之处参见方法部分说明即可。\n[0141] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说\n明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据\n本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不\n应理解为对本发明的限制。