IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

Zéro d'une fonction


Sujet :

MATLAB

  1. #1
    Nouveau membre du Club
    Homme Profil pro
    Enseignant Chercheur
    Inscrit en
    Novembre 2008
    Messages
    115
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Gironde (Aquitaine)

    Informations professionnelles :
    Activité : Enseignant Chercheur

    Informations forums :
    Inscription : Novembre 2008
    Messages : 115
    Points : 34
    Points
    34
    Par défaut Zéro d'une fonction
    Bonjour,

    Je cherche à obtenir le zéro d'une fonction. Pour cela j'utilise la fonction 'fzero', or celle-ci nécessite de connaitre autour de quel point on cherche ce zéro et ce n'est pas le cas pour moi. Par défaut j'y vais un peu en force brute en mettant un intervalle le plus étendu possible, mais cela n'est en rien rigoureux. Y aurait-il une méthode alternative à l'emploi de 'fzero' ?

    La fonction sur laquelle je cherche mon zéro est du type :

    f(x) = -(a*(sin(x))^2) + (b*cos(x)) --> f(x) = 0


    Merci.

  2. #2
    Rédacteur

    Homme Profil pro
    Comme retraité, des masses
    Inscrit en
    Avril 2007
    Messages
    2 978
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 83
    Localisation : Suisse

    Informations professionnelles :
    Activité : Comme retraité, des masses
    Secteur : Industrie

    Informations forums :
    Inscription : Avril 2007
    Messages : 2 978
    Points : 5 179
    Points
    5 179
    Par défaut
    Salut!
    Il y a une multitude de méthodes pour résoudre ton problème, par exemple:
    1. Méthode de la dichotomie
    2. Méthode de la sécante
    3. Regula falsi
    4. Méthode de Ridder
    5. Méthode de Brent
    6. Méthode de Newton-Raphson
    7. etc.

    A toi de choisir.
    Jean-Marc Blanc
    Calcul numérique de processus industriels
    Formation, conseil, développement

    Point n'est besoin d'espérer pour entreprendre, ni de réussir pour persévérer. (Guillaume le Taiseux)

  3. #3
    FLB
    FLB est déconnecté
    Modérateur
    Avatar de FLB
    Homme Profil pro
    Ing. Aérospatiale
    Inscrit en
    Juin 2003
    Messages
    770
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 35
    Localisation : France

    Informations professionnelles :
    Activité : Ing. Aérospatiale
    Secteur : Aéronautique - Marine - Espace - Armement

    Informations forums :
    Inscription : Juin 2003
    Messages : 770
    Points : 1 185
    Points
    1 185
    Par défaut
    Ta fonction possède une infinité de zéro, la grande majorité des algos (tous?) ne te donneront qu'un zéro, généralement proche du point de départ de la recherche.
    Dans ton cas tu pourrais essayer de trouver la période de ton signal et chercher à l'intérieur de cet intervalle.
    Tu ne peux pas résoudre analytiquement sinon?

  4. #4
    Modérateur
    Avatar de le fab
    Homme Profil pro
    Développeur informatique
    Inscrit en
    Mars 2005
    Messages
    1 881
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 48
    Localisation : France, Isère (Rhône Alpes)

    Informations professionnelles :
    Activité : Développeur informatique
    Secteur : Industrie

    Informations forums :
    Inscription : Mars 2005
    Messages : 1 881
    Points : 3 429
    Points
    3 429
    Par défaut
    sinon avec symbolic math toolbox :
    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
    >> syms a b x
    >> f = -(a*(sin(x))^2) + (b*cos(x))
     
    f =
     
    b*cos(x) - a*sin(x)^2
     
    >> solve(f,x)
     
    ans =
     
     pi + acos((b + (4*a^2 + b^2)^(1/2))/(2*a))
     pi + acos((b - (4*a^2 + b^2)^(1/2))/(2*a))
     pi - acos((b + (4*a^2 + b^2)^(1/2))/(2*a))
     pi - acos((b - (4*a^2 + b^2)^(1/2))/(2*a))

  5. #5
    Nouveau membre du Club
    Homme Profil pro
    Enseignant Chercheur
    Inscrit en
    Novembre 2008
    Messages
    115
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Gironde (Aquitaine)

    Informations professionnelles :
    Activité : Enseignant Chercheur

    Informations forums :
    Inscription : Novembre 2008
    Messages : 115
    Points : 34
    Points
    34
    Par défaut
    Citation Envoyé par FLB Voir le message
    Ta fonction possède une infinité de zéro, la grande majorité des algos (tous?) ne te donneront qu'un zéro, généralement proche du point de départ de la recherche.
    Dans ton cas tu pourrais essayer de trouver la période de ton signal et chercher à l'intérieur de cet intervalle.
    Tu ne peux pas résoudre analytiquement sinon?
    Bonjour,

    Analytiquement la fonction est assez compliquée à résoudre car b dépend aussi de x ici !

  6. #6
    Modérateur
    Avatar de le fab
    Homme Profil pro
    Développeur informatique
    Inscrit en
    Mars 2005
    Messages
    1 881
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 48
    Localisation : France, Isère (Rhône Alpes)

    Informations professionnelles :
    Activité : Développeur informatique
    Secteur : Industrie

    Informations forums :
    Inscription : Mars 2005
    Messages : 1 881
    Points : 3 429
    Points
    3 429
    Par défaut
    si tu me donnes la relation entre b et x je te fais la résolution analytique avec symbolic math toolbox

  7. #7
    Nouveau membre du Club
    Homme Profil pro
    Enseignant Chercheur
    Inscrit en
    Novembre 2008
    Messages
    115
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Gironde (Aquitaine)

    Informations professionnelles :
    Activité : Enseignant Chercheur

    Informations forums :
    Inscription : Novembre 2008
    Messages : 115
    Points : 34
    Points
    34
    Par défaut
    Citation Envoyé par le fab Voir le message
    si tu me donnes la relation entre b et x je te fais la résolution analytique avec symbolic math toolbox
    Effectivement si je prend la solution réelle et positive de cette résolution cela devrait convenir.

    La relation est telle que : b=sqrt(c^2 - (a*sin(x))^2)

  8. #8
    Modérateur
    Avatar de le fab
    Homme Profil pro
    Développeur informatique
    Inscrit en
    Mars 2005
    Messages
    1 881
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 48
    Localisation : France, Isère (Rhône Alpes)

    Informations professionnelles :
    Activité : Développeur informatique
    Secteur : Industrie

    Informations forums :
    Inscription : Mars 2005
    Messages : 1 881
    Points : 3 429
    Points
    3 429
    Par défaut
    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
    >> syms a c x
    >> b=sqrt(c^2 - (a*sin(x))^2)
     
    b =
     
    (c^2 - a^2*sin(x)^2)^(1/2)
     
    >> f = -(a*(sin(x))^2) + (b*cos(x))
     
    f =
     
    cos(x)*(c^2 - a^2*sin(x)^2)^(1/2) - a*sin(x)^2
     
    >> solve(f,x)
     
    ans =
     
      asin(c*(1/(a^2 + c^2))^(1/2))
     -asin(c*(1/(a^2 + c^2))^(1/2))

Discussions similaires

  1. Trouver le zéro d'une fonction
    Par marion15 dans le forum C++
    Réponses: 10
    Dernier message: 12/06/2013, 16h15
  2. [Débutant] recherche du zéro d'une fonction par recherche numerique
    Par x2thezi dans le forum C#
    Réponses: 1
    Dernier message: 28/12/2011, 18h30
  3. Réponses: 4
    Dernier message: 04/04/2011, 13h15
  4. Trouver rapidement les zéros d'une fonction
    Par nahouto dans le forum MATLAB
    Réponses: 17
    Dernier message: 23/01/2009, 14h32
  5. Différentes méthodes de détermination du zéro d'une fonction
    Par Desperados dans le forum Mathématiques
    Réponses: 4
    Dernier message: 30/05/2006, 23h29

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