首页> 中国专利> 基于微分测度随机相遇不确定性的搜救方法及系统

基于微分测度随机相遇不确定性的搜救方法及系统

摘要

本发明公开了一种基于微分测度随机相遇不确定性的搜救方法及系统,其中方法包括以下步骤:步骤1、确定走失者D所走失的路径L,并确定走失者D最后出现在路径L上的位置点;步骤2、测量搜寻者C在路径L上的位置点;步骤3、计算搜寻者C与走失者D分布在路径L的概率密度函数p

著录项

  • 公开/公告号CN105095680A

    专利类型发明专利

  • 公开/公告日2015-11-25

    原文格式PDF

  • 申请/专利权人 武汉理工大学;

    申请/专利号CN201510603233.5

  • 发明设计人 尹章才;吴杨;胡立夫;

    申请日2015-09-21

  • 分类号G06F19/00;

  • 代理机构湖北武汉永嘉专利代理有限公司;

  • 代理人许美红

  • 地址 430070 湖北省武汉市洪山区珞狮路122号

  • 入库时间 2023-12-18 12:21:18

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2018-01-23

    授权

    授权

  • 2015-12-23

    实质审查的生效 IPC(主分类):G06F19/00 申请日:20150921

    实质审查的生效

  • 2015-11-25

    公开

    公开

说明书

技术领域

本发明涉及搜救领域,尤其涉及一种基于微分测度随机相遇不确定性的搜 救方法及系统。

背景技术

近年来,人员搜寻和救援案例时有发生。例如,户外活动中的人员失联与 搜救,老人小孩的走失与搜寻等。一种典型案例是这样的:已知走失者最后出 现在路径L的中间位置且在路径L上自由移动;一搜寻者在路径L上遇到走 失者的可能性有多大?搜寻者在路径L的哪个位置点遇见走失者的可能性最 大?

1.改进前方法工作的机理

传统的救援过程往往出于人的本能或人道主义,将大量的搜救资源投入到 走失者最后出现地点周边的一定区域,缺乏精准的搜救规划(刘钊,等,2014)。 其中,预知搜寻者能遇见走失者的可能性,是搜救方案规划和最大可能地成功 救援的前提。已有的概率时间地理学采用概率值定量化表达相遇可能性,并提 出了一种计算相遇概率的离散型方法(Winter,YIN,2011)。该方法规定:搜寻者 C与走失者D能相遇的条件是在离散型地理空间中C、D位于同一离散单元中。

设:搜寻者C与走失者D所在的路径L的长度为l。C、D分布在路径L 上的概率密度函数为c、d。相遇概率的离散型方法的计算步骤如下:

步骤1:将路径L均匀划分为n小段:L1、L2、…、Ln(图1(a))。

步骤2:记C、D分别位于任一小段Li的概率值ci、di,有:0≤ci≤1,0≤di≤1(图1(b)),Σi=1nci=1,Σi=1ndi=1,i=1,2,...,n.

步骤3:在搜寻者找到走失者之前,两个体的移动可视为独立的。这样, 个体C、D位于或相遇于任一单元Li的概率值为ci×di。相应地,相遇于整个 路径L的概率值为Σi=1nci×di,i=1,2,...,n(图1(c))。

为了简单起见,令c、d为均匀分布。

(1)当n=1时,则c、d分布在L1的概率值c1=1,d1=1,相应的相遇 概率:Σi=1nci×di=Σi=11ci×di=c1×d1=1.

(2)当n=2时,则c、d分布在L1的概率值c1=d1=0.5,c、d分布在L2的概率值c2=d2=0.5,相应的相遇概率:Σi=1nci×di=Σi=12ci×di=c1×d1+c2×d2=0.5.

(3)当n=10时,则c、d分布在Li的概率值ci=di=0.1,相应的相遇概 率:Σi=1nci×di=Σi=110ci×di=10×(0.1×0.1)=0.1.

