rdivide.m

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

function y = rdivide(X,Y)
%RDIVIDE (overloaded)

% Author Johan L鰂berg 
% $Id: rdivide.m,v 1.5 2005/07/18 15:01:30 joloef Exp $   

% Check dimensions
[nx,mx] = size(X);
[ny,my] = size(Y);
if ~((prod(size(X))==1) | (prod(size(Y))==1))
    if ~((nx==ny & (mx ==my)))
        error('Matrix dimensions must agree.')
    end
end

% Quick exit for simple case X/scalar
if isa(Y,'double') & prod(size(Y))==1
    y = X;
    y.basis = y.basis/Y;
    return
end

% FIX : SLOOOOW BUT SOMEWHAT ROBUST
[nx,mx] = size(X);
[ny,my] = size(Y);
if prod(size(X)) == 1 & prod(size(Y))~=1
    X = repmat(X,ny,my);
end;
if prod(size(Y)) == 1 & prod(size(X))~=1
    Y = repmat(Y,nx,mx);
end;
[nx,mx] = size(X);
y = [];
for i = 1:nx
    ytemp = [];
    for j = 1:mx
        dummy = struct('type','()','subs',{{i,j}});
        ytemp = [ytemp subsref(X,dummy)*subsref(Y,dummy)^-1];
    end
    y = [y;ytemp];
end