mulaw.asv

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

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

for i=1:length(x)                   
    if x(i)<1/255           %序列值位于第1折线
        ypcm(i)=255/8*x(i);
    elseif x(i)<3/255       %序列值位于第2折线
        ypcm(i)=255/16*x(i)+1/16;
    elseif x(i)<7/255      %序列值位于第3折线
        ypcm(i)=255/32*x(i)+5/32;
    elseif x(i)<15/255      %序列值位于第4折线
        ypcm(i)=255/64*x(i)+17/64;
    elseif x(i)<31/255     %序列值位于第5折线
        ypcm(i)=255/128*x(i)+49/128;
    elseif x(i)<63/255     %序列值位于第6折线
        ypcm(i)=255/256*x(i)+129/256;
    elseif x(i)<127/2    %序列值位于第7折线
        ypcm(i)=255/512*x(i)+321/512;
    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);
for i=1:size(ypcm)
    if ypcm(i)==128
        ypcm(i)=127;
    end
end
t=dec2bin(ypcm);
f=[s,t];