牛顿插值法，可用于二维曲线拟合求解方程
需积分: 0 84 浏览量
2023-03-02
14:40:01
上传
评论
收藏 53KB DOCX 举报
牛顿插值多项式的一般形式是:
其中 n 是多项式的阶数 是 k 次均差, 被定义为
_k_次均差还可以表示为:
function [A,y]= newton(X,Y,x)
% Newton 插值函数
% X 为已知数据点的 x 坐标
% Y 为已知数据点的 y 坐标
% x 为插值点的 x 坐标
% 函数返回 A 差商表
% y 为各插值点函数值
n=length(X); m=length(x);
for t=1:m
z=x(t); A=zeros(n,n);A(:,1)=Y';
s=0.0; y=0.0; c1=1.0;
for j=2:n
for i=j:n
A(i,j)=(A(i,j-1)- A(i-1,j-1))/(X(i)-X(i-j+1));
end
end
C=A(n,n);
for k=1:n
p=1.0;
for j=1:k-1
p=p*(z-X(j));
end
s=s+A(k,k)*p;
end
ss(t)=s;
end
y=ss;
A=[X',A];
end
资源评论
zhuixundianfeng
- 粉丝: 0
- 资源: 1
