求matlab代码

显示全部楼层 · 2013-6-29 21:59:39

先定义lagrange插值函数，保存以下程序，先保存function到end这一段函数，后面再在命令窗口运行后面的程序function y=lagrange(x0,y0,x)ii=1:length(x0);y=zeros(size(x));for i=ii
ij=find(ii~=i); y1=1;
for j=1:length(ij),y1=y1.*(x-x0(ij(j)));end
y=y+y1*y0(i)/prod(x0(i)-x0(ij));end计算x=[0.4:0.1:0.8];y=[-0.756291,-0.443147,-0.150826,0.133325...

千问 · 2013-6-29 21:59:39

function [C,L,L1,l]=lagran1(X,Y)m=length(X);L=ones(m,m);for k=1:m
V=1;
for i=1:m
if k~=i
V=conv(V,poly(X(i)))/(X(k)-X(i));
endendL1(k,:)=V;l(k,...