数字信号分析实验1.docVIP

实验一、地震子波波形显示及一维地震记录合成 一、实验目的 认识地震子波（以雷克子波为例），对子波有直观的认识。 利用线性褶积公式合成一维地震记录。 二、实验内容 雷克子波： （零相位子波）、 （最小相位子波）， 其中代表子波的中心频率，代表子波频带宽度,随着的增大，子波能量后移，当=7时，最小相位子波可视为混合相位子波，这里取 = 25 Hz，= 3; 根据公式编程实现零相位子波、最小相位子波的波形显示； 设计反射系数(n=500)，其中，，，，，其它为0; 应用褶积公式合成一维地震记录，并图形显示； 根据所学知识对实验结果进行分析。 三、实验要求 1、独立完成程序设计工作； 2、独立完成书面报告（A4纸打印）； 3、提交程序源代码。 实习报告一 姓名：王忠成 学号：050422011042 专业年级：地信2011级 当γ=3时，汇出最小相位子波和零相位子波图形，如下所示： 当γ=5和γ=7时的波形图分别为： 由雷克子波公式知：(2πfm/γ)2为一个指数衰减因子，当γ增大时，(2πfm/γ)2减小，故衰减变慢。 设计反射系数(n=500)，其中，，，，，其它为0。应用褶积公式合成一维地震记录，对零相位子波和最小相位子波分别与r(n)褶积，汇出其图形如下： 由图可知，两种波形的地震记录都有5个峰值，分别对应于5个分界面:r(100)、r(200)、r(300)、r(400)和r(500)，并且当r0时，峰值为正，当r0时，峰值为负。w和r的长度都是500，则褶积后的长度为500+500-1=999。 程序源代码： FM=25; %形成两种相位子波并用图形显示 gama=3; t=linspace(0,0.1,500); w1=exp(-(2*pi*FM.*t/gama).*(2*pi*FM.*t/gama)).*cos(2*pi*FM.*t); w2=exp(-(2*pi*FM.*t/gama).*(2*pi*FM.*t/gama)).*sin(2*pi*FM.*t); subplot(1,2,1),plot(t,w1),grid on,title(零相位子波) subplot(1,2,2),plot(t,w2),grid on,title(最小相位子波) for i=1:500 %对r赋初值 r(i)=0; end r(100)=1.0; r(200)=-0.7; r(300)=0.5; r(400)=0.4; r(500)=0.6; for i=1:999 %对f1初始化 f1(i)=0; end for i=1:500 %求w1和r的褶积, 记为f1 for j=1:i f1(i)=f1(i)+r(j)*w1(i+1-j); end end for i=501:999 for j=(i+1-500):500 f1(i)=f1(i)+r(j)*w1(i+1-j); end end for i=1:999 n(i)=i; end plot(n,f1), title(最小相位子波的地震记录) for i=1:999 %对f2初始化 f2(i)=0; end for i=1:500 %求w2和r的褶积, 记为f1 for j=1:i f2(i)=f2(i)+r(j)*w2(i+1-j); end end for i=501:999 for j=(i+1-500):500 f2(i)=f2(i)+r(j)*w2(i+1-j); end end plot(n,f2), title(零相位子波的地震记录)