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