主权项 |
基于最大共信息熵的定量小波基选择方法,其特征在于:包括以下步骤:(1)根据采集到的信号,得到信号序列x(i),i=1,2,...,N,i表示等时间间隔采样时间点,N为信号长度;(2)采用小波基库中预选出的各小波基作为小波包变换的小波基函数,按照公式(a)对采集到的信号x(i)进行4层小波包变换;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>α</mi><mrow><mi>j</mi><mo>+</mo><mn>1,2</mn><mi>k</mi></mrow></msub><mo>=</mo><munder><mi>Σ</mi><mi>m</mi></munder><mi>h</mi><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>2</mn><mi>n</mi><mo>)</mo></mrow><msub><mi>α</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow>]]></math><img file="FDA0000509835430000011.GIF" wi="516" he="109" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>α</mi><mrow><mi>j</mi><mo>+</mo><mn>1,2</mn><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><munder><mi>Σ</mi><mi>m</mi></munder><mi>g</mi><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>2</mn><mi>n</mi><mo>)</mo></mrow><mi></mi><msub><mi>α</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000509835430000012.GIF" wi="656" he="90" /></maths>获得各个频带的小波包系数α<sub>j,k,m</sub>,其中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>j</mi><mo>=</mo><mn>0,1,2,3</mn><mo>,</mo><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msup><mrow><mo>,</mo><mn>2</mn></mrow><mi>j</mi></msup><mo>-</mo><mn>1</mn><mo>,</mo><mi>m</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msup><mn>2</mn><mrow><msub><mi>n</mi><mn>0</mn></msub><mo>-</mo><mi>j</mi></mrow></msup><mo>-</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0000509835430000017.GIF" wi="925" he="78" /></maths>n<sub>0</sub>=log<sub>2</sub>N≥J,J为最大分解层数,n为滤波器阶数,h和g分别为小波包分解的低通滤波器和高通滤波器组系数,且满足公式g(n)=(‑1)<sup>n</sup>h(1‑n),α<sub>j,k</sub>={α<sub>j,k,1</sub>,α<sub>j,k,2</sub>,...,α<sub>j,k,m</sub>}表示小波包分解中第j层、第k个子空间的小波包系数,m为该子空间小波包系数的个数;(3)根据所述步骤(2)获得的各个频带的小波包系数,按照公式(b)对小波包节点[4,0]的小波包系数进行重构;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>y</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><munder><mi>Σ</mi><mi>n</mi></munder><mrow><mo>(</mo><msub><mi>h</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>2</mn><mi>n</mi><mo>)</mo></mrow><msub><mi>α</mi><mrow><mi>j</mi><mo>+</mo><mn>1,2</mn><mi>k</mi></mrow></msub><mo>+</mo><msub><mi>g</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>2</mn><mi>n</mi><mo>)</mo></mrow><msub><mi>α</mi><mrow><mi>j</mi><mo>+</mo><mn>1,2</mn><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000509835430000013.GIF" wi="1086" he="132" /></maths>得到重构信号y,其中,h<sub>1</sub>和g<sub>1</sub>分别为小波包重构的低通滤波器和高通滤波器组系数。(4)把原始信号x看作数据序列X,把数据序列X经4层小波包变换得到的小波包节点[4,0]的重构信号y看作数据序列Y,根据公式(c)(d)分别计算数据序列X和数据序列Y的香农熵H(X)和H(Y),根据公式(e)计算两个数据序列X、Y的联合熵H(X,Y);<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munder><mi>Σ</mi><mrow><mi>x</mi><mo>∈</mo><mi>X</mi></mrow></munder><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>log</mi><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000509835430000014.GIF" wi="850" he="116" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munder><mi>Σ</mi><mrow><mi>x</mi><mo>∈</mo><mi>Y</mi></mrow></munder><mi>p</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mi>log</mi><mi>p</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000509835430000015.GIF" wi="846" he="115" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munder><mi>Σ</mi><mrow><mi>x</mi><mo>∈</mo><mi>X</mi></mrow></munder><munder><mi>Σ</mi><mrow><mi>y</mi><mo>∈</mo><mi>Y</mi></mrow></munder><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>log</mi><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000509835430000016.GIF" wi="871" he="121" /></maths>其中,p(x)和p(y)则分别代表数据序列X和数据序列Y的概率密度,p(x,y)表示数据序列X和数据序列Y的联合概率密度。(5)按照公式(f)分别计算采用各小波基函数对信号进行分解时两个数据序列的共信息熵值;I(X;Y)=‑H(X,Y)+H(X)+H(Y) (f)(6)选择共信息熵值最大时所采用的小波基作为最适合分析该类信号的小波基。 |