subsasgn.m

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

function X = subsasgn(X,I,Y)
% SUBSASGN (overloaded)

% Author Johan L鰂berg
% $Id: subsasgn.m,v 1.2 2006/07/13 20:37:30 joloef Exp $

base = reshape(1:size(X.basis,1),X.dim);
base = subsref(base,I);

if isa(Y,'ndsdpvar')
    Y = sdpvar(Y);
elseif isa(Y,'double')
    Y = Y(:);
end

dim = X.dim;
X = sdpvar(X);
X(base) = Y;
X = reshape(X,dim);
X = clean(X);