samer@4: function Y=h_todbs(X,dBs)
samer@4: % todbs - Convert values to decibels 
samer@4: %
samer@4: % todbs :: 
samer@4: %    real ~'dimensionless value', 
samer@4: %    real ~'dynamic range in dB'
samer@4: % -> real ~'value in dBs'.
samer@4: 
samer@4: if nargin>1,
samer@4: 	mmax=max(X(:));
samer@4: 	X=max(X,10^(-dBs/10)*mmax);
samer@4: end
samer@4: Y=10*log10(X);
samer@4: