工程物探上有一个实验要求合成地震记录，还是比较简单的，就是反射系数序列和雷克子波卷积就可。卷积公式这儿就不写出来了。反射系数序列可以自己设置，200个点同时赋值为0，地下有几层，就把几个点设置为其他的值。对于雷克子波，有公式就可以做出来，对于工程物探，频率尽量大一些，这样地震记录的图要好看一些。 下边把代码奉献给大家。   #include<stdio.h> #include<math.h> #include<stdlib.h> //#define dz1…

 一、实验目的 通过编制和运行相关C或matlab程序，理解地震子波的概念和特点，掌握地震合成记录的制作流程。 认识子波，对子波的波形有直观的认识。（零相位子波，混合相位子波，最小相位子波；了解子波的分辨率与频宽的关系；） 利用褶积公式合成一维地震记录 二、主要内容和原理 1、理论子波的计算及显示         （1）零相位子波、最小相位子波、钟型子波的计算 零相位子波的表达式 最小相位子波的表达式为； n=m1/m2，为最大波峰m1和最大波谷m2之比，取2.0或1.5，本次试验采用2.0.；代表子波的中心频率，…