Discussion: Sauts dans les valeurs calculees [Débutant]

  1. #1
    Invité de passage
    Homme Profil pro
    Inscrit en
    février 2012
    Messages
    7
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Finlande

    Informations forums :
    Inscription : février 2012
    Messages : 7
    Points : 1
    Points
    1

    Par défaut Sauts dans les valeurs calculees

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

  2. #2
    Invité de passage
    Homme Profil pro
    Inscrit en
    février 2012
    Messages
    7
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Finlande

    Informations forums :
    Inscription : février 2012
    Messages : 7
    Points : 1
    Points
    1

    Par défaut

    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

  3. #3
    Responsable FAQ C++

    Avatar de Winjerome
    Homme Profil pro
    Inscrit en
    septembre 2009
    Messages
    8 325
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 28
    Localisation : France, Pyrénées Atlantiques (Aquitaine)

    Informations forums :
    Inscription : septembre 2009
    Messages : 8 325
    Points : 30 478
    Points
    30 478

    Par défaut

    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.
    Avant de poser votre question : FAQ, Tutoriels et recherche sur le forum
    Une erreur ? Messages d'erreur et avertissements
    "Ça ne marche pas" n'apporte aucune information utile permettant de vous aider. Expliquez clairement votre problème (erreurs entières, résultat souhaité vs obtenu...).

    En essayant continuellement on finit par réussir. Donc: plus ça rate, plus on a de chance que ça marche. - Jacques Rouxel
    L'expérience, c'est le nom que chacun donne à ses erreurs - Oscar Wilde
    Retrouvez-moi en soirée sur le chat de developpez.com - Pas de question technique par MP, Merci

  4. #4
    Invité de passage
    Homme Profil pro
    Inscrit en
    février 2012
    Messages
    7
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Finlande

    Informations forums :
    Inscription : février 2012
    Messages : 7
    Points : 1
    Points
    1

    Par défaut

    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.

  5. #5
    Responsable FAQ C++

    Avatar de Winjerome
    Homme Profil pro
    Inscrit en
    septembre 2009
    Messages
    8 325
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 28
    Localisation : France, Pyrénées Atlantiques (Aquitaine)

    Informations forums :
    Inscription : septembre 2009
    Messages : 8 325
    Points : 30 478
    Points
    30 478

    Par défaut

    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
    Avant de poser votre question : FAQ, Tutoriels et recherche sur le forum
    Une erreur ? Messages d'erreur et avertissements
    "Ça ne marche pas" n'apporte aucune information utile permettant de vous aider. Expliquez clairement votre problème (erreurs entières, résultat souhaité vs obtenu...).

    En essayant continuellement on finit par réussir. Donc: plus ça rate, plus on a de chance que ça marche. - Jacques Rouxel
    L'expérience, c'est le nom que chacun donne à ses erreurs - Oscar Wilde
    Retrouvez-moi en soirée sur le chat de developpez.com - Pas de question technique par MP, Merci

  6. #6
    Invité de passage
    Homme Profil pro
    Inscrit en
    février 2012
    Messages
    7
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Finlande

    Informations forums :
    Inscription : février 2012
    Messages : 7
    Points : 1
    Points
    1

    Par défaut

    Ok. J'ai compris.
    Merci.

+ Répondre à la discussion
Cette discussion est résolue.

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo