《数字信号处理课程设计报告-地震信号的干扰去除》由会员分享,可在线阅读,更多相关《数字信号处理课程设计报告-地震信号的干扰去除(13页珍藏版)》请在金锄头文库上搜索。
1、数字信号数字信号处处理理课程设计报告课设题目:课设题目: 地震信号的干扰去除学学 院:院: 信息科学与工程学院专专 业:业: 电子与信息工程班班 级:级: 电信二班姓姓 名:名: 学学 号:号: 指导教师:指导教师: 哈尔滨工业大学(威海)哈尔滨工业大学(威海)20102010 年年7 7 月月1414 日日哈尔滨工业大学(威海)课程设计报告- I -目 录一. 课程设计任务.1二. 课程设计原理及设计方案.2三. 课程设计的步骤和结果.3四. 课程设计总结.4五. 设计体会.5六. 参考文献.6哈尔滨工业大学(威海)课程设计报告- 1 -一. 课程设计任务1. 自行产生 10M 采样点数,截
2、至带宽为 40hz 的符合要求的信号,放在文档中(至少要产生三个) 。2.由于在地震信号的采集中往往掺杂工频信号,设计一个滤波器将工频信号滤去。3.设计一个低通 FIR 滤波器对信号进行处理,由于信号较大,为了加快分析的速度,要采用学过的分段方法进行分段分析。4.采用 GUI 页面,方便操作。5.写实验报告。哈尔滨工业大学(威海)课程设计报告- 2 -二. 课程设计原理及设计方案截至频率为40hz的地震信号可以由正弦信号的线性相加得到,然后利用matlab中计算出的和得到文本文档保存,即可得到“保存的地震信号” 。对工频信号的滤除可以用陷波滤波器。而陷波滤波器又分自适应、小波变换、模板匹配滤波
3、器.同时实现FIR设计陷波滤波器又有窗函数法、切比雪夫逼近法、频率抽样法。本设计采用窗函数法,窗函数简单实用,且在设计中要先构造出低通滤波器在后续中也可能用到。低通滤波器可以将40hz以上的频率滤除。为了比较各个滤波器的优劣,本设计采用了三个典型的滤波器。其中观察傅里叶变换所得频域的结果更能明显看出实验结果。傅利叶变换原理2210()( )( )(01)kkmNjjNNnnnmX eX kx m ekN卷积()-1()()()2jwjjnj wnX eW eeX ed 由于信号的时宽带宽积为一常数,所以主瓣宽度与窗口宽度成反比,N越大,的主瓣越窄。同时窗的形状也对傅里叶频谱有影响,比如吉布斯效
4、应。虽然频()jwW e率分辨率很高但由于第一旁瓣的衰减很小但是求得的与的偏差较大。()jwnX e()jwX e所以対窗的选择要根据实际需要选择合适的窗及其大小形状。而在设计中对大量信号的处理为了加快处理熟读我们可以用分段卷积:分段卷积可以分为重叠相加法和重叠保留法。重叠相加法:设h(n)的长度为M,X(n)的长度为N,将长序列X按如下方法分:根据卷积分配率:因此分段后,h与x 的卷积可以认为是x的各个子段与h的卷积之和。function y=ovrlpadd(x,h,Nfft)lenx=length(x);哈尔滨工业大学(威海)课程设计报告- 3 -M=length(h);M1=M-1;N
5、=Nfft-M+1;x=x,zeros(1,N-1);K=floor(Lenx+Nfft-1);for i=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卷积之和即为结果。哈尔滨工业大学(威海)课程设计报告- 4 -三. 课程设计的步骤和结果1.信号的产生根据傅里叶变换,任何函数都可以写作正弦函数的叠加
6、,所以地震信号可以用频率不同的正弦信号叠加得到。在 matlab 的 workplace 区有使用的所有变量的值,将所需要的值复制到文档里,比如 txt 文档,然后保存。由于变量在 matlab 中以矩阵的形式存在,所以保存在文档里的数据也是矩阵。 Fs=str2double(get(handles.fs,string);00.511.522.533.54-8-6-4-202468时域信号哈尔滨工业大学(威海)课程设计报告- 5 -频域信号2.陷波滤波器的设计由于要求的地震信号的截至频率为 40HZ,而在实际地震的信号的获取时,工频信号的存在在所难免因此要设计一个陷波滤波器滤除 50HZ 的信
7、号。利用 matlab 设计 FIR 滤波器的方法有窗函数法、频率抽样法和切比雪夫逼近法等,但窗函数法我们比较熟悉且简单实用。窗函数思想:根据要求选择适当的理想低通滤波器,由于理想低通滤波器得脉冲响应是非因果且无限长,用最优化窗函数来截取就可以达到目标。function hd=ideal(Wc,M)a=(M-1)/2;n=0:(M-1);m=n-a+0.;hd=sin(Wc*m)./(pi*m);以上为一个低通滤波器,Wc为归一化截止角频率,M为阶数。陷波滤波器可以看做是一个低通滤波器加上一个高通滤波器,而高通滤波器有可以由全通滤波器减去一个低通滤波器得到。去除工频信号后的波形3.低通滤波器的
8、设计FIR 滤波器的设计问题实际上就是用要求设计的 FIR 滤波器的频率响应 H(w)去逼近所要求的理想滤波器的响应 Hd(w) 。从单位取样序列来看,就是使所涉及的滤波器的 h(n)逼近单位取样响应序列 hd(n) 。既哈尔滨工业大学(威海)课程设计报告- 6 -由于理想低通滤波器 ,其单位取样响应 hd(n)是无限长的,是非因果序列。因此为了构造有限长度的线性滤波器需要将理想滤波器截取一段,这可以用理想低通滤波器与一个窗函数相乘得到。窗函数可以有矩形窗、汉宁窗、海明窗、布莱克曼窗等。窗函数的要求为:A.窗谱主瓣尽可能窄,可以获取较陡的过渡带B.尽量减小窗谱的最大旁瓣的相对幅度;即使能量集中
9、于主瓣,使肩峰和波纹减小,增大阻带的衰减。本设计采用以下三种: 凯赛窗:rp=1;rs=-100;fcuts=39.8 40.2;rp 为通带衰减,rs 为阻带衰减,fcuts为过渡带宽度。而对于凯赛窗可以直接调用 matlab 中的函数。n,wn,beta,ftype=kaiserord(fcuts,mags,devs,Fs);hh=fir1(n,wn,ftype,kaiser(n+1,beta);Y=fft(p1);x=fftfilt(hh,p1);凯赛窗频率响应哈尔滨工业大学(威海)课程设计报告- 7 -经过低通滤波器后的频谱布莱克曼: fp1=39;N=2*ceil(12*pi/(ws
10、1-wp1)/2);window=blackman(N+1);wn=(wp1+ws1)/2;bf1=fir1(N,wn/pi,low,window);output1=fftfilt(bf1,p1);X=fft(output1);布莱克曼窗函数频率响应哈尔滨工业大学(威海)课程设计报告- 8 -汉宁窗bf2=fir1(n-1,wc,low,hanning(n);output2=fftfilt(bf2,p1);X=fft(output2);汉宁窗频率响应 哈尔滨工业大学(威海)课程设计报告- 9 -4课程设计总结Matlab 是一款强大的软件:其基本数据单位是矩阵,比 C,FORTRAN 等语言完
11、成相同的事情简捷,迅速。同时 Matlab 的应用范围非常广,包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。很多情况下对于滤波器的选择都要根据实际情况、成本、物理环境等等其他因素综合考虑。比如有写滤波器可能精度高但是成本高,而在很多高科技领域(军事、航天等)可忽略成本。大多数滤波器的原理是相通的,比如由高通滤波器可由全通滤波器减去低通滤波器,陷波滤波器可以看做是一个低通滤波器加上一个高通滤波器。哈尔滨工业大学(威海)课程设计报告- 10 -五设计体会由于平时很少使用 matlab 进行这类复杂的设计,所以一开始的困难时在所难免得,但是感觉万事开
12、头难,当你找到方向时问题就会迎刃而解。在很多情况下真本事是很重要的,在此次课程设计中有很多东西都要从网上找,但是找到的东西要么有错误要么就不适合,需要自己修正。这就需要我们懂的基础知识,把别人的东西掌握并能活学活用。哈尔滨工业大学(威海)课程设计报告- 11 -六参考文献具体格式参考如下:1, 余成波 杨如民等编著 数字信号处理及 MATLAB 实现清华大学出版社出版2, 程佩青数字信号处理教程 清华大学出版社出版3 李杰 张猛 邢笑雪 编著 信号处理 matlab 实验教程 北京大学出版社 4 薛年喜 编著 matlab 在数字信号处理中的应用 清华大学出版社 5 董紹平 陈世耕 编著 数字信号处理基础 哈尔滨工业大学出版社