%例程7－15  PCM編碼
%首先，A律13折線量化的Matlab程序?qū)崿F(xiàn)如下： 

function y=pcmA13(x)
%A律13折線壓縮特性
x=x/max(x);                           % x歸一化
z=sign(x);                            %語音信號是雙極性信號，負(fù)方向?qū)ΨQ
x=abs(x);
for i=1:length(x)
    if (x(i)>=0)&(x(i)<1/64)          %第一、二段斜率均為16，視為一條線段
        y(i)=16*x(i);
    elseif (x(i)>=1/64)&(x(i)<1/32)
        y(i)=8*x(i)+1/8;
    elseif (x(i)>=1/32)&(x(i)<1/16)
        y(i)=4*x(i)+2/8;
    elseif (x(i)>=1/16)&(x(i)<1/8)
        y(i)=2*x(i)+3/8;
    elseif (x(i)>=1/8)&(x(i)<1/4)
        y(i)=x(i)+4/8;
    elseif (x(i)>=1/4)&(x(i)<1/2)
        y(i)=1/2*x(i)+5/8;
    elseif (x(i)>=1/2)&(x(i)<=1)
        y(i)=1/4*x(i)+6/8;
    end
end
y=z.*y;

%接下來8位二進(jìn)制編碼的MATLAB程序?qū)崿F(xiàn)如下：
function f=pcmcode(y)
%8位二進(jìn)制PCM編碼
f=zeros(length(y),8);
z=sign(y);
y=y*128;
y=fix(y);
y=abs(y);
for i=1:length(y)
    if y(i)==128
        y(i)=127.999;
    end
end
for i=1:length(y)
    for j=6:-1:0
        f(i,8-j)=fix(y(i)/(2^j));
        y(i)=mod(y(i),(2^j));
    end
end
for i=1:length(y)
    if z(i)==-1
        f(i,1)=0;
    else 
        f(i,1)=1;
    end
end