p
u
b
l
i
c
i
t
é
publicité

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
    Expert Confirmé Sénior
    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 544
    Points
    30 544

    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
    Ma version améliorée de FireDVP

  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
    Expert Confirmé Sénior
    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 544
    Points
    30 544

    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
    Ma version améliorée de FireDVP

  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.

Discussions similaires

  1. Recherche d'une String dans les valeurs d'une HashMap
    Par 4rocky4 dans le forum Général Java
    Réponses: 2
    Dernier message: 04/05/2011, 15h24
  2. [AC-2003] « Saut » dans les enregistrements
    Par sabredebois dans le forum Modélisation
    Réponses: 4
    Dernier message: 13/01/2010, 13h46
  3. Remplacement des espaces dans les valeur d'un tableau
    Par calimero642 dans le forum Langage
    Réponses: 10
    Dernier message: 12/05/2006, 18h56
  4. [LDAP] Problèmes d'accent dans les valeurs
    Par Jack_le_Boulet dans le forum Bibliothèques et frameworks
    Réponses: 4
    Dernier message: 28/04/2006, 20h56
  5. [CR] doublon dans les valeurs
    Par Gaëtane dans le forum SAP Crystal Reports
    Réponses: 5
    Dernier message: 02/06/2004, 15h25

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