Programa 5.2 (Polinoio optimo en minimos cuadrados - pagina 296).

%Programa 5.2 (Polinoio optimo en minimos cuadrados - pagina 296).
% Construccion del polinomio de grado M dado por
% Pm(x)=C1 + C2·X + C3·X^2 + ··· + Cn·X^M
% que mejor se ajusta en el sentido de los minimos cuadrados a los N datos (X1,Y1),...,(Xn,Yn).
%
%
%Sintaxis
% lspoly(X,Y,M)
%
%Entrada
% - X es el vector fila de abscisas 1 x n
% - Y es el vector fila de ordenadas 1 x n
% - M es el grado del polinomio optimo
%Resultados
% - C es el vector de coeficientes del polinomio
% en potencias decrecientes de x

n=length(X);
B=zeros(1:M+1);
F=zeros(n,M+1);

%Fill the columns of F with the powers of X

for k=1:M+1
F(:,k)=X'.^(k-1);
end

%Solve the linear system from (25)

A=F'*F;
B=F'*Y';
C=A\B;
C=flipud(C);