由上可知,当离散单元的数量n不断增大时,搜寻者C在路径L上遇见 走失者D的可能性不断减小,即成功搜寻的概率与n成反比;这一结论同样 适用于c、d为非均匀分布的情形,如正态分布、三角形分布等。

2.改进前方法存在的问题

综上,离散型方法依赖于离散单元的尺度(Winter,YIN,2011),相遇概率会 随离散单元数量的增加而减小;这样,离散型方法中尺度或离散单元数量设置 的人为性势必造成相遇概率的随意性。然而,作为客观规律的相遇概率具有稳 定性,在理论上与计算方法无关。

发明内容

本发明针对在路径L上搜寻者C能遇见走失者D的概率问题,离散型方 法不能给出具有稳定性或唯一性的相遇概率。提出一种根据搜寻者与走失者分 布在路径L上的非均匀概率分布,利用连续微分推算出具有唯一性的相遇概率 的基于微分测度随机相遇不确定性的搜救方法及系统。

本发明解决其技术问题所采用的技术方案是:

提供一种基于微分测度随机相遇不确定性的搜救方法,包括以下步骤:

步骤1、确定走失者D所走失的路径L,并确定走失者D最后出现在路径 L上的位置点;

步骤2、测量搜寻者C在路径L上的位置点;

步骤3、计算搜寻者C与走失者D分布在路径L的概率密度函数pc、pd, C位于路经L上的位置点x0时能与D相遇的概率相应计算路径L上其他点x1、x2、…、xk、…的微分相遇概率获 得{p(Emeet|x0),p(Emeet|x1),…,p(Emeet|xk),…},找到微分相遇概率最大值所对 应的位置点为xm

步骤4、当搜寻者C位于全部点{x|x=x0,x1,…,xk,…}时,搜救者C能与 走失者D的相遇概率p(Emeet):

p(Emeet)=p(Emeet|x0)+p(Emeet|x1)+...+p(Emeet|xk)+...=pc(x0)dx×|y-x0|md,yLpd(y)dy+pc(x1)dx×|y-x1|md,yLpd(y)dy+...+pc(xk)dx×|y-xk|md,yLpd(y)dy+...=ΣxLpc(x)dx·|y-x|md,yLpd(y)dy=xLpc(x)·|y-x|md,yLpd(y)dydx

则搜寻者C在路径L上能找到走失者D的概率:

其中md为搜寻者C与走失者D可相遇的 最大距离,变量x表示搜寻者C所在位置点与路径L的一个端点O之间的距 离;变量y表示走失者D所在位置点与端点O之间的距离;

步骤5、根据计算的与微分相遇概率最大值相对应的位置点,以及相遇概 率p(Emeet)安排搜寻者C在路径L上进行搜救。

本发明所述的搜救方法中,若走失者D最后出现在路径L的中间点且只 在路径L上作自由移动,则可合理假设走失者D分布在路径L的概率密度函 数pd呈三角形分布;若搜寻者C从路径L的中间点开始找寻,则也可合理假 设搜寻者C分布在路径L的概率密度函数pc呈三角形分布。

本发明还提供一种基于微分测度随机相遇不确定性的搜救系统,该系统包 括:

确认模块,用于确定走失者D所走失的路径L,并确定走失者D最后出现 在路径L上的位置点;

数据获取模块,用于测量的搜寻者C在路径L上的位置点;

概率密度函数计算模块,用于计算搜寻者C与走失者D分布在路径L的概 率密度函数pc、pd,C位于路经L上的位置点x0时能与D相遇的概率相应计算路径L上其他点x1、x2、…、xk、…的微分相 遇概率,获得{p(Emeet|x0),p(Emeet|x1),…,p(Emeet|xk),…},找到微分相遇概率 最大值所对应的位置点为xm

相遇概率计算模块,用于计算当搜寻者C位于全部点{x|x=x0,x1,…, xk,…}时,搜救者C能与走失者D的相遇概率p(Emeet):

