frsod.m

控制系统计算机辅助设计——MATLAB语言与应用(源代码）

function s=frsod(w,f)
%FRSOD  Row sum of abs(off-diagonals) of MVFR matrix
%       FRSOD(W,F) returns the row sum of the absolute
%       values of off diagonals of the component matrices of
%       the MVFR matrix, F, with one row per component matrix.
%       W is the associated frequency vector.
%       See also FCSOD,FRDOM,FCDOM

%       Dr M.P. Ford 4th August 1987
% Copyright (c) 1987 by GEC Engineering Research Centre & Cambridge Control Ltd

[m,n]=fsize(w,f);
lw=length(w);
if m~=n
  error('Not a square system');
end
s=ones(lw,m)*inf;   % Set up output matrix
k=1:m;             % vector of rows of each matrix in F
for i=1:lw          % for each frequency
  fm=abs(f(k+(i-1)*m,:))';  % column sum of the transpose
  s(i,:)=sum(fm)-diag(fm)';
end % for i=1:lw