一种基于DSP的快速滤波方法技术

本发明专利技术涉及一种基于DSP的快速滤波方法。属于地球物理探测技术领域。方法包括以下步骤：步骤1根据滤波参数计算滤波缓冲阵列参数；步骤2建立滤波缓冲阵列，长度为N1；步骤3将实时采集数据写入步骤2建立的阵列中；步骤4在阵列中添加随机头序列；步骤5对缓冲阵列进行快速傅里叶变换；步骤6根据滤波系数去除带外频点；步骤7对带外置零的缓冲阵列进行快速傅里叶反变换；步骤8提取时域信息。本发明专利技术解决传统滤波算法参数过多导致的运算量过大，无法应用于嵌入式算法的限制性的问题，同时去除传统滤波算法的通带内纹波，减少相应的频率失真。减少相应的频率失真。

【技术实现步骤摘要】
一种基于DSP的快速滤波方法


[0001]本专利技术涉及一种快速滤波方法，特别是涉及一种用于基于DSP的快速滤波方法。属于地球物理探测


技术介绍

[0002]对于地质勘探而言，其仪器需要自身快速处理数据，且对数据进行滤波等操作，现存的滤波技术手段存在如下缺点：
[0003]1.通带内纹波较高，造成通带不平坦，引起频率失真效应；
[0004]2.传统滤波算法造成滤波系数较多，会引起运算大幅增加，进而占用系统资源，不适用于地震采集仪器嵌入式算法；
[0005]3.传统滤波算法对于频带的高频或低频截止，不能做到很高，如若做到很高，当采用FIR滤波器时，必须用非常多的系数进行拟合，如若用IIR滤波器，则会引起相位失真，对于线性相位所要求的应用所不适用；
[0006]4.对于带通滤波器，传统滤波算法实现比较困难，往往采用滤波器级联方式来实施，但如若用此方式，则会引起运算量的成倍增加。

技术实现思路

[0007]本专利技术要解决的技术问题是提供一种基于DSP的快速滤波方法，解决传统滤波算法参数过多导本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于DSP的快速滤波方法，其特征在于：所述方法包括以下步骤：步骤1根据滤波参数计算滤波缓冲阵列参数；步骤2建立滤波缓冲阵列，长度为N1；步骤3将实时采集数据写入步骤2建立的阵列中；步骤4在阵列中添加随机头序列；步骤5对缓冲阵列进行快速傅里叶变换；步骤6根据滤波系数去除带外频点；步骤7对带外置零的缓冲阵列进行快速傅里叶反变换；步骤8提取时域信息。2.根据权利要求1所述的基于DSP的快速滤波方法，其特征在于：所述方法步骤1具体为：对于低通滤波，根据FFT，对于采样率为f
s
的信号，设滤波缓冲阵列的计算长度为N1，则FFT变换后，频率分辨率为设低通截止频率为ω
p
，则是的整数倍，即存在整数n，使对于高通滤波，根据FFT，对于采样率为f
s
的信号，设滤波缓冲阵列的计算长度为N1，则FFT变换后，频率分辨率为设高通截止频率为ω
p
，则ω
p
是的整数倍，即存在整数n，使对于带通滤波，根据FFT，对于采样率为f
s
的信号，设滤波缓冲阵列的计算长度为N1，则FFT变换后，频率分辨率为设高通截止频率为ω
p
，则存在整数同时满足公式一和公式二；由此可计算出滤波缓冲阵列的长度N1；DSP中FFT数据头失真长度为：N
r
＝10 公式三在缓冲阵列前端，添加随机头序列，序列长度为N
r
；后续存入的实时数据的阵列长度为N
s
。3.根据权利要求2所述的基于DSP的快速滤波方法，其特征在于：所述方法中步骤2具体为：建立滤波缓冲区阵列，长度为N1具体为在内存中开辟长度为N1的数组，接收长度为N
s
的数据流，以及长度为N
r
的随机阵列，并将此两类数据应用于后续的滤波运算：
N
s
+N
r
＝N
1 公式四。4.根据权利要求3所述的基于DSP的快速滤波方法，其特征在于：所述方法中步骤3具体为：将采集到的实时数据实时添加到步骤2所开辟的缓冲区...

【专利技术属性】
技术研发人员：许晴，钟玉林，钟沛江，梁曾，龙建民，
申请(专利权)人：中地装重庆地质仪器有限公司，
类型：发明
国别省市：