Jessi00的博客

原创 两种插值计算方法

拉格朗日插值计算 function fx=lagrange(n) m=linspace(-5,5,n); y=1./(m.^2+1); syms x; %定义一个变量x fx=zeros(size(x)); %设fx初始为一个x阶0矩阵 for k=1:n lx=ones(size(x)); %设lx初始为一个x阶1矩阵 for j=[1:k-1 k+1:n] lx=lx.(x-m(j))/(m(k)-m(j)) end fx=fx+y(k)lx;