matlab高手近来看下画图问题
发布网友
发布时间：2023-07-19 10:05
共1个回答
热心网友
时间：2024-11-26 08:10
%求lagrange插值多项式
x=[0.4 0.55 0.65 0.80 0.95 1.05];
y=[0.41075 0.57815 0.69675 0.90 1.0 1.25382];
plot(x,y,'o');hold on;
disp([x;y])
n=length(x);%x的长度
v=1;%u的长度，且v的向量
xs=zeros(6);
for k=1:n%一共要的循环
w=1;%w为空的，长度为u
e1=1;
for j=[1:k-1 k+1:n]%少了k的
e=[1 -x(j)];
e1=conv(e1,e);%多项式乘法
w=(x(k)-x(j)).*w; %lagrange系数
end
v=y(k)*e1./w;%总系数
for i=1:n
xs(k,i)=v(i);%系数矩阵
end
end
xs
t=sum(xs)%最后系数矩阵
x=[0.596 0.99];
y2 =t(1)*x.^5+t(2)*x.^4+t(3)*x.^3+t(4)*x.^2+t(5)*x+t(6)
x=0.35:0.01:1.1;
y1 =t(1)*x.^5+t(2)*x.^4+t(3)*x.^3+t(4)*x.^2+t(5)*x+t(6);
%y2=121.6264*x.^5-422.7503*x.^4+572.5667*x.^3-377.2549*x.^2+121.9718*x-15.08452;
plot(x,y1,'b');%画lagrange图
