coefmatrix4.m

来自「可以用这个m文件构造5个因数4次的完全多项式（包括1次项2次项和3次项和交叉项）」· M 代码 · 共 155 行

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
    xx=[xx x(i,:)'.^4];
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

for i=1:m-1
    for j=i+1:m
        last=[];
        for k=1:n
            last=[last x(i,k)*x(j,k).^3]
        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).^3*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).^2*x(j,k).^2]
        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).^2*x(k,l)];
        end
        last=last';
        xx=[xx last];
    end
    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).^2];
        end
        last=last';
        xx=[xx last];
    end
    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).^2*x(j,l)*x(k,l)];
        end
        last=last';
        xx=[xx last];
    end
    end
end

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