lzhhun 发表于 2009-4-18 18:04

画拉格朗日插值的图像

小弟初学此软件，编了个程序用来画出拉格朗日插值的图像，可是老是出现错误，在这里跪求大侠们的指教了～～～
x=-1:1:1;
for i=1:3
    x0(i)=-1+1*(i-1);
    y0=1/(1+25*x0^2);
    s=0;
    for k=1:3
      L=1;
      for j=1:3
            if j~=k
            L=L*(x-x0(j))/(x0(k)-x0(j));
            end
      end
      s=s+L*y0(k);
    end
    y(i)=s;
end
plot(x,y)

[ 本帖最后由 ChaChing 于 2009-4-18 20:42 编辑 ]

ChaChing 发表于 2009-4-18 20:40

回复 楼主 lzhhun 的帖子

当loop i=1时, x0仅为1*1, loop j怎样赋值至3?
还有建议楼主看下本版规则! 标题要明确!

[ 本帖最后由 ChaChing 于 2009-4-18 20:41 编辑 ]
页: [1]
查看完整版本: 画拉格朗日插值的图像