brqtize.m

这是数字信号处理——matlab释义与实现一书的源代码

function y=brqtize(x,N)
%《数字信号处理教程——MATLAB释义与实现》
% 相对量化子程序
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% function y=brqtize(x,N)
% -------------------------------------------
% 将x数组按各个元素x的相对误差量化为N位二进制
% x--任意实数
% N--量化为二进制的位数
%
ax=abs(x);                  % 去掉符号
m=ceil(log(ax)/log(2));	    % 确定x的幅值占二进制m位
deltax=2.^(m-N); 		    % 求出量化步长，它是一个与x同长的数组
xint=round(ax./deltax);	    %将x除以量化单位，再取整。
% 若将round 改成fix,则属于截断量化
y=sign(x).*xint.*deltax;    % 再将此整数乘以2^(m-N),恢复量化后的原值