function f=ecode(yn)
%本函數(shù)實(shí)現(xiàn)國(guó)際通用的PCM量化A律13特性近似編碼
x=yn;
s=sign(x);                              %求得每一序列值的符號(hào)
x=abs(x);                               %求序列的絕對(duì)值
ypcm=zeros(length(x),1);

for i=1:length(x)                      %直接將序列的絕對(duì)值量化
    if mod(i,10000)==0
        i
    end
    if x(i)<1/64           %序列值位于第1和第2折線
        ypcm(i)=16*x(i);
    elseif x(i)<1/32      %序列值位于第3折線
        ypcm(i)=8*x(i)+1/8;
    elseif x(i)<1/16      %序列值位于第4折線
        ypcm(i)=4*x(i)+2/8;
    elseif x(i)<1/8     %序列值位于第5折線
        ypcm(i)=2*x(i)+3/8;
    elseif x(i)<1/4     %序列值位于第6折線
        ypcm(i)=x(i)+4/8;
    elseif x(i)<1/2    %序列值位于第7折線
        ypcm(i)=1/2*x(i)+5/8;
    else               %序列值位于第8折線
        ypcm(i)=1/4*x(i)+6/8;
    end
end

for i=1:size(s)
    if s(i)==1
        s(i)=1;
    else
        s(i)=0;
    end
end
s=dec2bin(s);
ypcm=ypcm.*(2^7);
ypcm=floor(ypcm);
t=dec2bin(ypcm);
f=[s,t];