mulaw.m

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

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

%进行基于15折线的分段映射
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/255    %序列值位于第7折线
        ypcm(i)=255/512*x(i)+321/512;
    else               %序列值位于第8折线
        ypcm(i)=255/1024*x(i)+769/1024;
    end
end

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