新冠插值法matlab_数值计算方法实验之Newton 多项式插值(MATLAB代码)

一、实验目的

在己知f(x),x∈[a,b]的表达式，但函数值不便计算或不知f(x),x∈[a,b]而又需要给出其在[a,b]上的值时,按插值原则f(xi)=yi (i=0,1,……, n)求出简单函数P(x)(常是多项式)，使其在插值基点xi处成立(xi)= yi(i=0,1,……,n),而在[a,b]上的其它点处成立f(x)≈P(x).

二、实验原理

三、实验内容

求f(x)=x4在[0,2]上按5个等距节点确定的Lagrange插值多项式

四、实验程序

(1).m文件

%输入的量:X是n+1个节点(x_i,y_i)(i = 1,2, ... , n+1)横坐标，Y是纵坐标，

%x是以向量形式输入的m个插值点，M在[a,b]上满足｜f~(n+1)(x)｜≤M

%注：f~(n+1)(x)表示f(x)的n+1阶导数

%输出的量：向量y是向量x处的插值，误差限R，n次牛顿插值多项式L及其系数向量C，

%差商的矩阵A

function[y,R,A,C,L] = newton(X,Y,x,M)

n = length(X);

m = length(x);

for t = 1 : m

z = x(t);

A = zeros(n,n);

A(:,1) = Y';

s = 0.0; p = 1.0; q1 = 1.0; c1 = 1.0;<