lagrange

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

clc
clear
%interpolação po lagrange!!
xs=input('digite os valores de x=(x0 x1 x2 ...xn) : ');
ys=input('digite os valores de y=(y0 y1 y2 ...yn) : ');
x=input('digite o valor a ser interpolado :');
n=length(xs); % o tamanho do vetor
if length(ys)~=n, error('x e y devem ter a mesmo tamanho');
end
yi=0;
L=zeros(n,n);
 for i=1:n
 produto=ys(i);
 v=1;
 for j=1:n
 if j~=i,
 v=conv(v,poly(xs(j)));
 v=v/(xs(i)-(xs(j)));
 end
 end
 for j=1:n
 if i~=j
 produto=produto*(x-xs(j))/(xs(i)-xs(j));
 end
 end
 yi=yi+ produto;
 L(i,:)=v; 
 end
 fprintf('\n')
fprintf(' O valor da interpolação é: %f\n',yi) 
p=ys*L;
syms x
 y=0;
 m = length(p);
 p=p(m:-1:1);
 i = 0;
 j = 0;
 s = 0;
 for i=m:-1:1
 y=(p(1,i))*x^(i-1);
 s=s + y;
 end
 fprintf('\n')
fprintf('Pn(x)= \n'),pretty(s)