fpart.m

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

function fp = fpart(w,f,rows,cols)
%FPART  Partitions MVFR matrix into blocks
%
%       fp = fpart(w,f,rows,cols)
%
%       f is an MVFR matrix, and w the associated frequency vector.
%       fp is an MVFR matrix with only the rows and columns specified
%       in rows and cols retained.

%       Dr J.M.Maciejowski, 22 November 1987.
% Copyright (c) 1987 by GEC Engineering Research Centre & Cambridge Control Ltd

if min(size(rows))~=1 | min(size(cols))~=1
  error('Arguments 3 and 4 must be vectors. (FPART)')
end
%
if length(rows)==1 & length(cols)==1    % 1x1 block
  fp = fget(w,f,[rows,cols]);
else
  fp=[];
  lr=length(rows);
  r=zeros(lr,1);
  r(:)=rows;      % Make sure r is a column vector
  for i=1:length(cols)
    fp = [fp, fget(w,f,[r,cols(i)*ones(lr,1)])];
  end
  fp = shpf(fp,length(rows),length(cols));
end