发明名称 一种基于模糊自助理论的乏信息多传感器数据融合估计方法
摘要 本发明提供一种基于模糊自助理论的乏信息多传感器数据融合估计方法,其具体步骤为:(1)乏信息多传感器测量数据自助再抽样;(2)乏信息多传感器数据融合;(3)求解数据模糊隶属函数离散值;(4)乏信息多传感器测量数据均值估计;(5)求解数据模糊拟合函数;(6)乏信息多传感器测量数据区间估计。本发明能实现测量数据样本量小、分布规律不明确等乏信息特征的多传感器数据融合估计,获得最代表原始信息的融合数据,实现乏信息多传感器测量数据的均值与区间估计。特点是:(1)对测量数据所服从的概率分布没有要求;(2)允许趋势项的变化规律未知;(3)允许测量数据样本量小。
申请公布号 CN103020715B 申请公布日期 2016.06.29
申请号 CN201210498844.4 申请日期 2012.11.29
申请人 北京航空航天大学 发明人 王中宇;王倩;李强;王岩庆
分类号 G06N7/02(2006.01)I 主分类号 G06N7/02(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 杨学明
主权项 一种基于模糊自助理论的乏信息多传感器数据融合估计方法,其特征在于,该方法包括如下步骤:步骤(1)乏信息多传感器测量数据自助再抽样;步骤(1)中所述的乏信息多传感器测量数据自助再抽样,具体实现过程如下:在测量过程中,m个传感器获得的时间序列数据为原始时间序列,构成的矩阵为:Y={y<sub>q</sub>(k)};k=1,2,…n;q=1,2,…m   (1)式中:y<sub>q</sub>(k)为第q个传感器获得的第k个数据;k为时间序号;n为数据个数;在时刻k,m个传感器获得的数据是矩阵Y的第k列列向量:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>Y</mi><mi>k</mi><mi>T</mi></msubsup><mo>=</mo><mo>{</mo><msub><mi>y</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>y</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>y</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971664520000011.GIF" wi="1315" he="70" /></maths>将式(2)看作初始样本,根据自助法,从<img file="FDA0000971664520000012.GIF" wi="49" he="54" />中等概率地随机抽取m次,每次抽取一个数据,就得到第一个自助样本,它有m个数据,抽取数据是可以放回的,连续地重复B次,就得到B个自助仿真样本,用矩阵Y<sub>kb</sub>表示为:Y<sub>kb</sub>={y<sub>k,b</sub>(u)};u=1,2,…m;b=1,2,…B;k=1,2,…n   (3)式中:y<sub>k,b</sub>(u)为时刻k的第b个自助样本中的第u个数据;步骤(2)乏信息多传感器数据融合;步骤(2)中所述的乏信息多传感器数据融合,具体实现过程如下:求解矩阵Y<sub>kb</sub>中各列的均值序列,用向量表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>Y</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mi>b</mi></mrow></msub><mo>=</mo><mo>{</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>u</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>y</mi><mrow><mi>k</mi><mo>,</mo><mi>b</mi></mrow></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>}</mo><mo>;</mo><mi>b</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>B</mi><mo>;</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971664520000013.GIF" wi="1500" he="135" /></maths>在时刻k,将<img file="FDA0000971664520000014.GIF" wi="77" he="68" />中的B个数据由小到大排序,并按一定的间隔分成T组,建立直方图,其中第t组的频率P<sub>t</sub>可表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>t</mi></msub><mo>=</mo><mfrac><msub><mi>n</mi><mi>t</mi></msub><mi>B</mi></mfrac><mo>;</mo><mi>t</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>T</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971664520000015.GIF" wi="1022" he="111" /></maths>式中:n<sub>t</sub>为第t组的数据个数;用连续变量x代替式(2)中的离散值y<sub>q</sub>(k),根据最大熵原理,可得自助分布概率密度函数p(x)的表达式:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>m</mi><mi>A</mi></msub></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971664520000021.GIF" wi="1093" he="158" /></maths>式中:λ<sub>i</sub>为第i个拉格朗日乘子i=0,1,…m<sub>A</sub>,由数值算法得出;m<sub>A</sub>为原点矩的最高阶数,取m<sub>A</sub>=5~8;在时刻k,m个传感器数据的融合值用数学期望表示:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>y</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Integral;</mo><mi>R</mi></munder><mi>x</mi><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>d</mi><mi>x</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971664520000022.GIF" wi="1094" he="119" /></maths>式中:R为积分空间,j=1,2,...m;式(7)可以用离散型式表示为加权均值:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>y</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>Y</mi><mi>t</mi></msub><msub><mi>P</mi><mi>t</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000971664520000023.GIF" wi="1086" he="136" /></maths>式中:Y<sub>t</sub>为第t组的组中值;对于m个传感器数据,融合值y<sub>j</sub>(k)构成一个时间序列,称为自助融合序列Y<sub>j</sub>,可表示为:Y<sub>j</sub>={y<sub>j</sub>(k)};k=1,2,…n   (9)步骤(3)求解数据模糊隶属函数离散值;其中,步骤(3)所述的求解数据模糊隶属函数离散值,具体实现过程如下:利用线性排序估计法得到隶属函数的离散值,线性估计排序法将步骤(2)形成的新序列Y<sub>j</sub>先由小到大排序,然后逐项相减,差值形成新的序列,利用新形成的差值序列,可构造线性隶属函数,线性隶属函数的值即为步骤(3)所需求解的模糊隶属函数离散值;步骤(4)乏信息多传感器测量数据均值估计;步骤(4)中所述的乏信息多传感器测量数据均值估计,具体实现过程如下:在步骤(3)求解得到的隶属函数的基础上,根据隶属最大原则估计多传感器测量数据均值,取隶属函数值为1时的自变量值为多传感器测量数据估计均值;步骤(5)求解数据模糊拟合函数;步骤(5)中所述的求解数据模糊拟合函数,具体实现过程如下:采取三阶多项式的拟合方式得到数据模糊拟合函数,所述的三阶多项式拟合是采样数据中与步骤(3)所需求解的隶属函数离散值之间残差无穷范数最小方法确定多项式系数;步骤(6)乏信息多传感器测量数据区间估计;步骤(6)中所述的乏信息多传感器测量数据区间估计,具体实现过程如下:根据模糊集合理论意义上的最优水平确定相应的水平截集λ,从而得到隶属区间,实现乏信息多传感器测量数据区间估计;该基于模糊自助理论的乏信息多传感器数据融合估计方法解决了乏信息数据融合估计问题,并且计算简便;该基于模糊自助理论的乏信息多传感器数据融合估计方法能实现测量数据样本量小、分布规律不明确乏信息特征的多传感器数据融合估计,获得最代表原始信息的融合数据,实现乏信息多传感器测量数据均值与区间估计;该基于模糊自助理论的乏信息多传感器数据融合估计方法对概率分布无任何要求,允许趋势项的变化规律未知,也允许测量数据样本量小。
地址 100191 北京市海淀区学院路37号