Sauts dans les valeurs calculees

Version imprimable

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.)

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

18/01/2013, 14h17
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.
19/01/2013, 16h03
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.

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

19/01/2013, 16h41
Howyd

Ok. J'ai compris.
Merci.