[发明专利]一种基于量子鸽群优化的无人机自主空中加油目标检测方法在审

专利信息
申请号: 201510481163.0 申请日: 2015-08-03
公开(公告)号: CN105066998A 公开(公告)日: 2015-11-18
发明(设计)人: 段海滨;孙永斌 申请(专利权)人: 北京航空航天大学
主分类号: G01C21/20 分类号: G01C21/20;G01C21/00
代理公司: 北京慧泉知识产权代理有限公司 11232 代理人: 王顺荣;唐爱华
地址: 100191*** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明是一种基于量子鸽群优化的无人机自主空中加油目标检测方法,其实施步骤为:步骤一:获取处理图像;步骤二:初始化量子鸽群优化算法参数;步骤三:设计代价函数;步骤四:利用量子鸽群优化算法的地图罗盘算子进行寻优;步骤五:利用量子鸽群优化算法的地标算子进行寻优;步骤六:储存结果并验证。该方法可以对无人机的标志点进行准确地检测,为空中加油的实现提供有效的保障。
搜索关键词: 一种 基于 量子 鸽群 优化 无人机 自主 空中加油 目标 检测 方法
【主权项】:
一种基于量子鸽群优化的无人机自主空中加油目标检测方法,其特征在于:该方法的步骤如下:步骤一:获取处理图像用相机,获取无人机的相关图像;步骤二:初始化量子鸽群优化算法参数(1)初始化优化参数维数D用群体智能算法在二维灰度图像中寻找目标区域,所以D为2;(2)初始化收缩扩张系数ω(t)的参数初始化收缩扩张系数ω(t)的最大值ωmax和最小值ωmin,这两个值的确定影响控制算法的收敛速度;(3)初始化种群数量Np群智能优化算法的种群数量Np对优化效果影响很大,选择一个适当的种群数量,保证算法的准确性和快速性;(4)初始化种群位置和速度在搜索空间没设定群体的位置上限Pmax和位置下限Pmin,以及速度上限Vmax和速度下限Vmin;给种群中的每个粒子都初始化一个初始的位置Xi和初始的速度Vi<mrow><mo>{</mo><mrow><mtable><mtr><mtd><mrow><msub><mi>X</mi><mi>i</mi></msub><mo>=</mo><msub><mi>P</mi><mi>min</mi></msub><mo>+</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><mrow><msub><mi>P</mi><mi>max</mi></msub><mo>-</mo><msub><mi>P</mi><mi>min</mi></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>V</mi><mi>i</mi></msub><mo>=</mo><msub><mi>V</mi><mi>min</mi></msub><mo>+</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><mrow><msub><mi>V</mi><mi>max</mi></msub><mo>-</mo><msub><mi>V</mi><mi>min</mi></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></mrow>(5)设置算法代数鸽群优化算法有两个算子,分别是地图罗盘算子和地标算子,算法运算前需要分别设定两个算法运行的最大代数NC1max和NC2max;步骤三:设计代价函数代价函数的确定是智能优化算法的核心,决定目标检测的准确性;采用最简单有效的累加计数方法,判断整幅图像中的各分区域的像素值与目标图像像素的大小关系,然后累加计数,得到在整幅图像中得到最优的区域;步骤四:利用QPIO地图罗盘算子进行寻优利用初始化的群体位置和速度,根据初始的个体的代价函数值选取全局最优位置Xg;根据公式<mrow><msub><mi>m</mi><mrow><mi>b</mi><mi>e</mi><mi>s</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><mi>p</mi></msub></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>p</mi></msub></munderover><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>PiPg(t+1)=f(t+1)×Pi(t)+(1‑f(t+1))×Pg(t)  (2)<mrow><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&omega;</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>max</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mi>min</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mfrac><mi>t</mi><msub><mi>t</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>如果f(t+1)≥0.5,那么<mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mi>i</mi></msub><msub><mi>P</mi><mi>g</mi></msub><mo>+</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>&times;</mo><mo>|</mo><msub><mi>m</mi><mrow><mi>b</mi><mi>e</mi><mi>s</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>&times;</mo><mi>ln</mi><mfrac><mn>1</mn><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>如果f(t+1)<0.5,那么<mrow><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mi>i</mi></msub><msub><mi>P</mi><mi>g</mi></msub><mo>-</mo><mi>&omega;</mi><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>&times;</mo><mo>|</mo><msub><mi>m</mi><mrow><mi>b</mi><mi>e</mi><mi>s</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo><mo>&times;</mo><mi>ln</mi><mfrac><mn>1</mn><mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>更新每个个体的位置Xi,计算新生成鸽子的代价函数值,如果新鸽子的代价函数值比全局最优位置的代价函数值更高,则把新生成的鸽子位置定义为新的全局最优位置Xg;反复应用地图罗盘算子进行寻优,直到运行代数大于地图罗盘算子最大代数NC1max时停止;公式中,mbest(t+1)为t+1代所有鸽子最佳位置的平均值,Pi(t)、Pg(t)分别为第i只鸽子的最优位置和所有鸽子的全局最优位置,PiPg(t+1)为鸽子个体最优值Pi(t)和群体全局最优值Pg(t)之间的随机点,ω(t)为收缩扩张系数,调节它的值能控制算法的收敛速度,较大的ω(t)有利于算法跳出局部最优值,而较小的ω(t)有利于算法的收敛,f(t+1)、u(t+1)是一个从0到1之间随机产生的一个随机数;步骤五:利用QPIO地标算子进行寻优利用地图罗盘算子寻优的结果作为地标算子的初始群体,根据公式<mrow><mo>(</mo><msub><mi>N</mi><mi>p</mi></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>=</mo><mfrac><mrow><msub><mi>N</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>X</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>&Sigma;X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>f</mi><mi>i</mi><mi>t</mi><mi>n</mi><mi>e</mi><mi>s</mi><mi>s</mi><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><msub><mi>N</mi><mi>p</mi></msub><mo>&CenterDot;</mo><mi>&Sigma;</mi><mi>f</mi><mi>i</mi><mi>t</mi><mi>n</mi><mi>e</mi><mi>s</mi><mi>s</mi><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>Xi(t)=Xi(t‑1)+rand·(Xc(t)‑Xi(t‑1))  (8)更新每个个体的速度Vi和位置Xi,计算新生成鸽子的代价函数值,如果新鸽子的代价函数值比全局最优位置的代价函数值更高,则把新生成的鸽子位置定义为新的全局最优位置Xg;根据公式(6)计算新种群的群体数量,根据公式(6)计算的结果舍弃群体中代价函数较小的一部分个体,选择当前群体中较优的群体作为保留群体进行下一轮寻优,反复应用地标算子进行寻优,直到运行代数大于地标算子最大代数NC2max时停止;公式中,Np为鸽群的数目,fitness是鸽子位置信息的代价函数,为了求得代价函数的最小值,取fmin作为目标函数,Xc是鸽群的加权位置中心;步骤六:储存结果并验证地标算子优化的结果被视为最终的目标检测结果,将此结果保存,并保存下寻优过程的进化曲线。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京航空航天大学,未经北京航空航天大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201510481163.0/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top