p(Emeet)=p(Emeet|x0)+p(Emeet|x1)+...+p(Emeet|xk)+...=pc(x0)dx×|y-x0|md,yLpd(y)dy+pc(x1)dx×|y-x1|md,yLpd(y)dy+...+pc(xk)dx×|y-xk|md,yLpd(y)dy+...=ΣxLpc(x)dx·|y-x|md,yLpd(y)dy=xLpc(x)·|y-x|md,yLpd(y)dydx

则搜寻者C在路径L上能找到走失者D的概率为:

其中md为搜寻者C与走失者D可相遇的 最大距离,变量x表示搜寻者C所在位置点与路径L的一个端点O之间的距 离;变量y表示走失者D所在位置点与端点O之间的距离;

分析处理模块,用于根据计算的与微分相遇概率最大值相对应的位置点, 以及相遇概率p(Emeet)安排搜寻者C在路径L上进行搜救。

本发明所述的搜救系统中,所述概率密度函数计算模块具体用于在走失者 D最后出现在路径L的中间点且只在路径L上作自由移动时,则合理假设走失 者D分布在路径L的概率密度函数pd呈三角形分布;在搜寻者C从路径L的 中间点开始找寻时,则也合理假设搜寻者C分布在路径L的概率密度函数pc 呈三角形分布。

本发明产生的有益效果是:本发明利用相遇概率的微分方法,根据可相遇 的最大距离md和搜寻者C、走失者D的概率密度函数pc、pd,可解决C能找 到D的概率以及在何处找到的概率最大等问题。本发明的相遇概率的微分方 法,不同于离散型方法,本发明计算的相遇概率值具有稳定性和唯一性。

附图说明

下面将结合附图及实施例对本发明作进一步说明,附图中:

图1:传统相遇概率的离散型方法图,其中(a)路径的离散化;(b)个体位于 离散单元的概率;(c)相遇概率;

图2:搜寻者C能遇见走失者D的事件,其中(a)变量定义;(b)相遇事件示 例;

图3:搜寻者C在何处能遇见走失者D的概率最大的推理方法,其中(a) 微分概率;(b)C位于给定点的微分相遇概率;(c)C位于序列点的序列微分相 遇概率;

图4:本发明中搜寻者C能遇见走失者D的概率的推理方法;

图5:本发明的相遇概率微分法的技术路线图;

图6:本发明一个实施例中相遇概率最大值所在的空间位置;

图7:本发明另一个实施例中相遇概率最大值所在的空间位置。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实 施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅 用以解释本发明,并不用于限定本发明。

1.本发明所采取的主要技术方案

在现实环境中,两个个体之间的相遇主要受两者之间的空间距离(如可视 距离)制约。这里,将两个体可相遇的最大距离,记为md(meetingdistance)。 据此,相遇语义可定义为:当且仅当两个体的相距距离不超过md时就认为能 相遇。这样,md在一定程度上确定了相遇的尺度,从而为相遇概率的唯一性 提供了理论基础。

下面介绍本发明所采取的主要技术方案。

1)相遇事件

令,直线段L的一个端点为O;搜寻者C与走失者D同位于L上。

设:变量x表示搜寻者C所在位置点与端点O之间的距离;变量y表示 走失者D所在位置点与端点O之间的距离(图2a)。

根据相遇语义,相遇事件Emeet={搜寻者C所在的位置x与走失者D所在 的位置y之间在距离上不超过md的事件},即:

Emeet={(x,y)||y-x|≤md,x∈L,y∈L}(公式1)

或者,Emeet={(x,y)|x-md≤y≤x+md,x∈L,y∈L}

这样,在搜寻者C位于位置点x时,如果走失者D所在的位置y1满足公 式1,即|y1-x|≤md,则相遇事件Emeet发生;如果走失者D所在的位置y2不 能满足公式1,即|y1-x|>md,则相遇事件Emeet不可能发生(图2b)。

2)搜寻者C在何处能以最大概率找到走失者D

相遇概率p(Emeet)就是相遇事件Emeet发生的概率。令,搜寻者C和走失者 D分布在路径L上的概率密度函数分别为:pc、pd,且C与D的运动相互独立。

通常情况下,变量x的增量Δx被称为x的微分,记作dx,即dx=Δx。

