发明名称 基于遗传算法多目标优化的流域传感器覆盖网优化方法
摘要 本发明提供一种基于遗传算法多目标优化的流域传感器覆盖网优化方法,其使用遗传算法多目标优化将河流流域传感器覆盖网优化问题转化为0/1多目标规划问题,通过染色体的遗传、交叉及变异等遗传算法操作,并通过监测传感器覆盖网适应度值的比较最终实现监测传感器覆盖网的合理优化选取。另外,本发明的算法采用了多个评估指标对监测传感器覆盖网适应度进行加权评估,分别是“监测节点流域特征系数”,“监测节点传感器网络传输范围系数”,“监测节点传感器使用费用系数”和“监测节点环境干扰系数”,这可以供决策者根据河流流域传感器覆盖网布局的需要,灵活调整指标权重,从而提升算法的适应性。
申请公布号 CN104270773A 申请公布日期 2015.01.07
申请号 CN201410554863.3 申请日期 2014.10.17
申请人 长江水利委员会长江科学院 发明人 李喆;谭德宝;申邵洪;张穗;陈蓓青;文雄飞;向大享
分类号 H04W16/18(2009.01)I;H04W24/02(2009.01)I;H04W84/18(2009.01)I 主分类号 H04W16/18(2009.01)I
代理机构 武汉楚天专利事务所 42113 代理人 孔敏
主权项 一种基于遗传算法多目标优化的流域传感器覆盖网优化方法,其特征在于包括如下步骤:(1)优化问题描述,针对河流流域传感器覆盖网布局的需要,假设在流域内均匀分布N个水文数据传感器采集节点,构成流域传感器覆盖网;(2)遗传算法参数初始化,并对染色体进行第一代编码,这里采用二进制随机编码方式,0表示不使用水文数据传感器采集节点,1表示使用水文数据传感器采集节点,则对于N个采集节点生成长度为N的二进制串,如下式所示:X=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N</sub>]x<sub>i</sub>={0,1}(3)计算每个染色体所对应的监测覆盖网适应度,监测覆盖网适应度函数为:f(X)=w<sub>1</sub>f<sub>1</sub>(X)+w<sub>2</sub>f<sub>2</sub>(X)‑w<sub>3</sub>f<sub>3</sub>(X)‑w<sub>4</sub>f<sub>4</sub>(X)<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><msup><mi>S</mi><mo>*</mo></msup><mo>)</mo></mrow></mrow></mfrac><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><msup><mi>S</mi><mo>*</mo></msup></mrow></munder><msup><mrow><mo>(</mo><msub><mi>v</mi><msub><mi>S</mi><mi>j</mi></msub></msub><mo>-</mo><msub><mover><mi>v</mi><mo>&OverBar;</mo></mover><msup><mi>S</mi><mo>*</mo></msup></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><msup><mi>S</mi><mo>*</mo></msup><mo>)</mo></mrow></mrow></mfrac><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><msup><mi>S</mi><mo>*</mo></msup></mrow></munder><msup><mrow><mo>(</mo><msub><mi>h</mi><msub><mi>S</mi><mi>j</mi></msub></msub><mo>-</mo><msub><mover><mi>h</mi><mo>&OverBar;</mo></mover><msup><mi>S</mi><mo>*</mo></msup></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><msup><mi>S</mi><mo>*</mo></msup><mo>)</mo></mrow></mrow></mfrac><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><msup><mi>S</mi><mo>*</mo></msup></mrow></munder><msup><mrow><mo>(</mo><msub><mi>p</mi><msub><mi>S</mi><mi>j</mi></msub></msub><mo>-</mo><msub><mover><mi>p</mi><mo>&OverBar;</mo></mover><msup><mi>S</mi><mo>*</mo></msup></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000588720850000011.GIF" wi="1597" he="151" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><msup><mi>S</mi><mo>*</mo></msup></mrow></munder><msub><mi>s</mi><msub><mi>S</mi><mi>j</mi></msub></msub><mo>/</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><mi>S</mi></mrow></munder><msub><mi>s</mi><msub><mi>S</mi><mi>j</mi></msub></msub></mrow>]]></math><img file="FDA00005887208500000111.GIF" wi="509" he="126" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>f</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><msup><mi>S</mi><mo>*</mo></msup></mrow></munder><msub><mi>c</mi><msub><mi>S</mi><mi>j</mi></msub></msub><mo>/</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><mi>S</mi></mrow></munder><msub><mi>c</mi><msub><mi>S</mi><mi>j</mi></msub></msub></mrow>]]></math><img file="FDA0000588720850000013.GIF" wi="488" he="130" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>f</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><msup><mi>S</mi><mo>*</mo></msup></mrow></munder><msub><mi>&alpha;</mi><msub><mi>S</mi><mi>j</mi></msub></msub><mo>/</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>S</mi><mi>j</mi></msub><mo>&Element;</mo><mi>S</mi></mrow></munder><msub><mi>&alpha;</mi><msub><mi>S</mi><mi>j</mi></msub></msub></mrow>]]></math><img file="FDA0000588720850000014.GIF" wi="512" he="138" /></maths>其中,f<sub>1</sub>表示监测节点流域特征系数,f<sub>2</sub>表示监测节点传感器网络传输范围系数,f<sub>3</sub>表示监测节点传感器使用费用系数,f<sub>4</sub>表示监测节点环境干扰系数,S={S<sub>1</sub>,S<sub>2</sub>,…,S<sub>N</sub>}表示监测传感器节点集合,S<sup>*</sup>={S<sub>j</sub>,|x<sub>j</sub>=1}表示S的子集,N(S<sup>*</sup>)表示S<sup>*</sup>的大小,<img file="FDA0000588720850000015.GIF" wi="234" he="88" />分别表示对应监测节点位置的河流流速、水位和水质,<img file="FDA0000588720850000016.GIF" wi="235" he="88" />分别表示包括所有监测节点的河流平均流速、平均水位流速和平均水质,<img file="FDA0000588720850000017.GIF" wi="78" he="72" />水质包括水温、PH值、电导率、溶解氧、叶绿素浓度和浊度。<img file="FDA0000588720850000018.GIF" wi="64" he="74" />表示监测节点网络数据传输范围,<img file="FDA0000588720850000019.GIF" wi="67" he="72" />表示监测节点传感器使用费用,<img file="FDA00005887208500000110.GIF" wi="78" he="77" />表示监测节点环境干扰,w<sub>1</sub>、w<sub>2</sub>、w<sub>3</sub>和w<sub>4</sub>分别表示f<sub>1</sub>、f<sub>2</sub>、f<sub>3</sub>和f<sub>4</sub>的权重;(4)对计算后监测覆盖网适应度进行排序,并按比例选取最优解;(5)按照轮盘赌方法,选择染色体到下一代染色体群体;(6)对下一代染色体群体执行交叉、变异操作;(7)如果满足终止条件,则结束,否则回到步骤3)继续进行计算,当计算结束时,其所得到的最优解为长度为N的二进制串,该串中的1的位置序号表示经过分布优化计算后所确定的需要使用的监测节点序号,通过其即可获知对于河流局部区域监测所需要采用的监测节点。
地址 430010 湖北省武汉市黄浦大街23号