发明名称 基于脉博熵的多点脉搏信息融合方法
摘要 基于脉搏熵的多点脉搏信息融合方法,属医学信号处理领域。本发明由多路传感器模块拾取人体多点脉搏信号,经放大、滤波和采样的预处理得多路脉搏信号S<SUB>ij</SUB>,i=1,2…M,j=1,2…N,M为信号路数,还包括以下步骤:用脉搏熵的思想,通过对S<SUB>ij</SUB>进行脉搏熵融合处理,提出近一化规则和选取最大熵比的输出优化的脉搏波形数据供分析使用。由于采用脉搏熵的思想,使得多路脉搏信号中最大的特征信息量会被提取和融合进一路脉搏信号,从而减少了存储和传输的数据量,另由于在脉搏熵融合过程中提出的近一化规则和选取最大熵比的方法,使得处理后的脉搏信号保持了信号的真实值的同时,又突显脉搏波信号特征信息,便于后续处理中得到脉搏信号的重要特征参数,可应用于无损伤测量和病理诊断。
申请公布号 CN1631317A 申请公布日期 2005.06.29
申请号 CN200510002277.9 申请日期 2005.01.20
申请人 北京工业大学 发明人 阮晓钢;乔俊飞;王家庆
分类号 A61B5/021;G06F17/00;G06F19/00;//G06F159∶00 主分类号 A61B5/021
代理机构 北京思海天达知识产权代理有限公司 代理人 张慧
主权项 1、一种基于脉搏熵的多点脉搏信息融合方法,首先由多路传感器模块拾取人体多点脉搏信号,经过放大、滤波和采样的预处理后,得到多路脉搏信号S<sub>ij</sub>,其中S<sub>ij</sub>是由压力传感器把脉搏压力信号转变为电压信号,是各采样点对应的脉搏信号幅度,i=1,2…M,j=1,2…N,M为信号路数,N为每路信号的采样点数,然后多路脉搏信号S<sub>ij</sub>通过RS 232总线或USB接口等方式传入计算机。本发明的特征还包括以下步骤:1)、计算机读取经过上述预处理的多路脉搏信号S<sub>ij</sub>存入内存;2)、计算M路信号在同时刻各采样点的比重p<sub>ij</sub>:<maths num="001"><![CDATA[ <math><mrow><msub><mi>P</mi><mi>ij</mi></msub><mo>=</mo><mfrac><msub><mi>S</mi><mi>ij</mi></msub><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>S</mi><mi>ij</mi></msub></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>M</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>3)、计算M路信号在同时刻各采样点的脉搏熵E<sub>ij</sub>:<maths num="002"><![CDATA[ <math><mrow><msub><mi>E</mi><mi>ij</mi></msub><mo>=</mo><mo>-</mo><msub><mi>p</mi><mi>ij</mi></msub><mi>lo</mi><msubsup><mi>g</mi><mn>2</mn><msub><mi>p</mi><mi>ij</mi></msub></msubsup><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>M</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>4)、计算M路信号在同时刻采样点的均熵<img file="A2005100022770002C3.GIF" wi="94" he="77" /><maths num="003"><![CDATA[ <math><mrow><mover><msub><mi>E</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>E</mi><mi>ij</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>M</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>5)、运用近一化规则,计算M路信号各采样点的近一化熵比e<sub>ij</sub>:<maths num="004"><![CDATA[ <math><mrow><msub><mi>e</mi><mi>ij</mi></msub><mo>=</mo><mfrac><msub><mi>E</mi><mi>ij</mi></msub><mover><msub><mi>E</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover></mfrac><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>M</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>6)、取近一化熵比e<sub>ij</sub>中同时刻采样点的最大熵比e<sub>1j</sub>:e<sub>1j</sub>=MAX(e<sub>ij</sub>)(i=1,2…M,j=1,2…N)    ⑤其中在步骤2)到步骤6)的过程中,进行计算M路信号在同时刻采样点的均值<maths num="005"><![CDATA[ <math><mrow><mover><msub><mi>s</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>:</mo><mover><msub><mi>s</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>s</mi><mi>ij</mi></msub><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>M</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>7)、输出一路优化脉搏信号W<sub>o</sub>:<maths num="006"><![CDATA[ <math><mrow><msub><mi>W</mi><mi>o</mi></msub><mo>=</mo><mover><msub><mi>s</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><msub><mi>e</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>W<sub>o</sub>是近一化后的最大熵比融合脉搏信号数据,其保持与脉搏信号实际值相对应,又能提取多路脉搏信号在同时刻的最大特征信息,为后续处理提供较好的脉搏信息。8)、对融合后的脉搏信号W<sub>o</sub>进行现有技术的时频域分析处理,计算诊断需要的指标参数。
地址 100022北京市朝阳区平乐园100号
您可能感兴趣的专利