%---------------------------------------------------------
%%%%%%%%%%%%%%%?PCM?A率13折線編碼?%%%%%%%%%%%%%%%%%%%%%%%%%
%---------------------------------------------------------

clcclear;
%A=87.6??y=Ax/（1+lnA）?（0dx=0.01;
x=-1:dx:1;
A=87.6;
%A_law
for?i=1:length（x）
????if?abs（x（i））<1/A
????????ya（i）=A*x（i）/（1+log（A））;
????else
????????ya（i）=sign（x（i））.*（1+log（A*abs（x（i）））/（1+log（A）））;
????end
end
figure;
plot（xya‘k‘）;
title（‘A_law‘）;
xlabel（‘x‘）;
ylabel（‘y‘）;
grid?on;
hold?on
xx=[-1-1/2-1/4-1/8-1/16-1/32-1/64-1/1281/1281/641/321/161/81/41/21];
yy=[-1-7/8-6/8-5/8-4/8-3/8-2/8-1/81/82/83/84/85/86/87/81];
plot（xxyy‘r‘）;
stem（xxyy‘b-‘）;
legend（‘A_law‘）;







%%
%下面對語音信號編碼

load?handel
x=y（1:20000）;%取handel語音信號的前20000個采樣點
sound（xFs）;
figure;
plot（x）;
title（‘原始語音信號‘）;

%編碼
x1=x/0.8.*4096;
yy=pcm_encode（x1）;%得到PCM編碼輸出yy（調用編碼函數）
figure;
stem（yy（1:80））;
title（‘前10個采樣點的PCM編碼‘）;

%receive=yy;?%不加噪聲時接收到的碼直接等于編碼輸出yy

%高斯信道
receive?=?awgn（yy10）;%對PCM碼加入信噪比10dB的高斯噪聲
figure;
stem（receive（1:80））;
title（‘加噪聲后?前10個采樣點的PCM碼‘）;

%譯碼
for?i=1:160000??%對收到的碼信號進行判決，還原為PCM碼
????if?receive（i）<0.5
????????receive（i）=0;
????else?receive（i）=1;
????end
end
zz=pcm_decode（receive0.8）;%得到由receive信號譯碼的語音信號（調用譯碼函數）
figure;
plot（zz）;
title（‘譯碼得到的語音信號‘）;
sound（zzFs）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1533??2010-12-26?16:57??第六章仿真作業?01089002李濤\my.m

?????文件????????526??2010-12-26?11:12??第六章仿真作業?01089002李濤\pcm_decode.m

?????文件???????1474??2010-12-26?15:31??第六章仿真作業?01089002李濤\pcm_encode.m

?????文件?????138244??2010-12-26?17:05??第六章仿真作業?01089002李濤\第六章仿真作業.docx

?????目錄??????????0??2010-12-26?17:05??第六章仿真作業?01089002李濤

-----------?---------??----------?-----??----

???????????????141777????????????????????5