Création d'une ondelette dans Matlab

Question

J'ai quelques données de réflectivité que je vais convoluer avec une ondelette de chapeau Ricker/mexicaine pour obtenir une trace sismique. Mon problème est dû à la création de l'ondelette. Je voudrais que l'ondelette ait une fréquence dominante d'environ 70Hz et que le pas de temps soit de 0.19 ms, ce qui est le même pas de temps que dans mes données de réflectivité. J'ai essayé d'utiliser la fonction mexihat dans MatLab et j'ai réglé les paramètres lb, ub et n pour résoudre mon problème, mais je ne peux pas le comprendre. Existe-t-il un moyen plus simple de résoudre mon problème? Est-ce que quelqu'un sait une formule pour l'ondelette de Ricker où la fréquence dominante se produit?

Toute aide serait grandement appréciée!

Answer 1

De Large Data in MATLAB: A Seismic Data Processing Case Study:

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