wucha.m

模拟信号数字化的PCM量化以及哈夫曼编码

close all
clear all
echo off
i=1:500;
a=randn(1,500);%产生500个高斯随机数

[sqnr,a_quan,code]=u_pcm(a,8);%量化级数为8
f=a-a_quan;                 %量化误差为输入值和量化值之间的差
subplot(3,1,1);
axis([0 500 -1 1]);
plot(i,f,'-');
title('量化误差,从上到下量化级数分别为8,16,64')

[sqnr,a_quan,code]=u_pcm(a,16);%量化级数为16
f=a-a_quan;                 %量化误差为输入值和量化值之间的差
subplot(3,1,2);
axis([0 500 -1 1]);
plot(i,f,'-');


[sqnr,a_quan,code]=u_pcm(a,64);
f=a-a_quan;                 %量化误差为输入值和量化值之间的差
subplot(3,1,3);
axis([0 500 -1 1]);
plot(i,f,'-');