import numpy as np#调用随机数
import matplotlib.pyplot as plt#调用画图库
from math import *
n=60
wt=np.zeros(n)
f0=30#给定主频
dt=0.001
i=0
for k in range(int(-n/2),int(n/2)):
wt[i]=(1-2.0*(pi*f0*k*dt)**2)*exp(-1*(pi*f0*k*dt)**2)
i+=1
plt.plot(list(range(int(-n/2),int(n/2))),wt)
plt.title('Ricker-time')
plt.show()
雷克子波是地震子波中具有代表性的一种,用雷克子波合成一维、二维、三维地震记录是勘探开发领域计算机编程的基础