matlab-拉格朗日插值法(共4页).docx

matlab 拉格朗日插值法 function f = Language(x,y,x0)%求已知数据点的拉格朗日插值多项式%已知数据点的x坐标向量: x%已知数据点的y坐标向量: y%插值点的x坐标: x0%求得的拉格朗日插值多项式或在x0处的插值: fsyms t;if(length(x) = length(y)n = length(x);else disp(x和y的维数不相等！); return; %检错endf=0.0;for(i=1:n)l=y(i);for(j=1:i-1) l=l*(t-x(j)/(x(i)-x(j);end;for(j=i+1:n) l=l*(t-x(j)/(x(i)-x(j);end;f=f+1;simplify(f);if(i=n) if(nargin = 3) f = subs (f,t,x0); %计算插值点的函数值 else f=collect(f);