CN105787444B - 基于v系统的信号去噪方法 - Google Patents
基于v系统的信号去噪方法 Download PDFInfo
- Publication number
- CN105787444B CN105787444B CN201610099844.5A CN201610099844A CN105787444B CN 105787444 B CN105787444 B CN 105787444B CN 201610099844 A CN201610099844 A CN 201610099844A CN 105787444 B CN105787444 B CN 105787444B
- Authority
- CN
- China
- Prior art keywords
- signal
- resampling
- fourier
- length
- denoised
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000012952 Resampling Methods 0.000 claims abstract description 35
- 230000009466 transformation Effects 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 230000011218 segmentation Effects 0.000 claims description 15
- 241001463014 Chazara briseis Species 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract 1
- 230000008929 regeneration Effects 0.000 abstract 1
- 238000011069 regeneration method Methods 0.000 abstract 1
- 238000009987 spinning Methods 0.000 description 6
- 101150065184 sym-2 gene Proteins 0.000 description 4
- 101150004636 SYM1 gene Proteins 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及基于V系统的信号去噪方法,其包括:首先,对含噪信号进行小波变换去噪后,计算该含噪信号的间断点位置;然后,根据该间断点位置对原含噪信号进行分段,对每段信号进行插值重采样,使间断点刚好落在V系统的分段处;最后,对每段信号进行V系统的傅立叶-V级数展开,保留满足阈值要求的系数进行信号重构,对重构后的信号重采样后得到去噪信号。根据本发明方法不仅能够有效地提高了信号检测的鲁棒性,还能够充分利用了V系统的再生性质,得到很好的去噪效果,从而极大地减少了重构信号过程中出现的Gibbs现象。
Description
技术领域
本发明涉及基于V系统的信号去噪方法,具体涉及采用小波变换检测信号间断点并对信号按间断点重采样后进行进行V系统重构信号的去噪方法。
背景技术
信号的产生及测量过程中都不可避免地引入噪声,干扰了原始信号的分析与处理,因此涌现出很多信号去噪的方法。常见的信号去噪方法主要采用正交变换的有限项级数逼近信号的原理进行去噪,常用的正交变换有傅立叶变换、正弦/余弦变换和小波变换,这些变换大都是用连续的正交基来逼近信号,在信号的奇异点处误差较大,易导致Gibbs振荡现象的发生。
发明内容
鉴于现有技术的上述不足,本发明提供一种基于V系统的信号去噪方法,其目的在于利用V系统的再生性对信号进行去噪。为了实现上述目的,本发明采用如下技术方案:
基于V系统的信号去噪方法,其包括:
步骤一:对含噪信号进行间断点检测;
步骤二:按间断点分割对信号分段后对分段信号进行重采样;
步骤三:对重采样后的分段信号进行V系统展开去噪后重构信号。
步骤一进一步包括:
步骤11对长度为N的含噪信号S(t),t=1,2,...,N进行小波变换,然后采用Heursure阈值法对含噪信号S(t)进行初步去噪分析得到恢复信号S'(t);
步骤12对恢复信号S'(t)进行一层小波变换,使用小波变换后的高频部分重构信号得到信号S'(t)的细节信号D1(t),t=1,2,...,N,找出D1(t)的最大值记为maxd1,记录D1(t)信号中大于t1*maxd1阈值的位置,将得到的所有间断点位置存储在间断位置x(i),i=1,2,...,M,其中t1为经验阈值,在此取0.2,M为检测到的间断点个数;
步骤13按间断位置x(i),i=1,2,...,M,将含噪信号S(t)分割成M+1段,找出M+1段信号长度的最大值记为maxL,然后计算大于等于maxL的最接近的能表示成2n的整数Ns=2n。
步骤二是指,将M+1段信号都采用分段三次Hermit多项式插值算法重采样成长度为Ns=2n点的重采样信号,该重采样信号记为Si(t),t=1,2,...,Ns,并记录各分段段信号重采样前后的对应位置信息表Pi(t),t=1,2,...,Ni,其中Ni为第i段信号的长度。
步骤三是指,对重采样的分段信号按V系统展开成傅立叶-V级数,接着对每个分段信号的傅立叶-V级数进行阈值去噪后重构分段信号拼接生成最后的去噪信号。
步骤三进一步包括:
步骤31对重采样信号Si(t)进行V系统的傅立叶-V级数展开,得到傅立叶-V级数的系数,a(k),k=1,2,...,Ns
其中Vk(t)为V系统第k个基函数的第t个值;
步骤32保留大于阈值t2的傅立叶-V级数的系数a(k)的值,其它系数的值为0,其中阈值σ为估计信号噪声方差,Ns为信号的长度;
步骤33采用V系统逼近信号得到S′i(t)
步骤34根据各分段信号位置信息表Pi(t)重构对应的去噪信号,并将各分段去噪信号拼合得到最后的去噪信号记为S′i(t)。
本发明的有益效果:本发明方法先采用基于Heursure阈值法的小波变换信号去噪方法对信号进行去噪,然后再对去噪后的信号进行小波变换,去除低频系数保留高频系数,从而有效避免噪声信号对间断点检测的干扰。与基于小波循环平移不变性的信号去噪方法相比,本发明的去噪信号的信噪比(SNR)有显著提升,极大地减少了重构信号过程中出现的Gibbs现象。
附图说明
图1是本发明方法的算法流程图;
图2a是根据本发明方法的一个具体实施例中长度为1024的原始Block信号;
图2b是图2a的具体实施例中长度为1024含高斯白噪声的Block信号,其中噪声信号与原信号的信噪比SNR为22.22dB;
图2c示出了图2b的含噪Block信号中的间断点位置;
图2d是根据图2c中的间断点对含噪Block信号分段并对分段信号重采样后的信号;
图2e对图2e示出的重采样信号基于V系统重构后得到的去噪信号;
图2f是根据图2e映射后的去噪信号,其与原Block信号的信噪比SNR为44.37dB;
图2g是基于sym1小波循环平移不变去噪方法得到的去噪信号,与原Block信号的信噪比SNR为35.98dB;
图3a是长度为1024的原始Skyline信号;
图3b是长度为1024含高斯白噪声Skyline信号与原始Skyline信号的信噪比SNR为23.92dB;
图3c是V系统去噪后的Skyline信号,与原始Skyline信号的信噪比SNR为40.47dB;
图3d是基于sym2小波循环平移不变去噪方法得到的去噪信号,与原始Skyline信号的信噪比SNR为30.37dB;
图4a是长度为1024的Piecewise1原始信号;
图4b是是长度为1024含高斯白噪声Piecewise1信号,与原始Piecewise1信号的信噪比SNR=20.39dB;
图4c是基于V系统的去噪信号,与原始Piecewise1信号的信噪比SNR为37.62dB;
图4d是基于sym2小波循环平移不变去噪方法得到的去噪信号,与原始Piecewise1信号的信噪比SNR为28.85dB。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明揭示一种基于V系统的信号去噪方法,下面以方波信号为例,本发明方法包括如下步骤:
步骤一:对含噪信号进行间断点检测;
步骤二:按间断点分割对信号分段后对分段信号进行重采样;
步骤三:对重采样后的分段信号进行V系统展开去噪后重构信号。
其中,对于步骤一,为了避免噪声信号对间断点检测的干扰,先采用基于Heursure阈值法的小波变换信号去噪方法对信号进行去噪,然后再对去噪后的信号进行小波变换,去除低频系数保留高频系数后进行小波逆变换得到带有间断信息的细节信号,然后采用阈值法检测出信号的间断点。
对于步骤二,对信号分段后,对分段信号采用插值方式的重采样,是每段信号都有统一的长度并满足V系统的再生性要求,并记录分段重采样前后每个点的位置变化关系,以备步骤三所用。
进一步地,步骤三是指,对重采样的分段信号按V系统展开成傅立叶-V级数,接着对每个分段信号的傅立叶-V级数进行阈值去噪后重构分段信号拼接生成最后的去噪信号。
在一个具体实施例中,步骤三中的V系统是一类由正交分段多项式构成的完备正交系,再对信号进行傅立叶-V级数展开时对不同信号采用的正交分段多项式次数不同,以适应信号的光滑性要求。对分段信号的傅立叶-V级数的系数进行阈值化处理,保留大于等于阈值的系数,其它系数的值置为零,然后用保留的系数重构信号得到去噪后的分段信号,最后采用步骤二中保存的信号点的位置变化关系得到重采样前位置的去噪信号,将分段的去噪信号拼接成完整的去噪信号。
更具体地,步骤一进一步包括:
步骤11对长度为N的含噪信号S(t),t=1,2,...,N进行小波变换,然后,采用Heursure阈值法对含噪信号S(t)进行初步去噪分析得到恢复信号S'(t);
步骤12对恢复信号S'(t)进行一层小波变换,使用小波变换后的高频部分重构信号得到信号S'(t)的细节信号D1(t),t=1,2,...,N,找出D1(t)的最大值记为maxd1,记录D1(t)信号中大于t1*maxd1阈值的位置,将得到的所有间断点位置存储在间断位置x(i),i=1,2,...,M,其中t1为经验阈值,在此取0.2,M为检测到的间断点个数;
步骤13按间断位置x(i),i=1,2,...,M,将含噪信号S(t)分割成M+1段,找出M+1段信号长度的最大值记为maxL,然后计算大于等于maxL的最接近的能表示成2n的整数Ns=2n。
具体地,步骤二是指,将M+1段信号都采用分段三次Hermit多项式插值算法重采样成长度为Ns=2n点的重采样信号,该重采样信号记为Si(t),t=1,2,...,Ns,并记录各分段段信号重采样前后的对应位置信息表Pi(t),t=1,2,...,Ni,其中Ni为第i段信号的长度。
具体地,步骤三是指,对重采样的分段信号按V系统展开成傅立叶-V级数,接着对每个分段信号的傅立叶-V级数进行阈值去噪后重构分段信号拼接生成最后的去噪信号。
步骤三进一步包括:
步骤31对重采样信号Si(t)进行V系统的傅立叶-V级数展开,得到傅立叶-V级数的系数a(k),k=1,2,...,Ns,
其中Vk(t)为V系统第k个基函数的第t个值;
步骤32保留大于阈值t2的傅立叶-V级数的系数a(k)的值,其它系数的值为0,其中阈值σ为估计信号噪声方差,Ns为信号的长度;
步骤33采用V系统逼近信号得到S′i(t)
步骤34根据各分段信号位置信息表Pi(t)重构对应的去噪信号,并将各分段去噪信号拼合得到最后的去噪信号记为S′i(t)。
图2a-2f是根据本发明的信号去噪算法的一个具体实施例。其中,图2a是长度为1024的原始Block信号;图2b是长度为1024含高斯白噪声Block信号,噪声信号与原信号的信噪比SNR为22.22dB;图2c中的上半部分图是检测到含噪Block信号中的间断点位置;图2d是对分段信号重采样后的信号;图2e是基于V系统重构后的去噪信号;图2f是根据图2d信号重采样映射后的去噪信号,与原Block信号的信噪比SNR为44.37dB;图2g是基于sym1小波循环平移不变去噪方法得到的去噪信号,与原Block信号的信噪比SNR为35.98dB。
根据本发明方法的具体实施方式中,采用小波变换进行检测的间断点位置如图2c所示,其中间断点个数为15,分割的区域为16。图2d下半部信号是采用分段三次Hermit多项式进行重采样,得到每个区域采样点数目均为64的信号。
图3a-3d是根据本发明方法对Skyline信号进行去噪的又一个具体实施例。在该实施例中,原始信号是长度为1024Skyline信号如图3a所示。往原始信号Skyline中加入噪声得到含噪信号,图3b是长度为1024含高斯白噪声信号的含噪信号,其中高斯白噪声Skyline信号与原始Skyline信号的信噪比SNR为23.92dB。图3c是V系统去噪后的Skyline信号,与原始Skyline信号的信噪比SNR为40.47dB。图3d是基于sym2小波循环平移不变去噪方法得到的去噪信号,与原始Skyline信号的信噪比SNR为30.37dB;
图4a-4d示出了根据本发明方法对Piecewise1信号的去噪又一个具体实施例。其中,图4a是长度为1024的Piecewise1原始信号;图4b是是长度为1024含高斯白噪声Piecewise1信号,与原始Piecewise1信号的信噪比SNR=20.39dB;图4c是基于V系统的去噪信号,与原始Piecewise1信号的信噪比SNR为37.62dB;图4d是基于sym2小波循环平移不变去噪方法得到的去噪信号,与原始Piecewise1信号的信噪比SNR为28.85dB。
以上是对本发明的描述,在不脱离本发明设计精神的前提下,本领域普通工程技术人员对本发明的技术方案作出的各种变形和替换,均应落入本发明的权利要求确定的保护范围内。
Claims (4)
1.基于V系统的信号去噪方法,其特征在于,包括:
步骤一:对含噪信号进行间断点检测;
步骤二:按间断点分割对信号分段后对分段信号进行重采样;
步骤三:对重采样后的分段信号进行V系统展开去噪后重构信号;
其中,步骤一进一步包括:
步骤11对长度为N的含噪信号S(t),t=1,2,...,N进行小波变换,然后采用Heursure阈值法对含噪信号S(t)进行初步去噪分析得到恢复信号S'(t);
步骤12对恢复信号S'(t)进行一层小波变换,使用小波变换后的高频部分重构信号得到信号S'(t)的细节信号D1(t),t=1,2,...,N,找出D1(t)的最大值记为maxd1,记录D1(t)信号中大于t1*maxd1阈值的位置,将得到的所有间断点位置存储在间断位置x(i),i=1,2,...,M,其中t1为经验阈值,在此取0.2,M为检测到的间断点个数;
步骤13按间断位置x(i),i=1,2,...,M,将含噪信号S(t)分割成M+1段,找出M+1段信号长度的最大值记为maxL,然后计算大于等于maxL的最接近的能表示成2n的整数Ns=2n。
2.根据权利要求1所述的基于V系统的信号去噪方法,其特征在于,步骤二是指,将M+1段信号都采用分段三次Hermit多项式插值算法重采样成长度为Ns=2n点的重采样信号,该重采样信号记为Si(t),t=1,2,...,Ns,并记录各分段段信号重采样前后的对应位置信息表Pi(t),t=1,2,...,Ni,其中Ni为第i段信号的长度。
3.根据权利要求1所述的基于V系统的信号去噪方法,其特征在于,步骤三是指,对重采样的分段信号按V系统展开成傅立叶-V级数,接着对每个分段信号的傅立叶-V级数进行阈值去噪后重构分段信号拼接生成最后的去噪信号。
4.根据权利要求1所述的基于V系统的信号去噪方法,其特征在于,步骤三进一步包括:
步骤31对重采样信号Si(t)进行V系统的傅立叶-V级数展开,得到傅立叶-V级数的系数,a(k),k=1,2,...,Ns
其中Vk(t)为V系统第k个基函数的第t个值;
步骤32保留大于阈值t2的傅立叶-V级数的系数a(k)的值,其它系数的值为0,其中阈值σ为估计信号噪声方差,Ns为信号的长度;
步骤33采用V系统逼近信号得到S′i(t)
步骤34根据各分段信号位置信息表Pi(t)重构对应的去噪信号,并将各分段去噪信号拼合得到最后的去噪信号记为S′i(t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610099844.5A CN105787444B (zh) | 2016-02-24 | 2016-02-24 | 基于v系统的信号去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610099844.5A CN105787444B (zh) | 2016-02-24 | 2016-02-24 | 基于v系统的信号去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105787444A CN105787444A (zh) | 2016-07-20 |
CN105787444B true CN105787444B (zh) | 2019-03-22 |
Family
ID=56402866
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610099844.5A Expired - Fee Related CN105787444B (zh) | 2016-02-24 | 2016-02-24 | 基于v系统的信号去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105787444B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106296609B (zh) * | 2016-08-09 | 2023-09-12 | 广东盈动高科自动化有限公司 | 用于容栅的全波和数字滤波去噪方法和电路 |
CN109029765A (zh) * | 2018-05-24 | 2018-12-18 | 北京石油化工学院 | 一种信号处理方法和信号处理装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005002200A2 (en) * | 2003-06-13 | 2005-01-06 | Nielsen Media Research, Inc. | Methods and apparatus for embedding watermarks |
CN102664013A (zh) * | 2012-04-18 | 2012-09-12 | 南京邮电大学 | 基于能量选择的离散余弦变换域的音频数字水印方法 |
CN102944252B (zh) * | 2012-11-28 | 2014-12-31 | 重庆邮电大学 | 一种基于平移不变量小波的fbg信号处理方法 |
CN104810022B (zh) * | 2015-05-11 | 2018-06-15 | 东北师范大学 | 一种基于音频断点的时域数字音频水印方法 |
-
2016
- 2016-02-24 CN CN201610099844.5A patent/CN105787444B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105787444A (zh) | 2016-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109767783B (zh) | 语音增强方法、装置、设备及存储介质 | |
WO2017202292A1 (zh) | 一种回声时延跟踪方法及装置 | |
CN107068161B (zh) | 基于人工智能的语音降噪方法、装置和计算机设备 | |
CN107169979A (zh) | 一种改进Canny算子的图像边缘检测方法 | |
CN107863099B (zh) | 一种新型双麦克风语音检测和增强方法 | |
CN103269412A (zh) | 一种视频图像的降噪方法及装置 | |
CN105787444B (zh) | 基于v系统的信号去噪方法 | |
CN108647597B (zh) | 一种手腕识别方法、手势识别方法、装置和电子设备 | |
CN108231084A (zh) | 一种基于Teager能量算子的改进小波阈值函数去噪方法 | |
CN110702149A (zh) | 光纤分布式扰动传感系统扰动信号的定位方法和系统 | |
CN113553997A (zh) | 一种联合改进小波阈值的信号处理方法、存储介质和系统 | |
CN109102818B (zh) | 一种基于信号频率概率密度函数分布的去噪音频采样算法 | |
TWI597720B (zh) | 嬰兒哭聲偵測電路及相關的偵測方法 | |
CN117708747A (zh) | 一种列车运行中轨道振动数据异常的监测方法 | |
CN106997768B (zh) | 一种语音出现概率的计算方法、装置及电子设备 | |
CN110755069B (zh) | 一种跳跃突变噪声的动态心电信号基线漂移校正方法 | |
CN115994870B (zh) | 一种用于增强去噪的图像处理方法 | |
CN116659646A (zh) | 一种基于机器视觉的风机叶片振动检测方法及装置 | |
CN115293219A (zh) | 一种融合小波和峭度的脉冲信号去噪方法 | |
CN111723677B (zh) | 一种基于自适应阈值的小波去噪方法 | |
CN109188370B (zh) | 一种雷达设备lfm脉冲信号包络曲线拟合方法及系统 | |
CN105513587B (zh) | 一种mfcc提取方法及装置 | |
CN111652040A (zh) | 一种基于射频及信号自适应卷积神经网络进行人体活动识别的方法与系统 | |
CN106197480A (zh) | 一种低信噪比信号的处理系统 | |
Xie et al. | An adaptive median filter using local texture information in images |
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: 20190322 Termination date: 20210224 |
|
CF01 | Termination of patent right due to non-payment of annual fee |