CN101246469B - 一种对数字信号使用dft理想滤波器的滤波方法 - Google Patents
一种对数字信号使用dft理想滤波器的滤波方法 Download PDFInfo
- Publication number
- CN101246469B CN101246469B CN2007100639481A CN200710063948A CN101246469B CN 101246469 B CN101246469 B CN 101246469B CN 2007100639481 A CN2007100639481 A CN 2007100639481A CN 200710063948 A CN200710063948 A CN 200710063948A CN 101246469 B CN101246469 B CN 101246469B
- Authority
- CN
- China
- Prior art keywords
- filtering
- dft
- filter
- ideal
- digital signal
- 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
Landscapes
- Complex Calculations (AREA)
Abstract
本发明涉及数字信号处理,具体涉及对数字信号使用DFT理想滤波器滤波。包括先对原始信号记录道进行希尔伯特变换,计算瞬时包络,计算相位函数,对相位函数进行DFT理想滤波器滤波,最后计算输出函数。使用本方法滤波,可以使滤波性能更佳,同时基本上没有吉普斯效应。本发明用于数字信号滤波领域。
Description
技术领域
本发明涉及数字信号处理,具体涉及对数字信号使用DFT理想滤波器滤波。
背景技术
在数字信号处理领域中,经常需要将特定频率或窄频带信号从宽频带信号中分解出来,以突出目标信号的效果;或者将特定频率带信号从宽频带信号中去除。离散傅立叶变换(DFT)的特定频率分量是无限长的余弦或正弦波列,使用DFT理想滤波器会产生吉普斯效应,因此,在实际数字信号处理中使用DFT理想滤波器一直被认为是不切实际的。然而,要克服吉普斯效应,DFT实际滤波器需有足够的镶边,这就使得分频结果内掺杂了过多不需要的频率分量。
在石油工业出版社1999年12月出版的R.E.谢里夫所著的《勘探地震学》(第二版)第704-705页,描述了常规DFT滤波及理想滤波器的问题。在该书第407-408页描述了希尔伯特Hilbert变换的基本应用。
发明内容
本发明要解决的技术问题是:
针对现有技术的不足,本发明的目的是提供一种对数字信号使用DFT理想滤波器的滤波方法,能够使宽频带信号的滤波性能更佳,信噪比高,基本无吉普斯效应。
本发明的技术方案是:
一种对数字信号使用DFT理想滤波器的滤波方法,原始数据为:一道时间序列数字信号,记为xr(t);
对xr(t)进行滤波,包括但不限于下述步骤:
第一步:对xr(t)进行希尔伯特变换,得到xh(t);
希尔伯特变换是公知的常规方法,本领域的技术人员都可以准确无误地从xr(t)惟一得到xh(t)。
第二步,分为两个独立的分步骤:第二A步、第二B步,这两个分步骤之间没有先后顺序,无论先完成哪一步都不影响本技术方案的实施。
所述的第二A步包括:计算瞬时包络a(t):
所述的第二B步包括两小步,
第二B步第1小步是:计算相位函数cosθ(t),
当先进行第二B步时,(2)式中的a(t)直接用(1)式右边代入,可以在不计算瞬时包络的情况下,直接计算出相位函数。
对原始信号进行希尔伯特变换、计算瞬时包络、计算相位函数,都是公知的常规方法,在三瞬分析中经常使用,其变换、计算方法在公知教科书中都已揭示。
第二B步第2小步是:对相位函数cosθ(t)进行DFT理想滤波器滤波,得到滤波后的相位函数yθ(t);
DFT理想滤波器滤波本身是公知方法,可以有:理想带通滤波器、理想高通滤波器、理想低通滤波器、理想带阻滤波器。
而对相位函数进行DFT理想滤波器滤波,是本发明的创新。
第三步:计算输出函数yr(t),
yr(t)=yθ(t)·a(t) (3)
也即把两个离散序列逐点相乘,得到作为结果的离散序列。
第三步也是本发明的创新。
通过上述步骤,得到一道时间序列的经过滤波的数字信号yr(t)。
完成滤波过程。
除了上述基本步骤之外,在实际应用时还可以增加步骤:
在所述的第一步,进行希尔伯特变换之前,先对xr(t)进行预处理。
并且还可以:
在所述的第二A步,计算瞬时包络之后,再对a(t)进行修饰处理。
例如:
所述的对xr(t)进行预处理,包括下列处理方式中的一种或一种以上:
剔野值,均衡,常规滤波;
上述的常规滤波包括:常规带通滤波、常规高通滤波、常规低通滤波、常规带阻滤波、小波变换滤波、平滑滤波;
上述的平滑滤波包括:波数滤波、中值滤波、积分滤波。
所述的对a(t)进行修饰处理,包括下列处理方式中的一种或一种以上:
中值滤波、积分滤波、波数滤波。
这些预处理、修饰处理方式,都是滤波中常用的公知技术手段,在本发明的技术方案中使用这些技术手段,可以针对不同具体应用情况,更好地提高滤波效果。
例如可以将上述手段组合使用:
在所述的第一步,进行希尔伯特变换之前,先对xr(t)进行预处理,所述的预处理包括下列处理方式中的一种或一种以上:剔野值、均衡、常规带通滤波、常规高通滤波、常规低通滤波、常规带阻滤波、小波变换滤波、波数滤波、中值滤波、积分滤波;
在所述的第二A步,计算瞬时包络之后,再对a(t)进行修饰处理,所述的修饰处理包括下列处理方式中的一种或一种以上:中值滤波、积分滤波、波数滤波;
在所述的第二B步第2小步,所述的DFT理想滤波器是下列之一:理想带通滤波器、理想高通滤波器、理想低通滤波器、理想带阻滤波器。
本发明的有益效果是:
使滤波性能更佳,可以把不需要的频段信号尽量抑制。同时基本上没有吉普斯效应。
附图说明
图1是本发明基本步骤示意图。
图2是本发明扩展步骤示意图。
图3是一个实际记录道xr(t)的示意图。
图4是图3所示实际记录道的频谱图。
图5是理想带通滤波器对xr(t)实施滤波后的10.5Hz~11.4Hz记录道的示意图。
图6是xr(t)经Hilbter变换后获得的相位函数cosθ(t)的示意图。
图7是理想带通滤波器对cosθ(t)实施滤波得到的10.5Hz~11.4Hz滤波后的相位函数yθ(t)示意图。
图8是实信号的瞬时包络示意图。
图9是xr(t)的原始瞬时包络a(t)示意图。
图10是基本无吉普斯效应的10.5Hz~11.4Hz滤波输出记录道示意图。
图11是图10记录道的频谱示意图。
图12是有3Hz镶边的常规带通滤波器对xr(t)的滤波输出示意图。
图13是图12的频谱示意图。
具体实施方式
下面结合实施例进一步描述本发明。本发明的范围不受这些实施例的限制,本发明的范围在权利要求书中提出。
一道数字信号的实际记录道,记为xr(t)。见图3。该实际记录道的频谱见图4。
如果对该实际记录道,直接用带通频率为10.5Hz~11.4Hz的理想带通滤波器进行滤波,得到图5。可以看出,有严重的吉普斯效应。
使用本发明,先对该实际记录道进行希尔伯特变换,计算瞬时包络;
并计算该实际记录道的相位函数cosθ(t),见图6。
然后对相位函数用理想带通滤波器进行滤波,带通频率为10.5Hz~11.4Hz,滤波后的相位函数yθ(t)见图7。
该实际记录道的原始瞬时包络a(t)见图9。
最后计算输出函数,见图10,得到基本无吉普斯效应的10.5Hz~11.4Hz记录道,其频谱见图11。
作为对照,使用现有技术,带有3Hz镶边的常规带通滤波器对xr(t)的滤波输出,见图12,其频谱见图13。
从图11与图13对比可以看出,使用本发明方法进行滤波后的频谱,明显比现有技术的频谱更尖锐,即杂波更少。
Claims (5)
1.一种对数字信号使用离散傅立叶变换DFT理想滤波器的滤波方法,原始数据为:一道时间序列数字信号,记为xr(t);
对xr(t)进行滤波,包括但不限于下述步骤:
第一步:先对xr(t)进行预处理;
对xr(t)进行希尔伯特变换,得到xh(t);
第二步,分为两个独立的分步骤:第二A步、第二B步,这两个分步骤之间没有先后顺序,
所述的第二A步包括:计算瞬时包络a(t):
再对a(t)进行修饰处理;
所述的第二B步包括两小步,
第二B步第1小步是:计算相位函数cosθ(t),
第二B步第2小步是:对相位函数cosθ(t)进行DFT理想滤波器滤波,得到滤波后的相位函数yθ(t);
第三步:计算输出函数yr(t),
yr(t)=yθ(t)·a(t) (3)
通过上述步骤,得到一道时间序列的经过滤波的数字信号yr(t)。
2.根据权利要求1所述的一种对数字信号使用离散傅立叶变换DFT理想滤波器的滤波方法,其特征是:
所述的对xr(t)进行预处理,包括下列处理方式中的一种以上:
剔野值,均衡,常规滤波;
上述的常规滤波包括:常规带通滤波、常规高通滤波、常规低通滤波、常规带阻滤波、小波变换滤波、平滑滤波;
上述的平滑滤波包括:波数滤波、中值滤波、积分滤波。
3.根据权利要求1所述的一种对数字信号使用离散傅立叶变换DFT理想滤波器的滤波方法,其特征是:
所述的对a(t)进行修饰处理,包括下列处理方式中的一种以上:
中值滤波、积分滤波、波数滤波。
4.根据权利要求1至3之一所述的一种对数字信号使用离散傅立叶变换DFT理想滤波器的滤波方法,其特征是:
在所述的第二B步第2小步,所述的离散傅立叶变换DFT理想滤波器是下列之一:
理想带通滤波器、理想高通滤波器、理想低通滤波器、理想带阻滤波器。
5.根据权利要求1所述的一种对数字信号使用离散傅立叶变换DFT理想滤波器的滤波方法,其特征是:
在所述的第一步,进行希尔伯特变换之前,先对xr(t)进行预处理,所述的预处理包括下列处理方式中的一种以上:剔野值、均衡、常规带通滤波、常规高通滤波、常规低通滤波、常规带阻滤波、小波变换滤波、波数滤波、中值滤波、积分滤波;
在所述的第二A步,计算瞬时包络之后,再对a(t)进行修饰处理,所述的修饰处理包括下列处理方式中的一种以上:中值滤波、积分滤波、波数滤波;
在所述的第二B步第2小步,所述的离散傅立叶变换DFT理想滤波器是下列之一:理想带通滤波器、理想高通滤波器、理想低通滤波器、理想带阻滤波器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100639481A CN101246469B (zh) | 2007-02-15 | 2007-02-15 | 一种对数字信号使用dft理想滤波器的滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007100639481A CN101246469B (zh) | 2007-02-15 | 2007-02-15 | 一种对数字信号使用dft理想滤波器的滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101246469A CN101246469A (zh) | 2008-08-20 |
CN101246469B true CN101246469B (zh) | 2010-09-08 |
Family
ID=39946929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007100639481A Active CN101246469B (zh) | 2007-02-15 | 2007-02-15 | 一种对数字信号使用dft理想滤波器的滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101246469B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043168A (zh) * | 2010-10-15 | 2011-05-04 | 中国石油化工股份有限公司 | 一种对数字信号进行仿真加噪的处理方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102053265B (zh) * | 2009-10-30 | 2014-11-26 | 中国石油化工股份有限公司 | 一种对地震资料剔除野值的滤波方法 |
CN103792577B (zh) * | 2012-10-26 | 2016-07-06 | 中国石油化工股份有限公司 | 一种消除伪频谱的频谱分析方法 |
CN105022917B (zh) * | 2015-07-06 | 2018-02-16 | 西安理工大学 | 一种信号精确提取与处理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5831934A (en) * | 1995-09-28 | 1998-11-03 | Gill; Stephen P. | Signal processing method for improved acoustic formation logging system |
US5917738A (en) * | 1996-11-08 | 1999-06-29 | Pan; Cheh | Removing the gibbs phenomenon in fourier transform processing in digital filters or other spectral resolution devices |
CN1713098A (zh) * | 2004-06-24 | 2005-12-28 | C.R.F.阿西安尼顾问公司 | 控制工业生产过程质量的方法及其*** |
-
2007
- 2007-02-15 CN CN2007100639481A patent/CN101246469B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5831934A (en) * | 1995-09-28 | 1998-11-03 | Gill; Stephen P. | Signal processing method for improved acoustic formation logging system |
US5917738A (en) * | 1996-11-08 | 1999-06-29 | Pan; Cheh | Removing the gibbs phenomenon in fourier transform processing in digital filters or other spectral resolution devices |
CN1713098A (zh) * | 2004-06-24 | 2005-12-28 | C.R.F.阿西安尼顾问公司 | 控制工业生产过程质量的方法及其*** |
Non-Patent Citations (1)
Title |
---|
王琼等.消除Gibbs现象的仿真研究.模式识别与仿真技术.2004,23(8),21-24. * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102043168A (zh) * | 2010-10-15 | 2011-05-04 | 中国石油化工股份有限公司 | 一种对数字信号进行仿真加噪的处理方法 |
CN102043168B (zh) * | 2010-10-15 | 2012-11-07 | 中国石油化工股份有限公司 | 一种对数字信号进行仿真加噪的处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101246469A (zh) | 2008-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101246469B (zh) | 一种对数字信号使用dft理想滤波器的滤波方法 | |
Iacobucci | Spectral analysis for economic time series | |
CN102546499B (zh) | 一种实线性调频信号的分数阶信道化接收方法 | |
CN106199532B (zh) | 基于混合傅立叶-小波分析的探地雷达信号降噪方法 | |
RU2009107093A (ru) | Устройство и способ для обработки действительного сигнала поддиапазона для ослабления эффектов наложения спектров | |
US9949029B2 (en) | Audio filtering with virtual sample rate increases | |
Sharma et al. | An empirical wavelet transform-based approach for cross-terms-free Wigner–Ville distribution | |
CN105044769B (zh) | 提高地震信号的分辨率的方法 | |
CN101881821B (zh) | 一种分数阶傅里叶域信道化接收方法 | |
CN109239780A (zh) | 基于同步挤压小波变换去除面波的方法 | |
RU2012148250A (ru) | Представление информационного сигнала с использованием преобразования с перекрытием | |
Wu et al. | Gearbox fault diagnosis using adaptive zero phase time-varying filter based on multi-scale chirplet sparse signal decomposition | |
CN106125134A (zh) | 基于双曲时窗的地震数据信噪比计算方法 | |
CN109813660A (zh) | 一种分段式的矿物光谱噪声处理方法 | |
CN100349007C (zh) | 消除地震记录信号中单频干扰的方法 | |
Lobov et al. | Digital compensation for uneven frequency response of analog filters from the hybrid filter bank | |
CN108181486B (zh) | 加速度信号的处理方法及装置 | |
CN102043168B (zh) | 一种对数字信号进行仿真加噪的处理方法 | |
CN102721663A (zh) | 一种基于自适应滤波的近红外土壤光谱去噪方法 | |
CN109975867B (zh) | 一种频率域信号混叠的地震数据拓频方法 | |
CN106027454A (zh) | 基于频率补偿的qam信号的频偏估计方法 | |
CN104836547A (zh) | 一种短群延时数字滤波方法 | |
CN106303879A (zh) | 一种基于时域分析的检测装置及检测方法 | |
RU2691122C1 (ru) | Способ и устройство компандирования звуковых вещательных сигналов | |
Rakshit et al. | Fir filter design using an adjustable novel window and its applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |