数字信号处理课程设计-地震信号的干扰去除

《数字信号处理课程设计-地震信号的干扰去除》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、.DONGFANGCOLLEGE，FUJIANAGRICULTUREANDFORESTRYUNIVERSITY课程名称：数字信号处理课设题目：地震信号的干扰去除系别：计算机系年级专业：10电子信息工程2班学号：10503020XX姓　　名：凯爷任课教师：罗志聪成绩：2013年月日一、作品简介现有1年的某地震信号，一共有365个文件，每个文件长度为10M个采样点，采样率为200Hz，由于受到50hz的信号干扰，请设计一FIR滤波器，对其进行滤波，为了加快分析的速度，请采用学过的分段方法进行分段分析。如果信号截至带宽为40hz，请设计一

2、低通滤波器FIR滤波器对以上信号进行处理；可自行产生10M采样点数的复合要求的信号，放在某文档中即可，为了证明可以对365个文档中的任意一个文档可以处理，请至少产生3个文档。二、设计目的数字信号处理是一门理论和技术发展十分迅速、应用非常广泛的前沿性学科，在雷达、通信、医学、地震等众多领域都有广泛的应用。《数字信号处理》课程是电子、通信、信号处理等专业的专业基础课程，数字信号处理课程设计是完成数字信号处理相关理论的学习后进行的综合性训练课程，使学生能够对信号的采集、处理、传输、显示和存储等理论及技术实现有一个系统的掌握和理解。本课程设

3、计的目的为：1、巩固和加深学生对数字信号处理的基本概念、基本原理、基本分析方法的理解；2、培养学生根据设计题目的要求，综合利用数字信号处理技术的能力，提高学生综合运用所学知识，解决实际工程技术问题的能力；3、通过课程设计，培养学生查阅文献、资料，独立解决实际问题的能力。通过课程设计，要求学生综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用MATLAB或LabView作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念，使学生基本掌握以下内容：1、学会MATLAB或LabView的使用

4、，掌握其程序设计方法。2、掌握在Windows环境下语音信号采集的方法。3、掌握数字信号处理的基本概念、基本理论和基本方法。4、掌握设计FIR和IIR数字滤波器的方法。5、学会用MATLAB或LabView对信号进行分析和处理。三、设计技术和方法截至频率为40hz的地震信号可以由正弦信号的线性相加得到，然后利用matlab中计算出的和得到文本文档保存，即可得到“保存的地震信号”。对工频信号的滤除可以用陷波滤波器。而陷波滤波器又分自适应、小波变换、模板匹配滤波器......同时实现FIR设计陷波滤波器又有窗函数法、切比雪夫逼近法、频率

5、抽样法。本设计采用窗函数法，窗函数简单实用，且在设计中要先构造出低通滤波器在后续中也可能用到。低通滤波器可以将40hz以上的频率滤除。为了比较各个滤波器的优劣，本设计采用了三个典型的滤波器。其中观察傅里叶变换所得频域的结果更能明显看出实验结果。傅利叶变换原理卷积由于信号的时宽带宽积为一常数，所以主瓣宽度与窗口宽度成反比，N越大，的主瓣越窄。同时窗的形状也对傅里叶频谱有影响，比如吉布斯效应。虽然频率分辨率很高但由于第一旁瓣的衰减很小但是求得的与的偏差较大。所以対窗的选择要根据实际需要选择合适的窗及其大小形状。而在设计中对大量信号的处理

6、为了加快处理熟读我们可以用分段卷积：分段卷积可以分为重叠相加法和重叠保留法。重叠相加法：设h（n）的长度为M，X（n）的长度为N，将长序列X按如下方法分：根据卷积分配率：因此分段后，h与x的卷积可以认为是x的各个子段与h的卷积之和。functiony=ovrlpadd(x,h,Nfft)lenx=length(x);M=length(h);M1=M-1;N=Nfft-M+1;x=[x,zeros(1,N-1)];K=floor(Lenx+Nfft-1);fori=0:1:Kix=i*N;xseg=x(ix+1:ix+N);yseg=

7、circonv(xseg,h,Nfft);y(ix+1:ix+Nfft)=y(ix+1:ix+Nfft)+yseg(1:Nfft);y=y(1:Lenx+M-1);重叠保留法：同样将序列X（n）按如下方式分成长度为L的字段：即每个子段之间有M-1哥点重叠。将分段后的x与h卷积之和即为结果。四、课程设计的步骤和结果1.信号的产生根据傅里叶变换，任何函数都可以写作正弦函数的叠加，所以地震信号可以用频率不同的正弦信号叠加得到。在matlab的workplace区有使用的所有变量的值，将所需要的值复制到文档里，比如txt文档，然后保存。由于

8、变量在matlab中以矩阵的形式存在，所以保存在文档里的数据也是矩阵。Fs=str2double(get(handles.fs,'string'));时域信号频域信号陷波滤波器的设计由于要求的地震信号的截至频率为40HZ，而在实际地震