CN109902594A - 一种基于小波的脉冲信号去噪方法 - Google Patents
一种基于小波的脉冲信号去噪方法 Download PDFInfo
- Publication number
- CN109902594A CN109902594A CN201910097404.XA CN201910097404A CN109902594A CN 109902594 A CN109902594 A CN 109902594A CN 201910097404 A CN201910097404 A CN 201910097404A CN 109902594 A CN109902594 A CN 109902594A
- Authority
- CN
- China
- Prior art keywords
- pulse signal
- wavelet
- small echo
- coefficient
- wavelet coefficient
- 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.)
- Withdrawn
Links
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于小波的脉冲信号去噪方法，包括以下步骤：(1)获取包含噪声的脉冲信号s，选用小波母函数做l层的离散小波变换；(2)对小波系数做阈值为T的阈值化处理，得到修饰过的小波系数(3)用修饰过的小波系数重构得到去噪信号本发明一种基于小波的脉冲信号去噪方法，相比于傅里叶变换中使用的三角函数，小波变换中的小波母函数具有很好的时频局部化特性，更适合脉冲信号的去噪处理；能够很好的保留信号特征，能够完整保留时间信息；具有丰富的小波函数和阈值取值方法，适用于不同噪声强度的脉冲信号的去噪优化。
Description
技术领域
本发明涉及信号去噪优化领域，特别涉及一种基于小波的脉冲信号去噪方法。
背景技术
在工程实际应用中，脉冲形状的信号十分常见，例如局部放电信号、雷电流信号和汽车变速箱故障信号。脉冲信号代表其相关的事件具有突发性和时效性，有效的检测脉冲信号具有重大的实际应用意义。而实际的脉冲信号通常伴有强烈的噪声，有效的去除噪声并且尽可能地保留信号的原有特征是检测脉冲信号中至关重要的一步。
传统的脉冲信号去噪方法通常基于傅里叶变换，将信号从时域变换到频域处理。对于脉冲信号，其持续时间很短，时频变换会使信号中的时间信息完全丢失。
发明内容
本发明的目的是提供一种基于小波的脉冲信号去噪方法，能够很好的保留信号特征，并且适用于不同噪声强度的脉冲信号的去噪优化。
为达到上述目的，本发明采用的技术方案是：
一种基于小波的脉冲信号去噪方法，包括以下步骤：
(1)获取包含噪声的脉冲信号s，选用小波母函数做l层的离散小波变换；
(2)对小波系数做阈值为T的阈值化处理，得到修饰过的小波系数
(3)用修饰过的小波系数重构得到去噪信号
优选地，所述小波母函数为Haar小波或Db小波或Sym小波。
优选地，在步骤(1)中，信号s可以分解为：
在公式(a)中，l为设定分解的层数；和φ(t)分别为小波函数和尺度函数；cl,k和dl,k分别为l尺度下的尺度系数和小波系数，且递推计算式为：
在公式(b)中，h0(n)和h1(n)为展开系数，不随尺度l的变化而变化，具体关系为：
其中，尺度系数cl,k和小波系数dl,k分别对应信号的低频部分和高频部分；脉冲信号s中的噪声包含在高频部分，而脉冲信号s中的特征则包含在低频部分。
更优选地，在步骤(2)中，将小波系数dl,k选择合适的阈值进行阈值化处理；小波分解各层的阈值Tl为：
在公式(d)中，ml为第l层小波系数的中位数，n为带噪脉冲信号s的长度；在具体阈值化处理中，包括硬阈值化和软阈值化两种方法；
硬阈值化之后的小波系数为：
软阈值化之后的小波系数为：
更进一步优选地，在步骤(3)中，使用阈值化之后的小波系数和尺度系数cl,k重构得到去噪之后的信号
由于上述技术方案的运用，本发明与现有技术相比具有下列优点：本发明一种基于小波的脉冲信号去噪方法，相比于傅里叶变换中使用的三角函数，小波变换中的小波母函数具有很好的时频局部化特性，更适合脉冲信号的去噪处理；能够很好的保留信号特征，能够完整保留时间信息；具有丰富的小波函数和阈值取值方法，适用于不同噪声强度的脉冲信号的去噪优化。
附图说明
附图1为本发明方法对脉冲信号的去噪流程图；
附图2为本发明方法拟去噪的脉冲信号及其经过小波分解之后的尺度系数和小波系数；
附图3为本发明方法拟去噪的脉冲信号和去噪信号。
具体实施方式
下面结合附图来对本发明的技术方案作进一步的阐述。
一种基于小波的脉冲信号去噪方法，包括以下步骤：
(1)对图2所示噪声污染脉冲信号s选用小波母函数(具体为Db7小波)做l＝5层的离散小波变换。实际计算中，常使用Mallet算法进行快速计算。信号s可以分解为：
在公式(a)中，l为设定分解的层数；和φ(t)分别为小波函数和尺度函数；cl,k和dl,k分别为l尺度下的尺度系数和小波系数，且递推计算式为：
在公式(b)中，h0(n)和h1(n)为展开系数，不随尺度l的变化而变化，具体关系为：
其中，尺度系数cl,k和小波系数dl,k分别对应信号的低频部分和高频部分；脉冲信号s中的噪声包含在高频部分，而脉冲信号s中的特征则包含在低频部分。
使用公式(a)、(b)、(c)，如图2所示，脉冲信号s被分解为近似部分a5和细节部分d1、d2、d3、d4和d5。
(2)将步骤1得到的尺度系数a5和小波系数d1、d2、d3、d4、d5分别对应信号的低频部分和高频部分。在本实施例中，将小波系数dl,k选择固定阈值法进行软阈值化处理，可以对脉冲信号s实现去噪。针对脉冲信号s持续时间等具体信号特征，固定阈值法可以实现较佳的去噪效果。小波分解各层的阈值Tl为：
在公式(d)中，ml为第l层小波系数的中位数，n为带噪脉冲信号s的长度；具体各层的阈值为：
软阈值化之后的小波系数为：
(3)如图3所示，使用阈值化之后的小波系数和原有的尺度系数cl,k重构得到去噪之后的信号
上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人士能够了解本发明的内容并加以实施，并不能以此限制本发明的保护范围，凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围内。
Claims (5)
1.一种基于小波的脉冲信号去噪方法，其特征在于：包括以下步骤：
(1)获取包含噪声的脉冲信号s，选用小波母函数做l层的离散小波变换；
(2)对小波系数做阈值为T的阈值化处理，得到修饰过的小波系数
(3)用修饰过的小波系数重构得到去噪信号
2.根据权利要求1所述的一种基于小波的脉冲信号去噪方法，其特征在于：所述小波母函数为Haar小波或Db小波或Sym小波。
3.根据权利要求1所述的一种基于小波的脉冲信号去噪方法，其特征在于：在步骤(1)中，信号s可以分解为：
在公式(a)中，l为设定分解的层数；和φ(t)分别为小波函数和尺度函数；cl,k和dl,k分别为l尺度下的尺度系数和小波系数，且递推计算式为：
在公式(b)中，h0(n)和h1(n)为展开系数，不随尺度l的变化而变化，具体关系为：
其中，尺度系数cl,k和小波系数dl,k分别对应信号的低频部分和高频部分；脉冲信号s中的噪声包含在高频部分，而脉冲信号s中的特征则包含在低频部分。
4.根据权利要求3所述的一种基于小波的脉冲信号去噪方法，其特征在于：在步骤(2)中，将小波系数dl,k选择合适的阈值进行阈值化处理；小波分解各层的阈值Tl为：
在公式(d)中，ml为第l层小波系数的中位数，n为带噪脉冲信号s的长度；在具体阈值化处理中，包括硬阈值化和软阈值化两种方法；
硬阈值化之后的小波系数为：
软阈值化之后的小波系数为：
5.根据权利要求4所述的一种基于小波的脉冲信号去噪方法，其特征在于：在步骤(3)中，使用阈值化之后的小波系数和尺度系数cl,k重构得到去噪之后的信号
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910097404.XA CN109902594A (zh) | 2019-01-31 | 2019-01-31 | 一种基于小波的脉冲信号去噪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910097404.XA CN109902594A (zh) | 2019-01-31 | 2019-01-31 | 一种基于小波的脉冲信号去噪方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109902594A true CN109902594A (zh) | 2019-06-18 |
Family
ID=66944498
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910097404.XA Withdrawn CN109902594A (zh) | 2019-01-31 | 2019-01-31 | 一种基于小波的脉冲信号去噪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109902594A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110430340A (zh) * | 2019-06-28 | 2019-11-08 | 北京大学 | 一种脉冲阵列信号的降噪方法与系统 |
| CN112418174A (zh) * | 2020-12-08 | 2021-02-26 | 中国石油天然气集团有限公司 | 随钻泥浆随机噪声的去除方法 |
| CN114461975A (zh) * | 2022-01-27 | 2022-05-10 | 无锡科晟光子科技有限公司 | 基于小波变换去噪的光纤测音方法 |
-
2019
- 2019-01-31 CN CN201910097404.XA patent/CN109902594A/zh not_active Withdrawn
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110430340A (zh) * | 2019-06-28 | 2019-11-08 | 北京大学 | 一种脉冲阵列信号的降噪方法与系统 |
| CN110430340B (zh) * | 2019-06-28 | 2020-08-25 | 北京大学 | 一种脉冲阵列信号的降噪方法与系统 |
| CN112418174A (zh) * | 2020-12-08 | 2021-02-26 | 中国石油天然气集团有限公司 | 随钻泥浆随机噪声的去除方法 |
| CN114461975A (zh) * | 2022-01-27 | 2022-05-10 | 无锡科晟光子科技有限公司 | 基于小波变换去噪的光纤测音方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103630808B (zh) | 一种基于提升小波变换的局部放电信号去噪方法 | |
| CN109902594A (zh) | 一种基于小波的脉冲信号去噪方法 | |
| CN103295204B (zh) | 一种基于非下采样轮廓波变换的图像自适应增强方法 | |
| CN107516065A (zh) | 经验模态分解结合字典学习的复杂信号去噪方法 | |
| CN104635223A (zh) | 基于经验模式分解与分数阶傅里叶变换的激光回波去噪方法 | |
| CN112084845B (zh) | 基于多尺度小波系数自相关的低频1/f噪声消除方法 | |
| CN105182418A (zh) | 一种基于双树复小波域的地震信号降噪方法及系统 | |
| CN104732498B (zh) | 一种基于非下采样Contourlet变换的阈值化图像去噪方法 | |
| CN101876585A (zh) | 基于小波包估计噪声方差的ica收缩去噪方法 | |
| CN106296583B (zh) | 基于图像块组稀疏编码与成对映射的含噪高光谱图像超分辨率重构方法 | |
| CN104978716A (zh) | 一种基于线性最小均方误差估计的sar图像降噪方法 | |
| CN113589253A (zh) | 一种基于伪时域的小波变换算法对微弱回波信号进行检测的方法 | |
| Zhao et al. | Speech signal feature extraction based on wavelet transform | |
| CN109409281A (zh) | 一种基于改进的小波阈值函数的降噪方法 | |
| CN103915102B (zh) | 一种lfm水声多途信号的噪声抑制方法 | |
| CN103426145A (zh) | 一种基于多分辨分析的合成孔径声纳相干斑噪声抑制方法 | |
| CN109003247A (zh) | 彩色图像混合噪声的去除方法 | |
| CN102289793B (zh) | 一种面向游牧服务的多尺度图像处理方法 | |
| CN102509268B (zh) | 基于免疫克隆选择的非下采样轮廓波域图像去噪方法 | |
| CN103986539B (zh) | 一种基于稀疏去噪的认知无线电频谱感知方法 | |
| CN106024006A (zh) | 基于小波变换的鲸类叫声信号去噪方法和装置 | |
| CN104463325A (zh) | 一种极地探冰雷达原始数据噪声抑制方法 | |
| CN104282042B (zh) | 一种基于段移位的可逆3d多边形网格数据隐藏方法 | |
| CN106651933A (zh) | 一种基于Radon变换的改进型纹理特征提取算法 | |
| Kong et al. | Application and optimization of wavelet threshold denoising algorithm in signal processing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20190618 |
|
| WW01 | Invention patent application withdrawn after publication |