包含噪声的测量仪测量模型

\qquad 测距仪的测量模型可以近似激光测距仪和超声波传感器。这类模型大致有四类误差包括小的测量噪声、意外对象引起的误差、未检测到对象引起的误差和随机意外噪声。
\qquad 1. 局部测量噪声 误差受到测距传感器的有限分辨率、大气对测量信号的影响等引起的。测量噪声通常由一个窄的均值为 z t k ∗ z_{t}^{k*} ztk(真实距离)、标准偏差为 σ h i t \sigma_{hit} σhit的高斯建模,标准偏差是测量模型的一个固有的噪声参数。使用 K K K表示在一个测量 z t z_{t} zt内的测量值数目,既有 z t = { z t 1 , z t 2 , ⋯   , z t K } z_{t}=\{z_{t}^{1},z_{t}^{2},\cdots,z_{t}^{K}\} zt={zt1,zt2,,ztK}, z t k z_{t}^{k} ztk表示一个独立的测量,使用 p h i t p_{hit} phit表示带有噪声的测量概率密度函数。测距传感器的值局限于 [ 0 ; z m a x ] [0;z_{max}] [0;zmax],其中 z m a x z_{max} zmax表示最大的传感器距离。概率密度函数为:
p h i t ( z t k ∣ x t , m ) = { η N ( z t k ; z t k ∗ , σ h i t 2 ) 0 ≤ z t k ≤ z m a x 0 o t h e r 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(ztkxt,m)={ηN(ztk;ztk,σhit2)0ztkzmax0other \qquad 单变量正态分布
N ( z t k ; z t t ∗ , σ h i t 2 ) = 1 2 π σ h i t 2 e − 1 2 z t k − z t k ∗ σ h i t 2 \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πσhit2 1e21σhit2ztkztk \qquad 归一化因子 η \eta η
η = ( ∫ 0 z m a x N ( z t k ; z t k ∗ , σ h i t 2 ) d z t k ) − 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 式中 z t k ∗ z_{t}^{k*} ztk通过 x t x_{t} xt m m m经过射线投射计算。
在这里插入图片描述
\qquad 2. 意外对象 移动机器人的环境是动态的,典型的移动对象就是与机器人共享操作空间的人,处理这类对象的一种方法是将它们作为状态向量的一部分来对待并估计它们的位置,另一种方法是将它们作为传感噪声来测量
\qquad 作为传感器噪声来处理,未建模对象具有这样的特性,即它们会导致比 z t k ∗ z_{t}^{k*} ztk更短而不是更长的距离!!!!!,检测到意外对象的可能性随着距离而减少。这种情况测量距离的概率使用指数分布来表示,该分布的参数 λ s h o r t \lambda_{short} λshort是测量模型的固有参数,根据指数分布的定义得到概率密度函数为
p s h o r t ( z t k ∣ x t , m ) = { η λ s h o r t e − λ s h o r t z t k 0 ≤ z t k ≤ z t k ∗ 0 o t h e r 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(ztkxt,m)={ηλshorteλshortztk0ztkztk0other \qquad 推导得
η = 1 1 − e − λ s h o r t z t k ∗ \eta=\frac{1}{1-e^{-\lambda_{short}z_{t}^{k*}}} η=1eλshortztk1 \qquad 该密度在 z t k z_{t}^{k} ztk范围内指数减少。
在这里插入图片描述
\qquad 3. 检测失败 环境中障碍会被完全忽略。传感器检测失败的典型结果是最大测量距离,传感器返回它的最大允许值 z m a x z_{max} zmax,概率密度为:
p m a x ( z t k ∣ x t , m ) = I ( z = z m a x ) = { 1 z = z m a x 0    o t h e r 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(ztkxt,m)=I(z=zmax)={1z=zmax0  other式中, I I I为一个指示函数,技术上, p m a x p_{max} pmax不具有一个概率密度函数。 p m a x p_{max} pmax是一个离散分布,下图中 p m a x p_{max} pmax使用一个很窄的以 z m a x z_{max} zmax为中心的均匀分布,假装密度的存在。
在这里插入图片描述
\qquad 4. 随机测量 测距仪偶尔产生完全无法解释的测量,该测量模型使用均匀分布在传感器测量范围 [ 0 ; z m a x ] [0;z_{max}] [0;zmax]来建立模型。
p r a n d ( z t k ∣ x t , m ) = { 1 z m a x 0 ≤ z t k ≤ z m a x    0 o t h e r 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(ztkxt,m)={zmax10ztkzmax  0other \qquad 密度函数为:
在这里插入图片描述
\qquad 将4种不同的分布通过四个参数 z h i t 、 z s h o r t 、 z m a x z_{hit}、z_{short}、z_{max} zhitzshortzmax z r a n d z_{rand} zrand进行加权平均混合,并且 z h i t + z s h o r t + z m a x + z r a n d = 1 z_{hit}+z_{short}+z_{max}+z_{rand}=1 zhit+zshort+zmax+zrand=1,有
p ( z t k ∣ x t , m ) = [ z h i t z s h o r t z m a x z r a n d ] T ⋅ [ p h i t ( z t k ∣ x t , m ) p s h o r t ( z t k ∣ x t , m ) p m a x ( z t k ∣ x t , m ) p r a n d ( z t k ∣ x t , 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(ztkxt,m)=zhitzshortzmaxzrandTphit(ztkxt,m)pshort(ztkxt,m)pmax(ztkxt,m)prand(ztkxt,m) \qquad 加权平均混合密度函数为:
在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值