coefmatrix.m

可以用这个m文件构造5个因数4次的完全多项式（包括1次项2次项和3次项和交叉项）

function [xx]=coefmatrix(x);
[m n] =size(x);
xx=ones(n,1);

for i=1:m
    xx=[xx x(i,:)'];
end

for i=1:m
    xx=[xx x(i,:)'.^2];
end

for i=1:m
    xx=[xx x(i,:)'.^3];
end

for i=1:m-1
    for j=i+1:m
        last=[];
        for k=1:n
            last=[last x(i,k)'*x(j,k)'];
        end
        last=last';
        xx=[xx last];
    end
end

for i=1:m-1
    for j=i+1:m
        last=[];
        for k=1:n
            last=[last x(i,k)'*(x(j,k)'.^2)];
        end
        last=last'
    xx=[xx last];
    end
end

for i=1:m-1
    for j=i+1:m
        last=[];
        for k=1:n
            last=[last (x(i,k)'.^2)*x(j,k)'];
        end
        last=last'
    xx=[xx last];
    end
end

for i=1:m-2
    for j=i+1:m-1
    for k=j+1:m
        last=[];
        for l=1:n
            last=[last x(i,l)*x(j,l)'*x(k,l)'];
        end
        last=last';
        xx=[xx last];
    end
    end
end