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

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

1、数字信号处理课程设计报告课设题目：地震信号的干扰去除学院：信息科学与工程学院专业：电子与信息工程班级：电信二班姓名：学号：070250212指导教师：哈尔滨工业大学（威海）2010年7月14日不要删除行尾的分节符，此行不会被打印哈尔滨工业大学（威海）课程设计报告目录一.课程设计任务1二.课程设计原理及设计方案2三.课程设计的步骤和结果3四.课程设计总结4五.设计体会5六.参考文献6-11-哈尔滨工业大学（威海）课程设计报告一.课程设计任务1.自行产生10M采样点数，截至带宽为40hz的符合要求的信号，放在文档中（至

2、少要产生三个）。2.由于在地震信号的采集中往往掺杂工频信号，设计一个滤波器将工频信号滤去。3.设计一个低通FIR滤波器对信号进行处理，由于信号较大，为了加快分析的速度，要采用学过的分段方法进行分段分析。4.采用GUI页面，方便操作。5.写实验报告。-11-哈尔滨工业大学（威海）课程设计报告一.课程设计原理及设计方案截至频率为40hz的地震信号可以由正弦信号的线性相加得到，然后利用matlab中计算出的和得到文本文档保存，即可得到“保存的地震信号”。对工频信号的滤除可以用陷波滤波器。而陷波滤波器又分自适应、小波变换、

3、模板匹配滤波器......同时实现FIR设计陷波滤波器又有窗函数法、切比雪夫逼近法、频率抽样法。本设计采用窗函数法，窗函数简单实用，且在设计中要先构造出低通滤波器在后续中也可能用到。低通滤波器可以将40hz以上的频率滤除。为了比较各个滤波器的优劣，本设计采用了三个典型的滤波器。其中观察傅里叶变换所得频域的结果更能明显看出实验结果。傅利叶变换原理卷积由于信号的时宽带宽积为一常数，所以主瓣宽度与窗口宽度成反比，N越大，的主瓣越窄。同时窗的形状也对傅里叶频谱有影响，比如吉布斯效应。虽然频率分辨率很高但由于第一旁瓣的衰减很

4、小但是求得的与的偏差较大。所以対窗的选择要根据实际需要选择合适的窗及其大小形状。而在设计中对大量信号的处理为了加快处理熟读我们可以用分段卷积：分段卷积可以分为重叠相加法和重叠保留法。重叠相加法：设h（n）的长度为M，X（n）的长度为N，将长序列X按如下方法分：根据卷积分配率：因此分段后，h与x的卷积可以认为是x的各个子段与h的卷积之和。functiony=ovrlpadd(x,h,Nfft)lenx=length(x);-11-哈尔滨工业大学（威海）课程设计报告M=length(h);M1=M-1;N=Nfft-M

5、+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=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卷积之和即为结果。-11-哈尔滨工业大学（威海）课程设计报告一.课程设计的步骤和结

6、果1.信号的产生根据傅里叶变换，任何函数都可以写作正弦函数的叠加，所以地震信号可以用频率不同的正弦信号叠加得到。在matlab的workplace区有使用的所有变量的值，将所需要的值复制到文档里，比如txt文档，然后保存。由于变量在matlab中以矩阵的形式存在，所以保存在文档里的数据也是矩阵。Fs=str2double(get(handles.fs,'string'));时域信号-11-哈尔滨工业大学（威海）课程设计报告频域信号2.陷波滤波器的设计由于要求的地震信号的截至频率为40HZ，而在实际地震的信号的获取时

7、，工频信号的存在在所难免因此要设计一个陷波滤波器滤除50HZ的信号。利用matlab设计FIR滤波器的方法有窗函数法、频率抽样法和切比雪夫逼近法等，但窗函数法我们比较熟悉且简单实用。窗函数思想：根据要求选择适当的理想低通滤波器，由于理想低通滤波器得脉冲响应是非因果且无限长，用最优化窗函数来截取就可以达到目标。functionhd=ideal(Wc,M)a=(M-1)/2;n=0:[(M-1)];m=n-a+0.000001;hd=sin(Wc*m)./(pi*m);以上为一个低通滤波器，Wc为归一化截止角频率，M为

8、阶数。陷波滤波器可以看做是一个低通滤波器加上一个高通滤波器，而高通滤波器有可以由全通滤波器减去一个低通滤波器得到。去除工频信号后的波形3.低通滤波器的设计FIR滤波器的设计问题实际上就是用要求设计的FIR滤波器的频率响应H（w）去逼近所要求的理想滤波器的响应Hd（w）。从单位取样序列来看，就是使所涉及的滤波器的h（n）逼近单位取样响应序列hd（n）。既-11