fmulf.m

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

function fout=fmulf(w,f1,f2)
%FMULF  Multiply two MVFR matrices.
%       FMULF(W,F1,F2) returns an MVFR matrix whose
%       component matrices = F1m*F2m,
%       where F1m and F2m are the component matrices
%       of the MVFR matrices F1 and F2 respectively.
%       W is the associated frequency vector.
%       See also FMUL,FINV,FADD,FADDF

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

[m1,n1]=fsize(w,f1);
[m2,n2]=fsize(w,f2);
lw=length(w);
fout=zeros(lw*m1,n2);
k1=1:m1;                % vector of rows of each matrix in F1
k2=1:m2;                % vector of rows of each matrix in F2
for i=0:lw-1     % for each frequency
  fout(k1+i*m1,:)=f1(k1+i*m1,:)*f2(k2+i*m2,:);
end % for i=0:lw-1