getexponentbase.m

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

function [exponents,base]=getexponentbase(p,x)
%GETEXPONENTBASE Internal function used in SOS programs

% Author Johan L鰂berg
% $Id: getexponentbase.m,v 1.6 2005/02/12 13:24:39 johanl Exp $

if isempty(p)
    exponents=[];
    base=[];
else
    p_vars = getvariables(p);
    x_vars = getvariables(x);

    base = getbase(p);
    %[sqrList,CompressedList]=yalmip('nonlinearvariables');
    %exponents = spalloc(length(p_vars),length(x_vars),0);

    monom_table = yalmip('monomtable');
    exponents = monom_table(p_vars,x_vars);
    if base(1)~=0
        exponents = [spalloc(1,size(exponents,2),0);exponents];
    else
        base = base(2:end);
    end
end