[Matlab] Méthode des moindres carrés Résolu/Fermé

Question

Bonjour,
svt donner moi le programme de la méthode de moinde careé en matlab

Fee Fay · Accepted Answer

Bonjour mon ami !

Il s'agit encore d'une demande d'une clarté exemplaire. La méthodes des moindres carrés est utilisée pour moult choses et écrite de cette manière, ta demande n'a pas grande signification. Pour rester dans le ton du sujet, voici quelques pistes.

Utilise la fonction polyfit pour fitter des données au sens des moindres carrés à l'aide d'une fonction polynômiale. Tape help polyfit ou doc polyfit pou avoir plus d'informations sur cette fonction.

Pour résoudre un système d'équations linéaires A*X=B au sens des moindres carrés, tu as juste à utiliser la fonction mldivide qui s'utilise plus simplement avec sa notation en terme d'opérateur arithmétique : un backslash \.
Il y a donc juste à écrire B=A\X; pour résoudre ton système d'équations linéaires au sens des moindres carrés. help mldivide et doc mldivide te fourniront des informations supplémentaires si tu en as besoin.

Je te souhaite une bien bonne journée !

che0po · Answer

function polmc=mcarre(mp,degre)
  taillemp=size(mp) // matrice (5,2)
  mc=ones(degre+1,degre+1) //créer une matrice remplie de 1 de taille degre+1, degre+1
  b=ones(degre+1,1) //matrice vecteur remplie de 1
  mc(1,1)=taillemp(1,1)  
  //calcul premier ligne de mc
  for col=2: degre+1 do
    s=0
    for i=1: taillemp(1,1) do
      s=s+mp(i,1)^(col-1)
    end
    mc(1,col)=s
  end
  //calcul derniere colonne de mc
  for ligne=2: degre+1 do
    s=0
    for i=1: taillemp(1,1) do
      s=s+mp(i,1)^(degre+ligne-1)
    end
    mc(ligne,degre+1)=s
  end
  //propagation en bas a gauche
  for ligne=2: degre+1 do
    for col=1: degre do
      mc(ligne,col)=mc(ligne-1,col+1)
    end
  end  
  //calcul de B
  for ligne=1: degre+1 do
    s=0
    for i=1: taillemp(1,1) do
      s=s+mp(i,2)*mp(i,1)^(ligne-1)
    end
    b(ligne,1)=s
  end
  //resolution
  X=linsolve(mc,-b)
  //calcul d'erreur
  x=poly(0,'x')
  polmc=0
  for k=1: degre+1 do
    polmc=polmc+X(k,1)*x^(k-1)
  end
  erreur=0
  for i=1: taillemp(1,1) do
    erreur=erreur+(mp(i,2)-horner(polmc,mp(i,1)))^2
  end 
  printf("
 erreur = %f",erreur)
    
endfunction

[Matlab] Méthode des moindres carrés

2 réponses

Discussions similaires