matlab - 在 Matlab 中创建小波

Question

我有一些反射率数据，我将用 Ricker/墨西哥帽子小波对这些数据进行卷积以获得地震轨迹。我的问题是由于创建了小波。我希望小波的主频率约为 70Hz，时间步长为 0.19 ms，这与我的反射率数据中的时间步长相同。我曾尝试在 MatLab 中使用 mexihat 函数并调整 lb、ub 和 n 参数来解决我的问题，但我无法弄清楚。有没有更简单的方法来解决我的问题？有谁知道出现主频率的 Ricker 小波的公式？

任何帮助将不胜感激！

Answer 1

从MATLAB 中的大数据：地震数据处理案例研究：

% N : number of points you want to plot
[rw,t] = ricker(70,N,0.019);
plot(t,rw), xlabel('Time'), ylabel('Amplitude')

Answer 2

% using the equation below
% https://wiki.seg.org/wiki/Dictionary:Ricker_wavelet
ric = @(t,fm)(1-2*pi()^2*fm^2*t.^2)*exp(-1*pi()^2*fm^2*t.^2); % Ricker equation
t = 0:001:2;
plot(ric(t,20));
hold on;
plot(ric(t,5));

Answer 3

请注意，我尝试了链接到的函数 lucasg，并在计算小波的公式中发现了一个错字：

s = (1-tau.*tau*f^2*pi^2).*exp(-tau.^2*pi^2*f^2);
应替换为：
s = (1-2*tau.*tau*f^2*pi^2).*exp(-tau.^2*pi^2*f^2);

否则，你会得到一个小脚变种人