lsmq.m

lu 分解 b样条曲线

x=[1 3 4 5 6 7 8 9 10];
y=[10 5 4 2 1 1 2 3 4];
n=length(x);
m=input('cishu:');
p=m+1;
xx=zeros(p,p);
for i=1:p
    for j=1:p
        for k=1:n
            xx(i,j)=xx(i,j)+x(k)^(i+j-2);
        end
    end
end
xy=zeros(p,1);
for i=1:p
    for k=1:n
        xy(i,1)=xy(i,1)+y(k)*x(k)^(i-1);
    end
end
xy=inv(xx)*xy