frgersh.m

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

function c=frgersh(w,f,i,n)
%FRGERSH         Generate row Gershgorin circles.
%       FRGERSH(W,F,i) calculates the row Gershgorin circles for row i
%       of the MVFR matrix, F, and returns one column per frequency
%       which contains the points on the circle at that frequency.
%
%       FRGERSH(W,F,i,n) calls FCOMP(W,F,n) befor
%       calculating the row Gershgorin circles for row i.
%
%       See also FCOMP

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

if nargin==4
   [f,w]=fcomp(w,f,n);
end
cent=fdiag(w,f);
rad=frsod(w,f);
c=circ(cent(:,i),rad(:,i),31);
end