truncation.m

来自「胡广书DSP信号处理光盘C fortran和matlab算法实现」· M 代码 · 共 25 行

%---------------------------------------------------------------------------
% truction: 完成给定序列的量化，并计算量化误差的信噪比。量化采用截尾方式。
% 输入   a: 需量化的序列，
%        n: 浮点位数
% 输出 aeq: 量化后的序列
%---------------------------------------------------------------------------
function aeq=truncation(a,n)

[m1,m2]=size(a);
for i = 1:m1
    for j = 1:m2
        r=floor(log2(abs(a(i,j))+eps));
        r = max(0,min(r)+1);
        b = floor(abs(a(i,j))*2^(n-r));
        aeq(i,j) = sign(a(i,j))*b/2^(n-r);
     end
end