%judge函数
function [B,m1,n1]=judge(I,m,n)
if I>=m&&I<(m+n)/2
B=0;
m1=m;
n1=(m+n)/2;
else
B=1;
m1=(m+n)/2;
n1=n;
end
%主函数如下
clear all;
clc;
figure(7)
fc=2000;
t=0:1/2000:1;
m=2000*sin(100*pi*t);
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(C);
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