\qquad文中主要研究具有区域约束的机器人网络执行静态最优覆盖。给定密度函数描述事件发生的概率,执行函数(performance functionperformance\ functionperformance function)代表到定位点的代价。目标是通过合理放置传感器使最小化探测代价。此外,由于平衡负载约束,预先设定给每个机器人的分配的区域。
一般性泰森分区(Generalized Voronoi Partitions)
\qquad 对于凸多边形Q⊂RdQ\subset\mathbb{R}^dQ⊂Rd和P=(p1,p2,...,pn)∈QnP=(p_{1},p_{2},...,p_{n})\in Q^{n}P=(p1,p2,...,pn)∈Qn,泰森分区为V(P)={V1(P),...,Vn(P)}\mathcal{V}(P)=\{V_{1}(P),...,V_{n}(P)\}V(P)={V1(P),...,Vn(P)},定义Delaunay graph\color{#F00}{Delaunay\ graph}Delaunay graph,对于相邻端点pi和pjp_{i}和p_{j}pi和pj如果Vi(P)∩Vj(P)≠ϕV_{i}(P)\cap V_{j}(P)\neq\phiVi(P)∩Vj(P)̸=ϕ,给定性能函数f:R→Rf:\mathbb{R}\rightarrow\mathbb{R}f:R→R严格递增,给定权重ω=(ω1,...,ωn)∈Rn\omega=(\omega_{1},...,\omega_{n})\in\mathbb{R}^{n}ω=(ω1,...,ωn)∈Rn,得出一般性泰森分区公式为:
Vif(P,ω)={q∈Q∣f(∣∣q−pi∣∣)−ωi≥f(∣∣q−pj∣∣)−ωj,i≠j}
V_{i}^{f}(P,\omega)=\{q\in Q|f(||q-p_{i}||)-\omega_{i}\ge f(||q-p_{j}||)-\omega_{j},i\ne j\}
Vif(P,ω)={q∈Q∣f(∣∣q−pi∣∣)−ωi≥f(∣∣q−pj∣∣)−ωj,i̸=j}一般的,一般化泰森分区既不是凸区域也不是星形。
区域约束下的定位最优问题
\qquad考虑区域约束条件多中心最优化问题,即在多中心最优问题(minimize H(p1,...,pn,W1,...,Wn)minimize\ \mathcal{H}(p_{1},...,p_{n},W_{1},...,W_{n})minimize H(p1,...,pn,W1,...,Wn))中增加可行域条件约束(aia_{i}ai)。给定可行域集{a1,...,an}⊂R>0\{a_{1},...,a_{n}\}\subset\mathbb{R}_{>0}{a1,...,an}⊂R>0满足∑i=1nai=∫Qϕ(q)dq=areaϕ(Q)\sum_{i=1}^{n}a_{i}=\int_{Q}\phi(q)dq=area_{\phi}(Q)∑i=1nai=∫Qϕ(q)dq=areaϕ(Q)。即:
minimizeH(p1,...,pn,W1,...,Wn)subject to∫Wiϕ(q)dq=ai,i∈{1,...,n}
minimize\quad\mathcal{H}(p_{1},...,p_{n},W_{1},...,W_{n})\\
subject\ to\quad \int_{W_{i}}\phi(q)dq=a_{i},i\in\{1,...,n\}
minimizeH(p1,...,pn,W1,...,Wn)subject to∫Wiϕ(q)dq=ai,i∈{1,...,n}等分域情形为ai=1/n∫Qϕ(q)dq for i∈{1,...,n}a_{i}=1/n\int_{Q}\phi(q)dq\ for\ i\in\{1,...,n\}ai=1/n∫Qϕ(q)dq for i∈{1,...,n}。
基于一般化泰森分区的权重-区域映射
\qquad定义权重-区域映射为M:U⊂Rn→Rn\mathcal{M}:U\subset\mathbb{R}^{n}\rightarrow\mathbb{R}^{n}M:U⊂Rn→Rn:
M(ω)=(∫V1f(P,ω)ϕ(q)dq,...,∫Vnf(P,ω)ϕ(q)dq)
\mathcal{M}(\omega)=(\int_{V_{1}^{f}(P,\omega)}\phi(q)dq,...,\int_{V_{n}^{f}(P,\omega)}\phi(q)dq)
M(ω)=(∫V1f(P,ω)ϕ(q)dq,...,∫Vnf(P,ω)ϕ(q)dq)其中,P={p1,...,pn}P=\{p_{1},...,p_{n}\}P={p1,...,pn}。令M\mathcal{M}M为梯度,∇F=−M\nabla F=-\mathcal{M}∇F=−M(???\color{#F00}{???}???),其中F:Rn→RF:\mathbb{R}^{n}\rightarrow\mathbb{R}F:Rn→R
F(ω)=∑j=1n∫Vjf(P,ω)(f(∣∣q−pj∣∣)−ωj)ϕ(q)dq
F(\omega)=\sum_{j=1}^{n}\int_{V_{j}^{f}(P,\omega)}(f(||q-p_{j}||)-\omega_{j})\phi(q)dq
F(ω)=j=1∑n∫Vjf(P,ω)(f(∣∣q−pj∣∣)−ωj)ϕ(q)dq\qquad若通过配置权值ωi\omega_{i}ωi使得M(ωi)=ai\mathcal{M}(\omega_{i})=a_{i}M(ωi)=ai,几乎不能直接求解映射的解析解,使用的方法为:定义函数g(ω1,...,ωn)=M(ω1,...,ωn)−(a1,...,an)g(\omega_{1},...,\omega_{n})=\mathcal{M}(\omega_{1},...,\omega_{n})-(a_{1},...,a_{n})g(ω1,...,ωn)=M(ω1,...,ωn)−(a1,...,an),其原函数为:F:Rn→Rn,F=−F(ω)−∑i=1nωiai\mathcal{F}:\mathbb{R}^{n}\rightarrow\mathbb{R}^{n},\mathcal{F}=-F(\omega)-\sum_{i=1}^{n}\omega_{i}a_{i}F:Rn→Rn,F=−F(ω)−∑i=1nωiai。寻找ω∈U\omega\in Uω∈U,使得:
∇F(ω)=g(ω)=0(vector)n
\nabla\mathcal{F}(\omega)=g(\omega)=0(vector)_{n}
∇F(ω)=g(ω)=0(vector)n\qquad所寻找的ω\omegaω即优化函数F\mathcal{F}F(多中心最优化问题???\color{#F00}{???}???)又使得泰森区域满足约束条件。使用雅克比(Jacobi)迭代算法配置权重。\color{#F00}{使用雅克比(Jacobi)迭代算法配置权重。}使用雅克比(Jacobi)迭代算法配置权重。
\qquad雅克比算法(生成线性方程的JORJORJOR算法)
x(t+1)=x(t)−γ[D(x(t))]−1∇F(x(t))
x(t+1)=x(t)-\gamma[D(x(t))]^{-1}\nabla F(x(t))
x(t+1)=x(t)−γ[D(x(t))]−1∇F(x(t))\qquad其中,γ\gammaγ是步长,D(x)D(x)D(x)是对角矩阵,其第i个对角元素为∇ii2F(x)\nabla_{ii}^{2}F(x)∇ii2F(x),假设其对角元素都是非零元素。
\qquad应用雅克比算法迭代权值为:
ωk+1=ωk−γdiag(∂g1∂ω1(ωk),...,∂gn∂ωn(ωk))−1g(ωk)
\omega_{k+1}=\omega_{k}-\gamma diag(\frac{\partial g_{1}}{\partial \omega_{1}}(\omega_{k}),...,\frac{\partial g_{n}}{\partial \omega_{n}}(\omega_{k}))^{-1}g(\omega_{k})
ωk+1=ωk−γdiag(∂ω1∂g1(ωk),...,∂ωn∂gn(ωk))−1g(ωk)\qquad 使用质量守恒(conservation−of−mass)定理\color{#F00}{使用质量守恒(conservation-of-mass)定理}使用质量守恒(conservation−of−mass)定理求解:
∂gi∂ωj=∫∂Viϕ ni∂q∂ωjdq
\frac{\partial g_{i}}{\partial \omega_{j}}=\int_{\partial V_{i}}\phi\ n_{i}\frac{\partial q}{\partial\omega_{j}}dq
∂ωj∂gi=∫∂Viϕ ni∂ωj∂qdq