拉格朗日插值

functiony=lagrange(x0,y0,x);n=length(x0);m=length(x);fori=1:mz=x(i);s=0.0;fork=1:np=1.0;forj=1:nif=k

已经编译运行确认：#include#include#includetypedefstructdata{floatx;floaty;}Data;//变量x和函数值y的结构Datad[20];//最多二十

线性插值也叫两点插值,已知函数y=f(x)在给定互异点x0,x1上的值为y0=f(x0),y1=f(x1)线性插值就是构造一个一次多项式　　P1(x)=ax+b　　使它满足条件　　P1(x0)=y0P

对第一个问题进行解答反证法n+1个点（设为(X1,Y1)(X2,Y2)……(Xn+1,Yn+1)）确定一个最高次为n的多项式假设可以确定两个多项式为P(X),Q(X)且P(X)不等于Q(X)令F(X)

应该是给你拉格朗日函数的输入参数不对.再问：你是说x0?应该不是吧。。再答：clearall;clcx0=1:5;y0=sin(x0);x=1:0.2:2;y0=lagrange(x0,y0,x)命令

拉氏定理是导数的,牛顿是积分的

x0=-2.00;x1=0.00;x2=1.00;x3=2.00;y0=17.00;y1=1.00;y2=2.00;y3=17.00;x=0.6y=(x-x1).*(x-x2).*(x-x3)/((x

functionyy=Nlagrange(x,y,xx)yy=0;j=1;n=length(x);while(j

% 为了说明language高次差值的龙格现象n=11;      % n次language差值 lb=-5;

用拉格朗日插值做y=1./(1+x.^2);本来就有很大的偏差,你的插值函数没写错这是我做过的一个y=1/(1+25x^2)的图像,这里取得是11点,你的应该是5点的吧.你的做完是这种么.这样没错,插

这个是MATLAB版本问题,7.0没有收录simplify和double函数,而7.9有,所以出现问题,谢谢!

x0=[30;45;60],y0=[0.5;0.7071;0.866];>>x=40;>>y=Lagrange(x0,y0,x)y=0.6434>>(y-sin(40*pi/180))/sin(40*

这个函数大体是正确的,但有小错误两个1  L=L.*abs((z-x0(j)))/(x0(k)-x0(j));这行不需要abs.2  x是个数值的话 

m=length(x);n=length(y);ifm~=n,error('向量x与y的长度必须一致');ends=0;fori=1:nz=ones(1,length(xi));forj=1:nifj

#include#include#includetypedefstructdata{floatx;floaty;}Data;//变量x和函数值y的结构Datad[20];//最多二十组数据floatf

拉格朗日插值多项式,用于离散数据的拟合C/C++code#include#include#includefloatlagrange(float*x,float*y,floatxx,intn)/*拉格朗

#include#include#includetypedefstructdata{floatx;floaty;}Data;//变量x和函数值y的结构Datad[20];//最多二十组数据floatf

已发送,望采纳!

你能把拉格朗日函数的内容说一下吗?再问：我已经解决了，谢谢。

functionf=Language(x,y,x0)%求已知数据点的拉格朗日差值多项式%已知数据点的x坐标向量:x%已知数据点的y坐标向量:y%插值点的x坐标:x0%求得的拉格朗日插值多项式:f%x0