ratedist.m

Matlab时频分析工具箱,希望能对大家有所帮助啊

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
% 

% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
% 
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
% 
% You should have received a copy of the GNU General Public License
% along with this program.  If not, see <http://www.gnu.org/licenses/>.

% 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);