alaw.m

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

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

%进行基于13折线的分段映射
for i=1:length(x)                   
    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

%舍掉后八位
ypcm=ypcm.*(2^7);    
ypcm=floor(ypcm);
ypcm=ypcm.*s;
%写入编码压缩文件
fid=fopen([name,'.a'],'wb');
fwrite(fid,ypcm,'int8');
fclose(fid);