lloyd_sqnr.asv

u率pcm编码,还有误差分析,通信原理实验的一个例子

clear all
close all

m=randn(1,500);%产生(0,1)高斯随机数
[a,y]=lloydmax('normal',20,10,0.001,0,1);%生成最佳量化器,得到量化级数和量化边界
for i=1:length(a)-1   %用最佳量化器对生成的序列进行量化
  m_quan(find((a(i)<=m)&(m<=a(i+1))))=...
  y(i).*ones(1,length(find((a(i)<=m)&(m<=a(i+1)))));
end
sqnr=20*log10(norm(m)/norm(m-m_quan))%求量化信噪比