dbp.m

来自「高阶sigma-delta调制器设计matlab工具包, 半波带滤波器设计工具」· M 代码 · 共 9 行

function y=dbp(x)
% dbp(x) = 10*log10(x): the dB equivalent of the power x
y = -Inf*ones(size(x));
if isempty(x)
    return
end
nonzero = x~=0;
y(nonzero) = 10*log10(abs(x(nonzero)));