[讨论]

每周一题（四）：编写拉格朗日插值法的MATLAB程序

2011-5-8 22:26:44 5082

0 要求：编写拉格朗日插值法的matlab程序。 0
2011-5-8 22:26:44　　 1 评论淘帖0 举报相关推荐 • 如何用MATLAB拉格朗日插值法作图 15358 • labview实现拉格l朗日插值法 2094 • matlab应用源代码拉格朗日插值法函数版 3642 • 每周一题（一） 4086 • 每周一题(三) 3098 • 每周一题（二） 3123 • 每周一题（五）滤波器设计 2632 • 拉格朗日算法求救！！！ 2969 • 基于MPI并行环境下拉格朗日插值的求解 11 • 基于拉格朗日插值多项式的光伏电池I-V特性建模方法 22 1 条评论电子工程师 2011-5-8 23:49 标题上要写上题目呵。赞回复举报 3 个讨论

0 回复 A670521546 的帖子 %MatLab 自定义函数 lag_func %输入参数：xi,yi为插值节点值、x为待插值点的x轴值；输出参数：y为插值获得的x对应的值。 function y=lag_func(xi,yi,x) % 计算l(i,j) [r,n]=size(xi); [r,m]=size(x); l=zeros(n,m); for i=1:n for j=1:m a=1; for k=1:n if (k~=i) a=a(x(j)-xi(k))/(xi(i)-xi(k)); end end l(i,j)=a; end end y=yil;

2011-5-12 22:09:54 评论举报 1^# vvg

0 这不是书上的典型吗？直接复制

2011-6-10 10:51:26 评论举报 2^# 王涛

拉格朗日插值的MATLAB函数文件agui_lagrange.m如下。
function f=agui_lagrange(x0,y0,x)
%x0为节点向量，y0为节点上的函数值，x为插值点，f返回插值
n=length(x0);m=length(x);
format long
s==0.0;
for k=1:n
   p=1.0;
   for j=1:n
      if j~=k
         p=p*(x-x0(j))/(x0(k)-x0(j));
      end
   end
   s=p*y0(k)+s;
end
f=s;
end

MATLAB命令窗口中实现过程：
>>x0=[已知的x取值]
>>y0=[已知的y取值]
>>f=agui_lagrange(x0,y0,要求的x值)