lagran2cichazhi.m

一些matlab程序

function[c,l]=lagran(x,y)
clear all
x=[10 11 12 ]';
y=[2.3026 2.3979 2.4849];
w=length(x);
n=w-1;
l=zeros(w,w);
for k=1:w
    v=1;
    for j=1:w
        if k~=j
            v=conv(v,poly(x(j)))/(x(k)-x(j));
        end
    end
    l(k,:)=v;
end
c= y*l ;
x=11.75;
xk=[x^2 x 1]';
f=c*xk;
ff=2.4639;
ek=ff-f