牛顿插值Matlab编程计算，并验证输出

最新推荐文章于 2022-10-11 11:17:10 发布

逐梦之程

最新推荐文章于 2022-10-11 11:17:10 发布

阅读量1.4k

点赞数 2

分类专栏：算法文章标签：算法机器学习线性代数

本文链接：https://blog.csdn.net/qq_36625174/article/details/122442969

版权

算法专栏收录该内容

14 篇文章 7 订阅

订阅专栏

所用到的公式

1、差商

2、差商表格如下

$xk$	$f(xk)$	一阶均差	二阶均差	三阶均差	...
$x0$	$f(x0)$
$x1$	$f(x1)$	$f(x0,x1)$
$x2$	$f(x2)$	$f(x1,x2)$	$f(x0,x1,x2)$
$x3$	$f(x3)$	$f(x2,x3)$	$f(x1,x2,x3)$	$f(x0,x1,x2,x3)$
$x4$	$f(x4)$	$f(x3,x4)$	$f(x2,x3,x4)$	$f(x1,x2,x3,x4)$	...
...	...	...	...	...	...

3、计算数值

4、代码

%给x轴赋值
for i=1:1:10
    x{i}=i;
end
%读出X轴的已知点的数量
num = numel(x);

%给y轴赋值，为了便于验证使用方程
for i=1:1:num
    y{i} = x{i}^3-6*x{i}^2-x{i}*5+16;
end

%绘制图形的点云
xInit=x{1};%初始值
xStep=0.1;%绘制图形步长
xEnd=x{num};%最后一个值

f = cell(num,num);%预分配存储空间

%原点云
for i=1:1:num
    f{i,1}=y{i};
end

%差商计算，计算储存在主对角线上
for j=2:1:num
    for i=j:1:num
        f{i,j} = (f{i,j-1} - f{i-1,j-1})/(x{i} - x{i-j+1});
    end
end

%计算其值
m=1;
for k=xInit:xStep:xEnd
    xy{m,1}=k;
    %插值计算结果
    result =f{1,1};
    for j=2:1:10
        product =1;
        for i=1:1:j-1
           product = product *( xy{m,1}-x{i});
        end
        result  = result  + f{j,j}*product;
    end
    %计算完成
    xi(m)=k;
    yi(m)=result;
    xy{m,2}=result;
    m=m+1;
end

%图像输出
plot(xi,yi);

5、运行结果