matlab求拉格朗日插值多项式并画图我用这个matlab程序可以求出x0点的函数值,但我想要求出插值多项式怎么求啊,还有把多项式画图象怎么弄?nx=length(x);ny=length(y);if nx~=nywarning('矢量x与y的长度应

来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/17 00:53:04
matlab求拉格朗日插值多项式并画图我用这个matlab程序可以求出x0点的函数值,但我想要求出插值多项式怎么求啊,还有把多项式画图象怎么弄?nx=length(x);ny=length(y);if

matlab求拉格朗日插值多项式并画图我用这个matlab程序可以求出x0点的函数值,但我想要求出插值多项式怎么求啊,还有把多项式画图象怎么弄?nx=length(x);ny=length(y);if nx~=nywarning('矢量x与y的长度应
matlab求拉格朗日插值多项式并画图
我用这个matlab程序可以求出x0点的函数值,但我想要求出插值多项式怎么求啊,还有把多项式画图象怎么弄?
nx=length(x);
ny=length(y);
if nx~=ny
warning('矢量x与y的长度应该相等')
return
end
m=length(x0);
%按照公式,对要求的插值点矢量x0的每个元素进行计算
for i=1:m
t=0.0;
for j=1:nx
u=1.0;
for k=1:nx
if =j
u=u*(x0(i)-x(k))/(x(j)-x(k));
end
end
t=t+u*y(j);
end
s(m)=t;
end
return
我要求的是这样的
x=[0 100 200 300 500];
y=[0.276 1.67 1.97 2.2 2.53];
figure
plot(x,y,'*')
x0=250
hold on
r=Lagrangenew(x,y,x0);
disp('r=');disp(r);
plot(x0,r,'r-*')
plot(x,y)
title('二氧化硫与电压关系分析')
求出插值多项式并画图出来,

matlab求拉格朗日插值多项式并画图我用这个matlab程序可以求出x0点的函数值,但我想要求出插值多项式怎么求啊,还有把多项式画图象怎么弄?nx=length(x);ny=length(y);if nx~=nywarning('矢量x与y的长度应
你需要把上面的部分保存为Lagrangenew.m文件.
文件头应加  function s=Lagrangenew(x,y,x0)
后面是
nx=length(x);
ny=length(y);
.
然后,可以命令窗口执行下面部分的命令
结果如下