convertlorentz.m

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

function [F,changed] = convertlorentz(F)
%catsdp          Internal function: converts ratated Lorentz to SOCC
  
% Author Johan L鰂berg 
% $Id: convertlorentz.m,v 1.3 2005/02/04 10:10:26 johanl Exp $
  
changed = 0;
Counter = size(F.clauses,2);
for i = 1:Counter
    % Yep, Lorentz
    if  (F.clauses{i}.type==5)
        changed=1;
        xyz = F.clauses{i}.data;
        x = xyz(1);
        y = xyz(2);
        z = xyz(3:end);
        F.clauses{i}.data = [(x+y)/sqrt(2);(x-y)/sqrt(2);z];
        if isa(x+y,'sdpvar')
            F = F+lmi(x+y>0);
        end
        F.clauses{i}.type = 4;
    end
end