Accueil > > > RECHERCHE DE ZÉRO D'UNE FONCTION PAR DICHOTOMIE
RECHERCHE DE ZÉRO D'UNE FONCTION PAR DICHOTOMIE
Information sur la source
Description
Le titre est assez explicite, c'est juste 60 petites lignes pour trouver le zéro d'une fonction par dichotomie.
Il s'agit d'un processus itératif où l'intervalle de recherche du zéro est divisé par deux à chaque étape,
ce qui le rend assez efficace.
J'ai essayé d'expliquer le fonctionnement avec des commentaires, ça me parait assez compréhensible...
Source
- #include <stdio.h>
- #include <stdlib.h>
- #include <math.h> //pour cos, sin ou autre...
-
- int dichotomie ( double (*)(double), double, double, double*, double); //fonction
-
- int main (int argc, char *argv[])
- {
- double zero , //zero recherché
- a, b, //intervalle [a,b]
- precision ; //précision voulue
-
- printf("Exemple : \n");
-
- if ( dichotomie (sin, -1.0, 1.0, &zero, 1.0e-2) != -1 ) {// exemple avec sin dans [-1;1]
- printf("Zero de la fonction sin dans l'intervalle [-1;1] a 1e-2 pres : \n%le\n", zero);
- }
- else printf("Erreur.\n");
-
- if ( dichotomie (sin, -1.0, 2.0, &zero, 1.0e-2) != -1 ) {// exemple avec sin dans [-1;2]
- printf("Zero de la fonction sin dans l'intervalle [-1;2] a 1e-2 pres : \n%le\n", zero);
- }
- else printf("Erreur.\n");
-
- system("PAUSE");
- return 0;
- }
-
- int dichotomie ( double (*f)(double), double a, double b, double* zero, double precision)
- //f : fonction dont on veut trouver le zero
- //a et b : intervalle
- //zero : zero
- //precision : precision voulue
- {double m, //milieu de l'intervalle courant
- fm, //valeur de f(m)
- fa, fb ; //valeur de f(a) et de f(b)
-
- fa = (*f)(a);
- fb = (*f)(b);
- if (fa*fb>=0 || a>=b ) {printf("Intervalle incorrect.\n"); return -1;} //Intervalle incorrect
-
- while ( b-a > precision )
- {
- m = (b+a) / 2 ;
- fm = (*f)(m) ;
- if (fm==0) break; //zero trouvé
-
- if (fa*fm < 0) {
- b = m;
- fb = fm;
- }
- else {
- a = m;
- fa = fm;
- }
-
- }
-
- *zero = m;
- return 0;
- }
#include <stdio.h>
#include <stdlib.h>
#include <math.h> //pour cos, sin ou autre...
int dichotomie ( double (*)(double), double, double, double*, double); //fonction
int main (int argc, char *argv[])
{
double zero , //zero recherché
a, b, //intervalle [a,b]
precision ; //précision voulue
printf("Exemple : \n");
if ( dichotomie (sin, -1.0, 1.0, &zero, 1.0e-2) != -1 ) {// exemple avec sin dans [-1;1]
printf("Zero de la fonction sin dans l'intervalle [-1;1] a 1e-2 pres : \n%le\n", zero);
}
else printf("Erreur.\n");
if ( dichotomie (sin, -1.0, 2.0, &zero, 1.0e-2) != -1 ) {// exemple avec sin dans [-1;2]
printf("Zero de la fonction sin dans l'intervalle [-1;2] a 1e-2 pres : \n%le\n", zero);
}
else printf("Erreur.\n");
system("PAUSE");
return 0;
}
int dichotomie ( double (*f)(double), double a, double b, double* zero, double precision)
//f : fonction dont on veut trouver le zero
//a et b : intervalle
//zero : zero
//precision : precision voulue
{double m, //milieu de l'intervalle courant
fm, //valeur de f(m)
fa, fb ; //valeur de f(a) et de f(b)
fa = (*f)(a);
fb = (*f)(b);
if (fa*fb>=0 || a>=b ) {printf("Intervalle incorrect.\n"); return -1;} //Intervalle incorrect
while ( b-a > precision )
{
m = (b+a) / 2 ;
fm = (*f)(m) ;
if (fm==0) break; //zero trouvé
if (fa*fm < 0) {
b = m;
fb = fm;
}
else {
a = m;
fa = fm;
}
}
*zero = m;
return 0;
}
Historique
- 29 novembre 2005 18:43:10 :
- Mise à jour des mots clés pour le moteur de recherche comme proposé dans la lettre d'info...
- 09 février 2008 20:15:28 :
- /
- 20 juillet 2009 14:37:01 :
- Modification de la description
Sources du même auteur
Sources de la même categorie
Commentaires et avis
Discussions en rapport avec ce code source dans le forum
Recherche fou de C pour réaliser une petite dll [ par loskiller62 ]
BonjourProgrammeur VB, j'ai un besoin (relativement urgent) pour une dll. Il s'agit d'une fonction réalisable qu'en Cpp mais dont j'ai besoin pour mon
Recherche algo de la fonction ulltoa() [ par akboy7015 ]
Actuellement, je développe sur VC7 et je travail beaucoup avec des __int64. Le problème est que la fonction ulltoa() n'a pas l'air d'exister.Bien ente
Aide pour une fonction recherche avec suivant [ par oeil123 ]
je n'arrive pas a trouver la fonction que microsoft utilise pour faire le dialog recherche avec le fameux bouton suivant qui m'interesse.En fait je vo
recherche de fonction [ par jeromedu94 ]
je recherche une fonction qui dirait si il y a eu un appui sur une touche ou pas
Problème Recherche + Helpppppppppppp [ par emmanuel9 ]
Bonjour, En faite j'ai une fenêtre destiné à la recherche et on peut recherche une personne sur son nom, prénom, sur la date de
Recherche méthodes pour faire pivoter le bureau de 0, 90, 180, 270 degrés [ par scrumby ]
Bonjour à tous,comme indiqué dans le sujet je cherche désespérément une fonction ou une méthode qui me permettrait de fa
retour d'une fonction [ par noulii ]
J'ai utilisé la fonction suivante, venant de "bruno" que j'ai trouvé sur ce site :HANDLE hfind;WIN32_FIND_DATA wfd;// Initialise la recherchehfind = F
Fonction recherche en C++ [ par Ciol ]
Salut a tous J ai bien lu les differents messages deposés par les 3 habitués du Chat et je ne souhaite donc pas trouver une solution tte faite pour m
recherche de l’optimum global de la fonction de Rastringin [ par angelmimitta ]
bonjour à tous, enfait j'ai un petit projet à faire et j'ai besoin d'aide, l"nnoncé est le suivant: Implémenter le recuit simulé pour la recherche de
Exercice sur Pointeurs, classes, héritage... [ par snoozriton ]
Bonjour à tous, étant en train d'étudier le C++, j'ai quelques soucis avec un exercice. C'est un classique dont je vous donne l'énoncé : [b]1. Défini
|
Derniers Blogs
UNE JOLIE-HORLOGE ET PAS QU'UN PEU !UNE JOLIE-HORLOGE ET PAS QU'UN PEU ! par neodante
Pour les possesseurs d'iPhone, ça y est Bijin Tokei - qui se traduit littéralement en Français par " Jolie Horloge " - est arrivé et GRATUITEMENT s'il vous plaît ! Après la version Tokyo, Hokkaido, night club, racing, Gal, "pour les mademoiselles'", . voi...
Cliquez pour lire la suite de l'article par neodante TECHDAYS PARIS 2010 : CONNECTEZ VOS DONNéES à SHAREPOINT 2010 AVEC LES BUSINESS CONNECTIVITY SERVICESTECHDAYS PARIS 2010 : CONNECTEZ VOS DONNéES à SHAREPOINT 2010 AVEC LES BUSINESS CONNECTIVITY SERVICES par ROMELARD Fabrice
Animé par: Gaetan Bouveret et Julien Chomarat Business Connectivity Services (BCS) est dans SharePoint 2010 la version 2 de Business Data Catalog (BDC dans SharePoint 2007). Il s'agit de la solution permettant de visualiser des données provenan...
Cliquez pour lire la suite de l'article par ROMELARD Fabrice [DIVERS] SUIVRE VOS SéRIES PRéFéRéS SUR LA TOILE[DIVERS] SUIVRE VOS SéRIES PRéFéRéS SUR LA TOILE par orion
Comme de nombreux geek, je suis un grand amateur de série TV et je rate régulièrement des épisodes de mes séries préférés. Une solution s'offre à vous avec ce merveilleux site : Tv Gorge - www.tvgorge.com Moteur de recherche à l'appui, vous pouvez ...
Cliquez pour lire la suite de l'article par orion TECHDAYS PARIS 2010 : LA BI DANS SHAREPOINT 2010TECHDAYS PARIS 2010 : LA BI DANS SHAREPOINT 2010 par ROMELARD Fabrice
Animé par: Vincent Bellet et Baptiste Giraudier La BI dans SharePoint 2010, Les nouveaux services d'application dans SP2010 et SQL Server Reporting services 2008 R2. La BI dans SharePoint est généralisée pour tous afin de permettre à tous les coll...
Cliquez pour lire la suite de l'article par ROMELARD Fabrice
Forum
RE : WIN APIRE : WIN API par racpp
Cliquez pour lire la suite par racpp WIN APIWIN API par omarino_007
Cliquez pour lire la suite par omarino_007
Logiciels
DB-MAIN (9.1.0)DB-MAIN (9.1.0)DB-MAIN is a data-modeling and data-architecture tool. It is designed to help developers and anal... Cliquez pour télécharger DB-MAIN Xilisoft DPG Convertisseur (5.1.37.0120)XILISOFT DPG CONVERTISSEUR (5.1.37.0120)Xilisoft DPG Convertisseur offre aux fans de Nintendo DS une bonne solution leur permettant de dé... Cliquez pour télécharger Xilisoft DPG Convertisseur GraphicsGale (2.01.01)GRAPHICSGALE (2.01.01)GraphicsGale est un logiciel de PixelArt avec de nombreuse fonctionnalités permettant de réalisé ... Cliquez pour télécharger GraphicsGale Architecte 3D (Platinum 2010)ARCHITECTE 3D (PLATINUM 2010)Architecte 3D Platinium vous permet de concevoir facilement les plans votre future maison, de l'é... Cliquez pour télécharger Architecte 3D TeamViewer 5 (TeamViewer 5)TEAMVIEWER 5 (TEAMVIEWER 5)Dépanner un ami,expliquer une manipulation devient un jeu d'enfant.
Prise en main d'un autre ord... Cliquez pour télécharger TeamViewer 5
|