声音信号的A律13折线（格雷码）编码仿真

function code=PCMcoding(S) z=sign(S); MaxS=max(abs(S)); S=abs(S/MaxS); %对采样信号S进行归一化 Q=2048*S; %对采样信号S进行量化，分为2048个区间， %采样13折线律作为压缩律进行非均匀量化 code=zeros(length(S),8); %代码存储矩阵（全零） % 段落码判断程序 for i=1:length(S) if (Q(i)>=128)&&(Q(i)<=2048) code(i,2)=1; %对段位码的第一位进行判断 end if (Q(i)>32)&&(Q(i)<128)||(Q(i)>=512)&&(Q(i)<=2048) code(i,3)=1; %对段位码的第二位进行判断 end if (Q(i)>=16)&&(Q(i)<32)||(Q(i)>=64)&&(Q(i)<128)||(Q(i)>=256)&&(Q(i)<512)||(Q(i)>=1024)&&(Q(i)<=2048) code(i,4)=1; %对段位码的第三位进行判断 end end %段内码判断程序 N=zeros(length(S)); for i=1:length(S) N(i)=bin2dec(num2str(code(i,2:4)))+1; %判断该量化值所在的段位 end a=[0,16,32,64,128,256,512,1024]; b=[1,1,2,4,8,16,32,64]; for i=1:length(S) q=ceil((Q(i)-a(N(i)))/b(N(i))); %求段内码（10进制） if q==0 code(i,(5:8))=[0,0,0,0]; else k=num2str(dec2bin(q-1,4)); %进行段内码10进制与2进制之间的转换 code(i,5)=str2num(k(1)); code(i,6)=str2num(k(2)); code(i,7)=str2num(k(3)); code(i,8)=str2num(k(4)); end if z(i)>0 %判断编码的正负（8位码组的第1位） code(i,1)=1; elseif z(i)<0 code(i,1)=0; end end code = reshape(code', 1, []); end