对于matlab，大家应该比较熟悉。今天小哥就带大家来做均匀量化。
量化是将模拟信号转化为数字信号必不可少的一步。均匀量化相对比较简单，就是比如将信号在[-1,1]分成相同的很多段，要是某个值在其中某一段，就将其的值归为那一类。每一类都可以用一个唯一的二进制编码表示。
现在我们要做的就是对wav的音频信号进行量化。具体一点：在road.wav文件左声道取前200个点做4bit量化，并且用matlab作图比较量化前和量化后的图。那么第一步是什么呢？那就是获取样本点？不是的，我们第一步是获取我们需要量化的音频。可以用
step1:
我们可以用[y,fs,nbits]= wavread('F:\a.wav',1024);即可得到我们所需处理的样本点。
step2:对采集到的信号进行量化。
方法一：直接使用量化函数->quantizer(做一个量化器)，quantize(对某信号进行用量化器量化)
syms sample_point;
syms quantizion_bits;
sample_point = 200;
quantizion_bits = 5;
[y,fs,nbits]= wavread('road.wav',sample_point);
sample = y(1:sample_point);
n = 1:sample_point;
q_riser = quantizer('fixed','round','saturate',[10,quantizion_bits]); %定义一个量化器
q_tread = quantizer('fixed','ceil','saturate',[10,quantizion_bits]); %定义一个量化器
mid_riser = quantize(q_riser,sample);%对采样信号进行量化
mid_tread = quantize(q_tread,sample);%对采样信号进行量化
subplot(1,2,1); plot(n,sample,'r',n,mid_riser,'b'); title('mid riser');
grid;
subplot(1,2,2); plot(n,sample,'r',n,mid_tread,'b'); title('tread riser');
grid;
方法二：自己量化，更加灵活。
syms sample_point;
sample_point = 200; %所需采集的声音样本点的数量
[y,fs,nbits]= wavread('road.wav',sample_point); %提取出音频信号的前200个点
sample = y(1:sample_point); %提取出左声道的信号
n = 1:sample_point;
mid_riser = ones(sample_point,1);
mid_tread = ones(sample_point,1);
delta = 0.7/16; %采样间隔
for i = 1:sample_point %采用mid-riser的方法
inteval_rise = floor((sample(i) + 0.4)/delta); %计算有的采样间隔数
inteval_tread = floor((sample(i) + 0.4 + delta/2 )/delta); %
mid_riser(i) = delta* (2*inteval_rise+1)*0.5 - 0.4; %mid_riser的方法
mid_tread(i) = delta*inteval_tread - 0.4; %mid_tread的方法
end
subplot(1,2,1); plot(n,sample,'r',n,mid_riser,'b'); title('mid riser');
grid;
subplot(1,2,2); plot(n,sample,'r',n,mid_tread,'b'); title('tread riser');
grid;
参考函数：
wavread（wavread的具体使用，举例说明）：
例子一：[Y,F,b]=wavread('test.wav');
Y2=Y(((F*10+1):F*20),:);
wavwrite(Y2,F,b,'test_new.wav') ;%你在开始可能觉得无从入手，
%因为你要把绝对路径加进去
%以上做的就是将这个音频的10~20秒给Y2。加入我们听到
%一首歌需要对其中某一部分进行操作，就可以这样做。
例子二：[y,fs,nbits]= wavread(wavFile1);
%y就是音频信号；
%fs是采样频率，比如说16000就是每秒16000次；
%nbit是采样精度，比如说16就是指16位的二进制码去表示数据；可以表示精
wavread('F:\a.wav',1024)%读取该音频文件前1024个采样点
例子三：Y=wavread(FILE);读取文件FILE里面的数据，是wav结尾的，如果没有加wav，
默认是wav
例子四：[...]=wavread(FILE,N);返回的是每个声道的开始N个样本点
例子五：[...]=wavread(FILE,[N1 N2]) 返回每个声道的从N1~N2的样本点