clear
clc
a=1024; %a表示元素个数，可以修改码元为任意个数如1024
r=unidrnd(2,1,a-1); %随机产生离散均匀二值分布（1，2），共a-1个
rn=r-1; %每个元素减1变为（0，1）二值随机数列
rn=[1,rn]; %跟书上一致，从1开始随机产生离散均匀二值分布（1，2），共a个
ori=rn; %将原序列保存起来，便于后面解码后比较
ori
l=1;
for n=1:a
if rn(n)==1 ;
l=-l; %每检测到一个1变一次符号
rn(n)=l;
end
end
AMI=rn; %变为AMI码：每过一个1就使原来的1变一次号
AMI
m=-1; %符号初始化
for k=1:a-3
if abs(rn(k))==1
m=rn(k); %m记录可能存在的连续4个0前面非零码元的符号
end
if rn(k)==0
if rn(k+1)==0
if rn(k+2)==0
if rn(k+3)==0 %检测是否连续4个码元都是0
rn(k+3)=2*m; %用2表示书上的V
end
end
end
end
end
rn ; %插入V(2)后 checked
l=1; %l=1表示两个V之间有偶数个非零码，0是偶数，所以初始为1，l=-1表示两个V之间有奇数个非零码
for s=1:a
if abs(rn(s))==2
for d=s+1:a
if abs(rn(d))==1 %V之间的非零码只有+1，-1遇到一个一使l变一次号，表示偶数或奇数
l=-l;
end
if abs(rn(d))==2
if l~=1
l=1; %不需要插入B时若检测到V也应把计数器清零
%最初版本的升级处1
else %检测到下一个V时，若为偶数个，则插入B
rn(d-3)=3*(-(rn(d-4))/(abs(rn(d-4)))); %3代表B,插入不带符号的B
rn(d:a)=-rn(d:a);
end %V后的符号再交替
end
end
end
end
rn %到此处完成了插入不带符号的B
%for s=1:a
% if abs(rn(s))==3 %找到B
%rn(s)=rn(s)*(-(rn(s-1))/(abs(rn(s-1)))); %B（3）与前一个位置的带符号的归一值相乘再取反，实现符号B(3)的极性与前一非符号的相反
%B后面第三个就是接下来的V,从它开始非零码正负号交替变化
%end
%end
HDB=rn; %给B(3)添加了符号，并且实现了V后的符号再交替
HDB ;
%以上便实现了HDB3的编码
%下面进行解码
for k=1:a
if abs(rn(k))==2
rn(k-3)=0;rn(k)=0; %每个V（2）都变成0，V前面第三个有可能是B(3)有可能是0也恢复为0
end
if abs(rn(k))==1
rn(k)=1;
end
end
rn; %解码后的恢复序列
rn-ori %解码与原码比较全为0则解码正确
