Discussion :

[Howyd]

Bonjour,
je souhaite tracer la vitesse de groupe d 'un filtre analogique.
J'ai tape ce code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
function [] = gen_butter()
 
ordre=1;
 
    [b,a]=butter(ordre, 10, 's');
    [h,w]=freqs(b,a);
 
    for jj=2:length(w)-1
        tgd(jj)=-(unwrap(angle(h(jj+1)))-unwrap(angle(h(jj-1))))/(w(jj+1)-w(jj-1));
    end
 
 
    tgd(1)=tgd(2);
    tgd(length(w))=tgd(length(w)-1);
    semilogx(w,tgd)
 
end

Ca a l'air ok pour ordre=1 ou 2 mais des pics apparaissent pour les valeurs superieures.

Merci.

(Clavier qwerty, teste sur Matlab 2012a et 2010a.)

[Howyd]

Alors je ne sais pas pourquoi mais avec une variable intermediaire, ca marche (variable "un").
Si quelqu'un connait la raison,...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
function [] = gen_butter()
 
ordre=10;
w=1:0.1:100;
 
for ii=1:ordre
 
%     [b,a]=cheby1(ii,0.5, 10, 's');
    [b,a]=butter(ii,10, 's');
    [h,w]=freqs(b,a,w);
    un=unwrap(angle(h));
 
    for jj=2:length(w)-1
%         tgd(jj)=-(unwrap(angle(h(jj+1)))-unwrap(angle(h(jj-1))))/(w(jj+1)-w(jj-1));
        tgd(jj)=-(un(jj+1)-un(jj-1))/(w(jj+1)-w(jj-1));
    end
 
 
    tgd(1)=tgd(2);
    tgd(length(w))=tgd(length(w)-1);
 
 
    semilogx(w,tgd)
    hold all
    clear tgd
end
hold off
 
end

[Invité]

Bonjour,

La fonction unwrap effectue une comparaison des valeurs successives pour savoir s'il doit les réajuster.
Si tu ne fournis qu'une seule valeur (comme dans ton premier code), la même valeur sera retournée, il n'y aura donc aucune correction effectuée.

[Howyd]

J'avais dépouillé les matrices calculées pour unwrap(angle(h)) et w. Dans les deux toutes les valeurs "se suivaient", mais au moment de faire

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

tgd(jj)=-(unwrap(angle(h(jj+1)))-unwrap(angle(h(jj-1))))/(w(jj+1)-w(jj-1));

, pour certaines valeurs de jj (aléatoires), mon tgd avaient des pics.

Histoire de faire encore plus bizarre, ma seconde fonction laisse des pics avec un filtre elliptique

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

[b,a]=ellip(ii, 0.5, 60, 1, 's');

. J'ai ajusté les valeurs à la main et à la louche.

[Invité]

Comme tu ne semble pas comprendre, je te laisse exécuter et regarder ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
ordre=10;
w=1:0.1:100;
 
for ii=1:ordre
    [b,a] = butter(ii,10, 's');
    [h,w] = freqs(b,a,w);
    un = unwrap(angle(h));
    plot(w,angle(h), w, un , '-g')
    hold on
    for jj = 1:length(w)
        un(jj) = unwrap(angle(h(jj))); % une seule valeur => un(jj) == angle(jj)
    end
    plot(w,un,'r.')
    legend('angle original' , 'unwrap' , 'unwrap(jj)')
    pause
    hold off
end

[Howyd]

Ok. J'ai compris.
Merci.