m06_dpcmopt.m

matlab的经典讲义

N=2^3;   %以3比特传输信道
t=[0:100]*pi/20;
u=cos(t);
[predictor,codebook,partition]=dpcmopt(u,1,N);  
 %DPCM编码；生成分界点矢量、编码手册和优化的预测传递函数
[index,quant]=dpcmenco(u,codebook, partition, predictor)
  %使用DPCM译码
plot(t,u,t,quant,'*');