pcm编码解码matlab实现

对信号进行抽样、量化和A律PCM编码，经过传输后，接收端进行PCM译码。要求画出不同幅度下PCM编码、译码后的波形以及未编码波形。

clc;
clear;
t = 0:0.01:10;
vm1 = -70:1:0; %输入的正弦信号幅度不同
vm = 10.^（vm1/20）;
figure（1）
for k = 1:length（vm）
for m = 1:2
x = vm（k）*sin（2*pi*t+2*pi*rand（1））;
v = 1;
xx = x/v;%normalize
sxx = floor（xx*4096）;
y = pcm_encode（sxx）;
yy = pcm_decode（yv）;

nq（m） = sum（（x-yy）.*（x-yy））/length（x）;
sq（m） = mean（yy.^2）;
snr（m） = （sq（m）/nq（m））;

drawnow
subplot（211）
plot（tx）;
title（‘sample sequence‘）;
subplot（212）
plot（tyy）
title（‘pcm decode sequence‘）;
end
snrq（k）=10*log10（mean（snr））;
end

figure（2）
plot（vm1snrq）;
grid;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        557  2009-05-11 22:24  pcm\pcm.m

     文件        564  2009-05-11 22:13  pcm\pcm_decode.m

     文件       1109  2009-05-11 23:13  pcm\pcm_encode.m

     目录          0  2009-05-12 12:53  pcm

----------- ---------  ---------- -----  ----

                 2230                    4