Lagrange插值函数
2007-12-07 20:52
260 查看
function yi=lagint(x,y,xi)
% x,y为已知节点向量值;yi为xi处所求插值点
d_xi=xi-x;
n=length(x);
L=zeros(size(y));
L(1)=prod(d_xi(2:n))/prod(x(1)-x(2:n));
%或者 L(1)=prod(d_xi(2:n)./(x(1)-x(2:n)));
L(n)=prod(d_xi(1:n-1))/prod(x(n)-x(1:n-1));
%或者 L(n)=prod(d_xi(1:n-1)./(x(n)-x(1:n-1)));
for j=2:n-1
fenzi=prod(d_xi(1:j-1))*prod(d_xi(j+1:n));
fenmu=prod(x(j)-x(1:j-1))*prod(x(j)-x(j+1:n));
L(j)=fenzi/fenmu;
end
yi=sum(y.*L);
%fprintf('x=%f,y=%f',xi,yi);
% x,y为已知节点向量值;yi为xi处所求插值点
d_xi=xi-x;
n=length(x);
L=zeros(size(y));
L(1)=prod(d_xi(2:n))/prod(x(1)-x(2:n));
%或者 L(1)=prod(d_xi(2:n)./(x(1)-x(2:n)));
L(n)=prod(d_xi(1:n-1))/prod(x(n)-x(1:n-1));
%或者 L(n)=prod(d_xi(1:n-1)./(x(n)-x(1:n-1)));
for j=2:n-1
fenzi=prod(d_xi(1:j-1))*prod(d_xi(j+1:n));
fenmu=prod(x(j)-x(1:j-1))*prod(x(j)-x(j+1:n));
L(j)=fenzi/fenmu;
end
yi=sum(y.*L);
%fprintf('x=%f,y=%f',xi,yi);
相关文章推荐
- newton插值 lagrange 插值
- 插值函数
- poj 2042 Lagrange's Four-Square Theorem
- [物理学与PDEs]第2章第5节 一维流体力学方程组的 Lagrange 形式 5.4 一维粘性热传导流体力学方程组的 Lagrange 形式
- 样条之抛物线(一元三点)插值函数
- NYOJ 178 找规律(Lagrange插值公式)
- Android应用:Lagrange
- Lagrange插值法
- 拉格朗日插值法在数据分析中的应用——Python插值scimpy,lagrange
- Lagrange插值问题
- POJ2042 Lagrange's Four-Square Theorem
- matlab 插值函数
- [物理学与PDEs]第2章习题10 一维理想流体力学方程组的 Lagrange 形式
- 样条之连分式插值函数
- NYOJ 题目178 找规律(Lagrange插值)
- 样条之EHMT插值函数
- 小观matlab插值函数
- 优化问题及其Lagrange对偶问题
- [转载]Lagrange插值方法
- Java实现lagrange 插值法