fmul.m

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

function fout=fmul(w,f,c)
%FMUL   Multiply MVFR matrix by a matrix.
%       FMUL(W,F,C) returns F*C
%       FMUL(W,C,F) returns C*F
%       where F is the MVFR matrix
%       W is the associated frequency vector
%       and C is a matrix compatible with the
%       component matrices of F.
%       See also FMULF,FINV,FADD,FADDF

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

[mf,nf]=size(f);
[mc,nc]=size(c);
lw=length(w);
if (rem(mf,lw)==0) & (mf/lw > 1) &...
   (rem(mc,lw)==0) & (mc/lw > 1)
  crlf=setstr([13,10]);
 disp(['Fmul Warning:  Both matrix arguments are consistent MVFR matrices.',...
       crlf,'               Did you mean to use  fmulf ?'])
end
if mc/nf==lw   % then C is the MVFR matrix
  [m,n]=fsize(w,c);
  fout=zeros(lw*mf,nc);
  j=1:mf;
  k=1:m;             % vector of rows of each matrix in C
    for i=0:lw-1     % for each frequency
      fout(j+i*mf,:)=f*c(k+i*m,:);
    end % for i=0:lw-1
else           %  F is the MVFR matrix
  [m,n]=fsize(w,f);
  fout=f*c;
end  % if mc/nf==lw