一种对数字信号使用ＤＦＴ理想滤波器的滤波方法技术

本发明专利技术涉及数字信号处理，具体涉及对数字信号使用ＤＦＴ理想滤波器滤波。包括先对原始信号记录道进行希尔伯特变换，计算瞬时包络，计算相位函数，对相位函数进行ＤＦＴ理想滤波器滤波，最后计算输出函数。使用本方法滤波，可以使滤波性能更佳，同时基本上没有吉普斯效应。本发明专利技术用于数字信号滤波领域。

【技术实现步骤摘要】

本专利技术涉及数字信号处理，具体涉及对数字信号使用DFT理想滤波器滤波。
技术介绍
在数字信号处理领域中，经常需要将特定频率或窄频带信号从宽频带信号中分解出来，以突出目标信号的效果；或者将特定频率带信号从宽频带信号中去除。离散傅立叶变换(DFT)的特定频率分量是无限长的余弦或正弦波列，使用DFT理想滤波器会产生吉普斯效应，因此，在实际数字信号处理中使用DFT理想滤波器一直被认为是不切实际的。然而，要克服吉普斯效应，DFT实际滤波器需有足够的镶边，这就使得分频结果内掺杂了过多不需要的频率分量。在石油工业出版社1999年12月出版的R.E.谢里夫所著的《勘探地震学》(第二版)第704-705页，描述了常规DFT滤波及理想滤波器的问题。在该书第407-408页描述了希尔伯特Hilbert变换的基本应用。
技术实现思路
本专利技术要解决的技术问题是：针对现有技术的不足，本专利技术的目的是提供一种对数字信号使用DFT理想滤波器的滤波方法，能够使宽频带信号的滤波性能更佳，信噪比高，基本无吉普斯效应。-->本专利技术的技术方案是：一种对数字信号使用DFT理想滤波器的滤波方法，原始数据为：一道时间序列本文档来自技高网...

【技术保护点】
一种对数字信号使用ＤＦＴ理想滤波器的滤波方法，原始数据为：一道时间序列数字信号，记为ｘ↓［ｒ］（ｔ）；对ｘ↓［ｒ］（ｔ）进行滤波，包括但不限于下述步骤：第一步：对ｘ↓［ｒ］（ｔ）进行希尔伯特变换，得到ｘ↓［ｈ］（ｔ）；第二步，分为两个独立的分步骤：第二Ａ步、第二Ｂ步，这两个分步骤之间没有先后顺序，所述的第二Ａ步包括：计算瞬时包络ａ（ｔ）：ａ（ｔ）＝＊＊＊（１）所述的第二Ｂ步包括两小步，第二Ｂ步第１小步是：计算相位函数ｃｏｓθ（ｔ），ｃｏｓθ（ｔ）＝ｘ↓［ｒ］（ｔ）／ａ（ｔ）（２）第二Ｂ步第２小步是：对相位函数ｃｏｓθ（ｔ）进行ＤＦＴ理想滤波器滤波，得到滤波后的相位函数ｙ↓［θ］（ｔ）；第三步：...

【技术特征摘要】
1. 一种对数字信号使用DFT理想滤波器的滤波方法，原始数据为：一道时间序列数字信号，记为xr(t)；对xr(t)进行滤波，包括但不限于下述步骤：第一步：对xr(t)进行希尔伯特变换，得到xh(t)；第二步，分为两个独立的分步骤：第二A步、第二B步，这两个分步骤之间没有先后顺序，所述的第二A步包括：计算瞬时包络a(t)：a(t)=xr2(t)+xh2(t)---(1)]]>所述的第二B步包括两小步，第二B步第1小步是：计算相位函数cosθ(t)，cosθ(t)=xr(t)a(t)---(2)]]>第二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步，计算瞬时包络之后，再对a(t)进行修饰处理。4. 根据权利要求1至3之一所述的一种对数字信号使用DFT理想滤波器的滤波方法，其特征是：在所述的第二B步第2小步，...

【专利技术属性】
技术研发人员：刘志成，方伍宝，朱海波，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：11[中国|北京]