搜寻者C位于点x0处的概率pc(x0)可采用微分概率pc(x0)dx表示。其中, 微分概率pc(x0)dx,等于搜寻者C的概率密度函数pc分布在点x0处的概率pc(x0) 与微分线段长度dx的乘积(图3a阴影部分的面积)。在Δx→0时,根据微分 原理有pc(x0)dx→pc(x0);这样,微分概率pc(x0)dx能表示概率pc(x0)。

当搜寻者C位于一个点x0(x0∈L)时,根据公式(1)相遇事件Emeet的发 生意味着走失者D位于点x0附近,即|y-x0|≤md(图3b粗横线)。走失者D 位于点x0附近的概率为(图3b阴影部分的面积),它表示D的概 率pd分布在区域{|y-x0|≤md,y∈L}上的概率值。这样,C位于一个点x0时 能与D相遇的概率

p(Emeet|x0)=pc(x0)dx×|y-x0|md,yLpd(y)dy(公式2)

其中,pc(x0)dx表示C位于点x0的微分概率。类似地,我们可以获得其他微分 相遇概率:

C位于一个点x1时能与D相遇的概率

C位于一个点x2时能与D相遇的概率

……

C位于一个点xk时能与D相遇的概率

……

综上,能获得{p(Emeet|x0),p(Emeet|x1),…,p(Emeet|xk),…}的最大值(如图 3(c)所示),不失一般性可设该最大值所在的位置点为xm。这样,搜寻者C 在xm处能以最大概率找到走失者D。

3)搜寻者C在路径L上能找到走失者D的概率

当搜寻者C位于全部点{x|x=x0,x1,…,xk,…}时,C能与D相遇的概率:

p(Emeet)=p(Emeet|x0)+p(Emeet|x1)+...+p(Emeet|xk)+...=pc(x0)dx×|y-x0|md,yLpd(y)dy+pc(x1)dx×|y-x1|md,yLpd(y)dy+...+pc(xk)dx×|y-xk|md,yLpd(y)dy+...=ΣxLpc(x)dx·|y-x|md,yLpd(y)dy=xLpc(x)·|y-x|md,yLpd(y)dydx

这样,搜寻者C在路径L上能找到走失者D的概率:

p(Emeet)=xLpc(x)·|y-x|md,yLpd(y)dydx(公式3)

从积分的角度,p(Emeet)就是p(Emeet|x)与X轴之间的面积(图4)。

2.技术路线

上述相遇概率技术方案,搜寻者C找到走失者D的概率计算可分为三步(图 5)。

步骤1:获取搜寻基本数据,包括路径L的长度,可相遇的最大距离md, 以及搜寻者C、走失者D分布在L的概率密度函数pc、pd等。

步骤2:根据公式(2),计算搜寻者C位于L上各点的微分相遇概率,分 析搜寻者C在何处能以最大概率找到走失者D。

步骤3:根据公式(3),推理搜寻者C在路径L上能找到走失者D的概率。

本发明的相遇概率的微分方法,根据可相遇的最大距离md和搜寻者C、 走失者D的概率密度函数pc、pd,能回答C能找到D的概率以及在何处找到 的概率最大等问题。相遇概率的微分方法,一方面不同于离散型方法,前者计 算的相遇概率值具有稳定性和唯一性,后者则具有随意性和非唯一性;另一方 面也不同于积分法,对于C、D的概率密度函数为非均匀分布时,相遇概率的 积分法由于将相遇事件作为一个整体因而往往不能计算复杂函数的积分,或计 算的复杂度大;相遇概率的微分法则由于将相遇事件划分成若干子事件,因而 能将复杂函数的积分问题转换为单个相遇事件上简单函数的积分问题,从而有 助于简化积分及相遇概率的计算。因此,对于C、D的概率密度函数呈非线性 特征时,采用微分法可降低计算复杂度;对于C、D的概率密度函数为均匀分 布或呈线性特征时,可以采用积分法或等效的微分法。

