CN108985179B - 一种基于改进小波阈值函数的电能质量信号去噪方法 - Google Patents
一种基于改进小波阈值函数的电能质量信号去噪方法 Download PDFInfo
- Publication number
- CN108985179B CN108985179B CN201810647073.8A CN201810647073A CN108985179B CN 108985179 B CN108985179 B CN 108985179B CN 201810647073 A CN201810647073 A CN 201810647073A CN 108985179 B CN108985179 B CN 108985179B
- Authority
- CN
- China
- Prior art keywords
- wavelet
- signal
- power quality
- threshold
- denoising
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/003—Environmental or reliability tests
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Measurement Of Current Or Voltage (AREA)
- Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
Abstract
本发明涉及一种基于改进小波阈值函数的电能质量信号去噪方法,首先采集原始电能质量信号;对其进行预处理,并选择db3小波基函数对预处理后的一维染噪电能质量模拟信号进行3层小波分解,得到一组小波系数;然后利用改进的小波阈值函数与使用统一阈值法计算出的阈值对进行阈值量化处理,得到估计小波系数;最后由阈值量化处理后的第1层到第3层的高频小波系数和第3层的低频小波系数实行小波逆变换,进行信号重构,得到重构信号。本发明改进后的小波阈值函数解决了去噪与保留原始电能质量信号局部特征信息之间的矛盾问题。
Description
技术领域
本发明涉及电能质量信号去噪方法,特别是一种基于改进小波阈值函数的电能质量信号去噪方法。
背景技术
近年来,随着各种电力电子器件的广泛使用及非线性、冲击性、波动性负荷的不断增加,由此带来的电能质量问题越来越引起电力部门和电力用户的高度重视,那么如何准确、快速、有效地检测各种类型的电能质量事件显得尤为重要。
然而在实际应用中,利用监测与数据采集设备获取的电能质量信号会因设备运行内、外部条件变化等原因引入随机噪声;另外,在测量和信道传输数据时,电能质量信号也会被噪声污染。噪声的存在会恶化扰动检测和识别方法的有效性,有时甚至会使它们失效。因此近年来关于电能质量信号去噪的研究也引起高度重视。
电能质量信号去噪的难点是在有效去除噪声的同时,保留原始信号的局部特征信息。近年来,国内外对电能质量信号去噪问题已经展开广泛而深入的研究和探讨,涌现了大量的去噪方法,如线性滤波技术去噪、高斯滤波技术去噪、EMD分解技术去噪、小波阈值技术去噪、似然比判决准则去噪、数学形态法、奇异值分解去噪、S变换去噪法等,这些方法各有利弊,但均无法克服电能质量信号去噪的难点。线性滤波技术只针对高斯噪声有明显的平滑去噪作用,而对于常见的突变的信号特征如异于高斯分布的信号和高频分量,却无法达到很好的去噪效果,甚至会对突变信息模糊化从而使有效的信号发生失真。高斯滤波技术需要进行多次迭代运算。EMD分解技术主要用于处理非线性或非平稳信号,不适用于平稳信号的固定参数提取。似然比判决准则去噪,需要预先进行突变点判决找到突变点位置,然而在噪声强度较高的情况下,突变点位置的确定将十分困难,从而影响去噪效果。对于数学形态法,由于实际脉冲信号的宽度未知,其形态滤波器可能会将脉冲扰动误当成噪声滤掉。奇异值分解去噪法的平滑效果较差,重构信号中噪声含量偏大。
小波阈值去噪法是由D.L.Donoho于1995年提出的一种基于小波变换去噪法的一种改进去噪方法,这种方法可以在有效去噪的同时,保留原始信号突变点信息,得到了许多学者的认可、关注及进一步的研究。针对传统小波软阈值、硬阈值函数去噪的不足,提出了许多的改进阈值函数,但均存在如下问题,如去噪后信号在高频处仍然存在信号模糊的现象、计算量较大等,均无法解决在有效去噪的同时最大程度保留原始信号的局部特征信息。
发明内容
有鉴于此,本发明的目的是提出一种基于改进小波阈值函数的电能质量信号去噪方法,可以在有效去噪的同时最大程度保留原始信号的局部特征信息。
本发明采用以下方案实现:一种基于改进小波阈值函数的电能质量信号去噪方法,包括以下步骤:
步骤S1:采集原始电能质量信号;
步骤S2:对步骤S1采集的原始电能质量信号进行预处理,并加入高斯白噪声(由于电能质量信号中夹杂的噪声信号是符合高斯分布的高斯白噪声信号),得到一维染噪电能质量模拟信号f(t);
步骤S3:对步骤S2获得的一维染噪电能质量模拟信号f(t)进行3层小波分解,得到一组小波系数ωj,k;
进一步地,步骤S1中,所述原始电能质量信号包括电压暂降的电能质量信号、电压暂升的电能质量信号、电压中断的电能质量信号。本发明利用利用监测与数据采集设备获取电能质量信号。
进一步地,步骤S2中,所述预处理包括对获得的原始电能质量信号进行时间段的截取与幅值处理。
进一步地,步骤S3中采用db3小波基函数对一维染噪电能质量模拟信号f(t)进行3层小波分解。
进一步地,步骤S4具体包括以下步骤:
步骤S41:采用统一阈值法对步骤S2得到的一维染噪电能质量模拟信号f(t)进行计算,得到相应的阈值T,具体采用下式:
式中,N为信号的采样长度,σ为噪声信号的标准差;
式中,sgn为符号函数,ξ表示加权因子。
进一步地,步骤S42中,为了解决传统小波软阈值函数与硬阈值函数去噪后造成的小波系数偏差、“伪吉布斯”现象等问题,所述加权因子的取值为:
这样在去噪过程中随着原始信号频率的变化,小波变换后的系数也随之变化,改进的阈值函数始终介于软、硬阈值函数之间,可使去噪后的小波估计系数最大化地接近有效信号的小波系数。
与现有技术相比,本发明有以下有益效果:本发明提出的改进后的小波阈值函数克服了传统小波阈值函数处理造成的“伪吉布斯”现象、小波系数偏差等问题,弥补了小波软、硬阈值去噪的缺点,其适应性更强。本发明可以在有效去噪的同时最大程度保留原始信号的局部特征信息。
附图说明
图1为本发明实施例的方法流程示意图。
图2为本发明实施例的改进阈值函数的去噪特性曲线图。
图3为本发明实施例的使用软阈值、硬阈值、折中阈值及改进阈值函数的电压暂降信号去噪对比图。
图4为本发明实施例的使用软阈值、硬阈值、折中阈值及改进阈值函数的电压暂升信号去噪对比图。
图5为本发明实施例的使用软阈值、硬阈值、折中阈值及改进阈值函数的电压中断信号去噪对比图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
如图1所示,本实施例提供了一种基于改进小波阈值函数的电能质量信号去噪方法,包括以下步骤:
步骤S1:采集原始电能质量信号;
步骤S2:对步骤S1采集的原始电能质量信号进行预处理,并加入高斯白噪声(由于电能质量信号中夹杂的噪声信号是符合高斯分布的高斯白噪声信号),得到一维染噪电能质量模拟信号f(t);
步骤S3:对步骤S2获得的一维染噪电能质量模拟信号f(t)进行3层小波分解,得到一组小波系数ωj,k;
在本实施例中,步骤S1中,所述原始电能质量信号包括电压暂降的电能质量信号、电压暂升的电能质量信号、电压中断的电能质量信号。本发明利用利用监测与数据采集设备获取电能质量信号。
在本实施例中,步骤S2中,所述预处理包括对获得的原始电能质量信号进行时间段的截取与幅值处理。
在本实施例中,步骤S3中采用db3小波基函数对一维染噪电能质量模拟信号f(t)进行3层小波分解。
在本实施例中,步骤S4具体包括以下步骤:
步骤S41:采用统一阈值法对步骤S2得到的一维染噪电能质量模拟信号f(t)进行计算,得到相应的阈值T,具体采用下式:
式中,N为信号的采样长度,σ为噪声信号的标准差;
式中,sgn为符号函数,ξ表示加权因子。
在本实施例中,步骤S42中,为了解决传统小波软阈值函数与硬阈值函数去噪后造成的小波系数偏差、“伪吉布斯”现象等问题,所述加权因子的取值为:
这样在去噪过程中随着原始信号频率的变化,小波变换后的系数也随之变化,改进的阈值函数始终介于软、硬阈值函数之间,可使去噪后的小波估计系数最大化地接近有效信号的小波系数。
特别的,在本实施例中,如图2所示为改进阈值函数的去噪特性曲线图,如图3、4、5所示为三种典型电能质量信号(电压暂降、电压暂升,电压中断)分别采用软阈值、硬阈值、折中阈值及改进阈值函数的去噪结果对比图。由图2可以看出,在去噪过程中随着原始信号频率的变化,小波变换后的系数也随之变化,改进的阈值函数始终介于软、硬阈值函数之间,使去噪后的小波估计系数最大化地接近有效信号的小波系数。由图3、4、5可以看出,软阈值函数及硬阈值函数去噪法对含有高斯白噪声的电能质量信号去噪后无法完整地重构出原始信号,尤其是在幅值较大处,信号模糊更为明显;折中阈值函数较好地重构出了原始信号,但在信号波峰、波谷处出现了微小的畸变;改进的阈值函数对含有高斯白噪声的电能质量信号去噪后,重构信号较其他三种更为平滑,与原始信号相似度极高,较好地在有效去噪的同时保留了原始电能质量信号包含的信息。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。
Claims (6)
1.一种基于改进小波阈值函数的电能质量信号去噪方法,其特征在于:包括以下步骤:
步骤S1:采集原始电能质量信号;
步骤S2:对步骤S1采集的原始电能质量信号进行预处理,并加入高斯白噪声,得到一维染噪电能质量模拟信号f(t);
步骤S3:对步骤S2获得的一维染噪电能质量模拟信号f(t)进行3层小波分解,得到一组小波系数ωj,k;
步骤S4具体包括以下步骤:
步骤S41:采用统一阈值法对步骤S2得到的一维染噪电能质量模拟信号f(t)进行计算,得到相应的阈值T,具体采用下式:
式中,N为信号的采样长度,σ为噪声信号的标准差;
式中,sgn为符号函数,ξ表示加权因子。
2.根据权利要求1所述的一种基于改进小波阈值函数的电能质量信号去噪方法,其特征在于:步骤S1中,所述原始电能质量信号包括电压暂降的电能质量信号、电压暂升的电能质量信号、电压中断的电能质量信号。
3.根据权利要求1所述的一种基于改进小波阈值函数的电能质量信号去噪方法,其特征在于:步骤S2中,所述预处理包括对获得的原始电能质量信号进行时间段的截取与幅值处理。
4.根据权利要求1所述的一种基于改进小波阈值函数的电能质量信号去噪方法,其特征在于:步骤S3中采用db3小波基函数对一维染噪电能质量模拟信号f(t)进行3层小波分解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810647073.8A CN108985179B (zh) | 2018-06-22 | 2018-06-22 | 一种基于改进小波阈值函数的电能质量信号去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810647073.8A CN108985179B (zh) | 2018-06-22 | 2018-06-22 | 一种基于改进小波阈值函数的电能质量信号去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108985179A CN108985179A (zh) | 2018-12-11 |
CN108985179B true CN108985179B (zh) | 2022-01-25 |
Family
ID=64538077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810647073.8A Active CN108985179B (zh) | 2018-06-22 | 2018-06-22 | 一种基于改进小波阈值函数的电能质量信号去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108985179B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110220708A (zh) * | 2019-04-24 | 2019-09-10 | 邯郸钢铁集团有限责任公司 | 一种基于改进小波算法的轴承信号降噪方法 |
CN110329910A (zh) * | 2019-04-30 | 2019-10-15 | 上海冠东国际集装箱码头有限公司 | 自动化集装箱岸桥远程操控小车路径规划方法 |
CN110703089B (zh) * | 2019-08-29 | 2021-09-24 | 广西电网有限责任公司电力科学研究院 | 一种用于低频振荡Prony分析的小波阈值去噪方法 |
CN110584649B (zh) * | 2019-09-09 | 2021-01-15 | 武汉大学 | 一种动态单导心电图中异常心电信号识别方法及系统 |
CN110726539A (zh) * | 2019-09-27 | 2020-01-24 | 国网江苏省电力有限公司盐城供电分公司 | 一种基于多种信号分析的变压器分接开关机械故障检测系统及其检测方法 |
CN111336066A (zh) * | 2020-03-04 | 2020-06-26 | 湖南城市学院 | 一种用于产生能量的风力发电系统及控制方法 |
CN111537789A (zh) * | 2020-05-09 | 2020-08-14 | 湖南省计量检测研究院 | 一种基于信号分离和精准积分的直流电能计量装置及方法 |
CN112132021B (zh) * | 2020-09-22 | 2023-10-31 | 三峡大学 | 基于wmm与hatf的暂态电能质量扰动检测方法 |
CN112307959B (zh) * | 2020-10-30 | 2022-10-11 | 合肥工业大学 | 一种用于心电信号分析的小波去噪方法 |
CN112395992A (zh) * | 2020-11-18 | 2021-02-23 | 云南电网有限责任公司电力科学研究院 | 一种基于改进小波阈值的电力谐波信号去噪方法 |
CN112528853B (zh) * | 2020-12-09 | 2021-11-02 | 云南电网有限责任公司昭通供电局 | 改进型双树复小波变换去噪方法 |
CN113405747A (zh) * | 2021-04-26 | 2021-09-17 | 深圳华星智感科技有限公司 | 一种基于小波分析及阈值处理的监测信号滤波方法 |
CN113239751A (zh) * | 2021-04-27 | 2021-08-10 | 电子科技大学 | 一种基于加权因子的小波阈值去噪方法 |
CN113835062A (zh) * | 2021-09-10 | 2021-12-24 | 中国科学院上海微系统与信息技术研究所 | 基于小波去噪和music算法的波达方向定位方法 |
CN113553997A (zh) * | 2021-09-17 | 2021-10-26 | 南京信息工程大学 | 一种联合改进小波阈值的信号处理方法、存储介质和系统 |
CN114942065B (zh) * | 2022-05-17 | 2023-07-11 | 湖北工业大学 | 称重信号降噪方法、装置、电子设备及计算机存储介质 |
CN114841213B (zh) * | 2022-05-19 | 2023-04-07 | 东南大学 | 基于改进小波阈值函数的硅微谐振式加速度计降噪方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103471623A (zh) * | 2013-09-24 | 2013-12-25 | 重庆邮电大学 | 一种基于邻域小波系数的mzi信号去噪方法 |
CN103513094A (zh) * | 2013-09-29 | 2014-01-15 | 天津理工大学 | 一种消除电力系统检测信号噪声的装置及方法 |
CN105572517A (zh) * | 2016-02-26 | 2016-05-11 | 贵州电网有限责任公司铜仁供电局 | 一种避雷器在线监测智能带电接入保护装置及控制方法 |
CN106562782A (zh) * | 2016-05-20 | 2017-04-19 | 彭慧敏 | 一种儿科护理专用心电监护仪 |
CN107274908A (zh) * | 2017-06-13 | 2017-10-20 | 南京邮电大学 | 基于新阈值函数的小波语音去噪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7526130B2 (en) * | 2005-10-18 | 2009-04-28 | Optoelectronics Co., Ltd. | Signal processing in barcode reading with a wavelet transformation |
-
2018
- 2018-06-22 CN CN201810647073.8A patent/CN108985179B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103471623A (zh) * | 2013-09-24 | 2013-12-25 | 重庆邮电大学 | 一种基于邻域小波系数的mzi信号去噪方法 |
CN103513094A (zh) * | 2013-09-29 | 2014-01-15 | 天津理工大学 | 一种消除电力系统检测信号噪声的装置及方法 |
CN105572517A (zh) * | 2016-02-26 | 2016-05-11 | 贵州电网有限责任公司铜仁供电局 | 一种避雷器在线监测智能带电接入保护装置及控制方法 |
CN106562782A (zh) * | 2016-05-20 | 2017-04-19 | 彭慧敏 | 一种儿科护理专用心电监护仪 |
CN107274908A (zh) * | 2017-06-13 | 2017-10-20 | 南京邮电大学 | 基于新阈值函数的小波语音去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108985179A (zh) | 2018-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108985179B (zh) | 一种基于改进小波阈值函数的电能质量信号去噪方法 | |
CN108158573B (zh) | 基于自适应阈值小波变换的心电信号降噪方法 | |
WO2021056727A1 (zh) | 一种基于变分模态分解和排列熵的联合降噪方法 | |
CN112395992A (zh) | 一种基于改进小波阈值的电力谐波信号去噪方法 | |
Lin et al. | A new threshold function for signal denoising based on wavelet transform | |
CN106570843A (zh) | 一种自适应小波阈值函数图像噪声抑制方法 | |
Kopsinis et al. | Empirical mode decomposition based soft-thresholding | |
CN113378661A (zh) | 一种基于改进小波阈值和相关检测的直流电能信号去噪方法 | |
Zhang et al. | Seismic exploration desert noise suppression based on complete ensemble empirical mode decomposition with adaptive noise | |
CN109409281A (zh) | 一种基于改进的小波阈值函数的降噪方法 | |
Zhang et al. | A novel ECG signal denoising method based on Hilbert-Huang Transform | |
CN110136086A (zh) | 基于bemd的区间阈值图像去噪方法 | |
CN110032968A (zh) | 基于双树复小波和自适应半软阈值法的去噪方法 | |
Khadtare et al. | Removal of 50Hz PLI using discrete wavelet transform for quality diagnosis of biomedical ECG signal | |
CN114358040A (zh) | 一种基于自适应估算阈值法的扰动信号去噪方法 | |
Zebin et al. | Research on ecg denoising method based on empirical mode decomposition and wavelet transform | |
Cui et al. | Research and application of new threshold de-noising algorithm for monitoring data analysis in nuclear power plant | |
CN113505688A (zh) | 一种改进的自适应小波阈值信号去噪方法 | |
CN113487505A (zh) | 一种基于噪声识别的红外图像混合降噪方法 | |
Ce et al. | An improved wavelet threshold de-noising method and its application | |
Wang et al. | Image edge detection algorithm based onwavelet fractional differential theory | |
Qin et al. | an image denoising method based on improved wavelet thresholding | |
CN109271971B (zh) | 一种用于时序金融数据的降噪方法 | |
Wang et al. | Comparative Study of Tongue Image Denoising Methods. | |
Lin et al. | Image denoising based on non-local means with Wiener filtering in wavelet domain |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |