minus.m

来自「求解线性矩阵不等式简单方便--与LMI工具箱相比」· M 代码 · 共 20 行

function sys = minus(X,Y)
%MINUS (overloaded)

% Author Johan L鰂berg
% $Id: minus.m,v 1.1 2004/12/07 08:46:53 johanl Exp $

if isempty(X)
    sys = set([]);
end

if isempty(Y)
    sys = X;
end

idX = getlmiid(X);
idY = getlmiid(Y);
YinX = find(~ismember(idX,idY));
% Get the correct subsref...
sys = subsref(X,struct('type','()','subs',{{YinX}}));