finv.m

MFD－多变量系统频域设计工具

function fout=finv(w,f)
%FINV   Invert an MVFR matrix
%       FINV(W,F) inverts the MVFR matrix, F.
%       If the component matrices of F are square 
%       INV is used, otherwise PINV is used.
%       W is the associated frequency vector.
%       See PINV,INV

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


[m,n]=fsize(w,f);
lw=length(w);
fout=zeros(n*lw,m);   % Set up output matrix
k=1:m;               % vector of rows of each matrix in F
p=1:n;
if m==n           % square use INV
  for i=0:lw-1     % for each frequency
    fout(p+i*n,:)=inv(f(k+i*m,:));
  end % for i=0:lw-1
else              % use Pinv
  for i=0:lw-1     % for each frequency
    fout(p+i*n,:)=pinv(f(k+i*m,:));
  end % for i=0:lw-1
end % if m=n