再lagrange 插值 matlab的实现

[复制链接]

字体大小: 正常放大

ZCR

4 主题	3 听众	142 积分

升级 21%

TA的每日心情

	开心 2014-11-30 18:43

签到天数: 48 天

[LV.5]常住居民I

群组: 数学建摸协会

群组: 学术交流A

电梯直达

1^#

发表于 2012-10-16 17:52 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

M文件：
function L = LagrangeInterp(x,y,u)
syms t l;
if(length(x) == length(y))
n = length(x);
else
disp('x和y的维数不相等！');
return;
end
L=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;
L=L+l;
end
simplify(L);
if(nargin == 3)
L = subs (L,'t',u);
else
L=collect(L);
L = vpa(L,6);
end
其中l=l*(t-x(j))/(x(i)-x(j));怎么解释……看不太懂！

zan

matlab