pcmcode.m

来自「用MATLAB做得fft谱分析」· M 代码 · 共 22 行

function f=pcmcode(y)
f=zeros(length(y),8);
z=sign(y);
y=y*128;
y=fix(y);
y=abs(y);
for i=1:length(y)
    if(y(i)==128)
        y(i)=127.999;
    end;end
for i=1:length(y)
    for j=6:-1:0
        f(i,8-j)=fix(y(i)/(2^j));
        y(i)=mod(y(i),(2^j));
    end;end
for i=1:length(y)
    if(z(i)==-1)
        f(i,1)=0;
    else
        f(i,1)=1;
    end;end