genmonoms.m

来自「国外专家做的求解LMI鲁棒控制的工具箱,可以相对高效的解决LMI问题」· M 代码 · 共 24 行

function monom = genmonoms(newton_m,x);
%GENMONOMS Internal function used in SOS programs

% Author Johan L鰂berg
% $Id: genmonoms.m,v 1.4 2004/10/01 11:48:10 johanl Exp $


precalc = [];
monom = [];
for i = 1:size(newton_m,2)
    powers = unique(newton_m(:,i));
    for j = 1:length(powers)         
        precalc{1+powers(j),i} = x(i)^powers(j);
    end   
end

monom = [];
for i = 1:size(newton_m,1)
    temp = 1;
    for j = 1:size(newton_m,2)
        temp = temp*precalc{newton_m(i,j)+1,j};%x(j)^newton_m(i,j);
    end
    monom = [monom;temp];
end