[Débutant] Sauts dans les valeurs calculees [Résolu]

[Howyd]

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

Code :

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 :

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

[Winjerome]

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 :

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 :

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

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

[Winjerome]

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

Code :

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.