%Programa 5.1 (Recta de regresion - pagina 280).
% Construccion de la recta de regresion y = Ax + B que mejor ajusta en el sentido de los
% minimos cuadrados a los N datos (X1,Y1),...,(Xn,Yn).
%
% Se resuelve el sig sistema para obtener los coef. A y B de la recta que mejor ajusta.
% suma(Xk^2) · A + suma(Xk) · B = suma(Xk·Yk)
% suma(Xk) · A + N · B = suma(Yk)
%
%Sintaxis
% lsline(X,Y)
%
%Entrada
% - X es el vector fila de abscisas 1 x n
% - Y es el vector fila de ordenadas 1 x n
%Resultados
% - A es el coeficiente de x en Ax + B
% - B es el termino independiente en Ax + B
xmean=mean(X);
ymean=mean(Y);
sumx2=(X-xmean)*(X-xmean)';
sumxy=(Y-ymean)*(X-xmean)';
A=sumxy/sumx2;
B=ymean-A*xmean;
viernes, 4 de diciembre de 2009
Suscribirse a:
Enviar comentarios (Atom)
No hay comentarios:
Publicar un comentario