tfr_ar.m

linear time－frequency toolbox

function [c,Ls]=tfr_ar(self,f,L)
%TFR_AR  Analysis transform, rectangular layout
%   Usage:  c=tfr_ar(tfr,f)
%           c=tfr_ar(tfr,f,L)
%           [c,Ls]=tfr_ar(tfr,f);
%           [c,Ls]=tfr_ar(tfr,f,L);
%
%   Input parameters:
%         tfr   : Time/frequency representation handle.
%         f     : Input data
%         L     : Length of transform to do.
%   Output parameters:
%         c     : Array of coefficients.
%         Ls    : Length of input signal.
%
%   TFR_AR(tfr,f) computes the analysis transform of the signal f as
%   speficied by tfr. The coefficients are returned in a rectangular layout.
%
%   The length of the transform will be the smallest possible that is
%   larger than the signal. f will be zero-extended to the length of the 
%   transform. If f is a matrix, the transformation is applied to each column.
%
%   The length of the transform done can be obtained by L=size(c,2)*a;   
%
%   TFR_AR(tfr,f,L) computes the analysis transform as above, but does
%   a transform of length L. f will be cut or zero-extended to length L
%   before the transform is done.
%
%   [c,Ls]=TFR_AR(tfr,f) or [c,Ls]=TFR_AR(tfr,f,L) additionally returns the
%   length of the input signal f. This is handy for reconstruction:
% 
%                [c,Ls]=tfr_ar(tfr,f)
%                fr=tfr_sr(c,Ls);
% 
%   will reconstuct the signal f no matter what the length of f is, provided
%   that the analysis and synthesis window of tfr are dual windows.
%
%   Not all transforms support this method. The supported transforms are:
%     'dgt', 'dgtii', 'dgtiii', 'dgtiv'
%     'rdgt', 'rdgtiii'
%     'edgtii'
%
%   SEE ALSO:  TFR_CREATE, TFR_SR, TFR_CLEAR

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

% Let comp_a do all the work.
if nargin==2
  [c,Ls]=comp_a(self,f,1,'TFR_AR',-1);
else
  [c,Ls]=comp_a(self,f,1,'TFR_AR',L);
end;