function f=fdly(w,f,dly)
%FDLY   Adds delays to MVFR matrix.
%       FDLY(W,F,DLY) returns an MVFR matrix
%       in which each component matrix of the MVFR matrix
%       F has been multiplied element-wise by exp(-j*w*dly).
%       The matrix DLY specifies the delay in seconds
%       for each element of F.

% Copyright (c) 1987 by GEC Engineering Research Centre & Cambridge Control Ltd

[m,n]=fsize(w,f);
[md,nd]=size(dly);
if (m~=md)|(n~=nd)
   error('Size of delay matrix not equal to component matrices of MVFR.');
end
j=sqrt(-1);
k=1:m;
for i=1:length(w)
    f((i-1)*m+k,:)=f((i-1)*m+k,:).*exp(dly.*(-j*w(i)));
end