CN105490516B - Atem发射机发射电流pwm噪声的去噪方法 - Google Patents
Atem发射机发射电流pwm噪声的去噪方法 Download PDFInfo
- Publication number
- CN105490516B CN105490516B CN201510917176.8A CN201510917176A CN105490516B CN 105490516 B CN105490516 B CN 105490516B CN 201510917176 A CN201510917176 A CN 201510917176A CN 105490516 B CN105490516 B CN 105490516B
- Authority
- CN
- China
- Prior art keywords
- threshold
- denoising
- soft
- threshold value
- hard
- 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.)
- Expired - Fee Related
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/44—Circuits or arrangements for compensating for electromagnetic interference in converters or inverters
Abstract
本发明涉及一种ATEM发射机发射电流PWM噪声的去噪方法,利用小波变换对待去噪电压信号进行小波分解,并得到细节系数与尺度系数;应用Stein无偏似然估计阈值选取方法对采样波形的PWM噪声进行阈值选取;应用上一步骤所选取的阈值,对细节系数与尺度系数进行改进软硬阈值方法结合的去噪处理,然后重构得到不含噪的信号;输出不含噪的重构信号。用Stein无偏似然估计阈值选取方法和软硬阈值结合的处理方法进行小波系数处理,既保留了Stein无偏似然估计阈值选取方法保留的有效信号成分更多的优点,使用具有改进软硬结合阈值处理方法对阈值选取方法进行了补偿,采取两种方法结合的方式,最大限度的保留了原始信号,能使去噪效果更为明显。
Description
技术领域
本发明涉及一种地球物理勘探设备的去噪方法,尤其是ATEM发射机发射电流PWM噪声的去噪方法。
背景技术
近年来,在电力电子领域,PWM控制技术的应用越来越广泛。尤其是在大功率发射机的应用方面,使用PWM控制策略对电压信号进行调制,可以达到很好的控制效果。但是,在输出波形达到预期效果的同时,在完整的频率范围内,尤其是在高频部分,PWM也会产生较大的噪声,并且几乎覆盖了原始的信号,从而不能观察到不失真的输出电压信号。
采用高通或低通滤波器的方法来去除噪声是最传统也是最简单的方法,由于PWM噪声大都集中在高频部分,因此可以使用LC低通滤波器来滤波,但是低通滤波器会把高频的有效信号部分同时滤除,会造成信号的失真。因此,更需要一种可以在全频范围内既可以滤除噪声信号又可以保留有效信号的滤波方法。在去噪的过程中,阈值的选取非常重要,好的阈值选取,可以有效降低失真度,从而获得更好的去噪效果。Stein无偏似然估计阈值选取方法是一种基于Stein无偏似然估计原理的自适应阈值获取方法,其特点是不易丢失混合信号中的有用信号,从而降低失真度。获得了阈值之后,对于阈值的处理同样是很重要,现有的阈值处理方法主要有软阈值与硬阈值两种方法。硬阈值是把小波系数的绝对值与阈值相比较,小于或等于阈值的小波系数设置为零,大于阈值的小波系数不变;软阈值处理则是把大于阈值的点变为该点与阈值的差值。硬阈值的处理方法容易出现间断,而软阈值虽然较为平滑,但去噪的效果不好。
ATEM发射机是一种发射双极性脉冲电流的发射装置,以低阻抗高感抗发射线圈为负载,产生大功率正负交变电流,激发出强烈的一次电磁场。但是由于ATEM发射机采用PWM脉冲宽度调制方式的控制策略,因此在测试发射波形时会产生干扰,采样波形如图2所示。因此本文采用电感很低的采样电阻接在发射机的输出端,最初用示波器测量输出电压波形,但由于精度较低,测量结果几乎全部都是噪声。
发明内容
本发明的目的就在于针对上述现有技术的不足,提供一种基于Stein无偏似然估计阈值选取方法和改进软硬阈值结合阈值处理方法结合的ATEM发射机发射电流PWM噪声去噪方法。
本发明的目的是通过以下技术方案实现的:
ATEM发射机发射电流PWM噪声的去噪方法,是基于Stein无偏似然估计阈值选取方法和改进软硬阈值结合阈值处理方法的结合,包括以下步骤:
A、利用小波变换对待去噪电压信号进行多层小波变换处理,并得到各层的细节系数与尺度系数;
B、应用Stein无偏似然估计阈值选取方法对待去噪电压波形的PWM噪声进行阈值选取;
C、应用步骤B所选取的针对PWM噪声的阈值,对细节系数与尺度系数进行改进的软、硬阈值结合的阈值处理;
D、重构得到去噪后的信号序列。
步骤B所述的阈值选取步骤为:
a、将步骤2所得的小波系数的平方按从小打到大的顺序排列,即:
pα={p1,p2,…,pn},且p1≤p2≤…≤pn;
b、定义风险向量R,其元素为:
c、以R元素中的最小rmin作为风险值,找出rmin的下标变量min所对应的α,进而获得对应的pα,求出阈值T,即
其中,N为含噪信号在所有尺度上的小波分解系数的个数总和,σ为噪声信号的偏差。
步骤C所述的改进的软、硬阈值结合的阈值处理函数表达式为:
其中:γ是调节因子,介于0和1之间。在|wj,k|>λ时,可以通过调节α使新阈值函数在软硬阈值之间灵活选择,从而获得最优的去噪效果。
有益效果:本发明是针对由于PWM噪声存在于全频范围内,去噪失真度高的问题,在进行阈值处理的过程中,采用Stein无偏似然估计阈值选取方法和软硬阈值结合的与之处理方法相结合的方式进行小波系数处理,既保留了Stein无偏似然估计阈值选取方法保留的有效信号成分更多的优点,同时使用具有特定的改进软硬结合阈值处理方法对阈值选取方法进行了补偿,使去噪效果更为明显,最终在全频范围内既能最大限度的减小失真,又提高了去噪效果。
附图说明
图1为ATEM发射机发射电流PWM噪声的去噪方法流程图
图2为去噪之前采样得到的原始信号图
图3为去噪之后得到的梯形波信号图
具体实施方式
以下结合附图和实施例对本发明做进一步的详细描述:
ATEM发射机是一种发射双极性脉冲电流的发射装置,以低阻抗高感抗发射线圈为负载,产生大功率正负交变电流,激发出强烈的一次电磁场。但是由于ATEM发射机采用PWM脉冲宽度调制方式的控制策略,因此在测试发射波形时会产生电磁干扰,采样波形如图2所示。因此本文采用电感很低的采样电阻接在发射机的输出端,使用高精度的NI采集卡采集,在下降段可以看到有效的不含噪的信号波形,因此只需对上升段和平顶段做去噪即可。
基于Stein无偏似然估计阈值选取方法和软硬阈值结合阈值处理方法结合的ATEM发射机发射电流PWM噪声去噪方法,包括以下步骤:
步骤A、利用小波变换对待去噪电压信号进行多层小波分解,并得到各层的细节系数与尺度系数,包含以下两个步骤:
(1).本实施例选取使用采集卡采集出的数据中的第1至2019个数据,将其单独列入一个.mat文件中,采用sym6小波,对该组数据进行一维多尺度小波分解,得到细节系数数组CA3,及尺度系数数组CD1、CD2和CD3。
(2)将步骤一分解获得的四个数组合并成一个数组,其长度为2052。
步骤B、应用Stein无偏似然估计阈值选取方法对采样电压波形的特定的PWM噪声进行阈值选取。
Stein无偏似然估计阈值选取方法的原理是基于Stein无偏似然估计的一种自适应阈值选取方法,其在小波消噪的过程中,可以大约保留3%的小波系数,这对于在信号的高频信息有很少一部分在噪声范围内时,这种方法非常有用,可以将弱小的信号提取出来。
本实施例进行小波分解后,对于小波的阈值选取,采用Stein无偏似然估计阈值选取方法选取阈值,包括以下三个步骤:
(1)将步骤2所得的小波系数的平方按从小打到大的顺序排列。
p={p1,p2,…,pn}且p1≤p2≤…≤pn
(2)定义风险向量R,其元素为:
(3)以R元素中的最小值rmin作为风险值,找出rmin的下标变量min所对应的α,进而获得对应的pα,求出阈值T,即
其中,N为含噪信号在所有尺度上的小波分解系数的个数总和,σ为噪声信号的偏差。在本实施例中,最终求得的阈值T=0.3638.
步骤C、应用步骤B所选取的针对PWM噪声的阈值,对细节系数进行软硬阈值方法结合的去噪处理,然后重构得到不含噪的电压信号。
小波系数的阈值处理方法目前比较常用的有两种方法,分别是硬阈值处理与软阈值处理。其中硬阈值方程为:
软阈值方程为:
式中sgn(x)为符号函数;λ为阈值。
上述的两种与之处理方法各有其优缺点,硬阈值是将信号的绝对值与阈值进行比较,将小于或等于阈值的信号设为零,大于阈值的信号保持不变,它可以很好地保留信号的局部信息,但硬阈值的性质决定了它的不连续性,去噪后的信号仍有较明显的噪声。软阈值是将信号的绝对值与阈值进行比较,把绝对值小于或等于阈值的信号时设为零,对于绝对值大于阈值的信号,将其设为自身与阈值的差,这样信号就会向零收缩。软阈值的性质决定了它具有比硬阈值更好的连续性,但它的导数却是不连续的,估计小波系数与带噪小波系数存在恒定偏差,而且对大于阈值的系数进行定值压缩与噪声随着小波系数增大而减少的事实不符合,所以软阈值法的处理结果虽然相对平滑,但可能会使信号出现失真。
综合上述两种方法的优缺点,本发明提出采用如下的改进软硬阈值结合的方法进行阈值处理。
其中:γ是调节因子,介于0和1之间。在|wj,k|>λ时,可以通过调节α使新阈值函数在软硬阈值之间灵活选择,从而获得最优的去噪效果在本实施例中,由于要求补偿Stein无偏似然估计阈值选取方法的的缺点,因此本实施例选取α值为25。
在经过阈值处理后的小波系数中,大部分的小于阈值的小波系数被置零,大于阈值的小波系数被重新置换,获得了2052个处理完成的小波系数。
步骤D:应用步骤C得到的经过阈值处理后的小波系数,得到去噪后的信号的尺度系数与细节系数,然后进行多尺度一维重构,得到去噪后的信号序列。通过以上四个步骤实现了一种基于Stein无偏似然估计阈值选取方法和改进软硬阈值结合阈值处理方法结合的ATEM发射机发射电流PWM噪声去噪方法,成功去除了噪声并使失真度达到最小,如图3所示。
通过图2和图3对比可以看到:在图2中,由于受到环境因素的影响,在PWM噪声以外,还采集到了其他因素所产生的噪声,如在连接计算机时的USB接口方面的噪声,这种噪声在下降段可以看到的比较明显;及环境干扰所产生的噪声,本实施例所采集的环境为在闹市区的房间内,因此受到包括电磁信号等的干扰较大,在PWM噪声以外还有其他的噪声也占很大的成分,如在上升段与平顶段。图3为去噪结束后经过小波重构所产生的波形,可以看到经过去噪处理之后,基本呈现出了梯形波的变化趋势,包括上升段、平顶段、下降段及过冲。经对比,本发明所采用的方法能有效的去除PWM噪声并且失真很小。
Claims (1)
1.一种ATEM发射机发射电流PWM噪声的去噪方法,是基于Stein无偏似然估计阈值选取方法和改进软硬阈值结合阈值处理方法的结合,比单一的方法其去噪效果更为明显,括以下步骤:
A、利用小波变换对待去噪电压信号进行多层小波变换处理,并得到各层的细节系数与尺度系数;
B、应用Stein无偏似然估计阈值选取方法对待去噪电压波形的PWM噪声进行阈值选取;
C、应用步骤B所选取的针对PWM噪声的阈值,对细节系数与尺度系数进行改进的软、硬阈值结合的阈值处理;
D、重构得到去噪后的信号序列;
其特征在于:步骤B所述的阈值选取步骤为:
a、将步骤A所得的小波系数的平方按从小到大的顺序排列,即:
pα={p1,p2,...,pn},且p1≤p2≤…≤pn
b、定义风险向量R,其元素为:
c、以R元素中的最小值rmin作为风险值,找出rmin的下标变量min所对应的α,进而获得对应的pα,求出阈值T,即
其中,n为步骤A所求小波系数的个数,1≤i≤n,p1≤pi≤pn,p1≤pk≤pn,N为含噪信号在所有尺度上的小波分解系数的个数总和,σ为噪声信号的偏差;
步骤C所述的改进的软、硬阈值结合的阈值处理函数表达式为:
其中:Wj,k为步骤A所得到的小波系数,λ为步骤B所求的阈值T,γ是调节因子,介于0和1之间,在|wj,k|>λ时,通过调节α使新阈值函数在软硬阈值之间灵活选择,从而获得最优的去噪效果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510917176.8A CN105490516B (zh) | 2015-12-11 | 2015-12-11 | Atem发射机发射电流pwm噪声的去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510917176.8A CN105490516B (zh) | 2015-12-11 | 2015-12-11 | Atem发射机发射电流pwm噪声的去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105490516A CN105490516A (zh) | 2016-04-13 |
CN105490516B true CN105490516B (zh) | 2018-09-04 |
Family
ID=55677313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510917176.8A Expired - Fee Related CN105490516B (zh) | 2015-12-11 | 2015-12-11 | Atem发射机发射电流pwm噪声的去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105490516B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109035538A (zh) * | 2018-10-16 | 2018-12-18 | 深圳美特优科技有限公司 | 一种基于人脸识别的来访人员登记查验装置 |
CN112975980A (zh) * | 2021-03-09 | 2021-06-18 | 扬州哈工科创机器人研究院有限公司 | 小波变换在六轴机械臂去抖的应用 |
CN113074807A (zh) * | 2021-03-18 | 2021-07-06 | 中国水产科学研究院黄海水产研究所 | 养殖围栏设施结构振动及形变实时监测系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508076A (zh) * | 2011-11-09 | 2012-06-20 | 东北大学 | 一种基于多智能体系统和小波分析的故障诊断装置及方法 |
CN103163050A (zh) * | 2012-12-26 | 2013-06-19 | 大连理工大学 | 一种基于电磁感应信号的润滑油系统金属磨粒检测方法 |
CN104715461A (zh) * | 2015-04-02 | 2015-06-17 | 哈尔滨理工大学 | 图像去噪方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7508545B2 (en) * | 2004-09-27 | 2009-03-24 | Eastman Kodak Company | Color contour detection and correction |
US20100010780A1 (en) * | 2008-07-10 | 2010-01-14 | The Hong Kong Polytechnic University | Method for signal denoising using continuous wavelet transform |
-
2015
- 2015-12-11 CN CN201510917176.8A patent/CN105490516B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508076A (zh) * | 2011-11-09 | 2012-06-20 | 东北大学 | 一种基于多智能体系统和小波分析的故障诊断装置及方法 |
CN103163050A (zh) * | 2012-12-26 | 2013-06-19 | 大连理工大学 | 一种基于电磁感应信号的润滑油系统金属磨粒检测方法 |
CN104715461A (zh) * | 2015-04-02 | 2015-06-17 | 哈尔滨理工大学 | 图像去噪方法 |
Non-Patent Citations (2)
Title |
---|
图像噪声和小波阈值去噪;刘海霞;《内蒙古水利》;20060325;同上 * |
图像噪声和小波阈值去噪;刘海霞;《内蒙古水利》;20060325;第111-112页,正文第2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN105490516A (zh) | 2016-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ur Rehman et al. | Filter bank property of multivariate empirical mode decomposition | |
CN105490516B (zh) | Atem发射机发射电流pwm噪声的去噪方法 | |
DE102012206008A1 (de) | Reduzierung von Kopplungseffekten zwischen Spulenelementen einer Magnetresonanzspulenanordnung | |
CN105139371B (zh) | 一种基于pcnn与lp变换的多聚焦图像融合方法 | |
CN104618277B (zh) | 一种在窄带干扰下的短波ofdm信道估计方法 | |
CN102176312B (zh) | 一种通过小波陷波来降低突发噪音的系统及方法 | |
CN110081967A (zh) | 基于谱图小波变换的机械振动信号阈值降噪方法 | |
CN109359506A (zh) | 一种基于小波变换的心磁信号降噪方法 | |
Kumar et al. | Biosignal denoising via wavelet thresholds | |
CN105741305A (zh) | 一种基于平稳小波变换滤除肌电干扰的方法和系统 | |
CN108231084A (zh) | 一种基于Teager能量算子的改进小波阈值函数去噪方法 | |
CN112084845B (zh) | 基于多尺度小波系数自相关的低频1/f噪声消除方法 | |
CN107864020B (zh) | 水下小目标单分量声散射回波的变换域提取方法 | |
El Bouny et al. | ECG signal denoising based on ensemble emd thresholding and higher order statistics | |
CN109409281A (zh) | 一种基于改进的小波阈值函数的降噪方法 | |
CN109460614B (zh) | 基于瞬时带宽的信号时间-频率分解方法 | |
CN109724693A (zh) | 一种基于平稳小波的融合光谱去噪方法 | |
CN106097264B (zh) | 基于双树复小波和形态学的卫星遥测数据滤波方法 | |
Dotsinsky | An approach to successful power-line interference suppression in ECG signals | |
CN109586763B (zh) | 一种电子通信系统中扩频信号的去噪方法及其去噪系统 | |
CN107171665B (zh) | 带通信号的双通道tiadc非线性系统参数估计方法 | |
CN105662391A (zh) | 一种胃磁信号的特征提取与分类系统及方法 | |
Reddy et al. | Comparison of FIR and IIR Filters using ECG Signal with Different Sampling Frequencies | |
CN105763226B (zh) | 一种基于信噪比的最优小波基中压载波去噪方法 | |
CN109636729A (zh) | 一种基于小波变换的昆虫刺吸电位(epg)信号去噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180904 Termination date: 20191211 |
|
CF01 | Termination of patent right due to non-payment of annual fee |