资源描述：
《数字信号分析大作业一》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验一、地震子波波形显示及一维地震记录合成一、实验内容1、认识地震子波（以雷克型子波为例），对子波的波形有直观的认识。2、利用褶积公式合成一维地震记录。二、实验步骤1、雷克型子波：（零相位子波）、（最小相位子波），其中代表子波的中心频率，代表子波宽度,随着的增大，子波能量后移，这里取=25Hz，=3;2、根据公式编程实现零相位子波、最小相位子波的波形显示；3、设计反射系数(n=500)，其中，，，，，其它为0;4、应用褶积公式合成一维地震记录，并图形显示；5、根据所学知识对实验结果进行分析。三、实验要求1、独立完成程序设计工作；2、独立完成书面报告（A4纸打
2、印）；3、提交程序源代码。一、零相位与最小相位雷克子波的构建注：由采样定理可得t<=1/(2*fm)，故t<=0.02s，实验中采样间隔应小于0.02s，本次实验中dt取0.001s，满足采样定理。零相位子波：dt=0.001;fm=25;y=3;fori=1:300w(i)=exp(-(2*pi*fm/y)^2*(dt*i)^2)*cos(2*pi*fm*dt*i);endplot(w,'LineWidth',5));xlabel('t');ylabel('w（t）');title('零相位子波')；最小相位子波：dt=0.001;fm=25;y=3;fo
3、ri=1:300w(i)=exp(-(2*pi*fm/y)^2*(dt*i)^2)*sin(2*pi*fm*dt*i);endplot(w,'LineWidth',5);xlabel('t');ylabel('w（t）');title('最小相位子波');一、应用褶积公式合成一维地震记录注：序列总长度为:L=n+m-1=300+500-1=800，褶积公式中判定条件为：0。clearalldt=0.001;fm=25;y=3;fori=1:300wcos(i)=exp(-(2*pi*fm/y)^2*(dt*i)^2)*cos(2*pi*fm
4、*dt*i);wsin(i)=exp(-(2*pi*fm/y)^2*(dt*i)^2)*sin(2*pi*fm*dt*i);endsubplot(2,2,1);plot(wcos,'LineWidth',2);xlabel('t');ylabel('wcos（t）');title('零相位子波');subplot(2,2,2);plot(wsin,'LineWidth',2);xlabel('t');ylabel('wsin（t）');title('最小相位子波');forn=1:500r(n)=0;endr(100)=1.0;r(200)=-0.7;r(3
5、00)=0.5;r(400)=0.4;r(500)=0.6;forn=1:800fcos(n)=0;fsin(n)=0;endg(1)=0;forn=1:799form=1:500if(n-m>0&&n-m<300)fcos(n)=fcos(n)+r(m)*wcos(n-m);fsin(n)=fsin(n)+r(m)*wsin(n-m);endendendsubplot(2,2,3);plot(fcos,'LineWidth',2);title('零相位子波的一维地震记录');subplot(2,2,4);plot(fsin,'LineWidth',2);t
6、itle('最小位子波的一维地震记录')一、实验分析1、在进行实验设计时应保证采样频率满足采样定理t<=1/(2*fm)。2、应用褶积公式进行计算时，要注意序列的长度L=n+m-1，在进行循环判断时要控制好判断的条件，使其满足褶积公式。3、通过观察零相位子波的地震记录和最小相位的地震子波的地震记录，可以发现最小相位子波的地震记录振幅相对于零相位子波更小，同时图像的不像最小相位子波的地震记录那样有突变，而是经历了一个比较小的时间才达到最大值，相比之下使用零相位子波合成地震记录，开始检测到地震波，图像就会有突越，这样必须提高检波器的灵敏度才能准确的记录地震波，而
7、使用最小相位子波可以有效地弱化这个问题，提高地震检波的精度。记录地震波开始的时刻比记录地震波最大值的时刻要困难，所以最小相位子波合成地震记录是一种更好的选择。姓名：杨秋野年级：2014级勘察技术与工程学号：14040031047
