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