view general/numerical/scalar/todbs.m @ 61:eff6bddf82e3 tip

Finally implemented perceptual brightness thing.

function Y=h_todbs(X,dBs)
% todbs - Convert values to decibels 
%
% todbs :: 
%    real ~'dimensionless value', 
%    real ~'dynamic range in dB'
% -> real ~'value in dBs'.

if nargin>1,
	mmax=max(X(:));
	X=max(X,10^(-dBs/10)*mmax);
end
Y=10*log10(X);