ratedist.m

linear time－frequency toolbox

function RateDist(sig,a,M,R)
%RATEDIST  Rate distortion curve for Gabor coefficients.
%   Usage:  RateDist(sig,a,M,R);
%
%   Input parameters:
%      sig      : Input signal
%      a,M      : Parameters for dgt
%      R        : List of bitrates to be examined
%
%   RATEDIST(sig,a,M,R) computes rate distortion curves
%   for Gabor coefficients
%

% Author: Bruno Torresani

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

signorm = norm(sig);

gf=tfr_create('dgt',a,M,'pgauss','candual');
gc = tfr_ar(gf,sig);

SNR = zeros(max(R),1);

for Nbits = R,
  qgc = uquant(gc,Nbits);
  rsig = tfr_sr(gf,qgc);
  err = norm(rsig - sig);
  SNR(Nbits) = 20*log10(signorm/err);
end

plot(R,SNR(R));

tfr_clear(gf);