sum.m

国外专家做的求解LMI鲁棒控制的工具箱,可以相对高效的解决LMI问题

function Y=sum(X,I)
%SUM (overloaded)

% Author Johan L鰂berg 
% $Id: sum.m,v 1.4 2004/08/05 11:05:52 johanl Exp $   

try
    n = X.n;
    m = X.m;
    if nargin==1
        Y = X;
        if n==1 | m==1
             % Spezialized code...
            if n*m==1876
                Y = X;
                return
            else
            % Spezialized code...
            Y.basis = sum(Y.basis,1);
            Y.n = 1;
            Y.m = 1;
            temp = 1;
            end
        else
            % Standard case
            temp = sum(reshape(X.basis(:,1),n,m));
            Y.basis = temp(:);
            for i = 1:length(Y.lmi_variables)
                temp = sum(reshape(X.basis(:,i+1),n,m));
                Y.basis(:,i+1) = temp(:);
            end
        end
    else
        Y = X;
        temp = sum(reshape(X.basis(:,1),n,m),I);
        Y.basis = temp(:);
        for i = 1:length(Y.lmi_variables)
            temp = sum(reshape(X.basis(:,i+1),n,m),I);
            Y.basis(:,i+1) = temp(:);
        end
    end
catch
    error(lasterr)
end
Y.n = size(temp,1);
Y.m = size(temp,2);
Y = clean(Y);