fcomp.m

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

function [fc,wc]=fcomp(w,f,n)
%FCOMP  Compress MVFR matrix.
%       FCOMP(W,F,n) returns the first component matrix of
%       the MVFR matrix, F, and every Nth+1 component
%       matrix thereafter.
%       The input argument n is optional.
%       If n is missing the default of 10 is used.
%
%       [FC,WC]=FCOMP(W,F,n) also returns the
%       associated compressed frequency vector.
%       See also COMP,EXPD

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

nargs=nargin;
error(nargchk(2,3,nargs));
if nargs <3
      n = 10;
end
[mf,nf]=fsize(w,f);
if n==1
   fc=f
else
  k=1:mf;
  p = length(w);
  p = (1:p);
  p = rem(p,n);
  p = p==1;
  p=expd(p,mf);   % expand p to length(w)
  fc = f(p,:);
end % if n==1

if nargout==2
   wc=comp(w,n);
end