matlab 差分 小程序 求解释 谢谢

chinasnowlee

function[p,q]=cs(x,y)
n=length(x);
p(:,1)=x;
p(:,2)=y;
for j=3:n+1,
    p(1:n+2-j,j)=diff(p(1:n+3-j,j-1))./(x(j-1:n)-x(1:n+2-j));   
end
q=p(1,2:n+1)';

chinasnowlee

for j=3:n+1,
    p(1:n+2-j,j)=diff(p(1:n+3-j,j-1))./(x(j-1:n)-x(1:n+2-j));   
end
程序没问题，我运行过。这步没看明白？？多谢指导

ht465331

chinasnowlee 发表于 2011-11-6 22:10
for j=3:n+1,
    p(1:n+2-j,j)=diff(p(1:n+3-j,j-1))./(x(j-1:n)-x(1:n+2-j));   
end

这是求一个高次多项式系数的值的问题    这个程序很巧妙  你自己再摸索一下   详细的说  比如是7次  给出了8个点的值（x，y），然后写了这么一个程序来求x的7次方、6次方、...、1次方前面的系数值