%程序:lagrange.m
%功能:拉格朗日插值算法（MATLAB R2009a）
%說明:x自變量取值向量
% ? ? y對應(yīng)x的函數(shù)取值
% ? ? xh可以是一個(gè)向量，計(jì)算xh處插值結(jié)果
% ? ? yh對應(yīng)xh的結(jié)果
%實(shí)例:給定x=0:2*pi,y=sin(x),xh=0:0.1:2*pi
% ? >>x=0:2*pi
% ? >>lagrange(x,sin(x),0:0.1:2*pi);
function yh=lagrange(x,y,xh)
n=length(x);
m=length(xh);
yh=zeros(1,m);
c1=ones(n-1,1);
c2=ones(1,m);
for i=1:n
? ? xp=x([1:i-1 i+1:n]);
? ? yh=yh+y(i)*prod((c1*xh-xp'*c2)./(x(i)-xp'*c2));
end
plot(xh,yh,'b');?
hold on;
plot(x,y,'r');
title('Lagrange插值');

總結(jié)

以上是生活随笔為你收集整理的MATLAB拉格朗日插值的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。