Apprendre Maple Index du Forum Apprendre Maple
Site dédié au logiciel de calcul formel Maple
 
  Page d'accueilPage d'accueil   FAQFAQ    RechercherRechercher    Liste des MembresListe des Membres    Groupes d'utilisateursGroupes d'utilisateurs 
S'enregistrerS'enregistrer    ProfilProfil   Se connecter pour vérifier ses messages privésSe connecter pour vérifier ses messages privés   ConnexionConnexion 

Approximation des moindres carrés.

 
Poster un nouveau sujet   Répondre au sujet    Apprendre Maple Index du Forum -> Maple en analyse
Voir le sujet précédent :: Voir le sujet suivant  
Auteur Message
minitim



Inscrit le: 11 Avr 2012
Messages: 4

MessagePosté le: 28 Mai 2012 10:10    Sujet du message: Approximation des moindres carrés. Répondre en citant

Bonjour a tous!

Je souhaite approcher un ensemble de points par une hyperbole par la méthode des moindres carrés mais je ne sais pas comment m'y prendre..
Existe-t-il une commande pour se faire?

Après une petite recherche je sais que mon hyperbole est une fonction du type :

y=(1/(-a*x+b))-(1/b) avec a et b constantes positives

Avez tous mes remerciements pour votre aide!

Bonne journée Smile
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
zozo



Inscrit le: 03 Jan 2013
Messages: 125

MessagePosté le: 28 Mai 2012 17:34    Sujet du message: Répondre en citant

Hi,
Tu peux utiliser les outils du package CurveFittting, mais ce qui gène est que la fonction n'est pas linéaire.

> with(CurveFitting);

A plus.
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
minitim



Inscrit le: 11 Avr 2012
Messages: 4

MessagePosté le: 28 Mai 2012 21:57    Sujet du message: Répondre en citant

Bah du coup j'ai essayé mais rien n'est concluant (du tout)
Penses-tu qu'il existerait une autre manière de faire?
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
ALS



Inscrit le: 11 Sep 2006
Messages: 647

MessagePosté le: 29 Mai 2012 8:25    Sujet du message: Re Répondre en citant

Bonjour,
L'équation n'est pas linéaire en les coefficients, mais on peut utiliser une régression linéaire en écrivant l'équation y=(1/(-a*x+b))-(1/b) sous la forme:
x/y= (b-ax)b/a = -bx +b^2/a
On fait donc une régression linéaire de W=x/y en X grace à la fonction LinearFit du package Statistics.
En en déduit les coefficients -b et b^2/a, donc a et b.


Code:

> restart:
> with(Statistics):
> X := Vector([1, 2, 3, 4, 5, 6], datatype=float):
> Y := Vector([0.029, 0.073, 0.15, 0.31, 0.85, -4.5], datatype=float):
> W:=Vector([seq(X[i]/Y[i],i=1..6)],datatype=float);

                            [34.4827586206897 ]
                            [                 ]
                            [27.3972602739726 ]
                            [                 ]
                            [       20.       ]
                       W := [                 ]
                            [12.9032258064516 ]
                            [                 ]
                            [5.88235294117647 ]
                            [                 ]
                            [-1.33333333333333]

> P:=LinearFit([1,x], X, W, x);

              P := 41.6275729810314 - 7.16348445605863 x

> b:=-coeff(P,x,1); solve(b^2/a=coeff(P,x,0),a): a:=%;

                        b := 7.16348445605863


                           a := 1.232728835

> plot(1/(b-a*x)-1/b,x=1..6,y=-10..10);



Les valeurs exactes attendues étaient a=1.23 et b:=7.15.
@+
Revenir en haut de page
Voir le profil de l'utilisateur Envoyer un message privé
Montrer les messages depuis:   
Poster un nouveau sujet   Répondre au sujet    Apprendre Maple Index du Forum -> Maple en analyse Toutes les heures sont au format GMT + 2 Heures
Page 1 sur 1

 
Sauter vers:  
Vous ne pouvez pas poster de nouveaux sujets dans ce forum
Vous ne pouvez pas répondre aux sujets dans ce forum
Vous ne pouvez pas éditer vos messages dans ce forum
Vous ne pouvez pas supprimer vos messages dans ce forum
Vous ne pouvez pas voter dans les sondages de ce forum


phpBB

Développé par phpBB © 2001, 2006 phpBB Group
Traduction par : phpBB-fr.com


Apprendre Maple - ©  - Alain Le Stang - Navigation optimisée pour une résolution 1024 x 768.