通信原理PCM编译码仿真

通信原理MATLAB PCM编译码仿真，得到语音波形对比

clc;
clear all;
% 产生信号
load handel
x=y（1:20000）;%取前20000个采样点
sound（xFs）;
%PCM编码
x1=x/0.8.*2048;
yy=pcm_encode（x1）;
figure;
subplot（211）;
stem（yy（1:80）‘.‘）;
title（‘PCM编码后的波形‘）;
%加噪声
snr=10;
sp=mean（yy.^2）;
attn=sp./ 10^（snr/10）;
attn = sqrt（attn）;
noise=randn（1length（yy））.*attn;
np=mean（noise.^2）;
snr1=10*log10（sp/np）;
data=yy+noise;
% data=yy;%不加噪声
subplot（212）;
stem（data（1:80）‘.r‘）;
title（‘PCM加噪声后波形‘）;
%译码
demodata=data>0.5;
zz=pcm_decode（demodata0.8）;
figure;
subplot（211）;
plot（x）;
title（‘原始语音信号‘）;
subplot（212）;
plot（zz）;
title（‘译码后的语音信号‘）;
sound（zzFs）;
figure;
plot（x‘b‘）;
hold on
plot（zz‘r‘）;
legend（‘编译码前的语音信号‘‘编译码后的语音信号‘2）;
title（‘编，译码前后的语音信号‘）;

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

     文件        793  2012-12-28 09:29  PCM\pcm.m

     文件        820  2012-12-28 22:04  PCM\pcm_decode.m

     文件       1421  2012-12-16 15:45  PCM\pcm_encode.m

     目录          0  2012-12-16 16:05  PCM

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

                 3034                    4