本发明的基于微分测度随机相遇不确定性的搜救系统,用于实现上述搜救 方法,该系统包括:

确认模块,用于确定走失者D所走失的路径L,并确定走失者D最后出现 在路径L上的位置点;

数据获取模块,用于测量的搜寻者C在路径L上的位置点;

概率密度函数计算模块,用于计算搜寻者C与走失者D分布在路径L的概 率密度函数pc、pd,C位于路经L上的位置点x0时能与D相遇的概率相应计算路径L上其他点x1、x2、…、xk、…的微分相 遇概率,获得{p(Emeet|x0),p(Emeet|x1),…,p(Emeet|xk),…},找到微分相遇概率 最大值所对应的位置点为xm

相遇概率计算模块,用于计算当搜寻者C位于全部点{x|x=x0,x1,…, xk,…}时,搜救者C能与走失者D的相遇概率p(Emeet):

p(Emeet)=p(Emeet|x0)+p(Emeet|x1)+...+p(Emeet|xk)+...=pc(x0)dx×|y-x0|md,yLpd(y)dy+pc(x1)dx×|y-x1|md,yLpd(y)dy+...+pc(xk)dx×|y-xk|md,yLpd(y)dy+...=ΣxLpc(x)dx·|y-x|md,yLpd(y)dy=xLpc(x)·|y-x|md,yLpd(y)dydx

则搜寻者C在路径L上能找到走失者D的概率为:

其中md为搜寻者C与走失者D可相遇的 最大距离,变量x表示搜寻者C所在位置点与路径L的一个端点O之间的距 离;变量y表示走失者D所在位置点与端点O之间的距离。

分析处理模块,用于根据计算的与微分相遇概率最大值相对应的位置点, 以及相遇概率p(Emeet)安排搜寻者C在路径L上进行搜救。

其中,所述概率密度函数计算模块具体用于在走失者D最后出现在路径L 的中间点且只在路径L上作自由移动时,则合理假设走失者D分布在路径L 的概率密度函数pd呈三角形分布;在搜寻者C从路径L的中间点开始找寻时, 则也合理假设搜寻者C分布在路径L的概率密度函数pc呈三角形分布。

以下利用具体实施例来说明上述搜救方法的具体应用。

设:路径L的长度l=10,能相遇的最大距离md=2,分析搜寻者C能随 机找到走失者D的概率以及在何处找到的概率最大。

实例1:

步骤1:当仅仅知道走失者D最后出现在线路L的中间点且只在L上作自 由移动,可以合理假设D分布在L的概率密度函数pd为三角形分布,即 pd(y)=y25,0y510-y25,5<y10.又,知道搜寻者C在L上随机移动,可合理假设C 分布在L的概率密度函数pc为均匀分布,即pc(x)=0.1,x、y分别表示C、D 的位置点。

步骤2:根据公式(2),C位于一个点xk时能与D相遇的概率在xk={1,2,3,4,5,6,7,8,9}时有

p(Emeet|xk=1)=pc(xk)dx0yxk+2pd(y)dy=dx10×03y25dy=4.5250dx.

p(Emeet|xk=2)=pc(xk)dx0yxk+2pd(y)dy=dx10×04y25dy=8250dx.

p(Emeet|xk=3)=pc(xk)dxxk-2yxk+2pd(y)dy=dx10×15y25dy=12250dx.

p(Emeet|xk=4)=pc(xk)dxxk-2yxk+2pd(y)dy=dx10×(25y25dy+5610-y25dy)=15250dx.

p(Emeet|xk=5)=pc(xk)dxxk-2yxk+2pd(y)dy=dx10×(35y25dy+5710-y25dy)=16250dx.

p(Emeet|xk=6)=pc(xk)dxxk-2yxk+2pd(y)dy=dx10×(45y25dy+5810-y25dy)=15250dx.

p(Emeet|xk=7)=pc(xk)dxxk-2yxk+2pd(y)dy=dx10×5910-y25dy=12250dx.

p(Emeet|xk=8)=pc(xk)dxxk-2yxk+2pd(y)dy=dx10×61010-y25dy=8250dx.

