x=randn（11000）;%產生一正態分布序列
xf=fft（x256）;?%1024為數據點數，即對信號采樣數據為1024點的處理
d=x;
subplot（211）;%繪制正態分布序列時域與頻域圖
plot（x）;
title（‘長度為1000的標準正態分布的隨機信號時域波形圖‘）;
subplot（212）;
plot（abs（xf））;
xlabel（‘頻率/Hz‘）;
ylabel（‘振幅‘）;
title（‘正態分布的頻域圖‘）;
figure;?????%打開新的窗口，繪制編碼與解碼信號時域與頻域圖
u=255;
xx=0:0.01:1;
xy=log（1+u*xx）/log（1+u）;
plot（xxxy）;
title（‘μ律PCM編碼輸入輸出關系曲線‘）;

%編碼過程
c=x
for?i=1:100

x（i）=x（i）/4
x（i）=fix（x（i）*4079）;??
s=sign（x（i））;
if?s<0
pcm（（i-1）*8+1）=0;
else
pcm（（i-1）*8+1）=1;
end
x（i）=abs（x（i））;
p（i）=x（i）;

if??x（i）<=15.5

pcm（（i-1）*8+2:（i-1）*8+4）=[000]
else?if?x（i）<=47.5
pcm（（i-1）*8+2:（i-1）*8+4）=[001]
p（i）=floor（（x（i）-15.5）/2）;
else?if?x（i）<=111.5
pcm（（i-1）*8+2:（i-1）*8+4）=[010];
p（i）=floor（（x（i）-47.5）/4）;
else?if?x（i）<=239.5
pcm（（i-1）*8+2:（i-1）*8+4）=[011];
p（i）=floor（（x（i）-111.5）/8）;
else?if?x（i）<=495.5
pcm（（i-1）*8+2:（i-1）*8+4）=[100];
p（i）=floor（（x（i）-239.5）/16）;
else?if?x（i）<=1007.5
pcm（（i-1）*8+2:（i-1）*8+4）=[101];
p（i）=floor（（x（i）