13折线a律编码译码输出_matlab实现A律13折线的编码和译码以及量化误差的计算

%主函数如下

clear all;

clc;

figure(7)

fc=2000;

t=0:1/2000:1;

m=2000sin(100pit);

subplot(311);

plot(t,m);%抽样序列

axis([0 0.02 -2000 2000]);

title(‘抽样序列’);

xlabel(‘t(s)’);

for i=1:length(m)

if m(i)>0 %确定极性码

B1=1;

else

B1=0;

end

m(i)=abs(m(i));

C=[0 16 32 64 128 256 512 1024 2048];%A律13折现非均匀区间

for j=1:length©;

if m(i)>=C(j)&&m(i)<=C(j+1)

L=j-1;

L1=dec2bin(L,3);

B2=L1(1);

B2=str2num(B2);

B3=L1(2);

B3=str2num(B3);

B4=L1(3);

B4=str2num(B4);

end

end

a=C(L+1);%确定段内码

b=C(L+2);

[B5,a1,b1]=judge(m(i),a,b);

[B6,a2,b2]=judge(m(i),a1,b1);