brqt.m

《数字信号处理实验(MATLAB版)》这

function [b1,a1]=brqt(b,a,N)
%二进制相对量化子程序
bx=abs(b); ax=abs(a);       %去掉符号
bxm=max(bx);axm=max(ax);
maxba=max(bxm,axm);         %寻找系数中最大的数
m=ceil(log2(maxba));        %确定幅度值占二进制m位
deltax=2.^(m-N);            %求出量化步长
aint=round(ax./deltax);     %将a除以量化单位,再取整.如将round改为fix则为截断量化.
bint=round(bx./deltax);     %将b除以量化单位,再取整.
b1=sign(b).*bint.*deltax;   %恢复为量化后的10进制数
a1=sign(a).*aint.*deltax;   %恢复为量化后的10进制数