ishermitian.m

来自「国外专家做的求解LMI鲁棒控制的工具箱,可以相对高效的解决LMI问题」· M 代码 · 共 26 行

function issym=ishermitian(X)
%ISHERMITIAN Check if variable is Hermitian

% Author Johan L鰂berg 
% $Id: ishermitian.m,v 1.4 2005/04/27 15:47:18 joloef Exp $  

n = X.n;
m = X.m;
issym = 0;
if (n==m)
    % What are the linear indicies to the transposed matrices
    ind = reshape(reshape(1:n^2,n,n)',n^2,1);   
    
    if ~isreal(X.basis)
        residual = X.basis-conj(X.basis(ind,:));
    else      
        residual = X.basis-X.basis(ind,:);
    end

    if nnz(residual)>0
        issym = norm(residual,1)<1e-10;
    else
        issym = 1;
    end
end