Filtre numérique

[karika]

Salut tout le monde,voila ce que j'ai fait avec MATLAB:

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
>> %filtre butterworth
Wp =0.1 ; Ws =0.2 ;
Rp = 1; Rs = 30;
[n1,Wn1] = buttord(Wp,Ws,Rp,Rs)
[b,a] = butter(n1,Wn1);
figure(1);
freqz(b,a)
title('n=6 Butterworth Bandpass Filter')
%filtre elliptique
Wp =0.1 ; Ws = 0.2;
Rp = 1; Rs = 30;
[n2,Wp] = ellipord(Wp,Ws,Rp,Rs)
[b,a] = ellip(n2,Rp,Rs,Wp);
figure(2);
freqz(b,a);
title('n=3 Elliptic Bandpass Filter')
%filtre chebychev type 1
Wp =0.1 ; Ws =0.2 ;
Rp = 1; Rs = 30;
[n3,Wp] = cheb1ord(Wp,Ws,Rp,Rs)
[b,a] = cheby1(n3,Rp,Wp);
figure(3);
freqz(b,a);
title('n=4 Chebyshev Type 1 Bandpass Filter')

maintenant si vous pouvez m'aidez,montrez moi comment representer les modules de ces filtres en échelle logarithmique...
et merci d'avance.

[Winjerome]

Bonjour,

• Module: doc abs
• Échelle logarithmique: doc semilogy

[karika]

dacord,et merci beaucoup

[ol9245]

Envoyé par karika

dacord,et merci beaucoup

cliquer sur 'Résolu' ?

[karika]

Salut tout le monde,voila ce que j'ai fait comme programme, mais la représentation en échelle logarithmique ne s'affiche pas.
Aidez-moi s'il vous plaît et merci d'avance:

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
%filtre butterworth
Wp1 =0.1 ; Ws1 =0.2 ;
Rp1 = 1; Rs1 = 30;
[n1,Wn1] = buttord(Wp1,Ws1,Rp1,Rs1)
[b,a] = butter(n1,Wn1);
figure(1);
freqz(b,a);
title('n=6 Butterworth Bandpass Filter');
w=logspace(1,3,1000);
s=j*w;
Adb=20*log10(abs([b,a]));
figure(12);
semilogx(w,Adb);grid;
title('le module de butterworth en echelle logarythmique');
xlabel('pulsation');
ylabel('module du filtre');

[Winjerome]

Bonjour,

Quel est le but de cette ligne: Adb=20*log10(abs([b,a])); ?
Ne souhaiterais-tu pas plutôt représenter la réponse fréquentielle?

[karika]

bonsoir,le but c'est que j'ai voulue convertir le module en décibel

[Winjerome]

Le module des coefficients a et b concaténés du filtre Butterworth? Dans quel but?

[karika]

moi je veut convertir le module du filtre en décibel,alors comment on fait????svp

[Winjerome]

Eh bien en utilisant la réponse fréquentielle du filtre et non ses coefficients
>> doc freqz que tu utilises déjà

[karika]

j'ai essayé mais matlab m'affiche qu'il ya une erreur d'utilisation de ""semilogy"",mais je sais pas ou!!et merci beaucoup

[Winjerome]

Pourtant l'erreur te le dit: (traduction directe de l'anglais) les deux vecteurs rentrés doivent être de la même longueur.

[karika]

alor comment je fait???aidez moi svp

[Winjerome]

Tu ajustes tes deux vecteurs pour qu'ils aient la même taille
Cela ne te semble pas logique?

[karika]

c'est bon ça marcher,merci beaucoup

[karika]

maintenant j'ai eu un autre problème sur l'échelle logarithmique dans un autre programme qui est le suivant:

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
%filtre elliptique
Wp2 =0.1 ; Ws2 = 0.2;
Rp2 = 1; Rs2 = 30;
[n2,Wp2] = ellipord(Wp2,Ws2,Rp2,Rs2)
[b,a] = ellip(n2,Rp2,Rs2,Wp2);
A=20*log10(abs(freqz(b,a)));
w=logspace(1,3,512);
s=j*w;
subplot(2,1,1),freqz(b,a);
title('n=3 Elliptic Bandpass Filter');
subplot(2,1,2),semilogy(w,A);grid;
title('le module du filtre elliptique en echelle logarythmique');
xlabel('pulsation');
ylabel('module du filtre');

mais je sais pas ou il est ce problème....aidez moi svp

[Winjerome]

MATLAB te parle: écoute-le

Warning: Negative data ignored

Qu'est-ce que cela veut dire selon toi?

Remarque: FREQZ travaillant déjà sur 2 SUBPLOTs, préfère mettre le SEMILOGX sur une autre figure, sinon tu ne verras pas grand chose de la phase.

[karika]

ça veut dire que les données négatives ignoré,mais quelles données???