发明名称 基于改进的人工蜂群算法的漏磁检测缺陷重构方法
摘要 本发明涉及一种基于改进的人工蜂群算法的漏磁检测缺陷重构方法,采用径向基函数神经网络作为前向模型,以前向模型预测的漏磁信号与实测漏磁信号的误差平方和作为目标函数,改进人工蜂群算法,引入当前个体最优解和全局最优解用来加快算法收敛速度,将改进的人工蜂群算法作为求解重构问题的迭代算法,最终得到的全局最优解即为重构的缺陷轮廓。本发明提高了漏磁检测缺陷重构的速度与精度。
申请公布号 CN104965941A 申请公布日期 2015.10.07
申请号 CN201510295609.0 申请日期 2015.06.02
申请人 上海电力学院 发明人 韩文花;汪胜兵;王建;吴正阳
分类号 G06F17/50(2006.01)I;G06N3/00(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 上海申汇专利代理有限公司 31001 代理人 吴宝根
主权项 一种基于改进的人工蜂群算法的漏磁检测缺陷重构方法,其特征在于,具体包括如下步骤:1)初始化算法参数及蜜源位置,并设定最大迭代数,设置初始蜜源,随机产生F<sub>N</sub>个初始蜜源<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>=</mo><msubsup><mi>x</mi><mi>j</mi><mi>min</mi></msubsup><mo>+</mo><mi>rand</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>j</mi><mi>max</mi></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>j</mi><mi>min</mi></msubsup><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>F</mi><mi>N</mi></msub><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>D</mi><mo>,</mo></mrow>]]></math><img file="FDA0000729835230000011.GIF" wi="1306" he="106" /></maths>式中的rand为[0,1]的随机值,x<sub>ij</sub>为第i个蜜源的第j维的值,<img file="FDA0000729835230000012.GIF" wi="257" he="93" />分别为j维的最大值和最小值,D表示维度;2)设定当前迭代次数为1;3)计算各蜜源的适应度值,采用改进的人工蜂群算法对蜜源的位置进行更新,具体包括:301)对蜜源按下式计算对应的适应度fitness值:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>fitness</mi><mi>i</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn><mo>+</mo><mo>|</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>F</mi><mi>N</mi></msub></mrow>]]></math><img file="FDA0000729835230000013.GIF" wi="907" he="238" /></maths>其中f(x<sub>i</sub>)是以测量的漏磁信号与RBFNN预测的漏磁信号间的误差平方和为目标函数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>f</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>D</mi></munderover><msup><mrow><mo>(</mo><msub><mi>Z</mi><mi>j</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000729835230000014.GIF" wi="420" he="174" /></maths>式中D为漏磁信号的维度,Z=[Z<sub>1</sub>,Z<sub>2</sub>,...,Z<sub>D</sub>]是实测漏磁信号,Y=[Y<sub>1</sub>,Y<sub>2</sub>,...,Y<sub>D</sub>]是RBFNN预测的漏磁信号,Z<sub>j</sub>与Y<sub>j</sub>分别为第j维实测漏磁信号与预测漏磁信号;302)通过比较适应度值的大小,得到当前个体最优解pbest<sub>i</sub>和全局最优解gbest,当前个体最优解pbest<sub>i</sub>为第i个蜜源在迭代中适应度值最大的解,全局最优解gbest为所有蜜源在迭代中适应度值最大的解;303)采蜜蜂按改进的下式搜索蜜源,并计算fitness,若蜜源质量变好则更新蜜源当前位置、pbest<sub>i</sub>和gbest的值并将计数变量counter<sub>i</sub>置0,否则采蜜蜂当前位置不变将counter<sub>i</sub>加1,<img file="FDA0000729835230000015.GIF" wi="1183" he="103" />其中t为当前迭代次数;j∈{1,2,...,D},k∈{1,2,...F<sub>N</sub>},j、k为其值域内的随机取值,且k1i;x<sub>kj</sub>为第k个蜜源的第j维的值;<img file="FDA0000729835230000021.GIF" wi="211" he="85" />为服从均匀分布的随机数;304)计算P<sub>i</sub>,<img file="FDA0000729835230000022.GIF" wi="367" he="234" />若P<sub>i</sub>大于一个随机值,则观察蜂转化为采蜜蜂,按下式搜索蜜源并计算fitness,若蜜源质量变好,则更新蜜源当前位置、pbest<sub>i</sub>和gbest的值,并将计数变量counter<sub>i</sub>置0,否则蜜源的当前位置不变,并将counter<sub>i</sub>加1;<img file="FDA0000729835230000023.GIF" wi="1103" he="116" />305)若counter<sub>i</sub>&gt;limit,则抛弃该蜜源,采蜜蜂转换为侦查蜂,随机选取新蜜源,其中limit为允许开采的最大次数;4)判断是否达到最大迭代次数,若是,则全局最优解作为漏磁信号的重构轮廓,若否,则将迭代次数加1,并用当前蜜源位置作为粒子的初始位置,并返回步骤3)。
地址 200090 上海市杨浦区平凉路2103号