FORUM DE DISCUSSION SUR LE LANGAGE PANORAMIC

Développement d'applications avec le langage Panoramic
 
AccueilAccueil  FAQFAQ  RechercherRechercher  S'enregistrerS'enregistrer  MembresMembres  GroupesGroupes  Connexion  
Derniers sujets
» Problème de compilateur.
par Pedro Alvarez Hier à 20:01

» COMPILATEUR V 0.9 beta 8 du 18 aout 2017
par Pedro Alvarez Hier à 19:54

» Pb 17 (en analyse): ITEM_SELECT ne fonctionne pas
par Jack Jeu 17 Aoû 2017 - 19:26

» Compilateur FBPano
par Mike Jeu 17 Aoû 2017 - 13:52

» un nouveau editeur panobasic
par Jean Claude Jeu 17 Aoû 2017 - 10:18

» Le compilateur.
par Pedro Alvarez Jeu 17 Aoû 2017 - 8:36

» Pb 16 (en analyse): ON_CLOSE plante à l'exécution
par Jack Mer 16 Aoû 2017 - 20:00

» Pb 15 (en analyse): TIMER_ON plante à l'exécution
par Jack Mer 16 Aoû 2017 - 19:58

» KGF_dll - nouvelles versions
par Yannick Dim 13 Aoû 2017 - 17:35

» probleme d'outil
par Yannick Dim 13 Aoû 2017 - 17:32

» Carte de France des régions
par Yannick Sam 12 Aoû 2017 - 21:33

» Pb 14 (en analyse): PRINT_LOCATE plante à l'exécution
par Jack Ven 11 Aoû 2017 - 22:37

» Petit avertissement [Mots réservés]
par papydall Ven 11 Aoû 2017 - 13:45

» Distances sur plan
par JL35 Jeu 10 Aoû 2017 - 21:29

» Tracé : Triangle, Carrée, Dents de scie, Sinusoïde redressée
par papydall Jeu 10 Aoû 2017 - 14:52

Navigation
 Portail
 Index
 Membres
 Profil
 FAQ
 Rechercher
Rechercher
 
 

Résultats par :
 
Rechercher Recherche avancée
Août 2017
LunMarMerJeuVenSamDim
 123456
78910111213
14151617181920
21222324252627
28293031   
CalendrierCalendrier

Partagez | 
 

 AX^3 + BX^2 + CX + D = 0

Voir le sujet précédent Voir le sujet suivant Aller en bas 
AuteurMessage
papydall

avatar

Nombre de messages : 5508
Age : 67
Localisation : Moknine (Tunisie) Entre la chaise et le clavier
Date d'inscription : 03/03/2012

MessageSujet: AX^3 + BX^2 + CX + D = 0   Sam 16 Mar 2013 - 17:38

Après l’équation du second degré, voici celle du troisième.

Soit l’équation du 3ème degré suivante :
ax3 + bx2 + cx + d = 0 ( forme 1)

D’une manière générale, on ne sait pas comment résoudre directement une telle équation.
Mais on sait comment résoudre une équation, de la forme :
X3 + pX + q = 0 ( forme 2)

Or, toute équation du 3ème degré de la forme 1, peut être mise sous la forme 2.
On résout cette dernière et on remonte pour avoir la solution de la forme 1

Un peu de théorie
La méthode de Cardan
L’astuce pour résoudre une équation de la forme 1 consiste à faire sauter le terme du second degré en faisant un changement de variable :

Posons X = x + b/(3*a) c.a.d x = X-b/(3*a)
En remplaçant x dans la forme 1 par sa nouvelle valeur et en développant puis en divisant tous les termes par a (a est différent de zéro), on obtient :
X3 + pX + q = 0
avec : p = (c/a) – (b*b/(3*a*a)) et q = 2*b*b*b/(27*a*a*a) – (b*c/(3*a*a)) + d/a

Ensuite on calcule :
Delta = q*q/4 + p*p*p/27

Si Delta > 0 on a une racine réelle et 2 racines complexes conjuguées
Si delta = 0 on a 3 racines réelles ( soit une racine réelle simple et une racine réelle double, soir une seule racine réelle triple).
Si delta < 0 on a 3 racines réelles distinctes.

Les formules sont assez compliquées : des racines cubiques qui font intervenir des racines carrées, des nombres complexes, des fonctions trigonométriques et trigonométriques inverses.
Il est inutile de développer tout ça ici.

Le programme ci-après permet de résoudre toute équation du 3ème degré.
Il vous demande d’entrer les coefficients a,b,c,d (tous réels).
Le coefficient a doit être différent de 0.
Le programme contrôle ces entrées et refuse toute valeur incorrecte.
La solution s’affiche immédiatement.
Vous pouvez alors, soit faire un autre calcul, soit quitter.

Code:
' ******************************************************************************

'    3DEGRE.BAS par Papydall

