tril.m

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

function Y=tril(X,r)
%TRIL (overloaded)

% Author Johan L鰂berg 
% $Id: tril.m,v 1.4 2005/02/23 17:44:48 johanl Exp $   

if nargin==1
    r = 0;
end

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

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