CN103237345B - 一种基于二元量化数据的声源目标迭代定位方法 - Google Patents
一种基于二元量化数据的声源目标迭代定位方法 Download PDFInfo
- Publication number
- CN103237345B CN103237345B CN201310122120.4A CN201310122120A CN103237345B CN 103237345 B CN103237345 B CN 103237345B CN 201310122120 A CN201310122120 A CN 201310122120A CN 103237345 B CN103237345 B CN 103237345B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- target
- msubsup
- theta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000004807 localization Effects 0.000 title claims abstract description 10
- 239000002245 particle Substances 0.000 claims abstract description 26
- 239000000203 mixture Substances 0.000 claims abstract description 22
- 238000013139 quantization Methods 0.000 claims abstract description 15
- 238000009499 grossing Methods 0.000 claims abstract description 12
- 230000005236 sound signal Effects 0.000 claims abstract description 9
- 230000004927 fusion Effects 0.000 claims abstract description 5
- 238000009826 distribution Methods 0.000 claims description 28
- 238000005070 sampling Methods 0.000 claims description 17
- 238000005315 distribution function Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 8
- 238000012544 monitoring process Methods 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 5
- 238000013178 mathematical model Methods 0.000 claims description 5
- 230000000295 complement effect Effects 0.000 claims description 4
- 239000006185 dispersion Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000005265 energy consumption Methods 0.000 abstract description 2
- 230000002776 aggregation Effects 0.000 abstract 2
- 238000004220 aggregation Methods 0.000 abstract 2
- 238000001914 filtration Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 4
- 230000008602 contraction Effects 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于二元量化数据的声源目标迭代定位方法,其利用若干无线声传感器节点模块和汇聚节点构成目标定位无线传感器网络,基于声传感器节点接收到的目标声音信号强度二元量化数据,利用粒子滤波器迭代估计目标位置。本发明首先将传感器节点接收到的声源信号强度数据量化压缩为一位二进制数据,然后采用粒子滤波与高斯混合核平滑方法对目标位置进行迭代估计,汇聚节点每次只处理一个传感器节点的一位二进制量化数据,降低了对无线传感器网络带宽和频率的要求,减少了网络能量消耗。本发明通过多传感器信息融合获得较高的目标定位精度,在传感器网络能量和带宽等制约下,实现了实时目标定位。
Description
技术领域
本发明属于无线传感器网络技术领域,涉及一种声源目标迭代定位方法,尤其是一种基于二元量化数据的声源目标迭代定位方法。
背景技术
无线传感器网络近年来成为研究热点,其由大量低功耗的无线传感器节点组成,通过节点之间的协同可以完成许多复杂的任务,目标定位是无线传感器网络的一个重要应用,其应用领域包括物流跟踪、环境监测、网络控制等。精确的目标定位在实际中具有非常重要的意义,目前国内已有很多关于目标定位方法的研究。
国内专利号为CN201110435631.2的一种粒子状态估计方法,通过观测数据对目标状态进行估计,得到状态信息,将得到的状态信息用于粒子滤波器的重要性采样过程,对粒子权值进行修正,克服了由于滤波初始阶段先验信息不足而对粒子最初先验状态进行错误估计进而导致滤波不稳定甚至发散的问题,但是该方法在重采样阶段复制高权值的粒子,抛弃低权值的粒子,可能会丧失粒子的多样性,导致定位不准确。
国内专利号为CN200980150679.X的使用粒子滤波估计声源地点的方法,所述方法利用声学传递函数之间的相关性,该相关性可以被用作粒子滤波的测量函数的伪似然函数,确定粒子的权重,将音频信息权重与视频信息权重组合产生多模式的粒子滤波方法,该方法没有考虑到当前可以利用的观测信息,使粒子退化现象比较严重。
国内专利号为CN02150140.8的一种位置估计方法,所述方法利用实时获取的非直视(NLOS)信号误差的均值对到达时间(TOA)测量中的直视(LOS)信号误差进行零均值矫正,然后结合NLOS识别结果,利用由NLOS信号误差的概率密度函数和参数确定NLOS信号误差的方差,对位置估计方法中的加权矩阵进行调整,实现对NLOS信号误差的抑制,最后再对多个位置估计结果进行平均,进一步抑制NLOS信号误差,但是该方法不能有效的抑制测距噪声的影响,不能很好的确定位置估计的准确性。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供一种基于二元量化数据的声源目标迭代定位方法,该方法利用无线声传感器节点模块和汇聚节点构成目标定位无线传感器网络,将收到声源目标信号强度数据量化压缩为二进制数据,利用粒子滤波器迭代估计目标未知参数的概率分布,从而得到目标位置估计。该方法能够减少通信量,同时提高了目标定位的精度。
本发明的目的是通过以下技术方案来解决的:
这种基于二元量化数据的声源目标迭代定位方法,利用若干无线声传感器节点模块和汇聚节点构成目标定位无线传感器网络,将传感器节点接收到声源目标信号强度数据量化压缩为二进制数据,利用粒子滤波器迭代估计目标位置参数的概率分布,从而得到目标位置估计。
进一步,以上将传感器节点接收到声源目标信号强度数据量化压缩为一位二进制数据。
进一步,以上利用粒子滤波器迭代估计目标位置参数的概率分布的方法包括以下步骤:
(2.1)设未知参数θ=[P0,x,y]包含未知的声源信号强度参数P0和目标位置坐标(x,y);在初始迭代i=0,从区间[0,PM]上均匀概率分布函数中采样声源信号强度P0的样本,其中PM为声音信号的最大值,在整个无线监测区域内均匀概率分布函数中采样目标位置(x,y)的样本,得到M个目标未知参数样本θ0 (m),m=1,2,…,M,给每个粒子分配相等的权值得到初始样本集{θ0 (m),w0 (m)|m=1,…,M};
(2.2)在后续迭代步骤中,汇聚节点依次序贯激活每一个传感器节点;设在第i+1次迭代中,节点i+1被激活,将其观测到的目标声音强度数据量化压缩为di+1=0或者di+1=1,并传递给汇聚节点;
(2.3)基于高斯混合核平滑方法,汇聚节点利用蒙特卡洛方法从核密度函数的M个高斯分量N(μi (m),b2Vi)中分别采样得到M个新的目标参数样本其中:N(μ,V)代表均值为μ,协方差矩阵为V的高斯分布概率密度函数,为高斯混合分布中第m个核函数的均值,为样本集合{θi (m),wi (m)|m=1,…,M}的协方差矩阵,常数α和b分别表示高斯混合分布的收缩和分散程度,其取值依赖于折衷因子δ,b2=1-[(3δ-1)2δ]2,0.95≤δ≤0.99,为上次迭代获得的样本集合{θi (m),wi (m)|m=1,…,M}的均值;
(2.4)汇聚节点利用传感器节点i+1的量化数据di+1=0或者di+1=1计算每个样本的权值
(2.5)汇聚节点利用权值对所得样本集合进行重采样,得到一组加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M};
(2.6)汇聚节点利用重采样后的加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M},迭代计算目标参数估计值:
(2.7)重复上述步骤(2.1)至(2.6),直到遍历完所有传感器节点。
进一步,设在第i次迭代中,汇聚节点得到的加权样本集为{θi (m),wi (m)|m=1,…,M},利用高斯混合核平滑方法,第i次迭代中目标未知参数θ的概率分布表示为高斯混合分布:
进一步,所述目标参数样本权值的迭代计算公式为:
其中:代表二元量化观测数据di+1的似然函数。
其中参数[P0,x,y]=θi+1 (m),η为预先设定好的量化阈值,σ为传感器观测噪声的协方差,Q(·)为高斯分布的补分布函数,(sx(i+1),sy(i+1))为第i+1个传感器节点的位置坐标。
进一步,在目标定位前,首先建立声音信号强度随目标距离衰减的数学模型:
其中:ai代表第i个声传感器节点接收到的声源信号强度,Gi是第i个传感器节点的增益,P0'为在参考距离d0处声源的信号强度,di为目标和第i个传感器节点的距离,(sxi,syi)和(x,y)分别为第i个传感器节点和目标的位置坐标,n=2为声源信号强度衰减指数。
进一步,以上汇聚节点对每个节点的量化数据进行序贯处理,在每次迭代中,汇聚节点只处理一个传感器节点的二元量化数据,其他传感器节点处于休眠状态。
以上所述基于二元量化数据的声源目标迭代定位方法具体包括下列步骤:
步骤A,建立声音信号强度随目标距离衰减的数学模型:
步骤B,各传感器节点对实际接受到的目标信号强度进行量化压缩处理,只传递一个二进制数字信号到汇聚节点;每个传感器节点的观测数据与一个阈值η进行比较,如果高过阈值,传感器发送信号1到融合中心,否则传感器保持静默,不发送信号;在每次迭代中,汇聚节点序贯激活每个传感器节点,为了节约能量,其余节点处于休眠状态;
步骤C,基于高斯混合核平滑方法,利用蒙特卡洛方法从核密度的高斯分量中分别采样新的声源目标参数样本
步骤D,汇聚节点利用接收到的传感器节点二元量化观测数据,计算观测数据的似然函数作为样本的权值并归一化;
步骤E,汇聚节点利用样本权值对所得的样本集合进行重采样;
步骤F,汇聚节点利用所得到的加权样本集合,以序贯方式迭代计算目标的位置估计。
本发明具有以下有益效果:
本发明基于二元量化数据的声源目标迭代定位方法,其传感器观测模式简单,讯通量小,具有较小的计算量,对通讯带宽和频率要求低,具有较高的定位精度。利用粒子滤波算法以及高斯混合核平滑方法,能有效克服粒子的退化现象,加速滤波器收敛,能够完全满足传感器网络目标定位的实时性要求。
附图说明
图1为本发明基于量化数据的声源目标迭代定位方法流程图;
图2为监测区域内无线传感器节点的部署及目标的真实位置示意图;
图3为利用本发明的目标定位方法进行迭代定位的效果图。
图4为利用本发明的目标定位方法对目标位置估计的均方根误差(RMSE)图。
具体实施方式
本发明基于二元量化数据的声源目标迭代定位方法,是利用若干无线声传感器节点模块和汇聚节点构成目标定位无线传感器网络,将传感器节点接收到声源目标信号强度数据量化压缩为二进制数据,利用粒子滤波器迭代估计目标位置参数的概率分布,从而得到目标位置估计。其中,利用粒子滤波器迭代估计目标位置参数的概率分布的方法包括以下步骤:
(2.1)设未知参数q=[P0,x,y]包含未知的声源信号强度参数P0和目标位置坐标(x,y);在初始迭代i=0,从区间[0,PM]上均匀概率分布函数中采样声源信号强度P0的样本,其中PM为声音信号的最大值,在整个无线监测区域内均匀概率分布函数中采样目标位置(x,y)的样本,得到M个目标未知参数样本θ0 (m),m=1,2,…,M,给每个粒子分配相等的权值得到初始样本集{θ0 (m),w0 (m)|m=1,…,M};
(2.2)在后续迭代步骤中,汇聚节点依次序贯激活每一个传感器节点;设在第i+1次迭代中,节点i+1被激活,将其观测到的目标声音强度数据量化压缩为di+1=0或者di+1=1,并传递给汇聚节点;
(2.3)基于高斯混合核平滑方法,汇聚节点利用蒙特卡洛方法从核密度函数的M个高斯分量N(μi (m),b2Vi)中分别采样得到M个新的目标参数样本其中:N(μ,V)代表均值为μ,协方差矩阵为V的高斯分布概率密度函数,为高斯混合分布中第m个核函数的均值,为样本集合{θi (m),wi (m)|m=1,…,M}的协方差矩阵,常数α和b分别表示高斯混合分布的收缩和分散程度,其取值依赖于折衷因子δ,b2=1-[(3δ-1)2δ]2,0.95≤δ≤0.99,为上次迭代获得的样本集合{θi (m),wi (m)|m=1,…,M}的均值;
(2.4)汇聚节点利用传感器节点i+1的量化数据di+1=0或者di+1=1计算每个样本的权值
(2.5)汇聚节点利用权值对所得样本集合进行重采样,得到一组加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M};
(2.6)汇聚节点利用重采样后的加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M},迭代计算目标参数估计值:
(2.7)重复上述步骤(2.1)至(2.6),直到遍历完所有传感器节点。
以下结合附图对本发明进行详细说明:
本发明的基于二元量化数据的声源目标迭代定位方法的具体实施步骤如下:
步骤A,建立声音信号强度随目标距离衰减的数学模型,
其中ai是第i个声传感器节点接收到的声源信号幅度,Gi是第i个传感器节点的增益,P0'为在参考距离d0处声源的强度,di为目标和第i个传感器节点的距离,(sxi,syi)和(x,y)分别为第i个传感器节点和声源目标的坐标,n=2为信号衰减指数。为了简单起见,假设对所有的传感器Gi=G,且P0=GP0',d0=1。所以信号强度衰减模型可以简化为
考虑进建模误差和背景噪声,第i个传感器实际观测到的来自目标的声信号强度观测值为:
zi=ai+wi(3)
这里假设噪声wi对所有的传感器节点是独立同分布的,且服从高斯分布,即wi~N(0,σ2),σ为观测噪声的协方差。设目标未知参数θ=[P0,x,y],三个分量包含未知待估计的声源信号强度参数P0和目标位置坐标(x,y)。
步骤B,为了减少网络通讯量从而减少能量消耗和带宽,各传感器节点对实际接受到的目标信号强度进行量化压缩处理,只传递一个二进制数字信号到汇聚节点。每个传感器节点的观测数据与一个设定的阈值η进行比较,如果高过阈值,传感器发送信号1到融合中心,否则传感器保持静默,不发送信号。由观测模型可知,传感器的观测值在区间范围内,将量化阈值设为设传感器节点i的实际观测值为zi,则量化后的二元观测数据为:
在每次迭代中,汇聚节点序贯激活每个传感器节点,为了节约能量,其余节点处于休眠状态。
步骤C,在迭代开始时刻i=0,从区间[0,PM](PM为声源信号的最大值)上均匀概率分布函数中采样声源信号强度P0的样本,在整个无线监测区域内均匀概率分布函数中采样目标位置坐标(x,y)的样本,得到M个样本粒子θ0 (m),并且分配相等的权值得到初始样本集{θ0 (m),w0 (m)|m=1,…,M}。设在第i次迭代中得到的加权样本集合为在i+1次迭代中,汇聚节点基于高斯混合核平滑方法,利用蒙特卡洛方法从核密度的高斯分量中分别采样得到新的声源参数样本具体方法如下:
使用上一次迭代的样本计算高斯混合核函数的参数:
从核密度的高斯分量N(μi (m),b2Vi)中分别采样目标位置的样本m=1,2,…,M。其中:为高斯混合分布中第m个核函数的均值,为上次迭代获得的样本集合{θi (m),wi (m)|m=1,…,M}的均值,为样本集合{θi (m),wi (m)|m=1,…,M}的协方差矩阵;
步骤D,汇聚节点利用接收到的传感器节点i+1的二元量化数据di+1,迭代计算每个样本的权值并归一化
其中:代表二元量化观测数据di+1的似然函数。
其中参数[P0,x,y]=θi+1 (m),η为预先设定好的量化阈值,σ为传感器观测噪声的协方差,Q(·)为高斯分布的补分布函数,(sx(i+1),sy(i+1))为第i+1个传感器节点的位置坐标。
步骤E,汇聚节点利用样本权值对所得的样本集合进行重采样,重采样后的加权样本集合为
步骤F,汇聚节点利用所得到的重采样后的加权样本集合 迭代计算目标的参数估计值:
转入步骤B,进入下次迭代循环,直到遍历完所有传感器节点。
本发明采用粒子滤波器与高斯混合核平滑方法对目标位置进行迭代估计,汇聚节点每次只需处理一个传感器节点的一位二进制量化数据,降低了对无线传感器网络带宽和频率的要求,减少了网络能量消耗。与传统的目标定位方法相比,本发明通过多传感器信息融合获得较高的目标定位精度,在传感器网络能量和带宽等制约下,实现了实时目标定位。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明的基于量化数据的声源目标迭代定位方法进一步进行详细说明。
参见图1,为基于量化数据的声源目标迭代定位方法流程图,该方法主要包括以下步骤:
步骤S101,在初始迭代i=0,从区间[0,PM]上均匀概率分布函数中采样声源信号强度P0的样本,其中PM为声音信号的最大值,在整个无线监测区域内均匀概率分布函数中采样目标位置(x,y)的样本,得到M个目标未知参数样本θ0 (m),m=1,2,…,M,给每个粒子分配相等的权值得到初始样本集{θ0 (m),w0 (m)|m=1,…,M}。设在第i次迭代中,粒子滤波器得到的加权样本集为{θi (m),wi (m)|m=1,…,M},利用高斯混合核平滑方法,第i次迭代中目标未知参数θ的概率分布可以表示为高斯混合和分布:
其中:N(μ,V)代表均值为μ,协方差矩阵为V的高斯分布密度函数, 常数α和b分别表示高斯混合的收缩和分散程度,其取值依赖于一个折衷因子0.95≤δ≤0.99,并且b2=1-[(3δ-1)2δ]2,N(μi (m),b2Vi)为高斯混合分布中的高斯分量核函数。
步骤S102,在后续迭代步骤中,汇聚节点依次序贯激活每一个传感器节点;设在第i+1次迭代中,传感器节点i+1被激活,获得原始观测数据zi+1,设定量化阈值则量化后的二元观测数据为:
并将其二元量化观测数据传递给汇聚节点。
步骤S103,基于高斯混合核平滑方法,汇聚节点利用蒙特卡洛方法从核密度函数的M个高斯分量N(μi (m),b2Vi)中分别采样得到M个新的目标参数样本
具体如下:
使用上一次迭代的样本集计算高斯混合核函数的参数:
其中,为样本的均值,Vi为样本方差。从而得到M个高斯分量核密度函数N(μi (m),b2Vi),从每个高斯分量N(μi (m),b2Vi)中分别采样得到新的目标参数的样本
步骤S104,汇聚节点利用传感器节点i+1的量化数据di+1=0或者di+1=1计算目标参数的样本的权值并归一化:
其中:代表二元量化观测数据di+1的似然函数。
其中参数[P0,x,y]=θi+1 (m),η为预先设定好的量化阈值,σ为传感器观测噪声的协方差,Q(·)为高斯分布的补分布函数,(sx(i+1),sy(i+1))为第i+1个传感器节点的位置坐标。
步骤S105,汇聚节点利用权值对所得样本集合进行重采样,得到一组加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M};
步骤S106,汇聚节点利用重采样后的加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M},迭代计算目标参数估计值:
步骤S107,重复步骤S101~S106,迭代估计目标未知参数,直到遍历完所有传感器节点。
以下以一实例详细说明本发明的基于二元量化数据的声源目标迭代定位方法:
如图2所示,设25个麦克风声传感器节点均匀部署在100×100m2方形区域内,每个传感器节点的位置由自身携带的GPS获得,声源目标的真实位置坐标为(60,60),我们利用所有声传感器节点的量化数据对声源进行定位。
为了节约汇聚节点的能量、带宽和频率,延长网络使用寿命,对每个节点的原始观测量化压缩为一位二进制数据,汇聚节点对二元量化数据进行序贯处理,每次只激活一个传感器节点,其余节点处于休眠状态,直到遍历完所有传感器节点,迭代结束,完成对目标的定位。
设粒子滤波器的样本个数为M=2000,传感器节点观测噪声方差σ2=1。在初始迭代中,声源信号强度P0的初始样本来自于均匀分布U[0,500](U[a,b]代表区间[a,b]上的均匀概率分布函数),目标位置坐标(x,y)的初始样本来自于二维均匀分布U[(0,100),(0,100)],并均匀分配权值为w0 (m)=1/2000。按照本发明所述声源目标定位方法,在整个迭代定位过程中对目标位置估计的效果如图3所示。
经过50次蒙特卡洛仿真,目标位置估计的均方根误差(RMSE)如图4所示,从图中可以看出,随着迭代次数的增加,误差越来越小,实现了精确的目标定位。
以上所述,仅为本发明较佳的具体实施方式,并不用于对本发明进行限制,本发明的保护范围应以权利要求的保护范围为准。
Claims (6)
1.一种基于二元量化数据的声源目标迭代定位方法,其特征在于,利用若干无线声传感器节点模块和汇聚节点构成目标定位无线传感器网络,将传感器节点接收到声源目标信号强度数据量化压缩为二进制数据,利用粒子滤波器迭代估计目标位置参数的概率分布,从而得到目标位置估计;
利用粒子滤波器迭代估计目标位置参数的概率分布的方法包括以下步骤:
(2.1)设未知参数θ=[P0,x,y]包含未知的声源信号强度参数P0和目标位置坐标(x,y);在初始迭代i=0,从区间[0,PM]上均匀概率分布函数中采样声源信号强度P0的样本,其中PM为声音信号的最大值,在整个无线监测区域内均匀概率分布函数中采样目标位置(x,y)的样本,得到M个目标未知参数样本θ0 (m),m=1,2,…,M,给每个粒子分配相等的权值得到初始样本集{θ0 (m),w0 (m)|m=1,…,M};
(2.2)在后续迭代步骤中,汇聚节点依次序贯激活每一个传感器节点;设在第i+1次迭代中,节点i+1被激活,将其观测到的目标声音强度数据量化压缩为di+1=0或者di+1=1,并传递给汇聚节点;
(2.3)基于高斯混合核平滑方法,汇聚节点利用蒙特卡洛方法从核密度函数的M个高斯分量N(μi (m),b2Vi)中分别采样得到M个新的目标参数样本其中:N(μ,V)代表均值为μ,协方差矩阵为V的高斯分布概率密度函数,为高斯混合分布中第m个核函数的均值,为样本集合{θi (m),wi (m)|m=1,…,M}的协方差矩阵,常数α和b分别表示高斯混合分布的收缩和分散程度,其取值依赖于折衷因子δ,b2=1-[(3δ-1)/2δ]2,0.95≤δ≤0.99,为上次迭代获得的样本集合{θi (m),wi (m)|m=1,…,M}的均值;
(2.4)汇聚节点利用传感器节点i+1的量化数据di+1=0或者di+1=1计算每个样本的权值
(2.5)汇聚节点利用权值对所得样本集合进行重采样,得到一组加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M};
(2.6)汇聚节点利用重采样后的加权样本集合{θi+1 (m),wi+1 (m)|m=1,…,M},迭代计算目标参数估计值:
(2.7)重复上述步骤(2.1)至(2.6),直到遍历完所有传感器节点;
所述目标参数样本权值的迭代计算公式为:
其中:代表二元量化观测数据di+1的似然函数:
其中参数[P0,x,y]=θi+1 (m),η为预先设定好的量化阈值,σ为传感器观测噪声的协方差,Q(·)为高斯分布的补分布函数,(sx(i+1),sy(i+1))为第i+1个传感器节点的位置坐标。
2.根据权利要求1所述的基于二元量化数据的声源目标迭代定位方法,其特征在于,将传感器节点接收到声源目标信号强度数据量化压缩为一位二进制数据。
3.根据权利要求1或2所述的基于二元量化数据的声源目标迭代定位方法,其特征在于,设在第i次迭代中,汇聚节点得到的加权样本集为{θi (m),wi (m)|m=1,…,M},利用高斯混合核平滑方法,第i次迭代中目标未知参数θ的概率分布表示为高斯混合分布:
4.根据权利要求1或2所述的基于二元量化数据的声源目标迭代定位方法,其特征在于,在目标定位前,首先建立声音信号强度随目标距离衰减的数学模型:
其中:ai代表第i个声传感器节点接收到的声源信号强度,Gi是第i个传感器节点的增益,P′0为在参考距离d0处声源的信号强度,di为目标和第i个传感器节点的距离,(sxi,syi)和(x,y)分别为第i个传感器节点和目标的位置坐标,n=2为声源信号强度衰减指数。
5.根据权利要求1或2所述的一种基于二元量化数据的声源目标迭代定位方法,其特征在于,汇聚节点对每个节点的量化数据进行序贯处理,在每次迭代中,汇聚节点只处理一个传感器节点的二元量化数据,其他传感器节点处于休眠状态。
6.一种基于二元量化数据的声源目标迭代定位方法,其特征在于,所述方法具体包括下列步骤:
步骤A,建立声音信号强度随目标距离衰减的数学模型:其中ai是第i个声传感器节点接收到的声源信号幅度,Gi是第i个传感器节点的增益,P′0为在参考距离d0处声源的强度,di为目标和第i个传感器节点的距离,(sxi,syi)和(x,y)分别为第i个传感器节点和声源目标的坐标,n=2为信号衰减指数;
步骤B,各传感器节点对实际接受到的目标信号强度进行量化压缩处理,只传递一个二进制数字信号到汇聚节点;每个传感器节点的观测数据与一个阈值η进行比较,如果高过阈值,传感器发送信号1到融合中心,否则传感器保持静默,不发送信号;在每次迭代中,汇聚节点序贯激活每个传感器节点,为了节约能量,其余节点处于休眠状态;
步骤C,基于高斯混合核平滑方法,利用蒙特卡洛方法从核密度的高斯分量中分别采样新的声源目标参数样本
步骤D,汇聚节点利用接收到的传感器节点二元量化观测数据,计算观测数据的似然函数作为样本的权值并归一化;
步骤E,汇聚节点利用样本权值对所得的样本集合进行重采样;
步骤F,汇聚节点利用所得到的加权样本集合,以序贯方式迭代计算目标的位置估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310122120.4A CN103237345B (zh) | 2013-04-09 | 2013-04-09 | 一种基于二元量化数据的声源目标迭代定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310122120.4A CN103237345B (zh) | 2013-04-09 | 2013-04-09 | 一种基于二元量化数据的声源目标迭代定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103237345A CN103237345A (zh) | 2013-08-07 |
CN103237345B true CN103237345B (zh) | 2015-12-02 |
Family
ID=48885352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310122120.4A Expired - Fee Related CN103237345B (zh) | 2013-04-09 | 2013-04-09 | 一种基于二元量化数据的声源目标迭代定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103237345B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104994579B (zh) * | 2015-06-30 | 2018-08-17 | 浙江工业大学 | 一种适用于事件触发的分布式移动目标跟踪方法 |
CN107064878B (zh) * | 2017-06-28 | 2019-08-20 | 山东大学 | 一种基于高精度gps的声源定位方法及其实现系统 |
CN110376290B (zh) * | 2019-07-19 | 2020-08-04 | 中南大学 | 基于多维核密度估计的声发射源定位方法 |
CN110728000B (zh) * | 2019-09-29 | 2021-06-29 | 华中科技大学 | 一种无砟轨道弹性扣件损伤识别方法 |
CN113055815B (zh) * | 2021-03-10 | 2021-12-07 | 电子科技大学 | 基于最大似然估计的多传感器1比特直接定位方法 |
CN114900817B (zh) * | 2022-07-12 | 2022-09-30 | 香港中文大学(深圳) | 一种基于二元量化的分布式无线传感网络信息估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102186241A (zh) * | 2011-04-26 | 2011-09-14 | 江苏科技大学 | 基于并行分布式粒子滤波的无线传感器网络目标跟踪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5602231B2 (ja) * | 2009-08-20 | 2014-10-08 | パーデュー・リサーチ・ファウンデーション | イベント駆動ワイヤレス・センサ・ネットワークのための予測的デューティ・サイクル順応方法 |
-
2013
- 2013-04-09 CN CN201310122120.4A patent/CN103237345B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102186241A (zh) * | 2011-04-26 | 2011-09-14 | 江苏科技大学 | 基于并行分布式粒子滤波的无线传感器网络目标跟踪方法 |
Non-Patent Citations (2)
Title |
---|
WSN中基于粒子滤波的目标跟踪研究;匡兴红等;《传感器技术学报》;20090731;第1节至第2节 * |
无线多跳传感器网络下基于粒子滤波的信道容错的信道跟踪方法;杨小军,刑科义;《自动化学报》;20110430;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103237345A (zh) | 2013-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103237345B (zh) | 一种基于二元量化数据的声源目标迭代定位方法 | |
CN105044662B (zh) | 一种基于wifi信号强度的指纹聚类多点联合室内定位方法 | |
CN103152820B (zh) | 一种无线传感器网络声源目标迭代定位方法 | |
CN104066179B (zh) | 一种改进的自适应迭代ukf的wsn节点定位方法 | |
CN105607039B (zh) | 非视距环境下基于到达时间差的稳健最小二乘定位方法 | |
CN107247259B (zh) | 基于神经网络的k分布海杂波形状参数估计方法 | |
Shi | A new weighted centroid localization algorithm based on RSSI | |
CN102890265B (zh) | 一种基于水下声传感器网络的目标被动定位方法 | |
CN107453993B (zh) | 基于Fisher信息矩阵与SUKF的目标跟踪方法 | |
CN105334495B (zh) | 一种无线网络中基于信号到达时间的非视距稳健定位方法 | |
CN108051779A (zh) | 一种面向tdoa的定位节点优选方法 | |
Wang et al. | Research on APIT and Monte Carlo method of localization algorithm for wireless sensor networks | |
CN111157943B (zh) | 在异步网络中基于toa的传感器位置误差抑制方法 | |
CN104902567A (zh) | 基于最大似然估计的质心定位方法 | |
CN112182961B (zh) | 一种换流站无线网络信道大尺度衰落建模预测方法 | |
CN106353722A (zh) | 基于代价参考粒子滤波的rssi测距方法 | |
CN102186241A (zh) | 基于并行分布式粒子滤波的无线传感器网络目标跟踪方法 | |
CN104023390A (zh) | 基于pso和ukf组合的wsn节点定位方法 | |
CN113543026A (zh) | 一种基于径向基函数网络的多楼层室内定位方法 | |
CN110536410B (zh) | 非视距环境下基于rss和tdoa测量的定位方法 | |
CN104113911A (zh) | 基于mle和ukf组合的wsn节点定位方法 | |
CN114915931A (zh) | 一种基于rssi技术的室内定位方法 | |
CN108848447A (zh) | 一种采用未知节点修正的差分DV_Distance节点定位方法 | |
Cheng et al. | An improved recurrent neural network for radio propagation loss prediction | |
CN104125640A (zh) | 一种基于无线传感网加权阈值质心定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151202 Termination date: 20160409 |