发明名称 一种多机器人系统合作定位危险气味源方法
摘要 本发明涉及一种多机器人系统合作定位危险气味源方法。现有的定位方法成功率较低,系统耗能过多。本发明首先建立气味源的观测模型,然后获得气味源位置的估计值,气味源位置的先验概率分布;其次如果检测到气味,修正先验概率分布,得到后验概率分布,并对后验概率分布采样,获得该气味源位置的估计值,产生机器人的新位置;如果没有检测到气味,直接采样气味源位置的先验概率分布,获得该气味源位置的估计值,产生机器人的新位置;最后用一致性算法控制机器人向新位置移动。本发明方法弥补了传统方法的不足,并有效地提高气味源定位的精度,保证多机器人系统能量的较小消耗,同时满足实际中要求快速定位的要求。
申请公布号 CN102034030A 申请公布日期 2011.04.27
申请号 CN201010609631.5 申请日期 2010.12.28
申请人 杭州电子科技大学 发明人 吕强;谢小高;罗平
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1.一种多机器人系统合作定位危险气味源方法,其特征在于该方法包括以下步骤:(1)利用气味分子的运动学模型,建立气味源位置的观测模型,具体方法是:A.建立单个气味分子团的运动模型:<img file="684940DEST_PATH_IMAGE001.GIF" wi="109" he="33" />;其中<img file="385042DEST_PATH_IMAGE002.GIF" wi="29" he="22" />指气味分子团在<img file="753576DEST_PATH_IMAGE003.GIF" wi="9" he="17" />时刻的位置,<img file="208828DEST_PATH_IMAGE004.GIF" wi="29" he="33" />指气味分子团在<img file="46334DEST_PATH_IMAGE003.GIF" wi="9" he="17" />时刻位置的微分,<img file="300860DEST_PATH_IMAGE005.GIF" wi="30" he="22" />是在<img file="274632DEST_PATH_IMAGE003.GIF" wi="9" he="17" />时刻的均值风速度,<img file="900786DEST_PATH_IMAGE006.GIF" wi="30" he="22" />表示一个随机过程,该随机过程服从均值为零,方差为<img file="474855DEST_PATH_IMAGE007.GIF" wi="24" he="22" />的高斯分布;对气味分子团的运动模型在时间段[<img file="516761DEST_PATH_IMAGE008.GIF" wi="13" he="25" />,<img file="672935DEST_PATH_IMAGE009.GIF" wi="16" he="25" />]内进行积分,获得气味分子团在<img file="158406DEST_PATH_IMAGE008.GIF" wi="13" he="25" />时刻被气味源释放,在当前时刻<img file="767242DEST_PATH_IMAGE009.GIF" wi="16" he="25" />时的位置:<img file="878417DEST_PATH_IMAGE010.GIF" wi="236" he="53" />其中[<img file="76049DEST_PATH_IMAGE008.GIF" wi="13" he="25" />,<img file="44005DEST_PATH_IMAGE009.GIF" wi="16" he="25" />]指积分的时间段,<img file="77820DEST_PATH_IMAGE008.GIF" wi="13" he="25" />指气味源释放气味分子团的时间,<img file="523845DEST_PATH_IMAGE009.GIF" wi="16" he="25" />指当前时间,<img file="654612DEST_PATH_IMAGE011.GIF" wi="50" he="25" />指的是在当前时间<img file="481885DEST_PATH_IMAGE012.GIF" wi="16" he="25" />气味分子团的位置,<img file="2996DEST_PATH_IMAGE013.GIF" wi="37" he="25" />是在时刻<img file="252712DEST_PATH_IMAGE014.GIF" wi="14" he="25" />气味源的位置;B.通过定义<img file="424936DEST_PATH_IMAGE015.GIF" wi="191" he="53" />和<img file="734695DEST_PATH_IMAGE016.GIF" wi="118" he="53" />,将运动模型离散化;其中离散时刻<img file="743102DEST_PATH_IMAGE017.GIF" wi="9" he="20" />是气味源释放气味分子团时间<img file="484924DEST_PATH_IMAGE014.GIF" wi="14" he="25" />的整数值,离散时刻<img file="324704DEST_PATH_IMAGE018.GIF" wi="14" he="20" />是当前时间<img file="743047DEST_PATH_IMAGE012.GIF" wi="16" he="25" />的整数值,<img file="35488DEST_PATH_IMAGE019.GIF" wi="20" he="20" />是指时间的微小变化,气味分子团的离散模型为:<img file="79536DEST_PATH_IMAGE020.GIF" wi="202" he="25" />其中<img file="711506DEST_PATH_IMAGE021.GIF" wi="44" he="22" />表示气味源在<img file="363067DEST_PATH_IMAGE017.GIF" wi="9" he="20" />时刻释放气味分子团,在<img file="96799DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻的位置;<img file="492008DEST_PATH_IMAGE022.GIF" wi="45" he="22" />表示在时间段[<img file="978485DEST_PATH_IMAGE017.GIF" wi="9" he="20" />,<img file="66526DEST_PATH_IMAGE018.GIF" wi="14" he="20" />]内风推动气味分子团的移动距离;<img file="333560DEST_PATH_IMAGE023.GIF" wi="49" he="22" />表示随机过程,该随机过程服从均值为零,方差<img file="453831DEST_PATH_IMAGE024.GIF" wi="61" he="25" />的正态分布;令<img file="60393DEST_PATH_IMAGE025.GIF" wi="123" he="46" />,<img file="319336DEST_PATH_IMAGE026.GIF" wi="129" he="46" />则气味分子团的离散运动模型进一步表示如下:<img file="496502DEST_PATH_IMAGE027.GIF" wi="169" he="28" />其中<img file="436776DEST_PATH_IMAGE028.GIF" wi="33" he="22" />为<img file="960161DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻气味分子团的位置,<img file="311377DEST_PATH_IMAGE029.GIF" wi="37" he="25" />为<img file="490686DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻气味源的位置,<img file="31388DEST_PATH_IMAGE030.GIF" wi="34" he="26" />为<img file="85977DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻风推动气味分子团移动的距离,<img file="686722DEST_PATH_IMAGE031.GIF" wi="38" he="26" />是测量噪声并且服从零均值,<img file="87748DEST_PATH_IMAGE032.GIF" wi="97" he="46" />方差的正态分布;C.通过定义<img file="432142DEST_PATH_IMAGE033.GIF" wi="124" he="26" />,得到气味源位置的观测模型:<img file="117070DEST_PATH_IMAGE034.GIF" wi="132" he="28" />其中,<img file="560820DEST_PATH_IMAGE035.GIF" wi="38" he="25" />是第<img file="777038DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人在时刻<img file="347959DEST_PATH_IMAGE018.GIF" wi="14" he="20" />对气味源位置<img file="700443DEST_PATH_IMAGE029.GIF" wi="37" he="25" />的测量值;(2)基于最大熵粒子滤波原理,获得气味源位置的先验概率分布;如果在<img file="315095DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻,检测到气味,使用测量值修正气味源位置的先验概率分布,得到气味源位置的后验概率分布,通过采样气味源位置的后验概率分布产生气味源位置新的估计值;然后,获得第<img file="753030DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人的下一步的位置;如果在<img file="439226DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻,没有检测到气味,则采样气味源位置的先验概率分布产生气味源位置新的估计值,获得第<img file="833167DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人的下一步的位置,具体方法是:a.建立关于气味源位置的背景分布<img file="946617DEST_PATH_IMAGE037.GIF" wi="36" he="22" />;<img file="809530DEST_PATH_IMAGE038.GIF" wi="131" he="25" />其中,<img file="722254DEST_PATH_IMAGE039.GIF" wi="13" he="14" />表示气味源位置满足正态分布的随机变量,<img file="49330DEST_PATH_IMAGE040.GIF" wi="23" he="25" />是气味源位置的经验估计值,<img file="271364DEST_PATH_IMAGE041.GIF" wi="22" he="25" />是关于气味源位置的经验估计协方差矩阵,<img file="418311DEST_PATH_IMAGE042.GIF" wi="30" he="22" />表示正态分布;根据各个机器人对气味源位置的估计值,求解下式:<img file="633261DEST_PATH_IMAGE043.GIF" wi="148" he="48" /><img file="752527DEST_PATH_IMAGE044.GIF" wi="210" he="48" />其中<img file="896194DEST_PATH_IMAGE045.GIF" wi="21" he="22" />表示机器人还没有获得测量值时的<img file="733700DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻,<img file="299811DEST_PATH_IMAGE046.GIF" wi="57" he="26" />是第<img file="522851DEST_PATH_IMAGE047.GIF" wi="14" he="21" />个机器人对气味源位置的先验估计值,<img file="149004DEST_PATH_IMAGE048.GIF" wi="48" he="26" />表示的是测量预测值的先验平均值,并由第<img file="208227DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人获得;<img file="578028DEST_PATH_IMAGE049.GIF" wi="52" he="26" />表示测量预测值的先验矩阵,<img file="422619DEST_PATH_IMAGE050.GIF" wi="14" he="16" />是机器人的个数;b.基于最大熵原理,采用共轭梯度法求解下式:<img file="157357DEST_PATH_IMAGE051.GIF" wi="538" he="46" />其中:<img file="766192DEST_PATH_IMAGE052.GIF" wi="252" he="25" /><img file="861056DEST_PATH_IMAGE053.GIF" wi="459" he="103" /><img file="809421DEST_PATH_IMAGE054.GIF" wi="285" he="33" /><img file="777377DEST_PATH_IMAGE055.GIF" wi="442" he="30" />通过求解上式,获得拉格朗日乘子<img file="561924DEST_PATH_IMAGE056.GIF" wi="42" he="25" />和<img file="273528DEST_PATH_IMAGE057.GIF" wi="45" he="25" />,从而得到气味源位置的先验概率分布<img file="341978DEST_PATH_IMAGE058.GIF" wi="219" he="26" />:<img file="667786DEST_PATH_IMAGE059.GIF" wi="453" he="26" />其中<img file="985635DEST_PATH_IMAGE039.GIF" wi="13" he="14" />表示气味源位置满足正态分布的随机变量,<img file="438613DEST_PATH_IMAGE060.GIF" wi="114" he="26" />是气味源位置的均值,<img file="423887DEST_PATH_IMAGE061.GIF" wi="73" he="26" />是关于气味源位置的协方差矩阵,<img file="422061DEST_PATH_IMAGE042.GIF" wi="30" he="22" />表示正态分布;c.根据第<img file="227206DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人是否检测到气味,执行不同的步骤;如果检测到气味,则执行c-1,否则执行c-2;c-1. 检测到气味,执行以下两个步骤:①使用测量值更新拉格朗日乘子<img file="218296DEST_PATH_IMAGE056.GIF" wi="42" he="25" />和<img file="58076DEST_PATH_IMAGE057.GIF" wi="45" he="25" />,得:<img file="538736DEST_PATH_IMAGE062.GIF" wi="184" he="25" /><img file="283707DEST_PATH_IMAGE063.GIF" wi="158" he="25" />其中<img file="812908DEST_PATH_IMAGE064.GIF" wi="21" he="22" />表示机器人获得测量值时的<img file="772774DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻,<img file="112751DEST_PATH_IMAGE065.GIF" wi="38" he="22" />是测量噪声协方差矩阵,这样就获得了气味源位置的后验概率分布<img file="892488DEST_PATH_IMAGE066.GIF" wi="221" he="26" />:<img file="225380DEST_PATH_IMAGE067.GIF" wi="453" he="26" />采样气味源位置的后验概率分布<img file="226703DEST_PATH_IMAGE066.GIF" wi="221" he="26" />,得到气味源位置的最新估计值<img file="49166DEST_PATH_IMAGE068.GIF" wi="56" he="26" />;②基于气味源位置的最新估计值<img file="253882DEST_PATH_IMAGE068.GIF" wi="56" he="26" />,则机器人下一步的位置为气味源位置的最新估计值,即<img file="452782DEST_PATH_IMAGE069.GIF" wi="126" he="26" />,其中<img file="821795DEST_PATH_IMAGE070.GIF" wi="61" he="25" />是<img file="815159DEST_PATH_IMAGE071.GIF" wi="34" he="20" />时刻第<img file="507171DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人的位置;c-2.没有检测到气味,则执行以下两个步骤:③采样气味源位置的先验概率分布<img file="509763DEST_PATH_IMAGE058.GIF" wi="219" he="26" />,得到气味源位置的最新估计值<img file="220098DEST_PATH_IMAGE072.GIF" wi="56" he="26" />;④基于气味源位置的最新估计值<img file="384364DEST_PATH_IMAGE072.GIF" wi="56" he="26" />和没有检测到气味的周期数,计算机器人下一步的位置,如果在两个采用周期内没有检测到气味,执行步骤ⅰ,如果在两个采用周期及两个采用周期以上没有检测到气味,则执行ⅱ;ⅰ机器人的下一步的位置为:<img file="563672DEST_PATH_IMAGE073.GIF" wi="405" he="54" /><img file="104375DEST_PATH_IMAGE074.GIF" wi="178" he="26" />其中<img file="170682DEST_PATH_IMAGE075.GIF" wi="21" he="26" />和<img file="771428DEST_PATH_IMAGE076.GIF" wi="20" he="26" />分别是风在<img file="172453DEST_PATH_IMAGE077.GIF" wi="14" he="16" />轴和<img file="703798DEST_PATH_IMAGE078.GIF" wi="16" he="18" />轴方向的速度,<img file="936196DEST_PATH_IMAGE079.GIF" wi="58" he="26" />是第<img file="645526DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人在<img file="284580DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻对于气味源位置最新估计值在<img file="167085DEST_PATH_IMAGE077.GIF" wi="14" he="16" />轴方向的坐标,<img file="519569DEST_PATH_IMAGE080.GIF" wi="45" he="25" />第<img file="399800DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人在<img file="24686DEST_PATH_IMAGE018.GIF" wi="14" he="20" />时刻位置在<img file="710882DEST_PATH_IMAGE077.GIF" wi="14" he="16" />轴方向的坐标,<img file="917872DEST_PATH_IMAGE081.GIF" wi="65" he="25" />和<img file="969005DEST_PATH_IMAGE082.GIF" wi="65" he="26" />第<img file="582651DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人在<img file="806959DEST_PATH_IMAGE071.GIF" wi="34" he="20" />时刻位置在<img file="134035DEST_PATH_IMAGE077.GIF" wi="14" he="16" />轴和<img file="90490DEST_PATH_IMAGE078.GIF" wi="16" he="18" />轴方向的坐标;ⅱ机器人下一步位置是基于气味源位置最新估计值<img file="503017DEST_PATH_IMAGE046.GIF" wi="57" he="26" />产生,具体方法是:<img file="717966DEST_PATH_IMAGE083.GIF" wi="281" he="47" /><img file="899549DEST_PATH_IMAGE084.GIF" wi="277" he="48" />其中<img file="292484DEST_PATH_IMAGE085.GIF" wi="17" he="22" />和<img file="926728DEST_PATH_IMAGE086.GIF" wi="14" he="18" />分别是控制<img file="181254DEST_PATH_IMAGE077.GIF" wi="14" he="16" />轴和<img file="482922DEST_PATH_IMAGE078.GIF" wi="16" he="18" />轴方向搜索范围控制参数,<img file="46759DEST_PATH_IMAGE087.GIF" wi="49" he="22" />在<img file="168299DEST_PATH_IMAGE088.GIF" wi="32" he="22" />范围内产生一个均匀分布的随机数,<img file="725051DEST_PATH_IMAGE050.GIF" wi="14" he="16" />是机器人的个数;d.当第<img file="881226DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人到达新的位置后,重复步骤a至步骤c,直到群体机器人中的某一个机器人检测到的气味浓度达到预先设定值,或通过视觉传感器判断出气味源,然后通过无线网络通知第<img file="615963DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人,则第<img file="224799DEST_PATH_IMAGE036.GIF" wi="9" he="18" />个机器人停止搜索,并提供气味源的估计位置<img file="821128DEST_PATH_IMAGE068.GIF" wi="56" he="26" />。
地址 310018 浙江省杭州市下沙高教园区2号大街