拉格朗日函数

拉格朗日插值函数

function y=Lagrangian(x1,y1,x)

m=length(x1);

n=length(x);

for i=1:n;

    z=x(i);

    s=0;

    for k=1:m;

        p=1;

        for j=1:m;

            if j~=k

                p=p*(z-x1(j))/(x1(k)-x1(j));

            end

        end

        s=p*y1(k)+s;

    end

    y(i)=s;

end

 

 

A=[1  2.74

   2  2.93

   3  3.40

   4  2.87

   5  2.27

   6  1.45

   7  1.20

   8  1.34

   9  1.55

   10 1.71

   11 2.43

   12 2.51];

x=A(:,1);

y=A(:,2);

plot(x,y,'r*-');

title('09年茄子平均价格随时间变化曲线');

xlabel('月份');

ylabel('元/kg');