Bonjour,

Je dois comprendre un code matlab pour pouvoir l'utiliser et je bloque sur la partie où on veut calculer le filtre passe-bas dans l'espace de fourier pour des données de déplacements.

qmax=nr2/(pi*min_feature_size);
%min_feature_size: spatial resolution of the stress measurement in units of the grid spacing.
%nr2 = number of rows and columns across field (must be square)

% Get distance from of a grid point from the centre of the array
y=repmat((1:nr2)'-nr2/2,1,nr2);
x=y';
q=sqrt(x.^2+y.^2);

% Make the filter
qmsk=exp(-(q./qmax).^2);
qmsk=ifftshift(qmsk);

D'abord, je ne suis pas sure de quel filtre a été utilisé : "exponential low pass filter in Fourier space" correspond-il à un filtrage dans le domaine fréquentiel comme décrit dans ce cours : http://www.isima.fr/~vbarra/IMG/pdf/TP3.pdf ? ou est-ce un filtre gaussien?

Ensuite je ne comprends pas trop la démarche du code matlab, ni à quoi correspond la fonction ifftshift dans ce cas

Merci pour votre aide!!!
Aude