Publicité
+ Répondre à la discussion
Affichage des résultats 1 à 6 sur 6
  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 :
    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 :
    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 Jérôme
    Inscrit en
    septembre 2009
    Messages
    7 979
    Détails du profil
    Informations personnelles :
    Nom : Homme Jérôme
    Âge : 27
    Localisation : France, Pyrénées Atlantiques (Aquitaine)

    Informations forums :
    Inscription : septembre 2009
    Messages : 7 979
    Points : 29 372
    Points
    29 372

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

  5. #5
    Responsable FAQ C++

    Avatar de Winjerome
    Homme Profil pro Jérôme
    Inscrit en
    septembre 2009
    Messages
    7 979
    Détails du profil
    Informations personnelles :
    Nom : Homme Jérôme
    Âge : 27
    Localisation : France, Pyrénées Atlantiques (Aquitaine)

    Informations forums :
    Inscription : septembre 2009
    Messages : 7 979
    Points : 29 372
    Points
    29 372

    Par défaut

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

Liens sociaux

Règles de messages

  • Vous ne pouvez pas créer de nouvelles discussions
  • Vous ne pouvez pas envoyer des réponses
  • Vous ne pouvez pas envoyer des pièces jointes
  • Vous ne pouvez pas modifier vos messages
  •