Fréquence de coupure filtre binomial gaussien-matlab

[inviteb8477130]

Bonjour;

Je souhaite calculer la fréquences de coupure horizontale et verticale d'un filtre binomiale-gaussien 2D. J'ai écris le petit programme ci-dessous, mais j'ai un doute que ce soit correcte, car j'obtiens 4 fréquences correspondant à ~-3db en réponse fréquentielle(4 freq horizontales et 4 verticales). Je ne sais pas laquelle correspond à la réelle fréquence de coupure...j'ai supposé que c'est la fréquence maximale parmi les 4...Mais est-ce juste?? Si vous pouvez m'aider sur ce point svp!!
Voici mon programme:

h_bg=1/16*[1 2 1;2 4 2;1 2 1];
[H,f1,f2]=freqz2(h_bg,32,32);
Hdb=20*log10(H);
k1=[];
k2=[];
Hdb2=[];
for i=1:32
for j=1:32
if (Hdb(i,j)>(-3.2) && Hdb(i,j)<(-2.94))
k1=[k1 f1(i)];
frec1=max(k1) %fréq de coupure horizontale ??
k2=[k2 f2(i)];
frec2=max(k2) % freq de coupure verticale ??
break

end
end
end

Merci d'avance!!!
-----

[acx01b]

bonjour, ton filtre 2D peut être décomposé comme le produit de convolution de deux filtres << 1D >>

[inviteb8477130]

Certes...mais ça ne répond pas au problème que j'ai posé...autrement dit, existe t-il une fonction matlab qui donne directement la fréquence de coupure exacte d'un filtre 1D...car si je refais ma méthode (même pour le filtre 1D) j'aurais toujours le même pb (ie que j'ai 4 fréquences pour lesquelles le module est à ~-3db, comment choisir parmi les 4 fréquences la bonne??)

Merci bcp!!!

[acx01b]

si tu fais freqz([1 2 1],1) ?

ou plutôt freqz([1 2 1],[0,1]) pour avoir la bonne phase

[A voir en vidéo sur Futura]