triu.m

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

function Y=triu(X,r)
%TRIU (overloaded)

% Author Johan L鰂berg 
% $Id: triu.m,v 1.3 2005/02/02 17:10:22 johanl Exp $   


if nargin==1
    r = 0;
end

Y = X;
x_lmi_variables = X.lmi_variables;
lmi_variables = [];
n = X.n;
m = X.m;
triuX = triu(reshape(X.basis(:,1),n,m),r);
Y.basis = triuX(:);

j = 1;
for i = 1:length(x_lmi_variables)
	triuX = triu(reshape(X.basis(:,i+1),n,m),r);
	if (norm(triuX,inf)>0)
		Y.basis(:,j+1) = triuX(:);
		lmi_variables = [lmi_variables x_lmi_variables(i)];
		j = j+1;
	end
end   
Y.n = size(triuX,1);
Y.m = size(triuX,2);
Y.lmi_variables = lmi_variables;
Y = clean(Y);