数字信号分析实验1

《数字信号分析实验1》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验一、地震子波波形显示及一维地震记录合成一、实验目的1、认识地震子波(以雷克子波为例)，对子波有直观的认识。2、利用线性褶积公式合成一维地震记录。二、实验内容1、雷克子波：w(r)=£(2欢肪2cos2^”r(零相位子波)、w(/)=£・(2"少片m2对第(最小相位子波),其中.九代表子波的中心频率，/代表子波频带宽度，随着丫的增大，子波能量后移，当了=7时，最小相位子波可视为混合相位子波，这里取fm=25Hz,/=3;2、根据公式编程实现零相位子波、最小相位子波的波形显示；3、设计反射系数r(n)(n=

2、500),其中"100)=1.0,r(200)=-0.7,r(300)=0.5,r(400)=0.4,r(500)=0.6,其它为0；4、应用褶积公式f(x)=r(n)*w(n)=r(/?i)vv(n-m)合成一维地震记录，并m=l图形显示；5、根据所学知识对实验结果进行分析。三、实验要求1、独立完成程序设计工作；2、独立完成书面报告(A4纸打印)；3、提交程序源代码。实习报告一姓名：王忠成学号：050422011042专业年级：地信2011级当Y=3吋，汇出最小相位子波和零相位子波图形，如下所示：当Y=5

3、和Y=7时的波形图分别为丫令时的零相位子波丫巧时的最小相位子波山雷克子波公式W(f)=£一(如"尸cos2

4、=1位子波分别与r(n)褶积，汇出其图形如下：

5、最小相位子波的地震记录由图可知，两种波形的地震记录都育5个峰值，分别对应于5个分界面:r(100).r(200)、r(300)、r(400)和r(500),并且当「＞0时，峰值为正，当rvO时，峰值为负。w和r的长度都是500,则褶积后的长度为500+500-1=999o程序源代码：FM=25;%形成两种相位子波并用图形显示gama=3;t=linspace((),().l,500);%对「赋初值wl=exp(-(2*pi*FM.*t/gama).*(2*pi*FM.*t/gama)).*cos(2*pi*

6、FM.*t);w2=exp(-(2*pi*FM.*t/gama).*(2*pi*FM.*t/gama)).*sin(2*pi*FM.*t);subplot(i,2,l),plot(t,w1),gridon,title('零相位子波)subplot(1,2,2),plot(t,w2),gridon,title(,最小相位子波Jfori=1:500r(i)=0;endr(100)=1.0;r(200)=-0.7;r(300)=0.5;r(400)=0.4;r(500)=0.6;fori=1:999%对fl初始化

7、fl(i)=0;endfori=1:500%求wl和i•的褶积,记为flforj=l:ifl(i)=fl(i)+r(j)*wl(i+l-j);endendfori=501:999forj=(i+1-500):500fl⑴=fl(i)+r(j)*wl(i+l-j);endendfori=l:999n(i)=i；endplot(n,fl),titleC最小相位子波的地震记录)fori=1:999%对f2初始化f2(i)=0;endfori=l:500%求w2和r的褶积,记为flforj=l:if2(i)=f2(

8、i)+r(j)*w2(i+1-j);endendfori=501:999forj=(i+1-500):500f2(i)=f2(i)+r(j)*w2(i+l-j);endendplol(n,⑵，lillc(零相位子波的地震记录')