' Résolution de l'équation du troisième dégré ax^3 + bx^2 + cx + d = 0

' ******************************************************************************

' D'après le théorème fondamental de l'algèbre des polynômes : tout polynôme de
' dégré N admet exactement N solutions (on dit aussi racines ou zéros). Ces
' racines sont réelles ou complexes, simples ou multiples et si elles sont
' complexes, elles sont conjuguées deux-à-deux.
' Ce programme permet de résoudre un polynôme du troisième dégré

' L'équation du troisième dégré admet donc trois racines : soit réelles simples
' ou multiples , soit complexes conjuguées.

' RAPPEL : un nombre complexe s'écrit sous forme z = x + i*y , avec x et y réels
' et i nombre imaginaire dont le carré vaut -1 (i² = -1).
' x s'appelle partie réelle du nombre complexe z ; y s'appelle partie imaginaire
' du nombre complexe z.
' Si y = 0 , le nombre est réel. Si x = 0 , le nombre est imaginaire pur.
' Deux nombres complexes conjugués sont deux nombres complexes dont leurs
' parties imaginaires sont de signes contraires.

' ******************************************************************************
Run()
end
' ******************************************************************************
SUB Run()
    label suite    : ' sous-programme pour faire un autre calcul
    label quit    : ' sous-programme pour quitter l'application
    dim a,b,c,d    : ' les coefficients de l'équation
    dim x1,x2,x3  : ' les racines réelles
    dim re,im      : ' Les racines complexes (partie réelle et partie imaginaire)
    dim complexe%  : ' flag pour indiquer si les racines sont complexes (=1) ou réelles (=0)
    color 0,0,0,0  : font_color 0,0,255,0 : font_size 0,12 : font_bold 0
    width 0,width(0)*2
    caption 0,"Résolution de l'équation du troisième dégré : AX^3 + BX^2 + CX + D = 0"
    button 1 : top 1,250 : left 1,50  : width 1, 100 : on_click 1,suite
    button 2 : top 2,250 : left 2,300 : width 2, 80  : on_click 2,quit
    caption 1,"Autre calcul" : caption 2,"Quitter"
    Demarrer()
END_SUB
' ******************************************************************************
SUB demarrer()
    cls  : hide 1 : hide 2
    saisie() : resol_3degre(a,b,c,d) : Format() : Affiche()
    show 1 : show 2
END_SUB
' ******************************************************************************
' Procédure de saisie des coefficients a,b,c,d
SUB saisie()
    DIM_LOCAL r$
    ' Saisie de coefficient A
    repeat      : ' Cette boucle assure que le coefficient A soit différent de 0
        repeat  : ' Cette boucle assure que le coefficient soit numérique
              r$ = message_input$("Entrer le paramètre A (différent de zéro)","A = ","")
        until numeric(r$) = 1
    until val(r$) <> 0
    A = val(r$)
    ' Saisie du coefficient B
    repeat
      r$ = message_input$("Entrer le paramètre B","B = ","")
    until numeric(r$) = 1
    B = val(r$)
    ' Saisie du coefficient C
    repeat
      r$ = message_input$("Entrer le paramètre C","C = ","")
    until numeric(r$) = 1
    c = val(r$)
    ' Saisie du coefficient D
    repeat
      r$ = message_input$("Entrer le paramètre D","D = ","")
    until numeric(r$) = 1
    d = val(r$)
END_SUB
' ******************************************************************************
' Procédure de calcul de la solution
SUB resol_3degre(AX3,BX2,CX,D)
    DIM_LOCAL p,q,delta,radic1,radic2,phi,pi
    pi = 4*atn(1) : complexe% = 0
    p = (c/a)-(power(b,2)/(3*power(a,2)))
    q = (d/a)-(b*c)/(3*(power(a,2)))+(2*(power(b,3))/(27*power(a,3)))
    delta = power(q,2)/4 + power(p,3)/27
    if delta > 0  : ' Une racine réelle et 2 racines complexes conjuguées
      complexe% = 1
      radic1 = (0-q/2)+ sqr(delta) : radic2 = (0-q/2)- sqr(delta)
      x1 = sgn(radic1)*power((abs(radic1)),1/3)+ sgn(radic2)*power(abs(radic2),1/3)
      x1 = x1 - b/(3*a)
      re = 0-0.5*sgn(radic1)*(power(abs(radic1),1/3)+ sgn(radic2)*power(abs(radic2),1/3))
      re = re-b/(3*a)
      im = 0.5*sqr(3)*sgn(radic1)*(power(abs(radic1),1/3)- sgn(radic2)*power(abs(radic2),1/3))
      exit_sub : ' les calculs sont faits, on quitte la procédure
    end_if
    if delta = 0  : ' Trois racines réelles
      if  p = 0  : ' une racine réelle triple
          x1 = sgn(d/a)* power(abs(d/a),1/3) : x2 = x1 : x3 = x1
      else  : ' une racine simple et une racine double
          x1 = 3*q/p : x1 = x1 - b/(3*a)
          x2 = 0-3*q/(2*p) : x2 = x2- b/(3*a)
          x3 = x2
      end_if
      exit_sub : ' les calculs sont faits, on quitte la procédure
    end_if
    ' ici delta est forcément < 0 et l'équation admet 3 racines réelles distinctes
    phi = acos(((3*q)/(2*p))*sqr(0-3/p))
    x1 = 2*sqr(0-p/3)*cos(phi/3) : x1 = x1 - b/(3*a)
    x2 = 2*sqr(0-p/3)*cos((phi+2*pi)/3) : x2 = x2 - b/(3*a)
    x3 = 2*sqr(0-p/3)*cos((phi+4*pi)/3) : x3 = x3 - b/(3*a)

END_SUB
' ******************************************************************************
' Procédure de formatage de la solution
SUB Format()
    dim_local a$,b$,c$,d$
    dim_local epsilon : epsilon = power(10,0-10)
    if a = 1  then a$ = "x^3 " : else : a$ = str$(a)+"x^3 " : ' Pour ne pas afficher 1x^3
    if a = -1 then a$ = "-x^3" : ' pour ne pas afficher -1x^3
    if b = 1  then b$ = "+x^2 " : else : b$ = str$(b)+"x^2 " : ' Pour ne pas afficher 1x^2
    if b = -1 then b$ = "-x^2 " : ' pour ne pas afficher -1x^2
    if b > 1  then b$ = "+" + str$(b)+"x^2 "
    if b = 0  then b$ = ""
    if c = 1  then c$ = "+x " : else : c$ = str$(c)+"x " : ' Pour ne pas afficher 1x
    if c = -1 then c$ = "-x " : ' pour ne pas afficher -1x
    if c > 1  then c$ = "+" + str$(c)+"x "
    if c = 0  then c$ = ""
    if d = 0  then d$ = " = 0" : else : d$ = str$(d)+" = 0 :"
    if d > 0  then d$ = "+" + str$(d)+" = 0"

    x1 = x1 + sgn(x1)*epsilon : x1 = int(x1*100000+0.5)/100000
    x2 = x2 + sgn(x2)*epsilon : x2 = int(x2*100000+0.5)/100000
    x3 = x3 + sgn(x3)*epsilon : x3 = int(x3*100000+0.5)/100000

    re = re + sgn(re)*epsilon : re = int(re*100000+0.5)/100000
    im = im + sgn(im)*epsilon : im = int(im*100000+0.5)/100000
    print_locate 50,20 : print "La solution de : "+a$+b$+c$+d$
END_SUB
' ******************************************************************************
' Procédure d'affichage de la solution
SUB Affiche()
    print_locate 50,20  : ' on positionne le curseur pour l'affichage
    if complexe% = 0  : ' Les racines sont réelles
      print_locate 50,50
      if x1 = x2    : ' une racine triple
          print "Une racine triple"
          print_locate 50,80 : print "X1 = X2 = X3 = ";x1
      else
          if x2 = x3    : ' une racine double
            print "Une racine réelle simple et une racine réelle double :"
            print_locate 50,80  : print "X1 = "; x1
            print_locate 50,100 : print "X2 = X3 = " ; x2
          else : ' les trois racines sont réelles et distinctes
            print "Trois racines réelles distinctes :"
            print_locate 50,80  : print "X1 = " ; x1
            print_locate 50,100 : print "X2 = " ; x2
            print_locate 50,120 : print "X3 = " ; x3
          end_if
      end_if
    else  : ' une racine réelle et deux racines complexes conjuguées
      print_locate 50,50  : print "une racine réelle et deux racines complexes conjuguées"
      print_locate 50,80  : print "X1 = ";x1
      print_locate 50,100 : print "X2 = "; re;" + ";ABS(im); "i"
      print_locate 50,120 : print "X3 = "; re;" -  ";abs(im); "i"
    end_if
END_SUB

' ******************************************************************************
' sous-programme pour un autre calcul
suite:
  demarrer()
return
' ******************************************************************************
' sous-programme pour quitter l'application
quit:
  if message_confirmation_yes_no("Vous voulez vraiment quitter ?") = 1 then terminate
return
' ************************** FIN ***********************************************
Revenir en haut Aller en bas
Voir le profil de l'utilisateur http://papydall-panoramic.forumarabia.com/
 
AX^3 + BX^2 + CX + D = 0
Voir le sujet précédent Voir le sujet suivant Revenir en haut 
Page 1 sur 1

Permission de ce forum:Vous ne pouvez pas répondre aux sujets dans ce forum
FORUM DE DISCUSSION SUR LE LANGAGE PANORAMIC :: PANORAMIC 32 bits :: Vos sources, vos utilitaires à partager-
Sauter vers: