CN101900674A - 一种用于二氧化碳浓度检测信号静噪的小波滤波方法 - Google Patents
一种用于二氧化碳浓度检测信号静噪的小波滤波方法 Download PDFInfo
- Publication number
- CN101900674A CN101900674A CN 201010216322 CN201010216322A CN101900674A CN 101900674 A CN101900674 A CN 101900674A CN 201010216322 CN201010216322 CN 201010216322 CN 201010216322 A CN201010216322 A CN 201010216322A CN 101900674 A CN101900674 A CN 101900674A
- Authority
- CN
- China
- Prior art keywords
- wavelet
- signal
- carbon dioxide
- wavelet decomposition
- gradation
- 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.)
- Pending
Links
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明属生物医学工程领域,涉及一种用于二氧化碳浓度检测信号静噪的小波滤波方法,包括下列步骤:选择Mexican Hat作为小波基,并确定小波分解的层次为3,然后对二氧化碳监测仪检测到的二氧化碳浓度信号进行3层小波分解;设定某个固定阈值,对第1到第3层的每一高频系数,进行阈值量化处理;根据小波分解的第3层的低频系数和经过阈值量化处理后的第1层到第3层的高频系数,进行一维信号的小波重构。本发明能很好的实现信号中有用部分和噪声的有效分离,而不导致信号失真。
Description
技术领域
本发明属于生物医学工程技术领域,涉及一种用于二氧化碳浓度检测信号静噪的小波滤波方法。
背景技术
在医用临床上监测二氧化碳浓度以及呼吸末二氧化碳浓度值有着重要的意义,已成为重要的生命指标之一。这就使得如何准确的捕捉呼吸末潮时刻成为同时随着NDIR技术在气体测量应用中的成熟和普遍,主流式CO2浓度监测方法在医疗中被提出和应用。
气体分子对一定波长的红外线辐射的吸收服从Lambert-Beer吸收定律,所以二氧化碳可以吸收特定波长的红外辐射。如果气体进入检测室,并用传感器一侧的红外光源照射,而传感器的另一侧则可以测的红外辐射的衰减,这样,就可以通过测量信号的衰减得到二氧化碳浓度的变化。
由于呼吸气体是不同气体的混合,而每种入射光对特定气体的吸收系数是不同的,吸收光线和气体之间的关系可以用下式表示:
I=I0exp(-∑μiCiL)
其中I0:入射光强度;I:紧急光照强度;μ:气体吸收系数;C:气体浓度;L:气体厚度
常规呼吸气体分子的吸收系数是不同的。二氧化碳吸收峰值对应的光波长为4.26μm,如果只有近似4.26μm波长的光线,则气体对光的吸收刻印描述为:
根据这一原理可以实现在患者呼吸主通道内实时、连续的检测呼吸中二氧化碳的浓度信号。
一种典型的主流式二氧化碳测量装置如图1所示,红外光源1和测量端3分别位于主气体通道2的两侧。红外光源1发出的红外光经过4.26μm窄带滤波后的光线,被测量端的红外传感器所采集。
图1中的主流式二氧化碳浓度监测方法将CO2浓度监测模块直接放置在患者的呼吸气路导管中,通过传感器对通路中的CO2浓度直接进行测定。但是如果检测到的二氧化碳浓度信号存在各种干扰和失真,就不能很好的实现实时、连续的检测患者呼吸中二氧化碳浓度的意义,这就需要一种可行有效地滤波方法。
二氧化碳气体检测中,常常由于噪声的存在而不能从所记录的数据中提取有效的信息。二氧化碳光电传感器采集到的电信号属于微弱信号,很容易被噪声信号淹没。
光电传感器在光电转换过程中,不仅输出所需的信号电压或电流,同时还伴随着输出无用的随机性涨落或起伏量,即噪声。光电探测器的固有噪声主要表现为热噪声、散粒噪声、产生复合噪声、温度噪声及电流噪声.这些噪声污染了与二氧化碳有关的光电信号,给有用光电信号的检测判别及进一步的数据处理带来了不利因素,有必要对光电信号进行数字滤波,以获得真实的信号,提高检测灵敏度.
发明内容
本发明的目的是提供一种小波滤波方法,将二氧化碳监测仪检测到的二氧化碳浓度信号进行小波去噪处理,以实现埋于噪声中的光电信号的有效测量。本发明采用如下的技术方案:
一种用于二氧化碳浓度检测信号静噪的小波滤波方法,包括下列步骤:
第一步:信号的小波分解:选择Mexican Hat作为小波基,并确定小波分解的层次为3,然后对二氧化碳监测仪检测到的二氧化碳浓度信号进行3层小波分解;
第二步:小波分解高频系数的阈值化:设定某个固定阈值,对第1到第3层的每一高频系数,进行阈值量化处理;
第三步:一维小波的重构:根据小波分解的第3层的低频系数和经过阈值量化处理后的第1层到第3层的高频系数,进行一维信号的小波重构。
本发明针对CO2监测信号属于一种突发性的短过程并且是明显的非平稳性信号,能量主要分布在<10Hz的范围内的特点,选用“Mexican Hat”小波函数作为小波基进行小波分解,经过小波分解去噪后的信号在去除噪声之外还能将原始信号变化时的细节保留下来,而使用FFT进行去噪之后,得到的信号在细节方面有很大程度的损失,有些关键信号的突变在经过FFT消噪时候就丢失了。呼吸作用中CO2浓度值属于非平稳信号中由于夹带着瞬态的突变成分,使用本发明的小波分析的方法能很好的实现信号中有用部分和噪声的有效分离,而不导致信号失真。
附图说明
图1:一种典型的主流式二氧化碳浓度测量装置。
图2本发明采用的小波分解结构图。
图3本发明的消噪过程示意图。
具体实施方式
小波变换在时域、频域均有良好的局化性,在分析暂态和非平稳信号方面具有独特的优越性。小波变换的主要特征是在时域和频域上的双重局部性以及线性分解性,可以将信号线性分解为表征不同时间或频率的分量.小波多分辨率分析利用了小波变换对信号分频,各个频率成分在时间轴上的位置保持不变的性质.因此,根据光谱信号和噪声信号在频率上的区别,可以对其进行有效分离,从而实现从噪声中提取光电信号。
小波分析方法是窗口大小固定,但其形状可改变,时间窗和频率窗都可改变的时频局域化分析方法,即在低频部分具有较高的频率分辨率和较低的时间分辨率,在高频部分具有较高时间分辨率和较低频率分辨率,优于经典Fourier变换和短时Fourier的地方。
本发明采用“Mexican Hat”小波函数作为小波基进行小波分解,该小波函数具有简便快速,不需设置尺度等参数等优点,无论对高频、低频信号均有很好的除噪能力,是对CO2传感器端进行信号的处理的最优选择
Mexican Hat小波为:
由于Mexican Hat小波的中心频率f0和尺度a易于设定,由其构成的滤波器具有良好的窄带滤波能力,此外Mexican Hat具有线性位移的特性,因此可以利用Mexican Hat小波作为小波消噪的小波基。
呼吸作用过程中,CO2监测信号实际上是属于一种突发性的短过程并且是明显的非平稳性信号。对信号进行频谱分析可以知道能量主要分布在<10Hz的范围内。
设CO2信号表达为:
s(i)=f(i)+e(i),i=0,...,n-1
其中:
f(i)为真实信号
e(i)为噪声
s(i)为含噪声的信号。
本发明采用的小波分解结构如图2所示,首先对信号进行三层小波分解,分解过程如图2。噪声部分通常包含在cD1,cD2,cD3中,因而,可以通过门限阈值等形式对小波系数进行处理,然后对信号进行重构,可以达到消噪的目的。对信号s(i)消噪的目的就是要抑止信号中的噪声部分,从而在s(i)中恢复出真实信号f(i)。
信号的消噪过程可分为三个步骤进行,见图3:
第一步:信号的小波分解。选择Mexican Hat作为小波基,并确定小波分解的层次为3,然后对检测到的CO2信号进行3层小波分解。
第二步:小波分解高频系数的阈值化。对第1到第3层的每一高频系数,设定某个固定阈值进行量化处理。
第三步:一维小波的重构。根据小波分解的第3层的低频系数和经过量化处理后的第1层到第3层的高频系数,进行一维信号的小波重构。
Claims (1)
1.一种用于二氧化碳浓度检测信号静噪的小波滤波方法,包括下列步骤:
第一步:信号的小波分解:选择Mexican Hat作为小波基,并确定小波分解的层次为3,然后对二氧化碳监测仪检测到的二氧化碳浓度信号进行3层小波分解;
第二步:小波分解高频系数的阈值化:设定某个固定阈值,对第1到第3层的每一高频系数,进行阈值量化处理;
第三步:一维小波的重构:根据小波分解的第3层的低频系数和经过阈值量化处理后的第1层到第3层的高频系数,进行一维信号的小波重构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010216322 CN101900674A (zh) | 2010-07-02 | 2010-07-02 | 一种用于二氧化碳浓度检测信号静噪的小波滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010216322 CN101900674A (zh) | 2010-07-02 | 2010-07-02 | 一种用于二氧化碳浓度检测信号静噪的小波滤波方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101900674A true CN101900674A (zh) | 2010-12-01 |
Family
ID=43226393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010216322 Pending CN101900674A (zh) | 2010-07-02 | 2010-07-02 | 一种用于二氧化碳浓度检测信号静噪的小波滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101900674A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102818629A (zh) * | 2012-05-04 | 2012-12-12 | 浙江大学 | 基于平稳小波变换的微型光谱仪信号去噪方法 |
CN103513094A (zh) * | 2013-09-29 | 2014-01-15 | 天津理工大学 | 一种消除电力系统检测信号噪声的装置及方法 |
CN103674847A (zh) * | 2013-12-17 | 2014-03-26 | 重庆川仪自动化股份有限公司 | 红外传感器输出信号的滤波方法、装置及红外气体分析仪 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101153840A (zh) * | 2006-09-29 | 2008-04-02 | 深圳迈瑞生物医疗电子股份有限公司 | 一种提高气体分析仪测量精度的方法和装置 |
US20080166064A1 (en) * | 2007-01-05 | 2008-07-10 | Guoyi Fu | Method And Apparatus For Reducing Noise In An Image Using Wavelet Decomposition |
-
2010
- 2010-07-02 CN CN 201010216322 patent/CN101900674A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101153840A (zh) * | 2006-09-29 | 2008-04-02 | 深圳迈瑞生物医疗电子股份有限公司 | 一种提高气体分析仪测量精度的方法和装置 |
US20080166064A1 (en) * | 2007-01-05 | 2008-07-10 | Guoyi Fu | Method And Apparatus For Reducing Noise In An Image Using Wavelet Decomposition |
Non-Patent Citations (1)
Title |
---|
《硕士学位论文》 20070713 程涛 基于小波分析的上海市环境空气质量变化及与气象关系研究 第10-19页 1 , 2 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102818629A (zh) * | 2012-05-04 | 2012-12-12 | 浙江大学 | 基于平稳小波变换的微型光谱仪信号去噪方法 |
CN103513094A (zh) * | 2013-09-29 | 2014-01-15 | 天津理工大学 | 一种消除电力系统检测信号噪声的装置及方法 |
CN103674847A (zh) * | 2013-12-17 | 2014-03-26 | 重庆川仪自动化股份有限公司 | 红外传感器输出信号的滤波方法、装置及红外气体分析仪 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102028457B (zh) | 脉率测量方法及指环式脉率测量仪 | |
CN101972148B (zh) | 基于经验模态分解的近红外脑功能检测的扰动消除方法 | |
US8050730B2 (en) | Method and apparatus for eliminating interference in pulse oxygen measurement | |
CN108108889A (zh) | 一种水质监测数据在线处理方法及装置 | |
CN106130660B (zh) | 一种无线光通信盲检测接收机的方法 | |
CN103093759B (zh) | 一种基于移动终端的嗓音检测评估装置及方法 | |
CN106096579A (zh) | 一种心电信号预处理的方法 | |
CN102973279B (zh) | 独立成分分析联合最小二乘法的近红外脑机接口的信号检测方法 | |
CN101153840B (zh) | 一种提高气体分析仪测量精度的方法和装置 | |
CN102824172B (zh) | 一种脑电特征提取方法 | |
CN108458871A (zh) | 一种基于改进经验小波变换的齿轮箱故障识别方法 | |
CN106859648A (zh) | 基于非接触式检测的多目标人体呼吸信号监测方法及装置 | |
CN103616354A (zh) | 一种藻类浓度原位荧光检测装置 | |
CN107563969A (zh) | 基于变分模态分解的dspi相位滤波方法 | |
CN101900674A (zh) | 一种用于二氧化碳浓度检测信号静噪的小波滤波方法 | |
CN107693004A (zh) | 基于hilbert变换的胎儿心电提取与胎儿心率识别方法 | |
CN111027488A (zh) | 一种信号分类的方法和设备 | |
Schoen et al. | Performance of a holographic supernotch filter | |
CN103932687A (zh) | 一种脉象信号预处理方法和装置 | |
CN114707555A (zh) | 崩岸前兆次声波信号降噪处理方法、特征识别方法及装置 | |
Zhang et al. | Variable spectral segmentation empirical wavelet transform for noisy signal processing | |
CN115855957A (zh) | 一种基于光电信号的激光焊接质量在线监测系统和方法 | |
CN109712639A (zh) | 一种基于小波阈值滤波的音频采集系统及方法 | |
Wu et al. | An improved method for ECG signal feature point detection based on wavelet transform | |
Zheng et al. | A new noise reduction method based on re-weighted group sparse decomposition and its application in gear fault feature detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20101201 |