alaw.asv

pcm编码 有关线性、A律、mu律的pcm编码的matlab实现

function f=ecode(yn)
%本函数实现国际通用的PCM量化A律13特性近似编码
x=yn;
s=sign(x);                              %求得每一序列值的符号
x=abs(x);                               %求序列的绝对值
ypcm=zeros(length(x),1);

for i=1:length(x)                      %直接将序列的绝对值量化
    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];