地震记录合成代码

  工程物探上有一个实验要求合成地震记录,还是比较简单的,就是反射系数序列和雷克子波卷积就可。卷积公式这儿就不写出来了。反射系数序列可以自己设置,200个点同时赋值为0,地下有几层,就把几个点设置为其他的值。对于雷克子波,有公式就可以做出来,对于工程物探,频率尽量大一些,这样地震记录的图要…