clearclf hold on;X=-1021;Y=11+25X^2;CD=newpoly X Y;x=-10011;y=polyvalC x;plotxy X Y;grid on;xp=-1021;z=11+25xp^2;plotxpzrfunction CD=newpoly XYn=lengthX;D=zeros nnD1=Yfor j=2nfor k=jnDkj=Dkj-1-Dk-1j-1

% 清空图像窗口，清空命令窗口，将图像保持在同一画布上 clear, clf, hold on;

% 初始化X和Y的值 X=-1:0.2:1; Y=1./(1+25*X.^2);

% 调用函数newpoly计算多项式C和D [C,D]=newpoly (X, Y);

% 初始化x和y的值 x=-1:0.01:1; y=polyval(C, x);

% 绘制多项式拟合曲线和原始数据点图像 plot(x,y, X, Y,'.');

% 绘制用于计算多项式的节点 grid on;

% 初始化xp和z的值 xp=-1:0.2:1; z=1./(1+25*xp.^2);

% 绘制原始函数图像 plot(xp,z,'r')

% 定义函数newpoly function [C,D]=newpoly (X,Y) % 获取X的长度 n=length(X);

% 初始化D矩阵 D=zeros (n,n)

% 将Y值赋给D的第一列 D(:,1)=Y'

% 通过递推计算出D的所有元素 for j=2:n for k=j:n D(k,j)=(D(k,j-1)-D(k-1,j-1))/(X(k)-X(k-j+1)); end end

% 初始化多项式系数C C=D(n,n);

% 通过递推计算出多项式系数C for k=(n-1):-1:1 C=conv(C,poly (X(k))) m=length(C); C(m)=C(m)+D(k,k); end en