matlab编写拉格朗日插值代码函数

matlab编写拉格朗日插值代码函数

要求：根据拉格朗日多项式插值法原理，设计算法流程并且编写拉格朗日插值代码函数。

代码如下：

function[y]=lagrange(x0,y0,x)   %建立一个函数名为lagrange的函数，输入x0,y0为插值点的坐标，均为数组，x为要求的点的横坐标，此处为一个数组，长度为n，表示一次可以求n个点。 N=length(x0); n=length(x); y=zeros(1,n); for k=1:n %对n个待求点循环     l=ones(N); for i=1:N %对每个插值点循环     for j=1:N %计算基函数         if(j~=i)             l(i)=l(i)*(x(k)-x0(j))/(x0(i)-x0(j));         end     end     y(k)=y(k)+y0(i)*l(i);  %计算拉格朗日插值函数 end end scatter(x0,y0,'r','*'); %画出插值点的图 hold on  plot(x,y,'b');  %画出插值函数的图 end

在命令行窗口输入如下命令，x0,y0表示插值点的坐标分别为（-1，1），（0，1），（1，1）

现在要分别求横坐标为-3，-2，0，2，3的点的值。

 运行后可以看到如下结果

说明对应的纵坐标分别为 9，4，0，4，9