数值计算讨论组

欢迎来到数值计算讨论组。 本论坛旨在增强师生交流的时效性,主要讨论数值计算相关问题,包括初等数值分析、矩阵计算、非线性方程(组)数值解法、微分方程数值解法等方面的内容。


    [旧]Newton向前插值算例

    分享

    Admin
    Admin

    帖子数 : 46
    注册日期 : 10-05-13

    [旧]Newton向前插值算例

    帖子 由 Admin 于 周四 五月 13, 2010 10:23 pm

    % NewtonForward.m
    % P44 例4
    clc;
    format long
    x = 2.4:0.1:2.9;
    y = besselj(0,x);
    n = length(x);

    % 计算差分。通过下面的计算,将x0处的第i阶差分存入y(i+1)
    for i = 1:n-1
    for j = n:-1:i+1
    y(j) = y(j)-y(j-1);
    end
    end

    % 通过Newton向前插值公式计算2.45处的函数值
    f = y(1);
    t = 0.5;
    co = 1;
    for i = 1:n-1 % i 表示差分阶数
    % co 表示组合数(t,i)
    co = co * (t-i+1)/i;
    % 输出 利用从点2.4及后面i个节点得到插值多项式 在2.45处的函数值
    f = f + co*y(i+1)
    end
    % 2.45处精确函数值
    f_exact = besselj(0,2.45)

      目前的日期/时间是周五 十二月 14, 2018 2:32 am