均匀量化实验matlab,Matlab实现均匀量化的简单实例

教学任务是要读取一个声音文件，对前200个采样点进行4bit的均匀量化。自己在网上看了很多资料，博客，最后还是一脸茫然，最后还是自己一点一点写，

虽然也弄出了很多问题，不过最后还是依靠自己解决了，记录下这一次辛勤工作的结果。

PS:没想到matlab还可以读取声音文件，修改，和播放，再一次体验到matlab实在太强大了。

matlab实现代码如下：

clear all;

clc;

[y,fs,bits]=wavread('road.wav'); %fs得到声音信号的采样率

ft = y(:,1); %假设声音文件有两个声道，我们只分析第一个声道，如果要分析第二个声道可以改为：ft=y(:,2);

sigLength = length(ft); %获取声音长度

sound(y, fs, bits) ;%可以使用sound函数来播放声音

t=(0:sigLength-1)/fs;

figure;

plot(t(1:200), ft(1:200)),grid;

xlabel('Time(s)');

ylabel('Amplitude');

hold on;

min = 999;

max = -999;

M = 16;

for i = 1:200;

if(ft(i) < min)

min = ft(i);

end

if(ft(i) > max)

max = ft(i);

end

end

det = (max - min)/M;

m = zeros(1,M);

for i = 1:M

m(i) = min + i*det;%量化区间的端点

end

%m_i = [min m(1:M-1)]; %

%量化输出电平为量化间隔的中点：qi = (m(i)+m(i-1))/2

%q_i = (m + m_i)/2

qv = zeros(1,200);%存储200个点的均值量化后的数字化的值

node = [min-0.0001,m(1:M-1),max+0.0001];

display(node(1));

for i = 1 : 200; %声音信号的每个采样点

for j = 1:M;

if(node(j) < ft(i) && ft(i) <= node(j+1))

qv(i) = (node(j) + node(j+1))/2;

end

end

%fprintf('qv(%d) = %d\n',i,qv(i));

end

plot(t(1:200), qv(1:200),'r'),grid;

原始声音波形和量化波形效果：