《matlab 拉格朗日插值法.docx》由会员分享,可在线阅读,更多相关《matlab 拉格朗日插值法.docx(4页珍藏版)》请在金锄头文库上搜索。
1、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(
2、i)-x(j);end;f=f+1;simplify(f);if(i=n) if(nargin = 3) f = subs (f,t,x0); %计算插值点的函数值 else f=collect(f); f = vpa(f,6); %将插值多项式的系数化成6位精度的小数 end endend在matlab中输入x=18 31 66 68 70 72 70; y=23 33 52 51 43 40 46;f=Language(x,y)plot(x,y)x = 18 31 66 68 70 72 70出现错误? Function collect is not defined for values
3、of class double.Error in = Language at 32 f=collect(f);怎么解决,谢谢了,急!问题补充: 还是不行,? Error using = sym.mapleat offset 12, ) unexpectedError in = sym.collect at 34 r = maple(collect,s,x);Error in = Language at 32 f=collect(f);最佳答案 function f = Language(x,y,x0)%求已知数据点的拉格朗日插值多项式%已知数据点的x坐标向量: x%已知数据点的y坐标向量: y
4、%插值点的x坐标: x0%求得的拉格朗日插值多项式或在x0处的插值: fsyms t;if(length(x) = length(y)n = length(x);else disp(x和y的维数不相等!); return; %检错endf=0.0; %改为 f=sym(0);for(i=1:n)l=y(i); %改为 l=sym(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
5、 (f,t,x0); %计算插值点的函数值 else f=collect(f); f = vpa(f,6); %将插值多项式的系数化成6位精度的小数 end endend -下面的这个应该可以:function f = Language(x,y,x0)%求已知数据点的拉格朗日插值多项式%已知数据点的x坐标向量: x%已知数据点的y坐标向量: y%插值点的x坐标: x0%求得的拉格朗日插值多项式或在x0处的插值: fsyms t l;if(length(x) = length(y) n = length(x);else disp(x和y的维数不相等!); return; %检错endh=sym(0);for (i=1:n) l=sym(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; h=h+l;endsimplify(h);if(nargin = 3) f = subs (h,t,x0); %计算插值点的函数值else f=collect(h); f = vpa(f,6); %将插值多项式的系数化成6位精度的小数end