p(Emeet|xk=9)=pc(xk)dxxk-2yxk+2pd(y)dy=dx10×71010-y25dy=4.5250dx.

这样,搜寻者C位于点时x=5时能遇见走失者D的概率最大(图6)。

步骤3:由于pc、pd都关于L的中间点对称,因此有 x[0,5]pc(x)|y-x|md,yLpd(y)dydx=x[5,10]pc(x)·|y-x|md,yLpd(y)dydx;这样,根据公式(3),有p(Emeet)=2x[0,5]pc(x)·|y-x|md,yLpd(y)dydx.

其中,

综上,相遇概率p(Emeet)=2×0.19467=0.38933。

实例2:

步骤1:当仅仅知道走失者D最后出现在线路L的中间点且只在L上作自 由移动,可以合理假设D分布在L的概率密度函数pd为三角形分布,即 pd(y)=y25,0y510-y25,5<y10.又,知道搜寻者C从L的中间点开始找寻,也可合 理假设C分布在L的概率密度函数pc为三角形分布,即 pc(x)=x25,0x510-x25,5<x10,x、y分别表示C、D的位置点。

步骤2:根据公式(2),C位于一个点xk时能与D相遇的概率在xk={1,2,3,4,5,6,7,8,9}时有

p(Emeet|xk=1)=pc(xk)dx0yxk+2pd(y)dy=dx25×03y25dy=4.5625dx.

p(Emeet|xk=2)=pc(xk)dx0yxk+2pd(y)dy=2dx25×04y25dy=16625dx.

p(Emeet|xk=3)=pc(xk)dxxk-2yxk+2pd(y)dy=3dx25×15y25dy=36625dx.

p(Emeet|xk=4)=pc(xk)dxxk-2yxk+2pd(y)dy=4dx25×(25y25dy+5610-y25dy)=60625dx.

p(Emeet|xk=5)=pc(xk)dxxk-2yxk+2pd(y)dy=5dx25×(35y25dy+5710-y25dy)=80625dx.

p(Emeet|xk=6)=pc(xk)dxxk-2yxk+2pd(y)dy=10-625dx×(45y25dy+5810-y25dy)=60625dx.

p(Emeet|xk=7)=pc(xk)dxxk-2yxk+2pd(y)dy=10-725dx×5910-y25dy=36625dx.

p(Emeet|xk=8)=pc(xk)dxxk-2yxk+2pd(y)dy=10-825dx×61010-y25dy=16625dx.

p(Emeet|xk=9)=pc(xk)dxxk-2yxk+2pd(y)dy=10-925dx×71010-y25dy=4.5625dx.

这样,搜寻者C位于点时x=5时能遇见走失者D的概率最大(图7)。

步骤3:由于pc、pd都关于L的中间点对称,因此有 x[0,5]pc(x)·|y-x|md,yLpd(y)dydx=x[5,10]pc(x)·|y-x|md,yLpd(y)dydx;这样,根据公式(3),有p(Emeet)=2x[0,5]pc(x)·|y-x|md,yLpd(y)dydx.

其中,

综上,相遇概率p(Emeet)=2×(0.0181333+0.0405333+0.1898667)=0.49707。

由以上实例可知,基于微分方法的相遇概率具有如下特点:①相遇概率 p(Emeet)完全由可相遇距离md和移动对象本身的概率密度函数pc(x)、pd(y)决定, 与算法过程中的变量无关,因而具有稳定性和唯一性;②搜寻者C在路径L 上的不同位置点xk可遇见走失者D的概率p(Emeet|xk)是xk的函数,存在众数 xm或最大概率值对应的位置点xm,即C在点xk可遇见D的概率最大;③在走 失者D的概率密度函数不变时,搜寻者C采用均匀分布、三角形分布的不同 策略,将出现不同的搜寻成功率,前者为0.38933,后者为0.49707,后者较前 者高出21.68%,因此相遇概率是搜救方案的精准规划和优化的前提和基础。

应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进 或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号