matlab作图:找出函数f(x)=1/(x^2+1)的20阶插值多项式p(x),做出f(x...

n = 21;%%%function y=ex217去掉，运行即可

x = linspace(-5,5,n)';转置符号，英文输入

h = (5-(-5))/(n-1);

y = 1./(1+x.^2);

% form the differences table

for j = 2:n,

y(1:n+1-j,j) = diff(y(1:n+2-j,j-1))./(x(j:n)-x(1:n+1-j));

end

% newton coeff（以上为ex217.m里的内容）

y = y(1,:);

pz = [ ];

v = linspace(-5,5,80);

for t = v,

z = y(n);

for j = n-1:-1:1,

z = z * ( t - x(j) ) + y(j);

end

pz = [pz z];

end

plot(v,pz,’r+-‘,v,1./(1+v.^2),‘g--’)%%引号英文输入

function y = ex217
n = 21;
x = linspace(-5,5,n);
h = (5-(-5))/(n-1);
y = (1./(1+x.^2))';
% form the differences table
for j = 2:n,
y(1:n+1-j,j) = diff(y(1:n+2-j,j-1))./(x(j:n)-x(1:n+1-j))';
end
% newton coeff（以上为ex217.m里的内容）
y = y(1,:);
pz = [ ];
v = linspace(-5,5,80);
for t = v,
z = y(n);
for j = n-1:-1:1,
z = z * ( t - x(j) ) + y(j);
end
pz = [pz z];
end
plot(v,pz,'r+-',v,1./(1+v.^2),'g--')

本文如未解决您的问题请添加抖音号：51dongshi（抖音搜索懂视），直接咨询即可。