MAQ多比特自适应量化
%MAQ量化方法,m=1的情况
function key=MAQ(alpha_vec,m)
K=2^(m+2);
alpha_sort=sort(alpha_vec); %对采样值排序
for c=1:K
eta©=alpha_sort(length(alpha_sort)*c/K); %确定门限,等间隔
end
for i=1:length(alpha_vec)
for j=1:K
if alpha_vec(i)<=eta(j)
k(i)=j;
break
end
end
end
for c=1:K
if (mod(c,4))>=2
e©=1;
else
e©=0;
end
end
% 格雷码
d0=[0 0 1 1 1 1 0 0];
d1=[0 0 0 0 1 1 1 1];
for i=1:length(alpha_vec)
if e(k(i))==0
key(i)=d0(k(i));
else
key(i)=d1(k(i));
end
end