Je me suis fait une petite macro pour présenter la méthode des moindres carrés dans un tableau :
- Code: Tout sélectionner
TeXgraph#
{ stat2var.mac version 0.2 du 23/09/2011 }
Include
"variations.mod";
Mac
serie2varTab =
//serie2varTab(<Liste de complexe>,<nb de chiffres apres la virgule si necessaire>)
[
L:=%1, nbdeci:=%2,
NbLigne:=Nops(L),
usecomma:=1,
//Construction du tableau, des lignes et des colonnes
dimTab(NbLigne+2,6),
defLargeCol(0,3.75), defHautLig(1,1.25), Cadre(1,1,NbLigne+1,6),
for k from 1 to NbLigne do traitH(k) od,
for k from 1 to 6 do traitV(k,1,NbLigne+1) od,
traitH(NbLigne+2,1,2), traitH(NbLigne+2,5,6),
for k in [0,1,2,4,5,6,7] do traitV(k,NbLigne+1,NbLigne+2) od,
//Affichage des entetes des listes
InitLig1("x_i","y_i",-1,"x_i-\overline{x}",-1,"y_i-\overline{y}",-1,
"\left(x_i-\overline{x}\right)\left(y_i-\overline{y}\right)",-1,
"\left(x_i-\overline{x}\right)^2"),
//Creation des listes
//List1 = xi
List1:=[for k in L do Re(k) od],
//List2 = yi
List2:=[for k in L do Im(k) od],
xbar:=moy(List1), ybar:=moy(List2),
Label(pos(NbLigne+2,1),["$\overline{x}=",StrNum(xbar),"$"],
pos(NbLigne+2,2),["$\overline{y}=",StrNum(ybar),"$"]),
//List3 = xi-xbar
List3:=[for k from 1 to NbLigne do List1[k]-xbar od],
//List4 = yi-ybar
List4:=[for k from 1 to NbLigne do List2[k]-ybar od],
//List5 = (xi-xbar)(yi-ybar)
List5:=[for k from 1 to NbLigne do List3[k]*List4[k] od],
//List6 = (xi-xbar)^2
List6:=[for k from 1 to NbLigne do (List3[k])^2 od],
//Calcul et affichage des valeurs des listes
for k from 1 to NbLigne do
Label(pos(k+1,1),["$",StrNum(List1[k]),"$"],
pos(k+1,2),["$",StrNum(List2[k]),"$"],
pos(k+1,3),["$",StrNum(List1[k]-xbar),"$"],
pos(k+1,4),["$",StrNum(List4[k]),"$"],
pos(k+1,5),["$",StrNum(List5[k]),"$"],
pos(k+1,6),["$",StrNum(List6[k]),"$"]),
od,
Label(pos(NbLigne+2,4),["Somme :"],
pos(NbLigne+2,5),["$",StrNum(sum(List5)),"$"],
pos(NbLigne+2,6),["$",StrNum(sum(List6)),"$"]),
];
le code
[
L:=[19*i, 5+21*i,10+25*i,15+30*i, 20+37*i, 25+42*i, 30+50*i],
serie2varTab(L,2),
]
permet d'obtenir :
Bonne journée
mise à jour du 23/09/11 - 20:30