\qquad测距仪的测量模型可以近似激光测距仪和超声波传感器。这类模型大致有四类误差包括小的测量噪声、意外对象引起的误差、未检测到对象引起的误差和随机意外噪声。
\qquad1. 局部测量噪声 误差受到测距传感器的有限分辨率、大气对测量信号的影响等引起的。测量噪声通常由一个窄的均值为ztk∗z_{t}^{k*}ztk∗(真实距离)、标准偏差为σhit\sigma_{hit}σhit的高斯建模,标准偏差是测量模型的一个固有的噪声参数。使用KKK表示在一个测量ztz_{t}zt内的测量值数目,既有zt={zt1,zt2,⋯ ,ztK}z_{t}=\{z_{t}^{1},z_{t}^{2},\cdots,z_{t}^{K}\}zt={zt1,zt2,⋯,ztK},ztkz_{t}^{k}ztk表示一个独立的测量,使用phitp_{hit}phit表示带有噪声的测量概率密度函数。测距传感器的值局限于[0;zmax][0;z_{max}][0;zmax],其中zmaxz_{max}zmax表示最大的传感器距离。概率密度函数为:
phit(ztk∣xt,m)={ηN(ztk;ztk∗,σhit2)0≤ztk≤zmax0other
p_{hit}(z_{t}^{k}|x_{t},m)=
\begin{cases}
\eta\mathcal{N}(z_{t}^{k};z_{t}^{k*},\sigma_{hit}^{2})\qquad0\leq z_{t}^{k}\leq z_{max}\\
\qquad\quad0\qquad\quad\qquad\qquad other
\end{cases}
phit(ztk∣xt,m)={ηN(ztk;ztk∗,σhit2)0≤ztk≤zmax0other\qquad单变量正态分布
N(ztk;ztt∗,σhit2)=12πσhit2e−12ztk−ztk∗σhit2
\mathcal{N}(z_{t}^{k};z_{t}^{t*},\sigma_{hit}^{2})=\frac{1}{\sqrt{2\pi\sigma_{hit}^{2}}}e^{-\frac{1}{2}\frac{z_{t}^{k}-z_{t}^{k*}}{\sigma_{hit}^{2}}}
N(ztk;ztt∗,σhit2)=2πσhit21e−21σhit2ztk−ztk∗\qquad归一化因子η\etaη为
η=(∫0zmaxN(ztk;ztk∗,σhit2)dztk)−1
\eta=(\int_{0}^{z_{max}}\mathcal{N}(z_{t}^{k};z_{t}^{k*},\sigma_{hit}^{2})dz_{t}^{k
})^{-1}
η=(∫0zmaxN(ztk;ztk∗,σhit2)dztk)−1\qquad式中ztk∗z_{t}^{k*}ztk∗通过xtx_{t}xt和mmm经过射线投射计算。
\qquad 2. 意外对象 移动机器人的环境是动态的,典型的移动对象就是与机器人共享操作空间的人,处理这类对象的一种方法是将它们作为状态向量的一部分来对待并估计它们的位置,另一种方法是将它们作为传感噪声来测量。
\qquad作为传感器噪声来处理,未建模对象具有这样的特性,即它们会导致比ztk∗z_{t}^{k*}ztk∗更短而不是更长的距离!!!!!,检测到意外对象的可能性随着距离而减少。这种情况测量距离的概率使用指数分布来表示,该分布的参数λshort\lambda_{short}λshort是测量模型的固有参数,根据指数分布的定义得到概率密度函数为
pshort(ztk∣xt,m)={ηλshorte−λshortztk0≤ztk≤ztk∗0other
p_{short}(z_{t}^{k}|x_{t},m)=\begin{cases}
\eta\lambda_{short}e^{-\lambda_{short}z_{t}^{k}}\qquad0\leq z_{t}^{k}\leq z_{t}^{k*}\\
\qquad\quad0\qquad\qquad\qquad other
\end{cases}
pshort(ztk∣xt,m)={ηλshorte−λshortztk0≤ztk≤ztk∗0other\qquad推导得
η=11−e−λshortztk∗
\eta=\frac{1}{1-e^{-\lambda_{short}z_{t}^{k*}}}
η=1−e−λshortztk∗1\qquad该密度在ztkz_{t}^{k}ztk范围内指数减少。
\qquad 3. 检测失败 环境中障碍会被完全忽略。传感器检测失败的典型结果是最大测量距离,传感器返回它的最大允许值zmaxz_{max}zmax,概率密度为:
pmax(ztk∣xt,m)=I(z=zmax)={1z=zmax0 other
p_{max}(z_{t}^{k}|x_{t},m)=I(z=z_{max})=\begin{cases}
1\qquad z=z_{max}\\
0\qquad\ \ other
\end{cases}
pmax(ztk∣xt,m)=I(z=zmax)={1z=zmax0 other式中,III为一个指示函数,技术上,pmaxp_{max}pmax不具有一个概率密度函数。pmaxp_{max}pmax是一个离散分布,下图中pmaxp_{max}pmax使用一个很窄的以zmaxz_{max}zmax为中心的均匀分布,假装密度的存在。
\qquad4. 随机测量 测距仪偶尔产生完全无法解释的测量,该测量模型使用均匀分布在传感器测量范围[0;zmax][0;z_{max}][0;zmax]来建立模型。
prand(ztk∣xt,m)={1zmax0≤ztk≤zmax 0other
p_{rand}(z_{t}^{k}|x_{t},m)=\begin{cases}
\frac{1}{z_{max}}\qquad 0\leq z_{t}^{k}\leq z_{max}\\
\ \ 0 \qquad\qquad other
\end{cases}
prand(ztk∣xt,m)={zmax10≤ztk≤zmax 0other\qquad密度函数为:
\qquad将4种不同的分布通过四个参数zhit、zshort、zmaxz_{hit}、z_{short}、z_{max}zhit、zshort、zmax和zrandz_{rand}zrand进行加权平均混合,并且zhit+zshort+zmax+zrand=1z_{hit}+z_{short}+z_{max}+z_{rand}=1zhit+zshort+zmax+zrand=1,有
p(ztk∣xt,m)=[zhitzshortzmaxzrand]T⋅[phit(ztk∣xt,m)pshort(ztk∣xt,m)pmax(ztk∣xt,m)prand(ztk∣xt,m)]
p({z_{t}^{k}|x_{t},m})=\begin{bmatrix}
z_{hit} \\
z_{short} \\
z_{max} \\
z_{rand}
\end{bmatrix}^{T}\cdot
\begin{bmatrix}
p_{hit}(z_{t}^{k}|x_{t},m) \\
p_{short}(z_{t}^{k}|x_{t},m) \\
p_{max}(z_{t}^{k}|x_{t},m) \\
p_{rand}(z_{t}^{k}|x_{t},m)
\end{bmatrix}
p(ztk∣xt,m)=⎣⎢⎢⎡zhitzshortzmaxzrand⎦⎥⎥⎤T⋅⎣⎢⎢⎡phit(ztk∣xt,m)pshort(ztk∣xt,m)pmax(ztk∣xt,m)prand(ztk∣xt,m)⎦⎥⎥⎤\qquad加权平均混合密度函数为: