équation à 20 inconnue

Question

Bonjour,

Voilà mon souci : je ne suis ni programmeur ni mathématicien, et je cherche à résoudre une équation telle que celle-ci : 80 a + 148 b + 464 c + ... + 2498 m = 13378

Comme outils, je n'ai en tout et pour tout que excel... ^__^ Hem...

Pourquoi faire ? Je participe à un jeu de simulation en ligne, dans lequel on peut espionner les forces ennemies, ce qui nous donne un total de points. Et le but serait de pouvoir connaître la composition "précise" de son armée. J'essaye de faire ça "à la main", à tâton, mais ce n'est pas évident. 

J'y ai réfléchi un petit peu, et je pense que des solutions multiples sont envisageables pour chaque équation. Je pensais qu'on pouvait peut-être faire appel à un algorithme qui explorerait toutes les possibilités. A part ça, je n'ai pas la plus petite idée de la manière de procéder. J'ai lu sur une autre page de votre forum que l'on pouvait réduire l'équation de manière à arriver à une seule inconnue et que ce n'était pas facile. Est-ce possible sous excel ? Faisable par un débutant motivé ?

Merci de m'avoir lu.

Yann

mamiemando · Accepted Answer

Le truc c'est que mathématiquement parlant il faut au moins 20 points de mesures linéairement indépendants pour résoudre le système (soit 20 équations, 1 par relevé de point). Et ça suppose effectivement que ce soit bien une équation linéaire. Pour s'en assurer l'idéal serait de pouvoir relever des points en ne faisant varier qu'une inconnue et vérifier qu'ils sont bien alignés.

Pour résoudre ton système de nombreuses méthodes sont possibles, mais à la main ça risque d'être un peu galère. Tu peux par exemple utiliser un outil comme matlab ou mapple pour résoudre le système (version étudiante). 

Pour excel aucune idée.

Bon courage

pavarotti_39 · Answer

Faites un système de Kramer, avec ça vous résoudrez toutce que vous voudrez.

ADM · Answer

Salut,
Moi, je suis mathématicien ;) Et de toute évidence une seule équation à 20 inconnues donne des solutions multiples ! Tu ne trouveras jamais une solution unique.
Il te faudrait 20 équations à 20 inconnues pour espérer en avoir une... et encore il y a des exceptions si au moins 2 équations sont liées entre elles (càd dépendantes). 
Après, il est effectivement faisable de lier une seule inconnue aux 19 autres, c’est d’ailleurs la solution ’’hard’’ de passer de 20 équations/20 inconnues à 19équations/19 inconnues (...) Ainsi de suite jusqu’à ne trouver qu’1 équation avec 1 inconnue, puis de remonter la chaîne à l’envers pour les différents résultats.
Dans ton cas, avec 1 équation, il y auras donc des solutions multiples !

Cernunnos · Answer

Je vais jeter un coup d'oeil sur les logiciels proposés, en espérant que ça fonctionne.

Pour les solutions multiples, ce n'est pas dit : Les multiplicateurs sont "fixes" (ainsi l'équation demande combien d'unités à 80 points, combien d'unités à 148 points, etc...). Mais une valeur-cible de 160, par exemple, il n'y a qu'une solution à léquation : c'est 2 x 80.

Quand je disais "à la main", c'était sous excel, avec une procédure essai-erreur et approximation : je tapais des hypothèses et essayait de retrouver le résultat en enlevant 1 par ci, en ajoutant 2 par là... Je n'ai pas pensé à une résolution graphique et j'avoue ne pas réussir à me représenter la chose.

Quand au système de Kramer, je suis désolé pour vous (et pour lui ^^) mais je ne crois pas en avoir entendu parler pendant mes études (ou alors sous un autre nom)

Cernunnos · Answer

Ca veut dire que je vous embête avec mes questions ? ^__^
Nan, désolé...

Est-ce qu'un programme existe qui suive le même raisonnement que moi ? Qui fasse une approximation, puis qui l'affine en changeant les paramètres a+2 b-1 pour approcher la valeur-cible (x) ? Puis qui, une fois une solution trouvée pour l'équation, recommence en partant d'autres valeurs ?

Je ne sais pas si je m'exprime assez clairement...

Nassim5555 · Answer

salut , 
moi j'ai utilisé MATLAB ,il peut tout faire ce grand !

ADM · Answer

Salut !
Pour moi je reste persuadé qu'1 équation à 20 inconnues possède des solutions multiples !
Le système de Kramer fonctionne très bien, mais avec le meme nombre d'équations que d'inconnues, ce qui n'est pas le cas ici je pense... (à moins d'avoir mal compris).
Sinon, oui, j'appuis que Mapple est très bien (et facile à comprendre). Pour moi, Matlab est un peu plus délicat... mais c'est mon avis ! En tout cas, ils permettent de résoudre les Kramer, quand il en est question bien sur ;)

ADM · Answer

Mais Matlab est très grand aussi, bien sur Nassim5555 !!! c'est juste un peu plus difficile pour un débutant je pense...

Nassim5555 · Answer

oui c'est vrai faut connaitre son langage qui s'inspire du langage C ,comme moi je suis informaticien alors je me suis rapidement habitué ,mais quand on veut on peut ,faut juste un peu de volonté ;^) Bon courage !

Cernunnos · Answer

J'ai vu la version étudiante de maple à 99 $, à peu près 70 euros ? *aie*

Ca n'est pas possible. Je suis sans doute naïf, mais je pensais trouver un logiciel gratuit. Matlab est payant lui aussi. Et je ne sais pas si un clone pourrait faire l'affaire.

Dernière solution (?) : trouver une bonne âme disposant du dit logiciel et qui pourrait se lancer dans ce défi, pas très gratifiant, sans doute, mais oh combien glorieux ! Bien sûr je fournirai tous les détails et explications nécessaires... Peut-être un étudiant, voire un professeur ? En tant qu'exercice ? Je ne sais pas...

Un grand merci en tous cas à tou(te)s ceux et celles qui sont intervenus pour tenter de m'aider jusqu'ici.

Cernunnos · Answer

J'ai vu la version étudiante de maple à 99 $, à peu près 70 euros ? *aie*

Ca n'est pas possible. Je suis sans doute naïf, mais je pensais trouver un logiciel gratuit. Matlab est payant lui aussi. Et je ne sais pas si un clone pourrait faire l'affaire.

Dernière solution (?) : trouver une bonne âme disposant du dit logiciel et qui pourrait se lancer dans ce défi, pas très gratifiant, sans doute, mais oh combien glorieux ! Bien sûr je fournirai tous les détails et explications nécessaires... Peut-être un étudiant, voire un professeur ? En tant qu'exercice ? Je ne sais pas...

Un grand merci en tous cas à tou(te)s ceux et celles qui sont intervenus pour tenter de m'aider jusqu'ici.

mamiemando · Answer

Tu as un matlab gratuit : scilab

Faust0 · Answer

Bonjour

1)Les inconnues appartiennent elles à l'ensemble des entiers positifs?
ou non

2)Y a t il une borne haute et basse pour chaque inconnues?
par exemple "a"  est un entier qui ne varie qu'entre 100 et 150

3) les coefficients (80, 148, 464..) sont ils tous positifs?

Cernunnos · Answer

@ mamie : Je vais chercher le scilab.

@ Fausto 
1) oui
2) entre 0 et... pas de limite haute, mais 1000 me semble assez grand pour couvrir tous les cas de figure.
3) oui, tous positifs (et pairs, mais ça, c'est un détail)

Faust0 · Answer

ah jusqu'à 1000 
ben c'est pas trop compliqué
tu écris un petit programme en pascal, en fortran, en basic bref en ce que tu connais
tu fais 20 boucles imbriquées
tu testes pour voir si à un moment tu obtiens 13378 avec ta formule
dans ce cas tu imprimes les valeurs des indices de boucle
et tu continues pour les éventuelles valeurs suivantes
bestial mais efficace et un peu long mais c'est du temp machine c'est pour cela qu'il faudrait borner au mieux

quant à calculer le determinant d'un matrice de rang 20 je laisse ça aux matheux d'autant plus qu'on n'est pas dans un système de 20 équations à 20 inconnues.

Cernunnos · Answer

Alors, j'essaye d'imaginer ce programme

au milieu la fameuse formule, et un système de boucle ou on implémente les variables les unes après les autres de façon à couvrir tout le "spectre des possibilités" ? Ca ne fait pas 1000 exposant 20 solutions ? 10 exp 23 cas de figure ? Même en imaginant qu'on devrait trouver la solution dans le 1er 100ème de la recherche, ça fait quand même 10 puissance 19 calculs..., allez, on arrondit à 10 puissance 15 ^^

Sur la manière dont fonctionne le programme, je lui demande de faire a=1, b=0 c=0 puis de comparer le résultat à la cible. Si ce n'est pas bon, a=a+1 et on recommence. Quand a x 80 > 13378, on recommence avec a=0 et b=1 et ainsi de suite... Quand un résultat est bon, les valeurs a b c d... sont affichées.

C'est ça ? Je vais chercher un programme de basic, alors...

Faust0 · Answer

c'est bien cela et c'est d'ailleurs ce que tu essayes de faire à la main et au pifomètre.
mais comme tu l'as très bien vu le nombre d'iterations est ENORME c'est pour cela qu'il faudra bien borner tes variables
Si tes paramètres sont bien positifs comme tu l'as dit et que tu fais bien une somme.
"a" est compris entre 0 et 13378/80 soit  153
pour "m" l'amplitude est limitée entre 0 et 5
etc...

Cernunnos · Answer

En fait j'ai fait un peu de basic étant ado, je devrais pouvoir y arriver, je suppose.

Mais le bornage se fera de lui-même, étant donné que dès que le résultat de l'équation est supérieur à 13378, il recommence avec d'autres valeurs...

Sinon, pour scilab, apparemment il n'existe que sous Linux...

Faust0 · Answer

oui je pense que ça ne doit pas poser de problème de programmation
tu pourras utiliser just basic qui est gratuit
http://www.justbasic.com/
( je sens qu'on va se faire pourrir pour oser utiliser une pareille relique)

tu nous diras les resultats

cernunnos · Answer

Bé ?? Elle est où ma réponse ?
Zut, j'ai du faire une fausse manip, je recommence :-(

J'ai installé VisualBasic, mais je n'ai pas encore eu le temps de plonger dedans.

Ce qu'à fait JvDo est génial, c'est tout à fait ça. Par contre je suis surpris du nombre de solutions ! Peut-être que j'aurais du m'en douter...

Voilà les coefficients réels :
80 158 464 908 688 1612 1326 6042 5164 5180 8052 2724 1452 196 260 518 650 2040 1694 2498
Effectivement une pondération pourrait être faite pour classer les unités selon leur intérêt (certaines unités sont chères et de peu d'intérêt)

Une autre manière de limiter le nombre de variables est possible : dans le jeu on peut utiliser un espion qui nous informe du plus haut niveau d'unité chez l'ennemi.

Au niveau des calculs, y aurait-il un gain à simplifier tous les éléments par 2 (puisque tous les nombres sont pairs) ? Je suppose que ça doit être négligeable...

Vraiment, autant de résultats, ça risque effectivement d'être peu exploitable, sachant que le nombre donné, 13 378, était un exemple, et que ça peut aller très loin (100 000, 200 000...)

ADM · Answer

1 équation à 20 inconnues... votre résultat confirme ce que je disais sur une infinité de solutions ;) En ne gardant que les nombres entiers, effectivment il doit y en avoir un nombre fini, mais bonne chance quand même pour la quantité !

JvDo · Answer

Bonsoir,


pour tes coefficients, ça donne :

niveaux traités :___Nb_____durée___Nb solutions___Nb combinaisons
_____20 à 7_____14______8mn_______149________147 253   
_____20 à 6_____15_____31mn_______649 ________579 436   
_____20 à 5_____16___1h43mn______2,832_______2 343 533   



voici le code : 
(je précise que je n'ai pas fait dans la dentelle)
Sur une feuille, mettre en A1:T1 et A4:T4 les libellés I01 --> I20 pour rendre lisible l'affichage des résultats.
Puis, lancer la macro une_équation_20_inconnues()

Sub une_équation_20_inconnues()
  result = 13378
  C01 = 80: I01 = 0: C02 = 148: I02 = 0
  C03 = 196: I03 = 0: C04 = 260: I04 = 0
  C05 = 464: I05 = 0: C06 = 518: I06 = 0
  C07 = 650: I07 = 0: C08 = 688: I08 = 0
  C09 = 908: I09 = 0: C10 = 1326: I10 = 0
  C11 = 1452: I11 = 0: C12 = 1612: I12 = 0
  C13 = 1694: I13 = 0: C14 = 2040: I14 = 0
  C15 = 2498: I15 = 0: C16 = 2724: I16 = 0
  C17 = 5164: I17 = 0: C18 = 5180: I18 = 0
  C19 = 6042: I19 = 0: C20 = 8052: I20 = 0

  Som00 = 0: compteur = 0: [V3] = Now

  For I01 = 0 To Int((result - Som00) / C01)
      som01 = Som00 + I01 * C01: Cells(2, 1).Value = I01
      For I02 = 0 To Int((result - Som1) / C02)
          Som02 = Som1 + I02 * C02: Cells(2, 2).Value = I02
              For I03 = 0 To Int((result - Som02) / C03)
              Som03 = Som02 + I03 * C03: Cells(2, 3).Value = I03
                 For I04 = 0 To Int((result - Som03) / C04)
                 Som04 = Som03 + I04 * C04: Cells(2, 4).Value = I04
                    For I05 = 0 To Int((result - Som04) / C05)
                    Som05 = Som04 + I05 * C05: Cells(2, 5).Value = I05
                        For I06 = 0 To Int((result - Som05) / C06)
                        Som06 = Som05 + I06 * C06: Cells(2, 6).Value = I06
                            For I07 = 0 To Int((result - Som06) / C07)
                            Som07 = Som06 + I07 * C07: Cells(2, 7).Value = I07
                                For I08 = 0 To Int((result - Som07) / C08)
                                Som08 = Som07 + I08 * C08: Cells(2, 8).Value = I08
                                    For I09 = 0 To Int((result - Som08) / C09)
                                    Som09 = Som08 + I09 * C09: Cells(2, 9).Value = I09
                                        For I10 = 0 To Int((result - Som09) / C10)
                                        Som10 = Som09 + I10 * C10: Cells(2, 10).Value = I10
                                            For I11 = 0 To Int((result - Som10) / C11)
                                            Som11 = Som10 + I11 * C11: Cells(2, 11).Value = I11
                                                For I12 = 0 To Int((result - Som11) / C12)
                                                Som12 = Som11 + I12 * C12: Cells(2, 12).Value = I12
                                                    For I13 = 0 To Int((result - Som12) / C13)
                                                    Som13 = Som12 + I13 * C13: Cells(2, 13).Value = I13
                                                        For I14 = 0 To Int((result - Som13) / C14)
                                                        Som14 = Som13 + I14 * C14: Cells(2, 14).Value = I14
                                                            For I15 = 0 To Int((result - Som14) / C15)
                                                            Som15 = Som14 + I15 * C15: Cells(2, 15).Value = I15
                                                                For I16 = 0 To Int((result - Som15) / C16)
                                                                Som16 = Som15 + I16 * C16: Cells(2, 16).Value = I16
                                                                    For I17 = 0 To Int((result - Som16) / C17)
                                                                    Som17 = Som16 + I17 * C17: Cells(2, 17).Value = I17
                                                                        For I18 = 0 To Int((result - Som17) / C18)
                                                                        Som18 = Som17 + I18 * C18: Cells(2, 18).Value = I18
                                                                            For I19 = 0 To Int((result - Som18) / C19)
                                                                            Som19 = Som18 + I19 * C19: Cells(2, 19).Value = I19
                                                                                For I20 = 0 To Int((result - Som19) / C20)
                                                                                som20 = Som19 + I20 * C20: Cells(2, 20).Value = I20
                                                                                'If Abs(Som20 - result) / result < 0.0001 Then Call affiche(I01, I02, I03, I04, I05, I06, I07, I08, I09, I10, I11, I12, I13, I14, I15, I16, I17, I18, I19, I20,compteur)
                                                                                If som20 = result Then Call affiche(I01, I02, I03, I04, I05, I06, I07, I08, I09, I10, I11, I12, I13, I14, I15, I16, I17, I18, I19, I20, compteur)
                                                                                compteur = compteur + 1
                                                                                Next I20
                                                                            Next I19
                                                                        Next I18
                                                                    Next I17
                                                                Next I16
                                                            Next I15
                                                        Next I14
                                                    Next I13
                                                Next I12
                                            Next I11
                                        Next I10
                                    Next I09
                                Next I08
                            Next I07
                        Next I06
                    Next I05
                Next I04
            Next I03
    Next I02
Next I01
End Sub

Sub affiche(I01, I02, I03, I04, I05, I06, I07, I08, I09, I10, I11, I12, I13, I14, I15, I16, I17, I18, I19, I20, compteur)
  Set mapos = Range("A65000").End(xlUp).Offset(1, 0)
  mapos.Offset(0, 0).Value = I01: mapos.Offset(0, 1).Value = I02
  mapos.Offset(0, 2).Value = I03: mapos.Offset(0, 3).Value = I04
  mapos.Offset(0, 4).Value = I05: mapos.Offset(0, 5).Value = I06
  mapos.Offset(0, 6).Value = I07: mapos.Offset(0, 7).Value = I08
  mapos.Offset(0, 8).Value = I09: mapos.Offset(0, 9).Value = I10
  mapos.Offset(0, 10).Value = I11: mapos.Offset(0, 11).Value = I12
  mapos.Offset(0, 12).Value = I13: mapos.Offset(0, 13).Value = I14
  mapos.Offset(0, 14).Value = I15: mapos.Offset(0, 15).Value = I16
  mapos.Offset(0, 16).Value = I17: mapos.Offset(0, 17).Value = I18
  mapos.Offset(0, 18).Value = I19: mapos.Offset(0, 19).Value = I20
  mapos.Offset(0, 20).Value = compteur: mapos.Offset(0, 21).Value = Now
End Sub

cernunnos · Answer

Hé bien ça c'est sympa !
Je te remercie !
Je peux marquer le problème comme résolu !

JvDo · Answer

Bonsoir,

un complément pour ceux que ça intéresse : en C++ le temps de traitement devient raisonnable, environ 8mn pour toutes les solutions entières 

niveau 6 : Duree traitement =  0 m  0 s. Nb Solutions =       398  nb combinaisons =       365.389
niveau 5 : Duree traitement =  0 m  0 s. Nb Solutions =     1.060  nb combinaisons =       854.376
niveau 4 : Duree traitement =  0 m  0 s. Nb Solutions =     3.264  nb combinaisons =     2.392.336
niveau 3 : Duree traitement =  0 m  1 s. Nb Solutions =     9.673  nb combinaisons =     6.922.952
niveau 2 : Duree traitement =  0 m  3 s. Nb Solutions =    34.978  nb combinaisons =    23.202.170
niveau 1 : Duree traitement =  0 m 31 s. Nb Solutions =   331.138  nb combinaisons =   207.110.348
niveau 0 : Duree traitement =  8 m 25 s. Nb Solutions = 5.340.371  nb combinaisons = 3.158.498.324

Cordialement

équation à 20 inconnue

24 réponses

Votre réponse

Discussions similaires

Newsletters