tfr_mul.m

linear time－frequency toolbox

function h=tfr_mul(self,f,c,L)
%TFR_MUL  Multiplier
%   Usage:  h=tfr_mul(tfr,f,c)
%           h=tfr_mul(tfr,f,c,L)
%
%   TFR_MUL(tfr,f,c) will filter f by a multiplier using the symbol c.
%   The symbol c can be either a vector or a matrix.
%
%   TFR_MUL(tfr,f,c,L) do the same using a transform of length L.
%

error(nargchk(3,4,nargin));

if nargin<4
  L=-1;
else
  if  (prod(size(L))~=1 || ~isnumeric(L))
    error([callfun,': L must be a scalar']);
  end;
  if rem(L,1)~=0
    error([callfun,': L must be an integer']);
  end;
end;

c=c(:);

[c_f,Ls]=comp_a(self,f,0,'TFR_MUL',L);

if length(c)~=size(c_f,1)
  error('Symbol size mismatch with transform size.');
end;

c_f=c_f.*c;

% Let comp_s do all the work.
h=comp_s(self,c_f,0,'TFR_MUL',Ls);