《基于传感器网络的复杂环境下多目标信息融合方法.pdf》由会员分享,可在线阅读,更多相关《基于传感器网络的复杂环境下多目标信息融合方法.pdf(11页完整版)》请在专利查询网上搜索。
1、(10)申请公布号 CN 102509020 A (43)申请公布日 2012.06.20 C N 1 0 2 5 0 9 0 2 0 A *CN102509020A* (21)申请号 201110361602.6 (22)申请日 2011.11.15 G06F 19/00(2011.01) (71)申请人哈尔滨工程大学 地址 150001 黑龙江省哈尔滨市南岗区南通 大街145号哈尔滨工程大学科技处知 识产权办公室 (72)发明人林云 李靖超 李一兵 叶方 康健 葛娟 田雪宜 (54) 发明名称 基于传感器网络的复杂环境下多目标信息融 合方法 (57) 摘要 本发明的目的在于提供基于传感器网。
2、络的复 杂环境下多目标信息融合方法,分为以下步骤: 选用适当的小波基函数将量测数据分解在若干个 尺度上,在每个尺度上的低频子空间信号利用卡 尔曼算法进行滤波,得到新的滤波数据,将每个尺 度上的高频子空间中细节信号的极大值点去掉, 得到新的滤波数据,将每个尺度上处理后的数据 通过小波重构算法,得到原始量测数据在不同尺 度上滤波后的结果。本发明具有高灵敏度、高精 度、高准确率、高稳定性和高可靠性等特点。 (51)Int.Cl. 权利要求书1页 说明书6页 附图3页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 1 页 说明书 6 页 附图 3 页 1/1页 2 1.基于传。
3、感器网络的复杂环境下多目标信息融合方法,其特征是: (1)选用小波基函数将量测数据分解在尺度上; (2)在每个尺度上的低频子空间信号利用卡尔曼算法进行滤波: 离散系统在尺度i下的状态方程和量测方程为: 在尺度i下,X(i,k)为k时刻的n维状态矢量,也是被估计矢量;Z(i,k)为k时刻的 m维量测矢量;(i,k+1/k)为k到k+1时刻的一步转移矩阵(nn阶);w(i,k)为k时刻 的系统噪声;H(i,k)为k时刻的量测矩阵(mn阶);G(i,k)是系统噪声的加权噪声;v(i, k)为k时刻的m维量测噪声,w(i,k)和v(i,k)满足如下条件: 其中Q(i,k)和R(i,k)分别为系统噪声和。
4、量测噪声的方差矩阵,它们分别是已知非负 阵和正定阵; kj 是Kronecker函数, 卡尔曼滤波滤波包括时间更新和量测更新: 时间更新即预测估计: P(i,k/k-1)(i,k/k-1)P(i,k-1/k-1) T (i,k/k-1)+ G(i,k-1)Q(i,k-1)G T (i,k-1) 量测更新即滤波估计: K(i,k)P(i,k/k-1)H T (i,k)H(i,k)P(i,k/k-1)H T (i,k)+R(i,k) -1 P(i,k/k)(I-K(i,k)H(i,k)P(i,k/k-1) 式中K(i,k)为滤波增益阵; (3)将每个尺度上的高频子空间中细节信号的极大值点去掉,得到。
5、新的滤波数据; (4)将每个尺度上处理后的数据通过小波重构方法,得到原始量测数据在不同尺度上 滤波后的结果: 所述的小波重构方法为:通过一个脉冲响应为h(l)的低通滤波器可以从尺度i上获得 尺度i-1上的低频子空间信号x L (i-1,k),通过一个脉冲响应为g(l)的高通滤波器可以获 得尺度i-1上的高频子空间信号x H (i-1,k): x L (i-1,k) l h(l)x(i,2k-l) x H (i-1,k) l g(l)x(i,2k-l)。 权 利 要 求 书CN 102509020 A 1/6页 3 基于传感器网络的复杂环境下多目标信息融合方法 技术领域 0001 本发明涉及的是。
6、一种目标跟踪与识别领域的信息融合方法。 背景技术 0002 目标的跟踪与识别是传感器网络的主要用途之一,也是一个难点和关键问题,在 许多领域例如交通监控、机构安全和战场状况获取等方面具有广阔的应用前景。 0003 在信息时代,传感器网络所提供的信息量很大,电磁环境十分复杂,这就要求信息 融合装置必须具有高灵敏度(作用距离远)、高准确率、高稳定性和高可靠性。一个传感器 网络往往可以获得关于目标的多个特征信息,如目标的中心频率、脉宽、脉冲重复周期、脉 冲到达时间等,而一般的信息融合装置只利用了很少的一部分信息,从而造成了传感器信 息的大量浪费。另外,由于信息的利用率不足,判决依据单一,往往容易造成。
7、误判和错判,会 对整个跟踪和识别过程产生重大影响。因而综合利用传感器网络提供的多个特征信息是十 分必要的,这样既不浪费信息,又能保证信息融合装置的高灵敏度、高准确率、高稳定性和 高可靠性。 发明内容 0004 本发明的目的在于提供一种高灵敏度、高精度、高准确率、高稳定性和高可靠性的 基于传感器网络的复杂环境下多目标信息融合方法。 0005 本发明的目的是这样实现的: 0006 本发明基于传感器网络的复杂环境下多目标信息融合方法,其特征是: 0007 (1)选用小波基函数将量测数据分解在尺度上; 0008 (2)在每个尺度上的低频子空间信号利用卡尔曼算法进行滤波: 0009 离散系统在尺度i下的。
8、状态方程和量测方程为: 0010 0011 在尺度i下,X(i,k)为k时刻的n维状态矢量,也是被估计矢量;Z(i,k)为k时 刻的m维量测矢量;(i,k+1/k)为k到k+1时刻的一步转移矩阵(nn阶);w(i,k)为k 时刻的系统噪声;H(i,k)为k时刻的量测矩阵(mn阶);G(i,k)是系统噪声的加权噪声; v(i,k)为k时刻的m维量测噪声,w(i,k)和v(i,k)满足如下条件: 0012 0013 其中Q(i,k)和R(i,k)分别为系统噪声和量测噪声的方差矩阵,它们分别是已知 非负阵和正定阵; kj 是Kronecker函数, 0014 卡尔曼滤波滤波包括时间更新和量测更新: 。
9、说 明 书CN 102509020 A 2/6页 4 0015 时间更新即预测估计: 0016 0017 P(i,k/k-1)(i,k/k-1)P(i,k-1/k-1) T (i,k/k-1)+ 0018 G(i,k-1)Q(i,k-1)G T (i,k-1) 0019 量测更新即滤波估计: 0020 K(i,k)P(i,k/k-1)H T (i,k)H(i,k)P(i,k/k-1)H T (i,k)+R(i,k) -1 0021 0022 P(i,k/k)I-K(i,k)H(i,k)P(i,k/k-1) 0023 式中K(i,k)为滤波增益阵; 0024 (3)将每个尺度上的高频子空间中细节。
10、信号的极大值点去掉,得到新的滤波数 据; 0025 (4)将每个尺度上处理后的数据通过小波重构方法,得到原始量测数据在不同尺 度上滤波后的结果: 0026 所述的小波重构方法为:通过一个脉冲响应为h(l)的低通滤波器可以从尺度i上 获得尺度i-1上的低频子空间信号x L (i-1,k),通过一个脉冲响应为g(l)的高通滤波器可 以获得尺度i-1上的高频子空间信号x H (i-1,k): 0027 x L (i-1,k) l h(l)x(i,2k-l) 0028 x H (i-1,k) l g(l)x(i,2k-l)。 0029 本发明的优势在于:本发明具有高灵敏度、高精度、高准确率、高稳定性和。
11、高可靠 性等特点。 附图说明 0030 图1是本发明的多目标信息融合装置实现框图; 0031 图2是本发明流程图; 0032 图3是数据关联流程图; 0033 图4是融合决策流程图。 具体实施方式 0034 下面结合附图举例对本发明做更详细地描述: 0035 结合图14,本发明的目的是这样实现的:信息融合装置首先将传感器网络传送 过来的多个目标的特征信息和状态信息进行多尺度维上的数据预处理,将数据中的野值和 噪声去除。紧接着采用灰色关联算法对多个尺度维上的数据进行关联,将最终的关联结果 与信息融合装置中的目标识别数据库进行对比,利用灰色关联算法得出每个传感器所提供 的目标置信度,获得传感器网络。
12、报告,最后利用DS证据合成理论对传感器网络报告进行时 空域的二级判决,得出最终的综合判决结果。 0036 本发明是一个具有高灵敏度、高精度、高准确率、高稳定性和高可靠性的信息融合 装置,利用传感器提供的大量特征信息进行多尺度维上的数据预处理和数据关联,将所得 结果利用目标识别数据库形成传感器网络报告,进行基于DS证据合成理论的时空域二级 判决,最终得到精确的判决结果。 说 明 书CN 102509020 A 3/6页 5 0037 多目标信息融合装置由传感器网络(1)、串口芯片(2)、可编程逻辑器件FPGA(3)、 数字信号处理器DSP(4)和存储芯片SDRAM(5)等组成。 0038 当传感。
13、器网络(1)搜索到目标后,将获得的目标的多个特征信息及状态信息作为 数据0传送给串口芯片(2),串口芯片(2)将数据0量化后变成数据1传送给可编程逻辑芯 片FPGA(3)。可编程逻辑器件FPGA(3)对数据1进行多尺度维上的数据预处理和数据关联, 这里多尺度维上的数据预处理是指对数据1进行多尺度维上的滤波和野值剔除,多尺度维 是指利用小波分解算法将量测数据分解到多个尺度上。可编程逻辑器件FPGA(3)将处理后 的数据2传送给数字信号处理器DSP(4),数字信号处理器DSP(4)将数据2与目标识别数据 库对比,按照灰关联算法得出目标的置信度,并形成传感器网络报告,采用基于DS证据合 成理论的时空。
14、域二级判决模型,得到最终判决结果。数字信号处理器DSP(4)将最终的判决 结果作为数据3传送给存储器SDRAM(5)。 0039 图2是多尺度维数据预处理流程图。 0040 在目标跟踪与识别过程中,会不可避免的引入噪声和野值,导致后面的状态更新 才生错误,无法准确的跟踪和识别目标。因此需要首先在可编程逻辑器件FPGA(3)内对数 据1进行数据预处理,即数据滤波和野值剔除,提高量测数据的信噪比,本发明采用多尺度 维上的卡尔曼滤波算法。 0041 算法的基本思路: 0042 1根据实际情况,选用适当的小波基函数将量测数据分解在若干个尺度上。 0043 2在每个尺度上的低频子空间信号利用卡尔曼算法进。
15、行滤波,得到新的滤波数据。 0044 3将每个尺度上的高频子空间中细节信号的极大值点去掉,得到新的滤波数据。 0045 4将每个尺度上处理后的数据通过小波重构算法,得到原始量测数据在不同尺度 上滤波后的结果。 0046 卡尔曼滤波算法的基本原理: 0047 已知离散系统在尺度i下的状态方程和量测方程如下所示: 0048 0049 在尺度i下,X(i,k)为k时刻的n维状态矢量,也是被估计矢量;Z(i,k)为k时 刻的m维量测矢量;(i,k+1/k)为k到k+1时刻的一步转移矩阵(nn阶);w(i,k)为k 时刻的系统噪声;H(i,k)为k时刻的量测矩阵(mn阶);G(i,k)是系统噪声的加权噪。
16、声; v(i,k)为k时刻的m维量测噪声。w(i,k)和v(i,k)需要满足如下条件: 0050 0051 其中Q(i,k)和R(i,k)分别为系统噪声和量测噪声的方差矩阵,它们分别是已知 非负阵和正定阵; kj 是Kronecker函数。 0052 卡尔曼滤波滤波原理包括时间更新和量测更新。 0053 时间更新(即预测估计): 0054 说 明 书CN 102509020 A 4/6页 6 0055 P(i,k/k-1)(i,k/k-1)P(i,k-1/k-1) T (i,k/k-1)+ 0056 G(i,k-1)Q(i,k-1)G T (i,k-1) 0057 量测更新(即滤波估计): 0。
17、058 K(i,k)P(i,k/k-1)H T (i,k)H(i,k)P(i,k/k-1)H T (i,k)+R(i,k) -1 0059 0060 P(i,k/k)I-K(i,k)H(i,k)P(i,k/k-1) 0061 式中K(i,k)为滤波增益阵。 0062 小波算法的基本原理: 0063 由小波理论可知,通过一个脉冲响应为h(l)的低通滤波器可以从尺度i上获得尺 度i-1上的低频子空间信号x L (i-1,k),通过一个脉冲响应为g(l)的高通滤波器可以获得 尺度i-1上的高频子空间信号x H (i-1,k): 0064 x L (i-1,k) l h(l)x(i,2k-l) 006。
18、5 x H (i-1,k) l g(l)x(i,2k-l) 0066 图3是数据关联流程图。 0067 数据关联是指判别来自不同传感器的量测或相同传感器的量测数据是否代表同 一个目标的过程。本发明采用熵权灰色数据关联算法,关联概率的获取方法如下: 0068 传感器能够获取目标的多个特征信息,如状态信息(S)、中心频率(f)、脉冲重复 周期(T)、脉宽()、脉内特征(C)等,设:比较数列x i x i (j)|j1,2,M,11, 2m(k):表示k时刻有m(k)个量测数列。参考数列X 0 X 0 (j)|j1,2,M:为参考目 标的特征信息。 0069 1、计算指标绝对差 i (j) 0070。
19、 指标绝对差 i (j)X 0 (j)-X i (j):表示X 0 与X i 第j个指标的绝对差。 0071 2、计算关联系数 i (j) 0072 0073 其中,为所有量测所有指标差中的最小值,为所有量 测所有指标差中的最大值。称为分辨系数,取值范围为0,1,通常取0.5。 0074 3、计算熵权 0075 各个特征的权重反应了各个特征在整个观测过程中的重要程度。一般权重由专家 主观经验给出,没有科学可靠性。熵权可以自适应的根据特征信息给出特征的权重。其计 算过程如下: 0076 (1)确定好初始数据矩阵 0077 ( ij ) mn i1 i2 , in 0078 其中,i1,2,m;j。
20、1,2,n, ij 为第i个量测第j个特征的特征值。 0079 (2)计算第j项指标下第i单元的指标值比重p ij 0080 0081 (3)计算第j个特征值的输出熵 说 明 书CN 102509020 A 5/6页 7 0082 每个判别特征的熵值为: 0083 0084 用熵值的最大值为E max lnm对E j 进行归一化,得到表征第j个特征相对重要性 确定程度的熵值: 0085 e j E j /E max 0086 (4)计算第j个特征的变异度D j 0087 D j 1-e j (j1,2,n) 0088 (5)计算第j个特征的权重 j 0089 0090 4、计算灰关联度 i 0。
21、091 0092 图4是融合决策流程图。 0093 由灰关联度形成传感器报告的过程如下: 0094 定义目标识别数据库中所有目标的集合为传感器报告中各证据 形成如下: 0095 0096 其中m(R)为不确定性目标的基本概率赋值函数, i 为第i个量测与目标识别框 架的灰关联度。 0097 DS证据合成理论的时空域判决方法如下: 0098 1、单个传感器的时域判决 0099 根据单传感器报告,对于第s个传感器,s1,2,P,依据Q个测量周期的累积量 测,可得第i个命题A i (i1,2,N)的单个传感器时域判决结果为: 0100 0101 其中,z1,Q。 0102 不确定目标的时域判决结果为。
22、: 0103 说 明 书CN 102509020 A 6/6页 8 0104 2、多个传感器的空域融合判决 0105 对所有传感器的时域判决结果进行空域上的融合。第i个命题A i (i1,2,N) 的多传感器的空域判决结果为: 0106 0107 不确定目标的空域判决结果为: 0108 0109 3、判决原理 0110 判决规则:设有两焦元若满足: 0111 0112 则A 1 为判决结果,其中 1 , 2 为预先设定门限,为不确定集合,并且A 1 、A 2 满 足: 说 明 书CN 102509020 A 1/3页 9 图1 说 明 书 附 图CN 102509020 A 2/3页 10 图2 图3 说 明 书 附 图CN 102509020 A 10 3/3页 11 图4 说 明 书 附 图CN 102509020 A 11 。