MATLAB实现牛顿插值的源程序

最新推荐文章于 2024-08-20 11:44:39 发布

文宇肃然

最新推荐文章于 2024-08-20 11:44:39 发布

阅读量4w

点赞数 21

分类专栏： MATLAB深入理解高级教程（附源码）

本文链接：https://blog.csdn.net/wenyusuran/article/details/41725983

版权

MATLAB深入理解高级教程（附源码）专栏收录该内容

192 篇文章 1054 订阅 ¥29.90 ¥99.00

订阅专栏

这段代码提供了一个MATLAB函数，用于实现牛顿插值方法。给定一组插值点(x, y)，函数计算在任意点x=xi处的牛顿插值多项式的值。首先初始化差商表，然后通过迭代计算差商，最后根据差商表求得牛顿插值多项式并求解在指定点的值。" 133173876,19974269,CSS连接指南：如何将CSS应用于HTML,"['CSS', 'HTML', '前端开发', '编程']

摘要由CSDN通过智能技术生成

function yi=Newton(x,y,xi)
%Newton插值方法，给定一系列插值的点(x,y)，得到在x=xi处的，牛顿插值多项的值yi
n=length(x);
m=length(y);
if m~=n
    error('x,y的长度不一样,请重新输入!');
    return
end
A=zeros(n);    %定义差商表
A(:,1)=y;      %差商表第一列为y
for j=2:n               %j为列标
    for i=1:(n-j+1)     %i为行标
        A(i,j)=(A(i+1,j-1)-A(i,j-1))/(x(i+j-1)-x(i));   %计算差商表
    end
end
%根据差商表,求对应的牛顿插值多项式在x=xi处的值yi
N(1)=A(1,1);
for j=2:n
    T=1;
    for i=1:j-1
        T=T*(xi-x(i));
    end
    N(j)=A(1,j)*T;
end
yi=sum(N);   %将x=xi带入牛顿插值多项式，得到的yi的值
end